Home
last modified time | relevance | path

Searched refs:gas (Results 1 – 25 of 30) sorted by relevance

12

/honee/qfunctions/
H A Dnewtonian_state.h42 CEED_QFUNCTION_HELPER CeedScalar HeatCapacityRatio(NewtonianIGProperties gas) { return gas.cp / gas… in HeatCapacityRatio() argument
44 CEED_QFUNCTION_HELPER CeedScalar GasConstant(NewtonianIGProperties gas) { return gas.cp - gas.cv; } in GasConstant() argument
46 CEED_QFUNCTION_HELPER CeedScalar Prandtl(NewtonianIGProperties gas) { return gas.cp * gas.mu / gas.… in Prandtl() argument
48 …CeedScalar SoundSpeed(NewtonianIGProperties gas, CeedScalar T) { return sqrt(gas.cp * (HeatCapacit… in SoundSpeed() argument
50 …_HELPER CeedScalar Mach(NewtonianIGProperties gas, CeedScalar T, CeedScalar u) { return u / SoundS… in Mach() argument
52 CEED_QFUNCTION_HELPER CeedScalar TotalSpecificEnthalpy(NewtonianIGProperties gas, const State s) { in TotalSpecificEnthalpy() argument
54 CeedScalar e_internal = gas.cv * s.Y.temperature; in TotalSpecificEnthalpy()
58 CEED_QFUNCTION_HELPER CeedScalar TotalSpecificEnthalpy_fwd(NewtonianIGProperties gas, const State s… in TotalSpecificEnthalpy_fwd() argument
60 CeedScalar de_internal = gas.cv * ds.Y.temperature; in TotalSpecificEnthalpy_fwd()
64 CEED_QFUNCTION_HELPER StatePrimitive StatePrimitiveFromConservative(NewtonianIGProperties gas, Stat… in StatePrimitiveFromConservative() argument
[all …]
H A Dnewtonian.h20 NewtonianIGProperties gas = context->newt_ctx.gas; in ICsNewtonianIG() local
24 State s = StateFromPrimitive(gas, context->reference); in ICsNewtonianIG()
25 StateToQ(gas, s, q, state_var); in ICsNewtonianIG()
51 NewtonianIGProperties gas = context->gas; in MassFunction_Newtonian() local
56 const State s = StateFromQ(gas, qi, state_var); in MassFunction_Newtonian()
57 const State s_dot = StateFromQ(gas, qi_dot, state_var); in MassFunction_Newtonian()
70 Tau_diagPrim(context->tau_coeffs, gas, s, dXdx, context->dt, Tau_d); in MassFunction_Newtonian()
71 Stabilization(context->stabilization, gas, s, Tau_d, grad_s, U_dot, body_force, divFdiff, stab); in MassFunction_Newtonian()
88 CEED_QFUNCTION_HELPER void InternalDampingLayer_Residual(const NewtonianIGProperties gas, const Sta… in InternalDampingLayer_Residual() argument
91 State damp_s = StateFromY_fwd(gas, s, damp_Y); in InternalDampingLayer_Residual()
[all …]
H A Dbc_outflow.h45 const NewtonianIGProperties gas = outflow->newt_ctx.gas; in RiemannOutflow() local
53 const State s_int = StateFromQ(gas, qi, state_var); in RiemannOutflow()
63 State s_ext = StateFromPrimitive(gas, y_ext); in RiemannOutflow()
66 StatePhysicalGradientFromReference_Boundary(Q, i, gas, s_int, state_var, Grad_q, dXdx, grad_s); in RiemannOutflow()
70 NewtonianStress(gas, strain_rate, kmstress); in RiemannOutflow()
72 ViscousEnergyFlux(gas, s_int.Y, grad_s, stress, Fe); in RiemannOutflow()
74 StateConservative F_inviscid_normal = RiemannFlux_HLLC(gas, s_int, s_ext, normal); in RiemannOutflow()
114 const NewtonianIGProperties gas = outflow->newt_ctx.gas; in RiemannOutflow_Jacobian() local
127 State s_int = StateFromQ(gas, qi, state_var); in RiemannOutflow_Jacobian()
128 const State ds_int = StateFromQ_fwd(gas, s_int, dqi, state_var); in RiemannOutflow_Jacobian()
[all …]
H A Driemann_solver.h97 CEED_QFUNCTION_HELPER void ComputeHLLSpeeds_Roe(NewtonianIGProperties gas, State left, CeedScalar u… in ComputeHLLSpeeds_Roe() argument
99 const CeedScalar gamma = HeatCapacityRatio(gas); in ComputeHLLSpeeds_Roe()
107 CeedScalar H_left = TotalSpecificEnthalpy(gas, left); in ComputeHLLSpeeds_Roe()
108 CeedScalar H_right = TotalSpecificEnthalpy(gas, right); in ComputeHLLSpeeds_Roe()
117 CEED_QFUNCTION_HELPER void ComputeHLLSpeeds_Roe_fwd(NewtonianIGProperties gas, State left, State dl… in ComputeHLLSpeeds_Roe_fwd() argument
120 const CeedScalar gamma = HeatCapacityRatio(gas); in ComputeHLLSpeeds_Roe_fwd()
130 CeedScalar H_left = TotalSpecificEnthalpy(gas, left); in ComputeHLLSpeeds_Roe_fwd()
131 CeedScalar H_right = TotalSpecificEnthalpy(gas, right); in ComputeHLLSpeeds_Roe_fwd()
132 CeedScalar dH_left = TotalSpecificEnthalpy_fwd(gas, left, dleft); in ComputeHLLSpeeds_Roe_fwd()
133 CeedScalar dH_right = TotalSpecificEnthalpy_fwd(gas, right, dright); in ComputeHLLSpeeds_Roe_fwd()
[all …]
H A Dblasius.h70 CeedScalar mu = blasius->newt_ctx.gas.mu; in BlasiusSolution()
76 CeedScalar Rd = GasConstant(blasius->newt_ctx.gas); in BlasiusSolution()
90 return StateFromY(blasius->newt_ctx.gas, Y); in BlasiusSolution()
101 const NewtonianIGProperties gas = context->newt_ctx.gas; in ICsBlasius() local
102 const CeedScalar mu = context->newt_ctx.gas.mu; in ICsBlasius()
117 StateToQ(gas, s, q, context->newt_ctx.state_var); in ICsBlasius()
132 const NewtonianIGProperties gas = context->newt_ctx.gas; in Blasius_Inflow() local
136 …const CeedScalar x0 = U_infty * rho_0 / (gas.mu * 25 / Square(context->delta0)… in Blasius_Inflow()
149 State s_int = StateFromU(gas, qi); in Blasius_Inflow()
154 s.Y = StatePrimitiveFromConservative(gas, s.U); in Blasius_Inflow()
[all …]
H A Dstabilization.h24 CEED_QFUNCTION_HELPER void StabilizationMatrix(const NewtonianIGProperties gas, const State s, cons… in StabilizationMatrix() argument
31 State ds = StateFromY_fwd(gas, s, dY); in StabilizationMatrix()
32 FluxInviscid_fwd(gas, s, ds, dF); in StabilizationMatrix()
40 …oid Stabilization(StabilizationType stab_type, const NewtonianIGProperties gas, const State s, con… in Stabilization() argument
49 FluxInviscidStrong(gas, s, ds, strong_residual); in Stabilization()
52 FluxInviscidStrong(gas, s, ds, strong_residual); in Stabilization()
56 StabilizationMatrix(gas, s, Tau_d, strong_residual, stab); in Stabilization()
66 CEED_QFUNCTION_HELPER void Tau_diagPrim(TauDiagCoefficients tau_coeffs, NewtonianIGProperties gas, … in Tau_diagPrim() argument
74 const CeedScalar cv = gas.cv; in Tau_diagPrim()
75 const CeedScalar mu = gas.mu; in Tau_diagPrim()
H A Dchannel.h31 NewtonianIGProperties gas = context->newt_ctx.gas; in Exact_Channel() local
32 const CeedScalar cp = gas.cp; in Exact_Channel()
33 const CeedScalar mu = gas.mu; in Exact_Channel()
34 const CeedScalar k = gas.k; in Exact_Channel()
52 return StateFromY(gas, Y); in Exact_Channel()
63 NewtonianIGProperties gas = context->newt_ctx.gas; in ICsChannel() local
69 StateToQ(gas, s, q, context->newt_ctx.state_var); in ICsChannel()
86 NewtonianIGProperties gas = context->newt_ctx.gas; in Channel_Inflow() local
87 const CeedScalar gamma = HeatCapacityRatio(gas); in Channel_Inflow()
109 State s_inside = StateFromU(gas, q_inside); in Channel_Inflow()
[all …]
H A Dmonitor_totalkineticenergy.h15 const NewtonianIGProperties gas = newt_ctx->gas; in MonitorTotalKineticEnergy() local
18 const State s = StateFromQ(gas, qi, state_var); in MonitorTotalKineticEnergy()
23 StatePhysicalGradientFromReference(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s); in MonitorTotalKineticEnergy()
30 v[1][i] = wdetJ * -2 * gas.mu * DotN((CeedScalar *)strain_rate, (CeedScalar *)strain_rate, 9); in MonitorTotalKineticEnergy()
31 v[2][i] = wdetJ * -gas.lambda * gas.mu * Square(div_u); in MonitorTotalKineticEnergy()
35 v[4][i] = wdetJ * gas.mu * Dot3(vorticity, vorticity); in MonitorTotalKineticEnergy()
H A Dbc_slip.h18 const NewtonianIGProperties gas = newt_ctx->gas; in Slip() local
22 State s = StateFromQ(gas, qi, state_var); in Slip()
32 State s_reflect = StateFromY(gas, Y_reflect); in Slip()
34 StateConservative flux = RiemannFlux_HLLC(gas, s, s_reflect, normal); in Slip()
65 const NewtonianIGProperties gas = newt_ctx->gas; in Slip_Jacobian() local
75 State s = StateFromQ(gas, qi, state_var); in Slip_Jacobian()
76 State ds = StateFromQ_fwd(gas, s, dqi, state_var); in Slip_Jacobian()
82 State s_reflect = StateFromY(gas, Y_reflect); in Slip_Jacobian()
88 State ds_reflect = StateFromY_fwd(gas, s_reflect, dY_reflect); in Slip_Jacobian()
90 StateConservative dflux = RiemannFlux_HLLC_fwd(gas, s, ds, s_reflect, ds_reflect, normal); in Slip_Jacobian()
H A Ddensitycurrent.h89 NewtonianIGProperties gas = context->newt_ctx.gas; in Exact_DC() local
90 const CeedScalar cp = gas.cp; in Exact_DC()
91 const CeedScalar Rd = GasConstant(gas); in Exact_DC()
120 return StateFromY(gas, Y); in Exact_DC()
131 NewtonianIGProperties gas = context->newt_ctx.gas; in ICsDC() local
137 StateToQ(gas, s, q, context->newt_ctx.state_var); in ICsDC()
H A Dbc_freestream.h22 const NewtonianIGProperties gas = context->newt_ctx.gas; in Freestream() local
27 const State s = StateFromQ(gas, qi, state_var); in Freestream()
36 flux = RiemannFlux_HLL(gas, s, context->S_infty, normal); in Freestream()
39 flux = RiemannFlux_HLLC(gas, s, context->S_infty, normal); in Freestream()
84 const NewtonianIGProperties gas = context->newt_ctx.gas; in Freestream_Jacobian() local
96 State s = StateFromQ(gas, qi, state_var); in Freestream_Jacobian()
97 State ds = StateFromQ_fwd(gas, s, dqi, state_var); in Freestream_Jacobian()
102 dflux = RiemannFlux_HLL_fwd(gas, s, ds, context->S_infty, dS_infty, normal); in Freestream_Jacobian()
105 dflux = RiemannFlux_HLLC_fwd(gas, s, ds, context->S_infty, dS_infty, normal); in Freestream_Jacobian()
H A Dstg_shur14.h250 const CeedScalar mu = stg_ctx->newt_ctx.gas.mu; in StgShur14Preprocess()
253 const CeedScalar Rd = GasConstant(stg_ctx->newt_ctx.gas); in StgShur14Preprocess()
290 const NewtonianIGProperties gas = context->gas; in ICsStg() local
296 const CeedScalar rho = P0 / (GasConstant(gas) * theta0); in ICsStg()
297 const CeedScalar nu = gas.mu / rho; in ICsStg()
317 State s = StateFromY(gas, Y); in ICsStg()
318 StateToQ(gas, s, q, context->state_var); in ICsStg()
341 const NewtonianIGProperties gas = stg_ctx->newt_ctx.gas; in StgShur14Inflow() local
347 const CeedScalar mu = gas.mu; in StgShur14Inflow()
351 const CeedScalar cv = gas.cv; in StgShur14Inflow()
[all …]
H A Dtaylorgreen.h29 const NewtonianIGProperties gas = context->newt_ctx.gas; in ICsTaylorGreen() local
30 CeedScalar R = GasConstant(gas); in ICsTaylorGreen()
46 State s = StateFromY(gas, Y); in ICsTaylorGreen()
47 StateToQ(gas, s, q, context->newt_ctx.state_var); in ICsTaylorGreen()
H A Dgaussianwave.h27 const NewtonianIGProperties gas = newt_ctx->gas; in IC_GaussianWave() local
35 const CeedScalar gamma = HeatCapacityRatio(gas); in IC_GaussianWave()
52 State initCond = StateFromU(gas, U); in IC_GaussianWave()
53 StateToQ(gas, initCond, qi, state_var); in IC_GaussianWave()
H A Dvelocity_gradient_projection.h16 NewtonianIGProperties gas = ((NewtonianIdealGasContext)ctx)->gas; in VelocityGradientProjectionRHS() local
23 const State s = StateFromQ(gas, qi, state_var); in VelocityGradientProjectionRHS()
25 StatePhysicalGradientFromReference(Q, i, gas, s, state_var, (CeedScalar *)Grad_q, dXdx, grad_s); in VelocityGradientProjectionRHS()
H A Dsgs_dd_training.h33 const NewtonianIGProperties gas = sgsdd_ctx->newt_ctx.gas; in ComputeSGS_DDAnisotropicTrainingDataNodal() local
44 const State s = StateFromQ(gas, qi, state_var); in ComputeSGS_DDAnisotropicTrainingDataNodal()
48 …ComputeSgsDDInputs(grad_velo_aniso, km_A_ij, delta, gas.mu / s.U.density, eigenvectors, inputs, &g… in ComputeSGS_DDAnisotropicTrainingDataNodal()
H A Dsgs_dd_model.h76 const NewtonianIGProperties gas = sgsdd_ctx->newt_ctx.gas; in ComputeSgsDDNodal_Fused() local
87 const State s = StateFromQ(gas, qi, state_var); in ComputeSgsDDNodal_Fused()
90 ComputeSgsDD_Fused(grad_velo_aniso, km_A_ij, delta, gas.mu / s.U.density, km_sgs, sgsdd_ctx); in ComputeSgsDDNodal_Fused()
120 const NewtonianIGProperties gas = sgsdd_ctx->newt_ctx.gas; in ComputeSgsDDNodal_Sequential_Inputs() local
131 const State s = StateFromQ(gas, qi, state_var); in ComputeSgsDDNodal_Sequential_Inputs()
134 …ComputeSgsDDInputs(grad_velo_aniso, km_A_ij, delta, gas.mu / s.U.density, eigenvectors, model_inpu… in ComputeSgsDDNodal_Sequential_Inputs()
225 const State s = StateFromQ(newt_ctx->gas, qi, state_var); in IFunction_NodalSgs()
H A Dmonitor_cfl.h15 const NewtonianIGProperties gas = newt_ctx->gas; in MonitorCFL() local
19 const State s = StateFromQ(gas, qi, state_var); in MonitorCFL()
H A Dadvection.h210 …ysicalGradientFromReference_ND(CeedInt N, CeedInt Q, CeedInt i, NewtonianIGProperties gas, State s, in StatePhysicalGradientFromReference_ND() argument
220 grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference_ND()
223 grad_s[2] = StateFromU(gas, U); in StatePhysicalGradientFromReference_ND()
233 grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference_ND()
278 NewtonianIGProperties gas = {0}; in IFunction_AdvectionGeneric() local
282 const State s = StateFromU(gas, qi); in IFunction_AdvectionGeneric()
287 …StatePhysicalGradientFromReference_ND(dim, Q, i, gas, s, STATEVAR_CONSERVATIVE, grad_q, dXdx, grad… in IFunction_AdvectionGeneric()
356 NewtonianIGProperties gas = {0}; in MassFunction_AdvectionGeneric() local
360 const State s = StateFromU(gas, qi); in MassFunction_AdvectionGeneric()
413 NewtonianIGProperties gas = {0}; in RHSFunction_AdvectionGeneric() local
[all …]
/honee/problems/
H A Dbc_freestream.c17 static PetscErrorCode RiemannSolverUnitTests(NewtonianIGProperties gas, CeedScalar rtol);
141 State S_infty = StateFromPrimitive(newtonian_ig_ctx->gas, Y_inf); in FreestreamBCSetup()
172 if (run_unit_tests) PetscCall(RiemannSolverUnitTests(newtonian_ig_ctx->gas, 5e-7)); in FreestreamBCSetup()
217 static PetscErrorCode TestRiemannHLL_fwd(NewtonianIGProperties gas, CeedScalar rtol_0, CeedScalar r… in TestRiemannHLL_fwd() argument
222 const CeedScalar p = (HeatCapacityRatio(gas) - 1) * rho * gas.cv * T; in TestRiemannHLL_fwd()
230 State left0 = StateFromY(gas, Y0_left); in TestRiemannHLL_fwd()
231 State right0 = StateFromY(gas, Y0_right); in TestRiemannHLL_fwd()
245 State dleft0 = StateFromY_fwd(gas, left0, dY_left); in TestRiemannHLL_fwd()
246 State dright0 = StateFromY_fwd(gas, right0, dY_right); in TestRiemannHLL_fwd()
248 … StateConservative dFlux_state = RiemannFlux_HLL_fwd(gas, left0, dleft0, right0, dright0, normal); in TestRiemannHLL_fwd()
[all …]
H A Dnewtonian.c17 static PetscErrorCode UnitTests_Newtonian(Honee honee, NewtonianIGProperties gas);
389 .gas = in NS_NEWTONIAN_IG()
417 … "Heat capacity at constant volume", NULL, newtonian_ig_ctx->gas.cv, &newtonian_ig_ctx->gas.cv, NU… in NS_NEWTONIAN_IG()
418 …Heat capacity at constant pressure", NULL, newtonian_ig_ctx->gas.cp, &newtonian_ig_ctx->gas.cp, NU… in NS_NEWTONIAN_IG()
419 …ar("-lambda", "Stokes hypothesis second viscosity coefficient", NULL, newtonian_ig_ctx->gas.lambda, in NS_NEWTONIAN_IG()
420 &newtonian_ig_ctx->gas.lambda, NULL)); in NS_NEWTONIAN_IG()
421 …hear dynamic viscosity coefficient", NULL, newtonian_ig_ctx->gas.mu, &newtonian_ig_ctx->gas.mu, NU… in NS_NEWTONIAN_IG()
422 …sScalar("-k", "Thermal conductivity", NULL, newtonian_ig_ctx->gas.k, &newtonian_ig_ctx->gas.k, NUL… in NS_NEWTONIAN_IG()
477 newtonian_ig_ctx->gas.cv *= units->J_per_kg_K; in NS_NEWTONIAN_IG()
478 newtonian_ig_ctx->gas.cp *= units->J_per_kg_K; in NS_NEWTONIAN_IG()
[all …]
H A Dblasius.c24 NewtonianIGProperties gas = blasius->newt_ctx.gas; in CompressibleBlasiusResidual() local
27 …PetscScalar Ma = Mach(gas, S_infty.Y.temperature, U_infty), Pr = Prandtl(gas), gamma = HeatCapacit… in CompressibleBlasiusResidual()
190 State S_infty = StateFromPrimitive(newtonian_ig_ctx->gas, Y_inf); in NS_BLASIUS()
H A Dchannel.c92 channel_ctx->B = body_force_scale * 2 * umax * newtonian_ig_ctx->gas.mu / (H * H); in NS_CHANNEL()
96 CeedScalar cv = newtonian_ig_ctx->gas.cv, cp = newtonian_ig_ctx->gas.cp; in NS_CHANNEL()
/honee/qfunctions/spanstats/
H A Dcflpe.h25 const NewtonianIGProperties gas = context->newt_ctx.gas; in ChildStatsCollection_CflPe() local
30 const State s = StateFromQ(gas, qi, state_var); in ChildStatsCollection_CflPe()
H A Dturbulence.h47 NewtonianIGProperties gas = context->newt_ctx.gas; in ChildStatsCollection() local
54 const State s = StateFromQ(gas, qi, state_var); in ChildStatsCollection()

12