Home
last modified time | relevance | path

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

/libCEED/examples/fluids/qfunctions/
H A Dnewtonian_state.h39 for (int i = 0; i < 3; i++) Y[i + 1] = s.velocity[i]; in UnpackState_Y()
60 CeedScalar e_kinetic = 0.5 * Dot3(s.Y.velocity, s.Y.velocity); in TotalSpecificEnthalpy()
66 CeedScalar de_kinetic = Dot3(ds.Y.velocity, s.Y.velocity); in TotalSpecificEnthalpy_fwd()
73 for (CeedInt i = 0; i < 3; i++) Y.velocity[i] = U.momentum[i] / U.density; in StatePrimitiveFromConservative()
74 CeedScalar e_kinetic = .5 * Dot3(Y.velocity, Y.velocity); in StatePrimitiveFromConservative()
85 dY.velocity[i] = (dU.momentum[i] - s.Y.velocity[i] * dU.density) / s.U.density; in StatePrimitiveFromConservative_fwd()
87 CeedScalar e_kinetic = .5 * Dot3(s.Y.velocity, s.Y.velocity); in StatePrimitiveFromConservative_fwd()
88 CeedScalar de_kinetic = Dot3(dY.velocity, s.Y.velocity); in StatePrimitiveFromConservative_fwd()
104 const CeedScalar e_kinetic = 0.5 * Dot3(Y.velocity, Y.velocity); in StateEntropyFromPrimitive()
107 for (int i = 0; i < 3; i++) V.S_momentum[i] = rho_div_p * Y.velocity[i]; in StateEntropyFromPrimitive()
[all …]
H A Dturb_spanstats.h31 v[TURB_MEAN_PRESSURE_VELOCITY_X][i] = wdetJ * s.Y.pressure * s.Y.velocity[0]; in ChildStatsCollection()
32 v[TURB_MEAN_PRESSURE_VELOCITY_Y][i] = wdetJ * s.Y.pressure * s.Y.velocity[1]; in ChildStatsCollection()
33 v[TURB_MEAN_PRESSURE_VELOCITY_Z][i] = wdetJ * s.Y.pressure * s.Y.velocity[2]; in ChildStatsCollection()
35 …TURB_MEAN_DENSITY_TEMPERATURE_FLUX_X][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[0]; in ChildStatsCollection()
36 …TURB_MEAN_DENSITY_TEMPERATURE_FLUX_Y][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[1]; in ChildStatsCollection()
37 …TURB_MEAN_DENSITY_TEMPERATURE_FLUX_Z][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[2]; in ChildStatsCollection()
41 v[TURB_MEAN_MOMENTUMFLUX_XX][i] = wdetJ * s.U.momentum[0] * s.Y.velocity[0]; in ChildStatsCollection()
42 v[TURB_MEAN_MOMENTUMFLUX_YY][i] = wdetJ * s.U.momentum[1] * s.Y.velocity[1]; in ChildStatsCollection()
43 v[TURB_MEAN_MOMENTUMFLUX_ZZ][i] = wdetJ * s.U.momentum[2] * s.Y.velocity[2]; in ChildStatsCollection()
44 v[TURB_MEAN_MOMENTUMFLUX_YZ][i] = wdetJ * s.U.momentum[1] * s.Y.velocity[2]; in ChildStatsCollection()
[all …]
H A Ddifferential_filter.h48 v0[DIFF_FILTER_VELOCITY_X][i] = wdetJ * s.Y.velocity[0]; in DifferentialFilter_RHS()
49 v0[DIFF_FILTER_VELOCITY_Y][i] = wdetJ * s.Y.velocity[1]; in DifferentialFilter_RHS()
50 v0[DIFF_FILTER_VELOCITY_Z][i] = wdetJ * s.Y.velocity[2]; in DifferentialFilter_RHS()
52 v1[DIFF_FILTER_VELOCITY_SQUARED_XX][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[0]; in DifferentialFilter_RHS()
53 v1[DIFF_FILTER_VELOCITY_SQUARED_YY][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[1]; in DifferentialFilter_RHS()
54 v1[DIFF_FILTER_VELOCITY_SQUARED_ZZ][i] = wdetJ * s.Y.velocity[2] * s.Y.velocity[2]; in DifferentialFilter_RHS()
55 v1[DIFF_FILTER_VELOCITY_SQUARED_YZ][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[2]; in DifferentialFilter_RHS()
56 v1[DIFF_FILTER_VELOCITY_SQUARED_XZ][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[2]; in DifferentialFilter_RHS()
57 v1[DIFF_FILTER_VELOCITY_SQUARED_XY][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[1]; in DifferentialFilter_RHS()
H A Driemann_solver.h167 CeedScalar u_left = Dot3(left.Y.velocity, normal); in RiemannFlux_HLL()
168 CeedScalar u_right = Dot3(right.Y.velocity, normal); in RiemannFlux_HLL()
202 CeedScalar u_left = Dot3(left.Y.velocity, normal); in RiemannFlux_HLL_fwd()
203 CeedScalar u_right = Dot3(right.Y.velocity, normal); in RiemannFlux_HLL_fwd()
204 CeedScalar du_left = Dot3(dleft.Y.velocity, normal); in RiemannFlux_HLL_fwd()
205 CeedScalar du_right = Dot3(dright.Y.velocity, normal); in RiemannFlux_HLL_fwd()
227 side.Y.velocity[0] + (s_star - u_side) * normal[0], in RiemannFlux_HLLC_Star()
228 side.Y.velocity[1] + (s_star - u_side) * normal[1], in RiemannFlux_HLLC_Star()
229 side.Y.velocity[2] + (s_star - u_side) * normal[2], in RiemannFlux_HLLC_Star()
250 side.Y.velocity[0] + (s_star - u_side) * normal[0], in RiemannFlux_HLLC_Star_fwd()
[all …]
H A Dblasius.h80 CeedScalar U_infty = sqrt(Dot3(S_infty.Y.velocity, S_infty.Y.velocity)); in BlasiusSolution()
115 const CeedScalar U_infty = sqrt(Dot3(S_infty.Y.velocity, S_infty.Y.velocity)); in ICsBlasius()
143 const CeedScalar U_infty = sqrt(Dot3(S_infty.Y.velocity, S_infty.Y.velocity)); in Blasius_Inflow()
200 const CeedScalar U_infty = sqrt(Dot3(S_infty.Y.velocity, S_infty.Y.velocity)); in Blasius_Inflow_Jacobian()
219 CeedScalar dE_kinetic = .5 * drho * Dot3(s.Y.velocity, s.Y.velocity); in Blasius_Inflow_Jacobian()
229 const CeedScalar u_normal = Dot3(norm, s.Y.velocity); in Blasius_Inflow_Jacobian()
233 v[j + 1][i] = -wdetJb * (drho * u_normal * s.Y.velocity[j] + norm[j] * dP); in Blasius_Inflow_Jacobian()
H A Dgaussianwave.h48 …onst CeedScalar e_kinetic = 0.5 * S_infty.U.density * Dot3(S_infty.Y.velocity, S_infty.Y.velocity); in IC_GaussianWave()
53 U[1] = S_infty.Y.velocity[0] * U[0]; in IC_GaussianWave()
54 U[2] = S_infty.Y.velocity[1] * U[0]; in IC_GaussianWave()
55 U[3] = S_infty.Y.velocity[2] * U[0]; in IC_GaussianWave()
H A Dbc_slip.h31 const CeedScalar vel_normal = Dot3(s.Y.velocity, norm); in Slip()
32 for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * norm[j] * vel_normal; in Slip()
84 const CeedScalar vel_normal = Dot3(s.Y.velocity, norm); in Slip_Jacobian()
85 for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * norm[j] * vel_normal; in Slip_Jacobian()
90 const CeedScalar dvel_normal = Dot3(ds.Y.velocity, norm); in Slip_Jacobian()
91 for (CeedInt j = 0; j < 3; j++) dvel_reflect[j] = ds.Y.velocity[j] - 2. * norm[j] * dvel_normal; in Slip_Jacobian()
H A Dtaylorgreen.h26 const CeedScalar V0 = sqrt(Dot3(reference.velocity, reference.velocity)); in ICsTaylorGreen()
H A Dchannel.h123 const CeedScalar E_kinetic = .5 * rho_in * Dot3(s_exact.Y.velocity, s_exact.Y.velocity); in Channel_Inflow()
130 const CeedScalar u_normal = Dot3(norm, s_exact.Y.velocity); in Channel_Inflow()
137 …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.h262 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in Tau()
269 MatVecNM(gijd_mat, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, gij_uj); in Tau()
270 …return 1 / sqrt(Square(2 * context->Ctau_t / context->dt) + DotN(s.Y.velocity, gij_uj, dim) * cont… in Tau()
315 div_u += grad_s[k].Y.velocity[j]; in IFunction_AdvectionGeneric()
318 CeedScalar strong_conv = s.U.E_total * div_u + DotN(s.Y.velocity, Grad_E, dim); in IFunction_AdvectionGeneric()
324 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in IFunction_AdvectionGeneric()
393 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in MassFunction_AdvectionGeneric()
451 div_u += grad_s[k].Y.velocity[j]; in RHSFunction_AdvectionGeneric()
454 CeedScalar strong_conv = s.U.E_total * div_u + DotN(s.Y.velocity, Grad_E, dim); in RHSFunction_AdvectionGeneric()
457 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in RHSFunction_AdvectionGeneric()
H A Dstabilization.h94 MatVec3(gijd_mat, s.Y.velocity, CEED_NOTRANSPOSE, gij_uj); in Tau_diagPrim()
95 velocity_term = Dot3(s.Y.velocity, gij_uj); in Tau_diagPrim()
H A Dnewtonian_types.h51 CeedScalar velocity[3]; member
H A Dbc_freestream.h198 const CeedScalar u_normal = Dot3(y_ext.velocity, normal); in RiemannOutflow()
201 …y_ext.velocity[j] += normal[j] * proj; // (I - n n^T) projects into the plane tangent to the norm… in RiemannOutflow()
274 const CeedScalar u_normal = Dot3(s_int.Y.velocity, normal); in RiemannOutflow_Jacobian()
275 const CeedScalar du_normal = Dot3(ds_int.Y.velocity, normal); in RiemannOutflow_Jacobian()
279 y_ext.velocity[j] += normal[j] * proj; in RiemannOutflow_Jacobian()
280 dy_ext.velocity[j] += normal[j] * dproj; in RiemannOutflow_Jacobian()
H A Dstg_shur14.h433 CeedScalar velocity[3]; in StgShur14Inflow_Jacobian() local
434 for (CeedInt j = 0; j < 3; j++) velocity[j] = jac_data_sur[5 + j][i]; in StgShur14Inflow_Jacobian()
443 CeedScalar dE_kinetic = .5 * drho * Dot3(velocity, velocity); in StgShur14Inflow_Jacobian()
453 const CeedScalar u_normal = Dot3(norm, velocity); in StgShur14Inflow_Jacobian()
456 …for (int j = 0; j < 3; j++) v[j + 1][i] = -wdetJb * (drho * u_normal * velocity[j] + norm[j] * dP); in StgShur14Inflow_Jacobian()
/libCEED/examples/fluids/problems/
H A Dbc_freestream.c38 …StatePrimitive Y_inf = {.pressure = reference->pressure / Pascal, .velocity = {0}, .temperature = … in FreestreamBCSetup()
39 for (int i = 0; i < 3; i++) Y_inf.velocity[i] = reference->velocity[i] * second / meter; in FreestreamBCSetup()
46 …("-freestream_velocity", "Velocity at freestream condition", NULL, Y_inf.velocity, &narray, NULL)); in FreestreamBCSetup()
102 for (int i = 0; i < 3; i++) Y_inf.velocity[i] *= meter / second; in FreestreamBCSetup()
398 CeedScalar u_left0 = Dot3(left0.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
399 CeedScalar u_right0 = Dot3(right0.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
414 CeedScalar du_left = Dot3(dleft0.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
415 CeedScalar du_right = Dot3(dright0.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
432 CeedScalar u_left1 = Dot3(left1.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
433 CeedScalar u_right1 = Dot3(right1.Y.velocity, normal); in TestComputeHLLSpeeds_Roe_fwd()
H A Dblasius.c26 CeedScalar U_infty = sqrt(Dot3(S_infty.Y.velocity, S_infty.Y.velocity)); in CompressibleBlasiusResidual()
326 .pressure = P_inf, .velocity = {U_inf, 0, 0}, in NS_BLASIUS()
H A Dnewtonian.c259 StatePrimitive reference = {.pressure = 1.01e5, .velocity = {0}, .temperature = 288.15}; in NS_NEWTONIAN_IG()
346 …rArray("-reference_velocity", "Reference/initial velocity", NULL, reference.velocity, &dim, NULL)); in NS_NEWTONIAN_IG()
405 for (PetscInt i = 0; i < 3; i++) reference.velocity[i] *= meter / second; in NS_NEWTONIAN_IG()
/libCEED/examples/nek/boxes/
H A Db1e.rea32 0.00000 p028 TORDER: mesh velocity (0: p28=p27)
64 0.00000 p060 !=0 --> init. velocity to small nonzero
98 0.00000 p094 start projecting velocity after p94 step
/libCEED/examples/fluids/
H A DREADME.md250 If in a region where the flow velocity is known (e.g., away from viscous walls), use `bc_freestream…
255 These conditions may be either weak or strong, with the latter specifying velocity and temperature …
256 The strong approach gives sharper resolution of velocity structures.
450 - Background velocity vector
691 - Freestream velocity vector
745 - Freestream velocity vector
853 - Maximum/centerline velocity of the flow
897 - Freestream velocity
1006 - Convective velocity, $U_0$
1011 - Only impose the mean velocity (no fluctutations)
H A Dindex.md28 … density (defined as $\bm{U}=\rho \bm{u}$, where $\bm{u}$ is the vector velocity field), $E$ the t…
224 A velocity vector $\bm u$ can be pulled back to the reference element as $\bm u_{\bm X} = \nabla_{\…
226 So a small normal component of velocity will be amplified (by a factor of the aspect ratio $1/\epsi…
227 …bm X} \rVert$ is a covariant measure of (half) the element length in the direction of the velocity.
305 where $u_i = \bm u \cdot \hat{\bm n}_i$ is the velocity component in direction $i$ and $a = \sqrt{\…
318 {ref}`problem-advection`, the problem of the transport of energy in a uniform vector velocity field…
528 with $\bm{u}$ the vector velocity field and $\kappa$ the diffusion coefficient.
533 In this case, a uniform circular velocity field transports the blob of total energy.
538 …In this case, a background wind with a constant rectilinear velocity field, enters the domain and …
580 …r section and $P=0.1$, $\rho=0.125$ for the driven section. The initial velocity is zero in both s…
[all …]