Lines Matching refs:Y
21 StatePrimitive Y; member
30 CEED_QFUNCTION_HELPER void UnpackState_Y(StatePrimitive s, CeedScalar Y[5]) { in UnpackState_Y()
31 Y[0] = s.pressure; in UnpackState_Y()
32 for (int i = 0; i < 3; i++) Y[i + 1] = s.velocity[i]; in UnpackState_Y()
33 Y[4] = s.temperature; in UnpackState_Y()
53 CeedScalar e_kinetic = 0.5 * Dot3(s.Y.velocity, s.Y.velocity); in TotalSpecificEnthalpy()
54 CeedScalar e_internal = gas.cv * s.Y.temperature; in TotalSpecificEnthalpy()
55 return e_internal + e_kinetic + s.Y.pressure / s.U.density; in TotalSpecificEnthalpy()
59 CeedScalar de_kinetic = Dot3(ds.Y.velocity, s.Y.velocity); in TotalSpecificEnthalpy_fwd()
60 CeedScalar de_internal = gas.cv * ds.Y.temperature; in TotalSpecificEnthalpy_fwd()
61 …return de_internal + de_kinetic + ds.Y.pressure / s.U.density - s.Y.pressure / Square(s.U.density)… in TotalSpecificEnthalpy_fwd()
65 StatePrimitive Y; in StatePrimitiveFromConservative() local
66 for (CeedInt i = 0; i < 3; i++) Y.velocity[i] = U.momentum[i] / U.density; in StatePrimitiveFromConservative()
67 CeedScalar e_kinetic = .5 * Dot3(Y.velocity, Y.velocity); in StatePrimitiveFromConservative()
70 Y.temperature = e_internal / gas.cv; in StatePrimitiveFromConservative()
71 Y.pressure = (HeatCapacityRatio(gas) - 1) * U.density * e_internal; in StatePrimitiveFromConservative()
72 return Y; in StatePrimitiveFromConservative()
78 dY.velocity[i] = (dU.momentum[i] - s.Y.velocity[i] * dU.density) / s.U.density; in StatePrimitiveFromConservative_fwd()
80 CeedScalar e_kinetic = .5 * Dot3(s.Y.velocity, s.Y.velocity); in StatePrimitiveFromConservative_fwd()
81 CeedScalar de_kinetic = Dot3(dY.velocity, s.Y.velocity); in StatePrimitiveFromConservative_fwd()
91 …NCTION_HELPER StateEntropy StateEntropyFromPrimitive(NewtonianIGProperties gas, StatePrimitive Y) { in StateEntropyFromPrimitive() argument
94 const CeedScalar rho = Y.pressure / (GasConstant(gas) * Y.temperature); in StateEntropyFromPrimitive()
95 const CeedScalar entropy = log(Y.pressure) - gamma * log(rho); in StateEntropyFromPrimitive()
96 const CeedScalar rho_div_p = rho / Y.pressure; in StateEntropyFromPrimitive()
97 const CeedScalar e_kinetic = 0.5 * Dot3(Y.velocity, Y.velocity); in StateEntropyFromPrimitive()
100 for (int i = 0; i < 3; i++) V.S_momentum[i] = rho_div_p * Y.velocity[i]; in StateEntropyFromPrimitive()
108 … drho = (dY.pressure * s.Y.temperature - s.Y.pressure * dY.temperature) / (GasConstant(gas) … in StateEntropyFromPrimitive_fwd()
110 const CeedScalar e_kinetic = .5 * Dot3(s.Y.velocity, s.Y.velocity); in StateEntropyFromPrimitive_fwd()
111 const CeedScalar de_kinetic = Dot3(dY.velocity, s.Y.velocity); in StateEntropyFromPrimitive_fwd()
112 const CeedScalar rho_div_p = s.U.density / s.Y.pressure; in StateEntropyFromPrimitive_fwd()
113 …const CeedScalar drho_div_p = (drho * s.Y.pressure - s.U.density * dY.pressure) / Square(s.Y.press… in StateEntropyFromPrimitive_fwd()
115 CeedScalar dentropy = dY.pressure / s.Y.pressure - gamma * drho / s.U.density; in StateEntropyFromPrimitive_fwd()
118 …for (CeedInt i = 0; i < 3; i++) dV.S_momentum[i] = rho_div_p * dY.velocity[i] + drho_div_p * s.Y.v… in StateEntropyFromPrimitive_fwd()
124 StatePrimitive Y; in StatePrimitiveFromEntropy() local
125 for (int i = 0; i < 3; i++) Y.velocity[i] = -V.S_momentum[i] / V.S_energy; in StatePrimitiveFromEntropy()
126 Y.temperature = -1 / (GasConstant(gas) * V.S_energy); in StatePrimitiveFromEntropy()
128 const CeedScalar e_kinetic = 0.5 * Dot3(Y.velocity, Y.velocity); in StatePrimitiveFromEntropy()
131 Y.pressure = exp(log_P); in StatePrimitiveFromEntropy()
132 return Y; in StatePrimitiveFromEntropy()
137 StateEntropy V = StateEntropyFromPrimitive(gas, s.Y); in StatePrimitiveFromEntropy_fwd()
141 const CeedScalar e_kinetic = 0.5 * Dot3(s.Y.velocity, s.Y.velocity); in StatePrimitiveFromEntropy_fwd()
142 const CeedScalar de_kinetic = Dot3(dY.velocity, s.Y.velocity); in StatePrimitiveFromEntropy_fwd()
144 …dY.pressure = s.Y.pressure * (-dentropy - gamma * dV.S_energy / V.S_energy) / (gam… in StatePrimitiveFromEntropy_fwd()
148 …PER StateConservative StateConservativeFromPrimitive(NewtonianIGProperties gas, StatePrimitive Y) { in StateConservativeFromPrimitive() argument
150 U.density = Y.pressure / (GasConstant(gas) * Y.temperature); in StateConservativeFromPrimitive()
151 for (int i = 0; i < 3; i++) U.momentum[i] = U.density * Y.velocity[i]; in StateConservativeFromPrimitive()
152 CeedScalar e_internal = gas.cv * Y.temperature; in StateConservativeFromPrimitive()
153 CeedScalar e_kinetic = .5 * Dot3(Y.velocity, Y.velocity); in StateConservativeFromPrimitive()
161 …dU.density = (dY.pressure * s.Y.temperature - s.Y.pressure * dY.temperature) / (GasConstant(gas) *… in StateConservativeFromPrimitive_fwd()
163 dU.momentum[i] = dU.density * s.Y.velocity[i] + s.U.density * dY.velocity[i]; in StateConservativeFromPrimitive_fwd()
165 CeedScalar e_kinetic = .5 * Dot3(s.Y.velocity, s.Y.velocity); in StateConservativeFromPrimitive_fwd()
166 CeedScalar de_kinetic = Dot3(dY.velocity, s.Y.velocity); in StateConservativeFromPrimitive_fwd()
167 CeedScalar e_internal = gas.cv * s.Y.temperature; in StateConservativeFromPrimitive_fwd()
195 const CeedScalar p = s.Y.pressure; in StateEntropyFromConservative_fwd()
227 StateEntropy V = StateEntropyFromPrimitive(gas, s.Y); in StateConservativeFromEntropy_fwd()
228 …for (int i = 0; i < 3; i++) dvelocity[i] = (-dV.S_momentum[i] - s.Y.velocity[i] * dV.S_energy) / V… in StateConservativeFromEntropy_fwd()
230 const CeedScalar e_kinetic = 0.5 * Dot3(s.Y.velocity, s.Y.velocity); in StateConservativeFromEntropy_fwd()
231 const CeedScalar de_kinetic = Dot3(dvelocity, s.Y.velocity); in StateConservativeFromEntropy_fwd()
237 …for (int i = 0; i < 3; i++) dU.momentum[i] = dU.density * s.Y.velocity[i] + s.U.density * dvelocit… in StateConservativeFromEntropy_fwd()
246 CEED_QFUNCTION_HELPER State StateFromPrimitive(NewtonianIGProperties gas, StatePrimitive Y) { in StateFromPrimitive() argument
247 StateConservative U = StateConservativeFromPrimitive(gas, Y); in StateFromPrimitive()
250 s.Y = Y; in StateFromPrimitive()
258 ds.Y = dY; in StateFromPrimitive_fwd()
273 …rvativeAXPBYPCZ(CeedScalar a, StateConservative X, CeedScalar b, StateConservative Y, CeedScalar c, in StateConservativeAXPBYPCZ() argument
276 R.density = a * X.density + b * Y.density + c * Z.density; in StateConservativeAXPBYPCZ()
277 …for (int i = 0; i < 3; i++) R.momentum[i] = a * X.momentum[i] + b * Y.momentum[i] + c * Z.momentum… in StateConservativeAXPBYPCZ()
278 R.E_total = a * X.E_total + b * Y.E_total + c * Z.E_total; in StateConservativeAXPBYPCZ()
284 …teToY(NewtonianIGProperties gas, const State input, CeedScalar Y[5]) { UnpackState_Y(input.Y, Y); } in StateToY()
287 StateEntropy state_V = StateEntropyFromPrimitive(gas, input.Y); in StateToV()
317 dstate_v = StateEntropyFromPrimitive_fwd(gas, input, dinput.Y); in StateToQ_fwd()
330 s.Y = StatePrimitiveFromConservative(gas, s.U); in StateFromU()
341 ds.Y = StatePrimitiveFromConservative_fwd(gas, s, ds.U); in StateFromU_fwd()
345 CEED_QFUNCTION_HELPER State StateFromY(NewtonianIGProperties gas, const CeedScalar Y[5]) { in StateFromY()
347 s.Y.pressure = Y[0]; in StateFromY()
348 s.Y.velocity[0] = Y[1]; in StateFromY()
349 s.Y.velocity[1] = Y[2]; in StateFromY()
350 s.Y.velocity[2] = Y[3]; in StateFromY()
351 s.Y.temperature = Y[4]; in StateFromY()
352 s.U = StateConservativeFromPrimitive(gas, s.Y); in StateFromY()
358 ds.Y.pressure = dY[0]; in StateFromY_fwd()
359 ds.Y.velocity[0] = dY[1]; in StateFromY_fwd()
360 ds.Y.velocity[1] = dY[2]; in StateFromY_fwd()
361 ds.Y.velocity[2] = dY[3]; in StateFromY_fwd()
362 ds.Y.temperature = dY[4]; in StateFromY_fwd()
363 ds.U = StateConservativeFromPrimitive_fwd(gas, s, ds.Y); in StateFromY_fwd()
376 s.Y = StatePrimitiveFromEntropy(gas, state_V); in StateFromV()
389 ds.Y = StatePrimitiveFromEntropy_fwd(gas, s, state_dV); in StateFromV_fwd()
428 …for (CeedInt j = 0; j < 3; j++) Flux[i].momentum[j] = s.U.momentum[i] * s.Y.velocity[j] + s.Y.pres… in FluxInviscid()
429 Flux[i].E_total = (s.U.E_total + s.Y.pressure) * s.Y.velocity[i]; in FluxInviscid()
437 …dFlux[i].momentum[j] = ds.U.momentum[i] * s.Y.velocity[j] + s.U.momentum[i] * ds.Y.velocity[j] + d… in FluxInviscid_fwd()
439 …dFlux[i].E_total = (ds.U.E_total + ds.Y.pressure) * s.Y.velocity[i] + (s.U.E_total + s.Y.pressure)… in FluxInviscid_fwd()
512 grad_velocity[0][0] = grad_s[0].Y.velocity[0]; in VelocityGradient()
513 grad_velocity[0][1] = grad_s[1].Y.velocity[0]; in VelocityGradient()
514 grad_velocity[0][2] = grad_s[2].Y.velocity[0]; in VelocityGradient()
515 grad_velocity[1][0] = grad_s[0].Y.velocity[1]; in VelocityGradient()
516 grad_velocity[1][1] = grad_s[1].Y.velocity[1]; in VelocityGradient()
517 grad_velocity[1][2] = grad_s[2].Y.velocity[1]; in VelocityGradient()
518 grad_velocity[2][0] = grad_s[0].Y.velocity[2]; in VelocityGradient()
519 grad_velocity[2][1] = grad_s[1].Y.velocity[2]; in VelocityGradient()
520 grad_velocity[2][2] = grad_s[2].Y.velocity[2]; in VelocityGradient()
560 CEED_QFUNCTION_HELPER void ViscousEnergyFlux(NewtonianIGProperties gas, StatePrimitive Y, const Sta… in ViscousEnergyFlux() argument
563 …Fe[i] = -Y.velocity[0] * stress[0][i] - Y.velocity[1] * stress[1][i] - Y.velocity[2] * stress[2][i… in ViscousEnergyFlux()
567 CEED_QFUNCTION_HELPER void ViscousEnergyFlux_fwd(NewtonianIGProperties gas, StatePrimitive Y, State… in ViscousEnergyFlux_fwd() argument
570 …dFe[i] = -Y.velocity[0] * dstress[0][i] - dY.velocity[0] * stress[0][i] - Y.velocity[1] * dstress[… in ViscousEnergyFlux_fwd()
571 … Y.velocity[2] * dstress[2][i] - dY.velocity[2] * stress[2][i] - gas.k * grad_ds[i].Y.temperature; in ViscousEnergyFlux_fwd()