fluids: Add flag for turning on SGS model

fluids: Apply nodal SGS tensor to IFunction

IWYU fixes (#1182)

* iwyu - include fixes

* iwyu - silence some iwyu output

* minor - clearer macro names

* iwyu - fix suggestion of "ceed/ceed.h" externally

* iwyu - lighter petsc heade

IWYU fixes (#1182)

* iwyu - include fixes

* iwyu - silence some iwyu output

* minor - clearer macro names

* iwyu - fix suggestion of "ceed/ceed.h" externally

* iwyu - lighter petsc headers

* iwyu - ceed/ceed.h -> ceed.h

* iwyu - cuda/hip include fixes

show more ...

Merge pull request #1177 from CEED/jrwrigh/fluids-getcontextfield

fluids: Use ContextGetDouble for time and dt label updates

Merge pull request #1174 from CEED/jeremy/operator-context

CeedOperatorContext* -> CeedOperator*Context*

style - rename CeedOperatorContext functions for consistency

Fix for cubic elements in surface forces for lift/drag coefficients (#1158)

* update for surface forces to return dof

* update surface forces

* update on surface forces

* PetscCall section

Fix for cubic elements in surface forces for lift/drag coefficients (#1158)

* update for surface forces to return dof

* update surface forces

* update on surface forces

* PetscCall section get field components

* Update examples/fluids/src/setupts.c

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

* Update examples/fluids/src/setupts.c

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

* examples/fluids: lift query from inner loop

* examples/fluids: enable coarser meshes

* update yaml file

* fix path to creating figure

---------

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

show more ...

Merge pull request #1147 from CEED/jrwrigh/generic_l2_projection

fluids: Add generic L2 projection function

fluids: Add VecP2C, VecC2P, and friends

Lift, Drag, Heat Transfer (#1130)

* add functions for surface forces

* function for surface forces

* reaction forces

* add points to label

* WIP: examples/fluids: sketch surface force ca

Lift, Drag, Heat Transfer (#1130)

* add functions for surface forces

* function for surface forces

* reaction forces

* add points to label

* WIP: examples/fluids: sketch surface force calculations

* examples/fluids: fix for parallel, when there are no such points on a process

* add TSMonitor_FaceForce function

* update on TSMonitor FaceForce

* create .csv file

* cleanup TS Monitor FaceForce function

* include surface forces viewer to header file

* include ts monitor surface force viewer to cloptions

* update where TSMonitor FaceForce are called

* examples/fluids: update wall forces to write csv

mpiexec -n 5 build/fluids-navierstokes -options_file examples/fluids/vortexshedding.yaml -{ts,snes}_monitor -degree 1 -ts_monitor_wall_forces ascii:force.csv:ascii_csv

* update documentation

* reaction force data

* plotting script

* examples/fluids: seaborn plotting with twinx and colors

* delete data file

* tuck away csv file

* place python script in subdirectory

* examples/fluids/postprocess: put shedding period in title

* examples/fluids: document wall force and shedding postprocess

---------

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

show more ...

Merge pull request #1135 from CEED/jrwrigh/spanwise_stats

fluids: Add UpdateContextLabel utility

fluids: Stats collection weighted in QFunction

- Allows for only keeping one CeedVector for child statistics

fluids: Add turb spanstats test, change test option

- `-test` changed to `-test_type` enum to test different outputs

fluids: Change turb_spanstats CL option names

- Options controlling the actual TSMonitor are under
`ts_monitor_turbulence_spanstats_*`
- `stats_enable` removed in favor of
`ts_monitor_turbulence

fluids: Change turb_spanstats CL option names

- Options controlling the actual TSMonitor are under
`ts_monitor_turbulence_spanstats_*`
- `stats_enable` removed in favor of
`ts_monitor_turbulence_spanstats_viewer` being the statistics switch
- Related objects excluding the viewer (such as the dm and ksp) get the
prefix `turbulence_spanstats`
- Change the `stats_test` related functions and variables to use the
`mms` nomenclature, since they're method of manufactured solutions.

show more ...

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: 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

fluids: Implement statistics collection

fluids: Allow communication overlap in IFunction

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

Merge pull request #1134 from CEED/jed/fluids-ts-reason-no-error

fluids: report TS converged reason and exit cleanly

examples/fluids TS: report converged/diverged reason with clean exit

The final state does not get written by the monitor with the default of
-ts_error_if_step_fails 1. Now we exit cleanly regardless

examples/fluids TS: report converged/diverged reason with clean exit

The final state does not get written by the monitor with the default of
-ts_error_if_step_fails 1. Now we exit cleanly regardless of reason,
ensuring that the final state checkpoint and monitor has been called.

show more ...