| /honee/qfunctions/ |
| H A D | eulervortex.h | 81 CeedScalar rho, P, T, E, u[3] = {0.}; in Exact_Euler() local
90 rho = pow(T / S_vortex, 1 / (gamma - 1.)); in Exact_Euler()
91 P = rho * T; in Exact_Euler()
96 q[0] = rho; in Exact_Euler()
97 q[1] = rho * u[0]; in Exact_Euler()
98 q[2] = rho * u[1]; in Exact_Euler()
99 q[3] = rho * u[2]; in Exact_Euler()
100 q[4] = P / (gamma - 1.) + rho * (u[0] * u[0] + u[1] * u[1]) / 2.; in Exact_Euler()
103 rho = 1.; in Exact_Euler()
107 q[0] = rho; in Exact_Euler()
[all …]
|
| H A D | shocktube.h | 85 CeedScalar rho, P, u[3] = {0.}; in Exact_ShockTube() local
89 rho = rho_high; in Exact_ShockTube()
92 rho = rho_low; in Exact_ShockTube()
97 q[0] = rho; in Exact_ShockTube()
98 q[1] = rho * u[0]; in Exact_ShockTube()
99 q[2] = rho * u[1]; in Exact_ShockTube()
100 q[3] = rho * u[2]; in Exact_ShockTube()
101 q[4] = P / (gamma - 1.0) + rho * (u[0] * u[0]) / 2.; in Exact_ShockTube()
109 …void ConvectiveFluxJacobian_Euler(CeedScalar dF[3][5][5], const CeedScalar rho, const CeedScalar u… in ConvectiveFluxJacobian_Euler() argument
118 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
[all …]
|
| H A D | channel.h | 153 const CeedScalar rho = q[0][i]; in Channel_Outflow() local
154 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Channel_Outflow()
166 v[0][i] -= wdetJb * rho * u_normal; in Channel_Outflow()
169 for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + norm[j] * P); in Channel_Outflow()
|
| H A D | stg_shur14.h | 254 const CeedScalar rho = P0 / (Rd * theta0); in StgShur14Preprocess() local
255 const CeedScalar nu = mu / rho; in StgShur14Preprocess()
296 const CeedScalar rho = P0 / (GasConstant(gas) * theta0); in ICsStg() local
297 const CeedScalar nu = gas.mu / rho; in ICsStg()
356 const CeedScalar rho = prescribe_T ? q[0][i] : P0 / (Rd * theta0); in StgShur14Inflow() local
369 CalcSpectrum(X[1][i], eps, lt, h_node_sep, mu / rho, qn, stg_ctx); in StgShur14Inflow()
375 const CeedScalar E_kinetic = .5 * rho * Dot3(u, u); in StgShur14Inflow()
379 E_internal = rho * cv * theta0; in StgShur14Inflow()
381 P = rho * Rd * theta0; // interior rho with exterior T in StgShur14Inflow()
398 v[0][i] -= wdetJb * rho * u_normal; in StgShur14Inflow()
[all …]
|
| H A D | stabilization.h | 76 const CeedScalar rho = s.U.density; in Tau_diagPrim() local
93 …tau = Square(rho) * (4. * Square(dts) + velocity_term) + Ctau_v * Square(mu) * DotN((CeedScalar *)… in Tau_diagPrim()
97 Tau_d[0] = Ctau_C * fact / (rho * (gijd_mat[0][0] + gijd_mat[1][1] + gijd_mat[2][2])) * 0.125; in Tau_diagPrim()
|
| H A D | advection.h | 507 const CeedScalar rho = q[0][i]; in Advection_InOutFlowGeneric() local
508 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Advection_InOutFlowGeneric()
|
| H A D | newtonian_state.h | 94 const CeedScalar rho = Y.pressure / (GasConstant(gas) * Y.temperature); in StateEntropyFromPrimitive() local
95 const CeedScalar entropy = log(Y.pressure) - gamma * log(rho); in StateEntropyFromPrimitive()
96 const CeedScalar rho_div_p = rho / Y.pressure; in StateEntropyFromPrimitive()
235 const CeedScalar rho = exp(log_rho); in StateConservativeFromEntropy_fwd() local
236 dU.density = -rho / (gamma - 1) * (dentropy + dV.S_energy / V.S_energy); in StateConservativeFromEntropy_fwd()
|
| /honee/examples/postprocess/ |
| H A D | vortexshedding.py | 8 def coeff(force, rho=1, u=1, D=1, zspan=0.2): argument
10 return 2 * force / (rho * u**2 * S)
|
| /honee/problems/ |
| H A D | bc_freestream.c | 221 const CeedScalar rho = 1.2; in TestRiemannHLL_fwd() local
222 const CeedScalar p = (HeatCapacityRatio(gas) - 1) * rho * gas.cv * T; in TestRiemannHLL_fwd()
283 const CeedScalar rho = 1.2; in TestRiemannHLLC_fwd() local
284 const CeedScalar p = (HeatCapacityRatio(gas) - 1) * rho * gas.cv * T; in TestRiemannHLLC_fwd()
345 const CeedScalar rho = 1.2; in TestComputeHLLSpeeds_Roe_fwd() local
346 const CeedScalar p = (HeatCapacityRatio(gas) - 1) * rho * gas.cv * T; in TestComputeHLLSpeeds_Roe_fwd()
420 const CeedScalar rho = 1.2; in TestTotalSpecificEnthalpy_fwd() local
421 const CeedScalar p = (HeatCapacityRatio(gas) - 1) * rho * gas.cv * T; in TestTotalSpecificEnthalpy_fwd()
|
| H A D | newtonian.c | 740 const CeedScalar rho = 1.2 * kg / Cube(m); in UnitTests_Newtonian() local
741 const CeedScalar P = (HeatCapacityRatio(gas) - 1) * rho * gas.cv * T; in UnitTests_Newtonian()
748 const CeedScalar entropy = log(P) - gamma * log(rho); in UnitTests_Newtonian()
749 const CeedScalar rho_div_p = rho / P; in UnitTests_Newtonian()
751 const CeedScalar U0[5] = {rho, rho * u[0], rho * u[1], rho * u[2], rho * e_total}; in UnitTests_Newtonian()
|
| H A D | channel.c | 98 CeedScalar rho = P0 / (Rd * theta0); in NS_CHANNEL() local
99 CeedScalar g[] = {channel_ctx->B / rho, 0., 0.}; in NS_CHANNEL()
|
| /honee/doc/ |
| H A D | theory.md | 13 \frac{\partial \rho}{\partial t} + \nabla \cdot \bm{U} &= 0 \\
14 …t} + \nabla \cdot \left( \frac{\bm{U} \otimes \bm{U}}{\rho} + P \bm{I}_3 -\bm\sigma \right) - \rho…
15 …t} + \nabla \cdot \left( \frac{(E + P)\bm{U}}{\rho} -\bm{u} \cdot \bm{\sigma} - k \nabla T \right)…
20 …rho$ represents the volume mass density, $U$ the momentum density (defined as $\bm{U}=\rho \bm{u}$…
23 P = \left( {c_p}/{c_v} -1\right) \left( E - {\bm{U}\cdot\bm{U}}/{(2 \rho)} \right) \, ,
39 \rho \\
40 \bm{U} \equiv \rho \bm{ u }\\
41 E \equiv \rho e
57 {(\bm{U} \otimes \bm{U})}/{\rho} + P \bm{I}_3 \\
58 {(E + P)\bm{U}}/{\rho}
[all …]
|
| H A D | examples.md | 12 \rho &= \rho_\infty\left(1+A\exp\left(\frac{-(\bar{x}^2 + \bar{y}^2)}{2\sigma^2}\right)\right) \\
122 \begin{aligned} \rho &= \frac{P_0}{( c_p - c_v)\theta(\bm{x},t)} \pi(\bm{x},t)^{\frac{c_v}{ c_p - c…
389 \rho &= \frac{p}{R T_0} \\
423 \frac{\partial \rho}{\partial t} + \nabla \cdot \bm{U} &= 0 \\
424 \frac{\partial \bm{U}}{\partial t} + \nabla \cdot \left( \frac{\bm{U} \otimes \bm{U}}{\rho} + P \bm…
425 \frac{\partial E}{\partial t} + \nabla \cdot \left( \frac{(E + P)\bm{U}}{\rho} \right) &= 0 \, , \\
429 …te}`zhang2011verification`, the mean flow for this problem is $\rho=1$, $P=1$, $T=P/\rho= 1$ (Spec…
436 There is no perturbation in the entropy $S=P/\rho^\gamma$ ($\delta S=0)$.
477 …roblem. The default initial conditions are $P=1$, $\rho=1$ for the driver section and $P=0.1$, $\r…
494 \tau_{SHOCK} = \frac{h_{SHOCK}}{2u_{cha}} \left( \frac{ \,|\, \nabla \rho \,|\, h_{SHOCK}}{\rho_{re…
[all …]
|
| H A D | runtime_options.md | 479 …rho, \rho \bm{u}, \rho e$), `primitive` ($P, \bm{u}, T$), or `entropy` ($\frac{\gamma - s}{\gamma …
|
| H A D | auxiliary.md | 119 | $\mean{\rho}$ | MeanDensity |
123 | $\mean{\rho T}$ | MeanDensityTemperature |
124 | $\mean{\rho T u_i}$ | MeanDensityTemperatureFlux[$i$] |
125 | $\mean{\rho u_i}$ | MeanMomentum[$i$] |
126 | $\mean{\rho u_i u_j}$ | MeanMomentumFlux[$ij$] |
129 where [$i$] are suffixes to the labels. So $\mean{\rho u_x u_y}$ would correspond to MeanMomentumFl…
|