examples/fluids: use Roe speeds directly in HLL

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

fluids: Add newtonian_wave problem

- Mainly used to verify the HLL boundary condition.

fluids: Prevent Chebyshev calc if using STG

fluids: Style, refactor to PetscCall, etc.

fluids: Add STG fluctuations to initial condition

- Extrudes the STG fluctuations downstream of the inflow (using the STG
method's inherent dependence on streamwise distance)

fluids: Integrate freestream bc

fluids: Add SetupBCQFunctions to setuplibceed.c

- Condenses redundant code into a single function. Will be useful when
implementing the freestream bc

fluids: Create AddBCSubOperator in setuplibceed.c

- Condenses redundant code (and will help reduce redundancy when I add
in the freestream boundary conditions as well).

fluids: Create freestream bc qfunction

style(fluids): Visual cleanup of the fluids code

Also satisfy astyle for RiemannFlux_HLL commit

fluids: Add RiemannFlux_HLL helper function

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

fluids: Refactor STG wall distance variable

- Needless shortening.
- Confused me when debugging even though I wrote the code myself...

fluids: Fix STG Ek_tot calculation bug

When the spectrum is small enough, Ek_tot can be zero, causing a
divide-by-zero NaN, which then corrupts the calculation of the velocity
fluctuations later dow

fluids: Fix STG Ek_tot calculation bug

When the spectrum is small enough, Ek_tot can be zero, causing a
divide-by-zero NaN, which then corrupts the calculation of the velocity
fluctuations later down the line.

Originally found during OLCF Crusher Hackathon 4

show more ...

Fluids - Blasius: Some Minor Modifications (#1063)

* Fluids - blasius: default T_wall = T_inf = 288

* Fluids - blasius: break when Chebyshev solve diverges

* Fluids - blasius: add prefix for C

Fluids - Blasius: Some Minor Modifications (#1063)

* Fluids - blasius: default T_wall = T_inf = 288

* Fluids - blasius: break when Chebyshev solve diverges

* Fluids - blasius: add prefix for Chebyshev snes solve

* Fluids - blasius: N_Chebyshev -> n_chebyshev

* Update examples/fluids/tests-output/blasius_test.yaml

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

* Fluids - blasius: Modify the parameters in the regression test

* examples/fluids: Blasius solution is variable density

We use the ideal gas relation that density ~ 1/enthalpy when held at
constant pressure.

* examples/fluids: Blasius slower growth factor

* examples/fluids: fix Chebyshev on GPU

and close memory leak

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

show more ...

fluids: Fix time label for strong STG #1057

Fluids - fix blasius test

Fluids - Some fixes missing from #1039 (#1056)

* Fluids - Blasius: drop theta0, style, and update README

* Fluids - Add regression test for laminar blasius

fluids: Add STATIC and print-% to makefile

fluids: Set solution time label for strong stg

Fluids - Compressible Blasius boundary layer (#1039)

* Fluids - Compressible Blasius boundary layer solution

* Fluids - Blasius: default delta0 = 4.2e-3

* Fluids - Use the new helper functions

Fluids - Compressible Blasius boundary layer (#1039)

* Fluids - Compressible Blasius boundary layer solution

* Fluids - Blasius: default delta0 = 4.2e-3

* Fluids - Use the new helper functions in blasius

* Update examples/fluids/qfunctions/blasius_exact.h

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

* Update examples/fluids/qfunctions/blasius_exact.h

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

* Fluids - remove blasius_exact.h

* Fluids - N_Chebyshev: 10 -> 20

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

show more ...

doc: Update fluids state_var option

fluids: Switch to StateVar enum

fluids: Use newtonian with StateFromQi*_t func pointers

This to "unify" the primitive and conservative formulations of the
newtonian solver. By passing function pointers to a helper function, the
co

fluids: Use newtonian with StateFromQi*_t func pointers

This to "unify" the primitive and conservative formulations of the
newtonian solver. By passing function pointers to a helper function, the
compiler can optimize for a given state function(s).

show more ...

fluids: Use StateFromQi* for newt boundary QFs

fluids: Add FluxTotal_Boundary, Refactor newt BI QFs

- Refactoring the newtonian boundary integral QFunctions