deleted all the unused coordinate references that were throwing warnings

adding body force term to energy equation...for now the b vector is using what is put into gravity

more style edits

style cleanup

missed some places that BI restoration requires

Potential energy removed from all state helpers so no more x (or dx) dependence for them. Interior integrals still need x for IDL but boundary elements don't so removed. This touches lots of qfuncti

Potential energy removed from all state helpers so no more x (or dx) dependence for them. Interior integrals still need x for IDL but boundary elements don't so removed. This touches lots of qfunctions that make state calls and a few problems as well

show more ...

style errors induced by dropping coordinates where not used

deleted all the unused coordinate references that were throwing warnings

adding body force term to energy equation...for now the b vector is using what is put into gravity

more style edits

style cleanup

missed some places that BI restoration requires

Potential energy removed from all state helpers so no more x (or dx) dependence for them. Interior integrals still need x for IDL but boundary elements don't so removed. This touches lots of qfuncti

Potential energy removed from all state helpers so no more x (or dx) dependence for them. Interior integrals still need x for IDL but boundary elements don't so removed. This touches lots of qfunctions that make state calls and a few problems as well

show more ...

Merge pull request #1330 from CEED/jrwrigh/fluids_physical_grad

fluids: Add StatePhysicalGradientFromReference util

fluids: Add StatePhysicalGradientFromReference util

fluids: Add QF value qdata packing/unpacking

Merge pull request #1189 from CEED/jrwrigh/differential_filter

fluids: Add differential filtering capability

refactor(fluids): Make j, k consisitent

When doing double for loops, `j` should be components and `k` should be
coordinate direction, such that Grad_v[k][j][i] indexes the gradient
array coming from

refactor(fluids): Make j, k consisitent

When doing double for loops, `j` should be components and `k` should be
coordinate direction, such that Grad_v[k][j][i] indexes the gradient
array coming from libCEED

show more ...

Merge pull request #1225 from CEED/jrwrigh/fluids_state_enum

Fluids: Move to state enum for state functions

fluids: Move to state enum for state functions

Merge pull request #1149 from CEED/jrwrigh/newtonian_ddsgs

fluids: Data-driven SGS modeling

fluids: Enable TSSetErrorIfStepFails when testing

- Misc error fixes too
- And style change to match index sizes for similar loops
- Other similar nested for loops have `j<5` and `k<3`, so cha

fluids: Enable TSSetErrorIfStepFails when testing

- Misc error fixes too
- And style change to match index sizes for similar loops
- Other similar nested for loops have `j<5` and `k<3`, so changed
this one to match

show more ...

fluids: Add rotation rate, velocity gradient helpers

- Also change KMStrainRate to KMStrainRate_State, leaving KMStrainRate
to be a generic interface to the velocity gradient

Ken/idl (#1138)

* inflow sponge pushing p towards P0 compiles and runs but unknown yet whether it is effective or not

* looks like my first attempt at this is missing something as there is neglig

Ken/idl (#1138)

* inflow sponge pushing p towards P0 compiles and runs but unknown yet whether it is effective or not

* looks like my first attempt at this is missing something as there is negligible effect from Jacobian addition

* examples/fluids: IDL fixes and demo

* Fix IDL sign error and rename variables/clang-format
* Let newtonianwave.yaml use -reference_pressure and temperature
* Propagate P0 = reference.pressure

* Small tweak to add IDL to the 3 ramp parameters.

* fluids: Extract fluids binary vector load function

- Needs to be used in multiple places (for ICs and for testing data)

* fluids: Correct typo in previous commits

* fluids: Add IDL and newtonian wave test

* fluids: Rename Newtonian Wave --> Gaussian Wave

* fluids: Move RampCoefficient into utils

- Also change name to LinearRampCoefficient, as we may have different
kinds of ramps in the future

* fluids: Change IDL commandline options

* docs(fluids): Document internal damping layer

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

* fluids: Add ScaleN QF utility

* fluids: Extract InternalDampingLayer QFunction helper

---------

Co-authored-by: Jed Brown <jed@jedbrown.org>
Co-authored-by: KennethEJansen <jansenke@colorado.edu>
Co-authored-by: James Wright <james@jameswright.xyz>

show more ...

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