Home
last modified time | relevance | path

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

12

/honee/qfunctions/
H A Dnumerics.h15 …`, \f$ \overline{g}_{jk} \f$ is the element metric tensor `gij`, and \f$ u_j \f$ is the `velocity`.
18 @param[in] velocity Velocity
23 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()
38 …`, \f$ \overline{g}_{jk} \f$ is the element metric tensor `gij`, and \f$ u_j \f$ is the `velocity`.
41 @param[in] velocity Advection velocity
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()
[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 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 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 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 Dsgs_dd_utils.h17 // @brief Calculate Frobenius norm of velocity gradient from eigenframe quantities
50 * @param[in] grad_velo_aniso Gradient of velocity in physical (anisotropic) coordinates
54 * @param[out] eigenvectors Eigenvectors of the (anisotropic) velocity gradient
56 * @param[out] grad_velo_magnitude Frobenius norm of the velocity gradient
64 // -- Transform physical, anisotropic velocity gradient to isotropic in ComputeSgsDDInputs()
100 * @param[in] eigenvectors Eigenvectors of the (anisotropic) velocity gradient
102 * @param[in] grad_velo_magnitude Magnitude of the velocity gradient
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 Dchannel.h17 CeedScalar umax; // !< Centerline velocity
112 // Find inflow state using calculated P and prescribed velocity, theta0 in Channel_Inflow()
115 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()
163 const CeedScalar u_normal = Dot3(norm, u); // Normal velocity in Channel_Outflow()
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 Dstg_shur14.h30 * @param[out] ubar Mean velocity at wall_dist
146 * @param[in] ubar Mean velocity at X
149 * @param[out] u Velocity at X and t
183 * @param[in] ubar Mean velocity at X
190 * @param[out] u Velocity at X and t
330 * at each location, then calculate the actual velocity.
389 // Velocity normal to the boundary in StgShur14Inflow()
438 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()
440 // TODO This is almost certainly a bug. Velocity isn't stored here, only 0s. in StgShur14Inflow_Jacobian()
[all …]
H A Dstg_shur14_type.h17 CeedScalar u0; // !< Convective velocity
33 size_t ubar; // !< Mean velocity, [nprof, 3]
H A Dbc_outflow.h33 // The velocity in the exterior state has optional softplus regularization to
58 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 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()
/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/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 …]
/honee/examples/
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
/honee/src/
H A Dvelocity_gradient_projection.c4 /// Functions for setting up and projecting the velocity gradient
16 grad_velo_proj->num_comp = 9; // 9 velocity gradient in VelocityGradientProjectionCreateDM()
20 PetscCall(PetscObjectSetName((PetscObject)grad_velo_proj->dm, "Velocity Gradient Projection")); in VelocityGradientProjectionCreateDM()
90 …PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_rhs_assemble, "velocity gradient", grad_velo_proj->n… in VelocityGradientProjectionSetup()
97 …PetscCallCeed(ceed, CeedOperatorSetField(op_rhs_assemble, "velocity gradient", elem_restr_grad_vel… in VelocityGradientProjectionSetup()
/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 …PetscOptionsScalarArray("-freestream_velocity", "Velocity at freestream condition", NULL, Y_inf.ve… 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 Dsgs_dd_model.c54 #define GRAD_VELO_PROJ_KEY "Gradient of Velocity Projection"
133 { // Get velocity gradient information in SgsDDSetupNodalEvaluation_Fused()
135 …PetscCallCeed(ceed, CeedOperatorGetFieldByName(grad_velo_proj->l2_rhs_ctx->op, "velocity gradient"… in SgsDDSetupNodalEvaluation_Fused()
167 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_sgs_dd_nodal, "gradient velocity", num_comp_grad_velo… in SgsDDSetupNodalEvaluation_Fused()
175 …PetscCallCeed(ceed, CeedOperatorSetField(op_sgs_dd_nodal, "gradient velocity", elem_restr_grad_vel… in SgsDDSetupNodalEvaluation_Fused()
371 { // Get velocity gradient information in SgsDDSetupNodalEvaluation_Sequential()
373 …PetscCallCeed(ceed, CeedOperatorGetFieldByName(grad_velo_proj->l2_rhs_ctx->op, "velocity gradient"… in SgsDDSetupNodalEvaluation_Sequential()
414 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_sgs_dd_inputs, "gradient velocity", num_comp_grad_vel… in SgsDDSetupNodalEvaluation_Sequential()
422 …PetscCallCeed(ceed, CeedOperatorSetField(op_sgs_dd_inputs, "gradient velocity", elem_restr_grad_ve… in SgsDDSetupNodalEvaluation_Sequential()
661 …// TODO: Should probably move the elem_restr_q and basis_q creation to inside the velocity gradien… in SgsDDSetup()
/honee/src/smartsim/
H A Dsgs_dd_training.c110 { // -- Setup filtered velocity gradient projection in SetupTrainingDataCalculation()
119 // Get velocity gradient information in SetupTrainingDataCalculation()
120 …rGetFieldByName(sgs_dd_train->filtered_grad_velo_proj->l2_rhs_ctx->op, "velocity gradient", &op_fi… in SetupTrainingDataCalculation()
127 { // Get filtered velocity product information in SetupTrainingDataCalculation()
141 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_sgs_dd_train, "velocity product", num_comp_filtered_v… in SetupTrainingDataCalculation()
142 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_sgs_dd_train, "gradient velocity", num_comp_grad_velo… in SetupTrainingDataCalculation()
150 …PetscCallCeed(ceed, CeedOperatorSetField(op_sgs_dd_train, "velocity product", elem_restr_filtered_… in SetupTrainingDataCalculation()
151 …PetscCallCeed(ceed, CeedOperatorSetField(op_sgs_dd_train, "gradient velocity", elem_restr_grad_vel… in SetupTrainingDataCalculation()

12