examples/fluids: fix stats after restart

* Switch from reweighted average to running sum of child stats
* Initialize sp_stats->prev_time after restart
* Use TS reason instead of extra flag
* Skip

examples/fluids: fix stats after restart

* Switch from reweighted average to running sum of child stats
* Initialize sp_stats->prev_time after restart
* Use TS reason instead of extra flag
* Skip redundant monitor call; monitors always get called after last step (converged or not)
* Close memory leak

show more ...

fluids: Create UpdateBoundaryValues, fix statistics

- Statistics would fail when non-zero strong bcs were used due to Q_loc
not containing those values. Instead, use user->Q_loc and call
UpdateB

fluids: Create UpdateBoundaryValues, fix statistics

- Statistics would fail when non-zero strong bcs were used due to Q_loc
not containing those values. Instead, use user->Q_loc and call
UpdateBoundaryValues to verify that the boundaries are set.

show more ...

fluids: Add Stats Setup log stage

fluids: Misc cleanup to CreateStatsDM

fluids: Create Stats DM for periodic face

- Note this relies on the periodicity being isoperiodic rather than
topologicially periodic.

Co-authored-by: Jed Brown <jed@jedbrown.org>

fluids: Use persistent viewer for statistics

- Use prefix 'stats_viewer'
- Particularly useful for setting `stats_viewer_cgns_batch_size`

fluids: Remove StatsCollectFinalCall

- Effectively duplicate of TSMonitor_Statistics

doc(fluids): New spanwise statistics gathering flags

fluids: Destory statistics objects at the end

fluids: Add log stages to statistics

fluids: Check result of statistics testing

fluids: Implement statistics collection

fluids: Add stats collection QF

fluids: Add solution time to newtonian context

fluids: Allow communication overlap in IFunction

fluids: Add MatShell generating code

fluids: Add CreateElemRestrColloc for statistics

- These CeedElemRestrictions are for storing the quadrature point data
at the child points and for their implementation into the parent
quadratur

fluids: Add CreateElemRestrColloc for statistics

- These CeedElemRestrictions are for storing the quadrature point data
at the child points and for their implementation into the parent
quadrature points.

show more ...

fluids: Mass QF for different number of components

- In order to avoid "Function defined, but not used" warnings, the
ReturnMassQFunction is used to return the correctly sized mass QF

fluids: Create SF for child->parent communication

fluids: Add CreateStatsDM using DMPlexFilter

- Use DMPlexFilter to extract target plane from the base mesh

Fix -Wextra-semi -Wmissing-field-initializers

fluids: Use CeedElemRestrictionCreateVector for qdata

examples/fluids: document outflow and add to demos

examples/fluids/meshes: reparametrize wake region and fix Makefile

examples/fluids: Riemann outflow using pressure and temperature

* This allows stable sustained and transient recirculation.

* short-wake domain recirculation behind cylinder has sustained
rec

examples/fluids: Riemann outflow using pressure and temperature

* This allows stable sustained and transient recirculation.

* short-wake domain recirculation behind cylinder has sustained
recirculation. If the exterior temperature is different, the
recirculation will bring in the wrong temperature, but at least
it's stable.
* nwave with slow background causes transient flow reversal when
strong acoustic hits boundary.
* Clean outflow and lid for Blasius

* I tried the 2 P_ext - P_int strategy from Mengaldo (constant density)
and with specified exterior temperature, but they disrupted the interior
flow and were not significantly less reflective.

* The methods I tried for extrapolating interior temperature were not
stable. Simple extrapolation of velocity worked best, and did not show
signs of instability, perhaps because the neighboring flow and implicit
solve ties it together.

* This commit adds support for softmax control of recirculation, but the
parameters were finnicky, solver nonlinearity was a concern, and I'm not
convinced it's ever better than just allowing recirculation.

show more ...