Home
last modified time | relevance | path

Searched refs:velocity (Results 1 – 25 of 28) sorted by relevance

12

/honee/qfunctions/
H A Dnumerics.h23 CEED_QFUNCTION_HELPER CeedScalar CalculateCFL_3D(const CeedScalar velocity[3], CeedScalar timestep,… in CalculateCFL_3D()
26 MatVec3(gij, velocity, CEED_NOTRANSPOSE, gij_uj); in CalculateCFL_3D()
27 return sqrt(Dot3(velocity, gij_uj)) * timestep; in CalculateCFL_3D()
46 CEED_QFUNCTION_HELPER CeedScalar CalculateCFL_2D(const CeedScalar velocity[2], CeedScalar timestep,… in CalculateCFL_2D()
49 MatVec2(gij, velocity, CEED_NOTRANSPOSE, gij_uj); in CalculateCFL_2D()
50 return sqrt(Dot2(velocity, gij_uj)) * timestep; in CalculateCFL_2D()
69 CEED_QFUNCTION_HELPER CeedScalar CalculatePe_3D(const CeedScalar velocity[3], CeedScalar diffusion_… in CalculatePe_3D()
73 MatVec3(gij_inv_mat, velocity, CEED_NOTRANSPOSE, gij_uj); in CalculatePe_3D()
74 return sqrt(Dot3(velocity, gij_uj)) / diffusion_coeff; in CalculatePe_3D()
93 CEED_QFUNCTION_HELPER CeedScalar CalculatePe_2D(const CeedScalar velocity[2], CeedScalar diffusion_… in CalculatePe_2D()
[all …]
H A Dnewtonian_state.h32 for (int i = 0; i < 3; i++) Y[i + 1] = s.velocity[i]; in UnpackState_Y()
53 CeedScalar e_kinetic = 0.5 * Dot3(s.Y.velocity, s.Y.velocity); in TotalSpecificEnthalpy()
59 CeedScalar de_kinetic = Dot3(ds.Y.velocity, s.Y.velocity); in TotalSpecificEnthalpy_fwd()
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()
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()
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()
[all …]
H A Dsgs_dd_training.h52 …v[0 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XX][i] - Square(s.Y.velocity[0])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal()
53 …v[1 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_YY][i] - Square(s.Y.velocity[1])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal()
54 …v[2 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_ZZ][i] - Square(s.Y.velocity[2])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal()
55 …v[3 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_YZ][i] - s.Y.velocity[1] * s.Y.velocity[2]) … in ComputeSGS_DDAnisotropicTrainingDataNodal()
56 …v[4 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XZ][i] - s.Y.velocity[0] * s.Y.velocity[2]) … in ComputeSGS_DDAnisotropicTrainingDataNodal()
57 …v[5 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XY][i] - s.Y.velocity[0] * s.Y.velocity[1]) … in ComputeSGS_DDAnisotropicTrainingDataNodal()
H A Ddifferential_filter.h44 v0[DIFF_FILTER_VELOCITY_X][i] = wdetJ * s.Y.velocity[0]; in DifferentialFilter_RHS()
45 v0[DIFF_FILTER_VELOCITY_Y][i] = wdetJ * s.Y.velocity[1]; in DifferentialFilter_RHS()
46 v0[DIFF_FILTER_VELOCITY_Z][i] = wdetJ * s.Y.velocity[2]; in DifferentialFilter_RHS()
48 v1[DIFF_FILTER_VELOCITY_SQUARED_XX][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[0]; in DifferentialFilter_RHS()
49 v1[DIFF_FILTER_VELOCITY_SQUARED_YY][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[1]; in DifferentialFilter_RHS()
50 v1[DIFF_FILTER_VELOCITY_SQUARED_ZZ][i] = wdetJ * s.Y.velocity[2] * s.Y.velocity[2]; in DifferentialFilter_RHS()
51 v1[DIFF_FILTER_VELOCITY_SQUARED_YZ][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[2]; in DifferentialFilter_RHS()
52 v1[DIFF_FILTER_VELOCITY_SQUARED_XZ][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[2]; in DifferentialFilter_RHS()
53 v1[DIFF_FILTER_VELOCITY_SQUARED_XY][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[1]; in DifferentialFilter_RHS()
H A Driemann_solver.h163 CeedScalar u_left = Dot3(left.Y.velocity, normal); in RiemannFlux_HLL()
164 CeedScalar u_right = Dot3(right.Y.velocity, normal); in RiemannFlux_HLL()
198 CeedScalar u_left = Dot3(left.Y.velocity, normal); in RiemannFlux_HLL_fwd()
199 CeedScalar u_right = Dot3(right.Y.velocity, normal); in RiemannFlux_HLL_fwd()
200 CeedScalar du_left = Dot3(dleft.Y.velocity, normal); in RiemannFlux_HLL_fwd()
201 CeedScalar du_right = Dot3(dright.Y.velocity, normal); in RiemannFlux_HLL_fwd()
223 side.Y.velocity[0] + (s_star - u_side) * normal[0], in RiemannFlux_HLLC_Star()
224 side.Y.velocity[1] + (s_star - u_side) * normal[1], in RiemannFlux_HLLC_Star()
225 side.Y.velocity[2] + (s_star - u_side) * normal[2], in RiemannFlux_HLLC_Star()
246 side.Y.velocity[0] + (s_star - u_side) * normal[0], in RiemannFlux_HLLC_Star_fwd()
[all …]
H A Dgaussianwave.h42 …onst CeedScalar e_kinetic = 0.5 * S_infty.U.density * Dot3(S_infty.Y.velocity, S_infty.Y.velocity); in IC_GaussianWave()
47 U[1] = S_infty.Y.velocity[0] * U[0]; in IC_GaussianWave()
48 U[2] = S_infty.Y.velocity[1] * U[0]; in IC_GaussianWave()
49 U[3] = S_infty.Y.velocity[2] * U[0]; in IC_GaussianWave()
H A Dbc_slip.h29 const CeedScalar vel_normal = Dot3(s.Y.velocity, normal); in Slip()
30 for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * normal[j] * vel_normal; in Slip()
79 const CeedScalar vel_normal = Dot3(s.Y.velocity, normal); in Slip_Jacobian()
80 for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * normal[j] * vel_normal; in Slip_Jacobian()
85 const CeedScalar dvel_normal = Dot3(ds.Y.velocity, normal); in Slip_Jacobian()
86 … for (CeedInt j = 0; j < 3; j++) dvel_reflect[j] = ds.Y.velocity[j] - 2. * normal[j] * dvel_normal; in Slip_Jacobian()
H A Dblasius.h73 CeedScalar U_infty = Norm3(S_infty.Y.velocity); in BlasiusSolution()
108 const CeedScalar U_infty = Norm3(S_infty.Y.velocity); in ICsBlasius()
135 const CeedScalar U_infty = Norm3(S_infty.Y.velocity); in Blasius_Inflow()
190 const CeedScalar U_infty = Norm3(S_infty.Y.velocity); in Blasius_Inflow_Jacobian()
209 CeedScalar dE_kinetic = .5 * drho * Dot3(s.Y.velocity, s.Y.velocity); in Blasius_Inflow_Jacobian()
219 const CeedScalar u_normal = Dot3(normal, s.Y.velocity); in Blasius_Inflow_Jacobian()
223 v[j + 1][i] = -wdetJb * (drho * u_normal * s.Y.velocity[j] + normal[j] * dP); in Blasius_Inflow_Jacobian()
H A Dchannel.h115 const CeedScalar E_kinetic = .5 * rho_in * Dot3(s_exact.Y.velocity, s_exact.Y.velocity); in Channel_Inflow()
122 const CeedScalar u_normal = Dot3(norm, s_exact.Y.velocity); in Channel_Inflow()
129 …for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho_in * u_normal * s_exact.Y.velocity[j]… in Channel_Inflow()
H A Dadvection.h245 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in Tau()
254 MatVecNM(gijd_mat, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, gij_uj); in Tau()
255 …return 1 / sqrt(Square(2 * context->Ctau_t / context->dt) + DotN(s.Y.velocity, gij_uj, dim) * Squa… in Tau()
299 div_u += grad_s[k].Y.velocity[j]; in IFunction_AdvectionGeneric()
303 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in IFunction_AdvectionGeneric()
304 CeedScalar strong_conv = s.U.E_total * div_u + DotN(s.Y.velocity, Grad_E, dim); in IFunction_AdvectionGeneric()
374 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in MassFunction_AdvectionGeneric()
434 div_u += grad_s[k].Y.velocity[j]; in RHSFunction_AdvectionGeneric()
437 CeedScalar strong_conv = s.U.E_total * div_u + DotN(s.Y.velocity, Grad_E, dim); in RHSFunction_AdvectionGeneric()
440 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in RHSFunction_AdvectionGeneric()
H A Dbc_outflow.h58 const CeedScalar u_normal = Dot3(y_ext.velocity, normal); in RiemannOutflow()
61 …y_ext.velocity[j] += normal[j] * proj; // (I - n n^T) projects into the plane tangent to the norm… in RiemannOutflow()
134 const CeedScalar u_normal = Dot3(s_int.Y.velocity, normal); in RiemannOutflow_Jacobian()
135 const CeedScalar du_normal = Dot3(ds_int.Y.velocity, normal); in RiemannOutflow_Jacobian()
139 y_ext.velocity[j] += normal[j] * proj; in RiemannOutflow_Jacobian()
140 dy_ext.velocity[j] += normal[j] * dproj; in RiemannOutflow_Jacobian()
H A Dmonitor_totalkineticenergy.h25 v[0][i] = wdetJ * 0.5 * s.U.density * Dot3(s.Y.velocity, s.Y.velocity); in MonitorTotalKineticEnergy()
H A Dmonitor_cfl.h30 v[i] = CalculateCFL_2D(s.Y.velocity, 1, gijd_mat); in MonitorCFL()
40 v[i] = CalculateCFL_3D(s.Y.velocity, 1, gijd_mat); in MonitorCFL()
H A Dstabilization.h89 MatVec3(gijd_mat, s.Y.velocity, CEED_NOTRANSPOSE, gij_uj); in Tau_diagPrim()
90 velocity_term = Dot3(s.Y.velocity, gij_uj); in Tau_diagPrim()
H A Dstg_shur14.h438 CeedScalar velocity[3]; in StgShur14Inflow_Jacobian() local
439 for (CeedInt j = 0; j < 3; j++) velocity[j] = jac_data_sur[5 + j][i]; in StgShur14Inflow_Jacobian()
448 CeedScalar dE_kinetic = .5 * drho * Dot3(velocity, velocity); in StgShur14Inflow_Jacobian()
458 const CeedScalar u_normal = Dot3(normal, velocity); in StgShur14Inflow_Jacobian()
461 …for (int j = 0; j < 3; j++) v[j + 1][i] = -wdetJb * (drho * u_normal * velocity[j] + normal[j] * d… in StgShur14Inflow_Jacobian()
H A Dnewtonian_types.h51 CeedScalar velocity[3]; member
H A Dsgs_dd_model.h210 Flux[4][j] = Y.velocity[0] * sgs[0][j] + Y.velocity[1] * sgs[1][j] + Y.velocity[2] * sgs[2][j]; in FluxSubgridStress()
H A Dtaylorgreen.h32 const CeedScalar V0 = Norm3(reference.velocity); in ICsTaylorGreen()
/honee/qfunctions/spanstats/
H A Dturbulence.h59 v[TURB_MEAN_PRESSURE_VELOCITY_X][i] = wdetJ * s.Y.pressure * s.Y.velocity[0]; in ChildStatsCollection()
60 v[TURB_MEAN_PRESSURE_VELOCITY_Y][i] = wdetJ * s.Y.pressure * s.Y.velocity[1]; in ChildStatsCollection()
61 v[TURB_MEAN_PRESSURE_VELOCITY_Z][i] = wdetJ * s.Y.pressure * s.Y.velocity[2]; in ChildStatsCollection()
63 …TURB_MEAN_DENSITY_TEMPERATURE_FLUX_X][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[0]; in ChildStatsCollection()
64 …TURB_MEAN_DENSITY_TEMPERATURE_FLUX_Y][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[1]; in ChildStatsCollection()
65 …TURB_MEAN_DENSITY_TEMPERATURE_FLUX_Z][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[2]; in ChildStatsCollection()
69 v[TURB_MEAN_MOMENTUMFLUX_XX][i] = wdetJ * s.U.momentum[0] * s.Y.velocity[0]; in ChildStatsCollection()
70 v[TURB_MEAN_MOMENTUMFLUX_YY][i] = wdetJ * s.U.momentum[1] * s.Y.velocity[1]; in ChildStatsCollection()
71 v[TURB_MEAN_MOMENTUMFLUX_ZZ][i] = wdetJ * s.U.momentum[2] * s.Y.velocity[2]; in ChildStatsCollection()
72 v[TURB_MEAN_MOMENTUMFLUX_YZ][i] = wdetJ * s.U.momentum[1] * s.Y.velocity[2]; in ChildStatsCollection()
[all …]
H A Dcflpe.h44 cfl = CalculateCFL_2D(s.Y.velocity, context->timestep, gijd_mat); in ChildStatsCollection_CflPe()
45 Pe = CalculatePe_2D(s.Y.velocity, context->diffusion_coeff, gijd_mat); in ChildStatsCollection_CflPe()
57 cfl = CalculateCFL_3D(s.Y.velocity, context->timestep, gijd_mat); in ChildStatsCollection_CflPe()
58 Pe = CalculatePe_3D(s.Y.velocity, context->diffusion_coeff, gijd_mat); in ChildStatsCollection_CflPe()
/honee/problems/
H A Dbc_freestream.c125 …StatePrimitive Y_inf = {.pressure = reference->pressure / units->Pascal, .velocity = {0}, .tempera… in FreestreamBCSetup()
126 …for (int i = 0; i < 3; i++) Y_inf.velocity[i] = reference->velocity[i] * units->second / units->me… in FreestreamBCSetup()
133 …("-freestream_velocity", "Velocity at freestream condition", NULL, Y_inf.velocity, &narray, NULL)); in FreestreamBCSetup()
138 for (int i = 0; i < 3; i++) Y_inf.velocity[i] *= units->meter / units->second; in FreestreamBCSetup()
357 CeedScalar u_left0 = Dot3(left0.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
358 CeedScalar u_right0 = Dot3(right0.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
373 CeedScalar du_left = Dot3(dleft0.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
374 CeedScalar du_right = Dot3(dright0.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
391 CeedScalar u_left1 = Dot3(left1.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
392 CeedScalar u_right1 = Dot3(right1.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
H A Dblasius.c23 CeedScalar U_infty = Norm3(S_infty.Y.velocity); in CompressibleBlasiusResidual()
187 .pressure = P_inf, .velocity = {U_inf, 0, 0}, in NS_BLASIUS()
/honee/doc/
H A Dexamples.md19 … a strong acoustic wave and leaves behind a cold thermal bubble that advects at the fluid velocity.
22 …r constant-pressure wave that transports temperature and transverse momentum at the fluid velocity.
65 The domain is filled with an ideal gas at rest (zero velocity) with temperature 24.92 and pressure …
67 …and Riemann-type outflow on the right, with exterior reference state at velocity $(1, 0, 0)$ givin…
68 …s imposed at the top and bottom boundaries $(y = \pm 4.5)$ (zero normal velocity component, zero h…
85 where $\rho_\infty, u_\infty$ are the freestream (inflow) density and velocity respectively.
189 where $H$ is the channel half-height, $u_{\max}$ is the center velocity, $T_w$ is the temperature a…
205 - Maximum/centerline velocity of the flow
300 At the inflow, the velocity is prescribed by the Blasius soution profile, density is set constant, …
316 - Freestream velocity
[all …]
H A Dtheory.md20 … density (defined as $\bm{U}=\rho \bm{u}$, where $\bm{u}$ is the vector velocity field), $E$ the t…
228 A velocity vector $\bm u$ can be pulled back to the reference element as $\bm u_{\bm X} = \nabla_{\…
230 So a small normal component of velocity will be amplified (by a factor of the aspect ratio $1/\epsi…
231 …bm X} \rVert$ is a covariant measure of (half) the element length in the direction of the velocity.
324 where $u_i = \bm u \cdot \hat{\bm n}_i$ is the velocity component in direction $i$ and $a = \sqrt{\…
337 {ref}`problem-advection`, the problem of the transport of energy in a uniform vector velocity field…
432 The SGS tensor is calculated at nodes using an $L^2$ projection of the velocity gradient and grid a…
525 …Reynolds stress tensor $\bm{C}$ (such that $\bm{R} = \bm{CC}^T$ ), bulk velocity $U_0$, wavemode a…
560 …undary condition is identical to the blasius inflow; it weakly enforces velocity, with the option …
679 - Convective velocity, $U_0$
[all …]
H A Druntime_options.md346 If in a region where the flow velocity is known (e.g., away from viscous walls), use `bc_freestream…
351 These conditions may be either weak or strong, with the latter specifying velocity and temperature …
352 The strong approach gives sharper resolution of velocity structures.

12