Home
last modified time | relevance | path

Searched full:velocity (Results 1 – 25 of 33) sorted by relevance

12

/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 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 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 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 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 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 Dchannel.h25 CeedScalar umax; // !< Centerline velocity
120 // Find inflow state using calculated P and prescribed velocity, theta0 in Channel_Inflow()
123 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()
171 const CeedScalar u_normal = Dot3(norm, u); // Normal velocity in Channel_Outflow()
H A Dstg_shur14.h35 * @param[out] ubar Mean velocity at wall_dist
151 * @param[in] ubar Mean velocity at X
154 * @param[out] u Velocity at X and t
188 * @param[in] ubar Mean velocity at X
195 * @param[out] u Velocity at X and t
327 * at each location, then calculate the actual velocity.
385 // Velocity normal to the boundary in StgShur14Inflow()
433 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()
435 // TODO This is almost certainly a bug. Velocity isn't stored here, only 0s. in StgShur14Inflow_Jacobian()
[all …]
H A Dstg_shur14_type.h24 CeedScalar u0; // !< Convective velocity
40 size_t ubar; // !< Mean velocity, [nprof, 3]
H A Dtaylorgreen.h26 const CeedScalar V0 = sqrt(Dot3(reference.velocity, reference.velocity)); in ICsTaylorGreen()
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 Dbc_freestream.h173 // The velocity in the exterior state has optional softplus regularization to
198 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 Deulervortex.h43 // Velocity:
46 // Velocity/Momentum Density:
108 case 1: // Constant zero velocity, density constant, total energy constant in Exact_Euler()
119 case 2: // Constant nonzero velocity, density constant, total energy constant in Exact_Euler()
132 …case 3: // Velocity zero, pressure constant (so density and internal energy will be non-constant)… in Exact_Euler()
146 …case 4: // Constant nonzero velocity, pressure constant (so density and internal energy will be n… in Exact_Euler()
147 // It should be transported across the domain, but velocity stays constant in Exact_Euler()
184 CeedScalar u_sq = u[0] * u[0] + u[1] * u[1] + u[2] * u[2]; // Velocity square in ConvectiveFluxJacobian_Euler()
541 // For test cases 1 and 3 the background velocity is zero in TravelingVortex_Inflow()
621 … const CeedScalar u_normal = norm[0] * u[0] + norm[1] * u[1] + norm[2] * u[2]; // Normal velocity in Euler_Outflow()
H A Dbc_freestream_type.h22 …CeedScalar softplus_velocity; // "width" of the softplus velocity regulariz…
H A Dnewtonian_types.h51 CeedScalar velocity[3]; member
/libCEED/examples/fluids/
H A Dvortexshedding.yaml19 # Reference state is used for the initial condition, zero velocity by default.
22 # of 100. With velocity 1, this flow is Mach 0.01.
34 velocity: 1,0,0
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 …]
H A Dgaussianwave.yaml16 velocity: 2,2,0
/libCEED/examples/fluids/src/
H A Dvelocity_gradient_projection.c8 /// Functions for setting up and projecting the velocity gradient
20 grad_velo_proj->num_comp = 9; // 9 velocity gradient in VelocityGradientProjectionCreateDM()
23 PetscCall(PetscObjectSetName((PetscObject)grad_velo_proj->dm, "Velocity Gradient Projection")); in VelocityGradientProjectionCreateDM()
88 …PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_rhs_assemble, "velocity gradient", grad_velo_proj->n… in VelocityGradientProjectionSetup()
95 …PetscCallCeed(ceed, CeedOperatorSetField(op_rhs_assemble, "velocity gradient", elem_restr_grad_vel… in VelocityGradientProjectionSetup()
/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 …PetscOptionsScalarArray("-freestream_velocity", "Velocity at freestream condition", NULL, Y_inf.ve… in FreestreamBCSetup()
102 for (int i = 0; i < 3; i++) Y_inf.velocity[i] *= meter / second; in FreestreamBCSetup()
152 …Real("-outflow_recirc", "Fraction of recirculation to allow in exterior velocity state [0,1]", NUL… in OutflowBCSetup()
154 …PetscCall(PetscOptionsReal("-outflow_softplus_velocity", "Characteristic velocity of softplus regu… in OutflowBCSetup()
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()
[all …]
H A Dblasius.c26 CeedScalar U_infty = sqrt(Dot3(S_infty.Y.velocity, S_infty.Y.velocity)); in CompressibleBlasiusResidual()
276 …PetscCall(PetscOptionsScalar("-velocity_infinity", "Velocity at boundary layer edge", NULL, U_inf,… in NS_BLASIUS()
289 …PetscCall(PetscOptionsBoundedInt("-platemesh_Ndelta", "Velocity at boundary layer edge", NULL, mes… in NS_BLASIUS()
326 .pressure = P_inf, .velocity = {U_inf, 0, 0}, in NS_BLASIUS()
H A Deulervortex.c74 …PetscCall(PetscOptionsRealArray("-mean_velocity", "Background velocity vector", NULL, mean_velocit… in NS_EULER_VORTEX()
94 …PetscCall(PetscPrintf(comm, "Warning! Background velocity vector for -euler_test t1 and -euler_tes… in NS_EULER_VORTEX()
156 " Background Velocity : %f,%f,%f\n" in PRINT_EULER_VORTEX()
/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
186 T VELOCITY

12