Home
last modified time | relevance | path

Searched refs:u (Results 1 – 18 of 18) sorted by relevance

/honee/qfunctions/
H A Deulervortex.h81 CeedScalar rho, P, T, E, u[3] = {0.}; in Exact_Euler() local
92 u[0] = mean_velocity[0] - C * y0; in Exact_Euler()
93 u[1] = mean_velocity[1] + C * x0; 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()
108 q[1] = rho * u[0]; in Exact_Euler()
109 q[2] = rho * u[1]; in Exact_Euler()
110 q[3] = rho * u[2]; in Exact_Euler()
[all …]
H A Dshocktube.h85 CeedScalar rho, P, u[3] = {0.}; in Exact_ShockTube() local
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 …cobian_Euler(CeedScalar dF[3][5][5], const CeedScalar rho, const CeedScalar u[3], const CeedScalar… in ConvectiveFluxJacobian_Euler()
111 CeedScalar u_sq = u[0] * u[0] + u[1] * u[1] + u[2] * u[2]; // Velocity square in ConvectiveFluxJacobian_Euler()
114 dF[i][j + 1][0] = ((i == j) ? ((gamma - 1.) * (u_sq / 2.)) : 0.) - u[i] * u[j]; in ConvectiveFluxJacobian_Euler()
117 …dF[i][j + 1][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (… in ConvectiveFluxJacobian_Euler()
118 … = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.) * u[i] * u[k]; in ConvectiveFluxJacobian_Euler()
[all …]
H A Dutils.h27 CEED_QFUNCTION_HELPER void ScaleN(CeedScalar *u, const CeedScalar alpha, const CeedInt N) { in ScaleN() argument
28 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) u[i] *= alpha; in ScaleN()
32 CEED_QFUNCTION_HELPER void SetValueN(CeedScalar *u, const CeedScalar alpha, const CeedInt N) { in SetValueN() argument
33 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) u[i] = alpha; in SetValueN()
43 CEED_QFUNCTION_HELPER CeedScalar DotN(const CeedScalar *u, const CeedScalar *v, const CeedInt N) { in DotN() argument
45 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) output += u[i] * v[i]; in DotN()
55 …ION_HELPER CeedScalar Dot3(const CeedScalar *u, const CeedScalar *v) { return u[0] * v[0] + u[1] *… in Dot3() argument
58 …_QFUNCTION_HELPER CeedScalar Dot2(const CeedScalar *u, const CeedScalar *v) { return u[0] * v[0] +… in Dot2() argument
61 …UNCTION_HELPER CeedScalar Norm3(const CeedScalar *u) { return sqrt(u[0] * u[0] + u[1] * u[1] + u[2… in Norm3() argument
64 CEED_QFUNCTION_HELPER CeedScalar Norm2(const CeedScalar *u) { return sqrt(u[0] * u[0] + u[1] * u[1]… in Norm2() argument
[all …]
H A Dstg_shur14.h153 … const CeedScalar qn[], CeedScalar u[3], const StgShur14Context stg_ctx) { in StgShur14Calc()
173 u[0] = ubar[0] + cij[0] * vp[0]; in StgShur14Calc()
174 u[1] = ubar[1] + cij[3] * vp[0] + cij[1] * vp[1]; in StgShur14Calc()
175 u[2] = ubar[2] + cij[4] * vp[0] + cij[5] * vp[1] + cij[2] * vp[2]; in StgShur14Calc()
195 … const CeedScalar eps, const CeedScalar lt, const CeedScalar nu, CeedScalar u[3], in StgShur14Calc_PrecompEktot()
220 u[0] = ubar[0] + cij[0] * vp[0]; in StgShur14Calc_PrecompEktot()
221 u[1] = ubar[1] + cij[3] * vp[0] + cij[1] * vp[1]; in StgShur14Calc_PrecompEktot()
222 u[2] = ubar[2] + cij[4] * vp[0] + cij[5] * vp[1] + cij[2] * vp[2]; in StgShur14Calc_PrecompEktot()
291 CeedScalar qn[STG_NMODES_MAX], u[3], ubar[3], cij[6], eps, lt; in ICsStg() local
311 StgShur14Calc(x_i, time, ubar, cij, qn, u, stg_ctx); in ICsStg()
[all …]
H A Dtaylorgreen.h19 CeedScalar u[3]; member
41 Y[1] = context->u[0] + V0 * sin(x[0]) * cos(x[1]) * cos(x[2]); in ICsTaylorGreen()
42 Y[2] = context->u[1] - V0 * cos(x[0]) * sin(x[1]) * cos(x[2]); in ICsTaylorGreen()
43 Y[3] = context->u[2]; in ICsTaylorGreen()
H A Dmass.h20 const CeedScalar(*u)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in Mass_N() local
24 …nt i = 0; i < Q; i++) CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) v[j][i] = q_data[i] * u[j][i]; in Mass_N()
H A Dchannel.h154 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Channel_Outflow() local
163 const CeedScalar u_normal = Dot3(norm, u); // Normal velocity 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 Dadvection.h508 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Advection_InOutFlowGeneric() local
515 const CeedScalar u_normal = DotN(normal, u, dim); in Advection_InOutFlowGeneric()
H A Dnewtonian_state.h50 …R CeedScalar Mach(NewtonianIGProperties gas, CeedScalar T, CeedScalar u) { return u / SoundSpeed(g… in Mach() argument
/honee/problems/
H A Dbc_freestream.c224 const CeedScalar u[3] = {u_base, u_base * 1.1, u_base * 1.2}; in TestRiemannHLL_fwd() local
225 const CeedScalar Y0_left[5] = {p, u[0], u[1], u[2], T}; in TestRiemannHLL_fwd()
226 const CeedScalar Y0_right[5] = {1.2 * p, 1.2 * u[0], 1.2 * u[1], 1.2 * u[2], 1.2 * T}; in TestRiemannHLL_fwd()
286 const CeedScalar u[3] = {u_base, u_base * 1.1, u_base * 1.2}; in TestRiemannHLLC_fwd() local
287 const CeedScalar Y0_left[5] = {p, u[0], u[1], u[2], T}; in TestRiemannHLLC_fwd()
288 const CeedScalar Y0_right[5] = {1.2 * p, 1.2 * u[0], 1.2 * u[1], 1.2 * u[2], 1.2 * T}; in TestRiemannHLLC_fwd()
348 const CeedScalar u[3] = {u_base, u_base * 1.1, u_base * 1.2}; in TestComputeHLLSpeeds_Roe_fwd() local
349 const CeedScalar Y0_left[5] = {p, u[0], u[1], u[2], T}; in TestComputeHLLSpeeds_Roe_fwd()
350 const CeedScalar Y0_right[5] = {1.2 * p, 1.2 * u[0], 1.2 * u[1], 1.2 * u[2], 1.2 * T}; in TestComputeHLLSpeeds_Roe_fwd()
423 const CeedScalar u[3] = {u_base, u_base * 1.1, u_base * 1.2}; in TestTotalSpecificEnthalpy_fwd() local
[all …]
H A Dtaylorgreen.c34 PetscInt u_size = PETSC_STATIC_ARRAY_LENGTH(taylorgreen_ctx->u); in NS_TAYLOR_GREEN()
36 taylorgreen_ctx->u, &u_size, NULL)); in NS_TAYLOR_GREEN()
H A Dnewtonian.c743 const CeedScalar u[3] = {u_base, u_base * 1.1, u_base * 1.2}; in UnitTests_Newtonian() local
744 const CeedScalar e_kinetic = 0.5 * Dot3(u, u); in UnitTests_Newtonian()
750 const CeedScalar Y0[5] = {P, u[0], u[1], u[2], T}; in UnitTests_Newtonian()
751 const CeedScalar U0[5] = {rho, rho * u[0], rho * u[1], rho * u[2], rho * e_total}; in UnitTests_Newtonian()
752 …opy) / (gamma - 1) - rho_div_p * (e_kinetic), rho_div_p * u[0], rho_div_p * u[1], rho_div_p * u[2], in UnitTests_Newtonian()
/honee/doc/
H A Dtheory.md15 …dot \left( \frac{(E + P)\bm{U}}{\rho} -\bm{u} \cdot \bm{\sigma} - k \nabla T \right) - \rho \bm{b}…
19 where $\bm{\sigma} = \mu(\nabla \bm{u} + (\nabla \bm{u})^T + \lambda (\nabla \cdot \bm{u})\bm{I}_3)…
20 …mass density, $U$ the momentum density (defined as $\bm{U}=\rho \bm{u}$, where $\bm{u}$ is the vec…
40 \bm{U} \equiv \rho \bm{ u }\\
63 - \bm{u} \cdot \bm{\sigma} - k \nabla T
69 \rho \bm{b}\cdot \bm{u}
111 <!-- TODO: This should be reframed in terms of PETSc TS's F(t, u, \dot u) = G(t, u) rather than spe…
228 A velocity vector $\bm u$ can be pulled back to the reference element as $\bm u_{\bm X} = \nabla_{\…
231 The ratio $\lVert \bm u \rVert / \lVert \bm u_{\bm X} \rVert$ is a covariant measure of (half) the …
236 The cell Péclet number is classically defined by $\mathrm{Pe}_h = \lVert \bm u \rVert h / (2 \kappa…
[all …]
H A Dexamples.md122 …v}{ c_p - c_v}} \, , \\ e &= c_v \theta(\bm{x},t) \pi(\bm{x},t) + \bm{u}\cdot \bm{u} /2 + g z \, ,…
126 For this problem, we have used no-slip and non-penetration boundary conditions for $\bm{u}$, and no…
385 u &= V'_x + V \sin(\hat x) \cos(\hat y) \sin(\hat z) \\
429 …= 1$ (Specific Gas Constant, $R$, is 1), and $\bm{u}=(u_1,u_2,0)$ while the perturbation $\delta \…
545 …ection-diffusion, which holds density $\rho$ and momentum density $\rho \bm u$ constant while adve…
549 \frac{\partial E}{\partial t} + \nabla \cdot (\bm{u} E ) - \kappa \nabla E = 0 \, ,
552 with $\bm{u}$ the vector velocity field and $\kappa$ the diffusion coefficient.
555 There are three different definitions for $\bm{u}$:
560 …We have solved {eq}`eq-advection` applying zero energy density $E$, and no-flux for $\bm{u}$ on th…
604 …hat{\bm{n}} \,dS = \int_{\partial \Omega_{inflow}} \bm v \, E_{wind} \, \bm u \cdot \widehat{\bm{n…
[all …]
H A Druntime_options.md479 …o \bm{u}, \rho e$), `primitive` ($P, \bm{u}, T$), or `entropy` ($\frac{\gamma - s}{\gamma - 1} - \…
H A Dauxiliary.md38 The projection of a function $u$ onto the parent FEM space would look like:
41 \bm M u_N = \int_0^{L_x} \int_0^{L_y} u \psi^\mathrm{parent}_N \mathrm{d}y \mathrm{d}x
44 Substituting the spanwise average of $\phi$ for $u$, we get:
/honee/examples/postprocess/
H A Dvortexshedding.py8 def coeff(force, rho=1, u=1, D=1, zspan=0.2): argument
10 return 2 * force / (rho * u**2 * S)
/honee/src/
H A Dmisc.c229 Vec u, v; in HoneeCalculateDomainSize() local
244 PetscCall(CeedOperatorCreateLocalVecs(op_mass, DMReturnVecType(dm), PETSC_COMM_SELF, &u, NULL)); in HoneeCalculateDomainSize()
246 PetscCall(VecSet(u, 1.)); in HoneeCalculateDomainSize()
247 PetscCall(ApplyCeedOperatorLocalToGlobal(u, v, op_mass_ctx)); in HoneeCalculateDomainSize()
251 PetscCall(VecDestroy(&u)); in HoneeCalculateDomainSize()