Home
last modified time | relevance | path

Searched refs:gamma (Results 1 – 14 of 14) sorted by relevance

/libCEED/examples/fluids/qfunctions/
H A Dshocktube.h87 const CeedScalar gamma = 1.4; // ratio of specific heats in Exact_ShockTube() local
106 q[4] = P / (gamma - 1.0) + rho * (u[0] * u[0]) / 2.; in Exact_ShockTube()
115 const CeedScalar gamma) { in ConvectiveFluxJacobian_Euler() argument
119 dF[i][j + 1][0] = ((i == j) ? ((gamma - 1.) * (u_sq / 2.)) : 0.) - u[i] * u[j]; in ConvectiveFluxJacobian_Euler()
122 …][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (gamma - 1.); in ConvectiveFluxJacobian_Euler()
123 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
125 dF[i][j + 1][4] = ((i == j) ? (gamma - 1.) : 0.); in ConvectiveFluxJacobian_Euler()
127 dF[i][4][0] = u[i] * ((gamma - 1.) * u_sq - E * gamma / rho); in ConvectiveFluxJacobian_Euler()
128 dF[i][4][4] = u[i] * gamma; in ConvectiveFluxJacobian_Euler()
228 const CeedScalar gamma = 1.4; in EulerShockTube() local
[all …]
H A Deulervortex.h72 const CeedScalar gamma = 1.4; in Exact_Euler() local
84 …const CeedScalar delta_T = -(gamma - 1.) * vortex_strength * vortex_strength * exp(1 - r * r) / (… in Exact_Euler()
86 …const CeedScalar S_bubble = (gamma - 1.) * vortex_strength * vortex_strength / (8. * gamma * M_PI … in Exact_Euler()
96 rho = pow(T / S_vortex, 1 / (gamma - 1.)); in Exact_Euler()
106 q[4] = P / (gamma - 1.) + rho * (u[0] * u[0] + u[1] * u[1]) / 2.; in Exact_Euler()
183 const CeedScalar gamma) { in ConvectiveFluxJacobian_Euler() argument
187 dF[i][j + 1][0] = ((i == j) ? ((gamma - 1.) * (u_sq / 2.)) : 0.) - u[i] * u[j]; in ConvectiveFluxJacobian_Euler()
190 …][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (gamma - 1.); in ConvectiveFluxJacobian_Euler()
191 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
193 dF[i][j + 1][4] = ((i == j) ? (gamma - 1.) : 0.); in ConvectiveFluxJacobian_Euler()
[all …]
H A Dnewtonian_state.h100 const CeedScalar gamma = HeatCapacityRatio(gas); in StateEntropyFromPrimitive() local
102 const CeedScalar entropy = log(Y.pressure) - gamma * log(rho); in StateEntropyFromPrimitive()
106 V.S_density = (gamma - entropy) / (gamma - 1) - rho_div_p * e_kinetic; in StateEntropyFromPrimitive()
114 const CeedScalar gamma = HeatCapacityRatio(gas); in StateEntropyFromPrimitive_fwd() local
122 CeedScalar dentropy = dY.pressure / s.Y.pressure - gamma * drho / s.U.density; in StateEntropyFromPrimitive_fwd()
124 dV.S_density = -dentropy / (gamma - 1) - de_kinetic * rho_div_p - e_kinetic * drho_div_p; in StateEntropyFromPrimitive_fwd()
134 const CeedScalar gamma = HeatCapacityRatio(gas); in StatePrimitiveFromEntropy() local
136 const CeedScalar entropy = gamma - (gamma - 1) * (V.S_density - e_kinetic * V.S_energy); in StatePrimitiveFromEntropy()
137 const CeedScalar log_P = -(entropy + gamma * log(-V.S_energy)) / (gamma - 1); in StatePrimitiveFromEntropy()
147 const CeedScalar gamma = HeatCapacityRatio(gas); in StatePrimitiveFromEntropy_fwd() local
[all …]
H A Dgaussianwave.h41 const CeedScalar gamma = HeatCapacityRatio(newt_ctx); in IC_GaussianWave() local
56 U[4] = S_infty.Y.pressure / (gamma - 1) * perturbation + e_kinetic; in IC_GaussianWave()
H A Driemann_solver.h103 const CeedScalar gamma = HeatCapacityRatio(gas); in ComputeHLLSpeeds_Roe() local
114 CeedScalar a_roe = sqrt((gamma - 1) * (H_roe - 0.5 * Square(u_roe))); in ComputeHLLSpeeds_Roe()
124 const CeedScalar gamma = HeatCapacityRatio(gas); in ComputeHLLSpeeds_Roe_fwd() local
141 CeedScalar a_roe = sqrt((gamma - 1) * (H_roe - 0.5 * Square(u_roe))); in ComputeHLLSpeeds_Roe_fwd()
142 …CeedScalar da_roe = 0.5 * sqrt((gamma - 1) / (H_roe - 0.5 * Square(u_roe))) * dH_roe; // (da/dH) … in ComputeHLLSpeeds_Roe_fwd()
143 …da_roe -= 0.5 * sqrt(gamma - 1) * u_roe / sqrt(H_roe - 0.5 * Square(u_roe)) * du_roe; // (da/du) … in ComputeHLLSpeeds_Roe_fwd()
H A Dchannel.h95 const CeedScalar gamma = HeatCapacityRatio(&context->newtonian_ctx); in Channel_Inflow() local
122 const CeedScalar rho_in = P / ((gamma - 1) * e_internal); in Channel_Inflow()
H A Dblasius.h197 const CeedScalar gamma = HeatCapacityRatio(gas); in Blasius_Inflow_Jacobian() local
226 dP = dE * (gamma - 1.); in Blasius_Inflow_Jacobian()
H A Dstg_shur14.h349 const CeedScalar gamma = HeatCapacityRatio(&stg_ctx->newtonian_ctx); in StgShur14Inflow() local
380 P = E_internal * (gamma - 1.); in StgShur14Inflow()
419 const CeedScalar gamma = HeatCapacityRatio(&stg_ctx->newtonian_ctx); in StgShur14Inflow_Jacobian() local
449 dP = dE * (gamma - 1.); in StgShur14Inflow_Jacobian()
/libCEED/examples/fluids/problems/
H A Dnewtonian.c123 const CeedScalar gamma = HeatCapacityRatio(gas); in UnitTests_Newtonian() local
124 const CeedScalar entropy = log(P) - gamma * log(rho); in UnitTests_Newtonian()
128 …const CeedScalar V0[5] = {(gamma - entropy) / (gamma - 1) - rho_div_p * (e_kinetic), rho_div_… in UnitTests_Newtonian()
H A Dblasius.c30 gamma = HeatCapacityRatio(&blasius->newtonian_ctx); in CompressibleBlasiusResidual() local
59 …_rho_tilde[0] * h[2] + mu_rho_tilde[1] * h[1]) + Pr * f[0] * h[1] + Pr * (gamma - 1) * mu_rho_tild… in CompressibleBlasiusResidual()
/libCEED/examples/fluids/
H A DREADME.md613 …py` ($\frac{\gamma - s}{\gamma - 1} - \frac{\rho}{P} (e - c_v T),\ \frac{\rho}{P} \bm{u},\ -\frac{…
H A Dindex.md34 …pressure and $c_v$ is the specific heat at constant volume (that define $\gamma = c_p / c_v$, the …
305 …\hat{\bm n}_i$ is the velocity component in direction $i$ and $a = \sqrt{\gamma P/\rho}$ is the so…
313 Note that this wave speed is specific to ideal gases as $\gamma$ is an ideal gas parameter; other e…
570 …)} \, (-\bar{y}, \, \bar{x}) \, , \\ \delta T &= - \frac{(\gamma-1) \, \epsilon^2}{8 \, \gamma \, …
574 There is no perturbation in the entropy $S=P/\rho^\gamma$ ($\delta S=0)$.
625 E &= \frac{p_\infty}{\gamma -1}\left(1+A\exp\left(\frac{-(\bar{x}^2 + \bar{y}^2)}{2\sigma^2}\right)…
/libCEED/examples/solids/
H A Dindex.md90 …Phi(\bm \epsilon)$; it can only depend on *invariants*, scalar-valued functions $\gamma$ satisfying
93 \gamma(\bm \epsilon) = \gamma(Q \bm \epsilon Q^T)
/libCEED/
H A DDoxyfile1371 # The HTML_COLORSTYLE_GAMMA tag controls the gamma correction applied to the
1374 # darker. The value divided by 100 is the actual gamma applied, so 80 represents
1375 # a gamma of 0.8, The value 220 represents a gamma of 2.2, and 100 does not
1376 # change the gamma.