Home
last modified time | relevance | path

Searched full:k (Results 1 – 25 of 28) sorted by relevance

12

/honee/qfunctions/
H A Deulervortex.h182 for (CeedInt k = 0; k < 3; k++) { // Columns of each Jacobian matrix in ConvectiveFluxJacobian_Euler() local
183 dF[i][0][k + 1] = ((i == k) ? 1. : 0.); in ConvectiveFluxJacobian_Euler()
184 …dF[i][j + 1][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (… in ConvectiveFluxJacobian_Euler()
185 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
295 for (CeedInt k = 0; k < 3; k++) { in Euler() local
296 drhodx[j] += drho[k] * dXdx[k][j]; in Euler()
297 dEdx[j] += dE[k] * dXdx[k][j]; in Euler()
299 dUdx[j][k] += dU[j][l] * dXdx[l][k]; in Euler()
300 dXdxdXdxT[j][k] += dXdx[j][l] * dXdx[k][l]; // dXdx_j,k * dXdx_k,j in Euler()
312 for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0.; in Euler() local
[all …]
H A Dshocktube.h115 for (CeedInt k = 0; k < 3; k++) { // Columns of each Jacobian matrix in ConvectiveFluxJacobian_Euler() local
116 dF[i][0][k + 1] = ((i == k) ? 1. : 0.); 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 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
245 for (CeedInt k = 0; k < 3; k++) { in EulerShockTube() local
246 du[j][k] = (dU[j][k] - drho[k] * u[j]) / rho; in EulerShockTube()
247 drhodx[j] += drho[k] * dXdx[k][j]; in EulerShockTube()
248 dEdx[j] += dE[k] * dXdx[k][j]; in EulerShockTube()
250 dUdx[j][k] += dU[j][l] * dXdx[l][k]; in EulerShockTube()
251 dXdxdXdxT[j][k] += dXdx[j][l] * dXdx[k][l]; // dXdx_j,k * dXdx_k,j in EulerShockTube()
[all …]
H A Dadvection.h215 for (CeedInt k = 0; k < 2; k++) { in StatePhysicalGradientFromReference_ND() local
218 … = grad_q[(Q * 5) * 0 + Q * j + i] * dXdx[0 * N + k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXdx[1 * … in StatePhysicalGradientFromReference_ND()
220 grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference_ND()
227 for (CeedInt k = 0; k < 3; k++) { in StatePhysicalGradientFromReference_ND() local
230 … = grad_q[(Q * 5) * 0 + Q * j + i] * dXdx[0 * N + k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXdx[1 * … in StatePhysicalGradientFromReference_ND()
231 grad_q[(Q * 5) * 2 + Q * j + i] * dXdx[2 * N + k]; in StatePhysicalGradientFromReference_ND()
233 grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference_ND()
293 v[f][i] = wdetJ * q_dot[f][i]; // K Mass/transient term in IFunction_AdvectionGeneric()
298 for (CeedInt k = 0; k < dim; k++) { in IFunction_AdvectionGeneric() local
299 div_u += grad_s[k].Y.velocity[j]; in IFunction_AdvectionGeneric()
[all …]
H A Dsetupgeo_helpers.h146 // dxdX_k,j * dxdX_j,k in InvertBoundaryMappingJacobian_3D()
149 for (CeedInt k = 0; k < 2; k++) { in InvertBoundaryMappingJacobian_3D() local
150 for (CeedInt l = 0; l < 3; l++) dxdXTdxdX[j][k] += dxdX[l][j] * dxdX[l][k]; in InvertBoundaryMappingJacobian_3D()
165 for (CeedInt k = 0; k < 3; k++) { in InvertBoundaryMappingJacobian_3D() local
166 dXdx[j][k] = 0; in InvertBoundaryMappingJacobian_3D()
167 for (CeedInt l = 0; l < 2; l++) dXdx[j][k] += dxdXTdxdX_inv[l][j] * dxdX[k][l]; in InvertBoundaryMappingJacobian_3D()
H A Dnewtonian_state.h46 …_QFUNCTION_HELPER CeedScalar Prandtl(NewtonianIGProperties gas) { return gas.cp * gas.mu / gas.k; } in Prandtl()
481 for (CeedInt k = 0; k < 3; k++) Flux[k + 1][j] = F_inviscid[j].momentum[k] - stress[k][j]; in FluxTotal() local
491 for (CeedInt k = 0; k < 3; k++) { in FluxTotal_Boundary() local
492 Flux[k + 1] += (F_inviscid[j].momentum[k] - stress[k][j]) * normal[j]; in FluxTotal_Boundary()
501 for (CeedInt k = 0; k < 3; k++) Flux[k + 1] = F_inviscid_normal.momentum[k]; in FluxTotal_RiemannBoundary() local
504 for (CeedInt k = 0; k < 3; k++) { in FluxTotal_RiemannBoundary() local
505 Flux[k + 1] -= stress[k][j] * normal[j]; in FluxTotal_RiemannBoundary()
563 …s[0][i] - Y.velocity[1] * stress[1][i] - Y.velocity[2] * stress[2][i] - gas.k * grad_s[i].Y.temper… in ViscousEnergyFlux()
571 … Y.velocity[2] * dstress[2][i] - dY.velocity[2] * stress[2][i] - gas.k * grad_ds[i].Y.temperature; in ViscousEnergyFlux_fwd()
590 for (CeedInt k = 0; k < 5; k++) dqi[k] = dq[k][j]; in StatePhysicalGradientFromReference() local
[all …]
H A Dnewtonian.h75 for (CeedInt k = 0; k < 3; k++) { in MassFunction_Newtonian() local
76 …Grad_v[k][j][i] = wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] + stab[j][2] * dXdx[k… in MassFunction_Newtonian()
161 // Fe = u Fu + k grad( T )
178 // k , Thermal conductivity
184 // We require the product of the inverse of the Jacobian (dXdx_j,k) and its transpose (dXdx_k,j) to…
227 …for (CeedInt k = 0; k < 3; k++) Grad_v[k][j][i] = wdetJ * (dXdx[k][0] * Flux[j][0] + dXdx[k][1] * … in RHSFunction_Newtonian() local
251 …for (CeedInt k = 0; k < 3; k++) Grad_v[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] *… in RHSFunction_Newtonian() local
347 for (CeedInt k = 0; k < 3; k++) in IFunction_Newtonian() local
348 …grad_v[k][j][i] = wdetJ * (grad_v_i[j][0] * dXdx[k][0] + grad_v_i[j][1] * dXdx[k][1] + grad_v_i[j]… in IFunction_Newtonian()
458 …for (int k = 0; k < 3; k++) grad_v[k][j][i] = wdetJ * (grad_v_i[j][0] * dXdx[k][0] + grad_v_i[j][1… in IJacobian_Newtonian() local
[all …]
H A Dvelocity_gradient_projection.h31 for (CeedInt k = 0; k < 3; k++) { in VelocityGradientProjectionRHS() local
32 v[j * 3 + k][i] = wdetJ * grad_velocity[j][k]; in VelocityGradientProjectionRHS()
H A Dutils.h131 … CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[i * N + k] * B[k * N + j]; in MatMatN() local
138 … CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[i * N + k] * B[j * N + k]; in MatMatN() local
149 … CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[k * N + i] * B[k * N + j]; in MatMatN() local
156 … CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[k * N + i] * B[j * N + k]; in MatMatN() local
233 for (CeedInt k = 0; k < P; k++) mat_C[i * M + j] += mat_A[i * P + k] * mat_B[k * M + j]; in MatMatNM() local
260 // @brief Calculate metric tensor from mapping, g_{ij} = xi_{k,i} xi_{k,j} = dXdx^T dXdx
H A Ddifferential_filter.h128 for (int k = 0; k < 3; k++) { in DifferentialFilter_LHS_N() local
129 …dq[k] = Grad_q[0 * N + j][i] * dXdx[0][k] + Grad_q[1 * N + j][i] * dXdx[1][k] + Grad_q[2 * N + j][… in DifferentialFilter_LHS_N()
133 for (int k = 0; k < 3; k++) { in DifferentialFilter_LHS_N() local
134 Grad_v[k * N + j][i] = wdetJ * dq_dXdx_a[k]; in DifferentialFilter_LHS_N()
H A Dsgs_dd_model.h209 for (CeedInt k = 0; k < 3; k++) Flux[k + 1][j] = sgs[k][j]; in FluxSubgridStress() local
234 for (CeedInt k = 0; k < 3; k++) { in IFunction_NodalSgs() local
236 …Grad_v[k][j][i] = -wdetJ * (dXdx[k][0] * Flux[j][0] + dXdx[k][1] * Flux[j][1] + dXdx[k][2] * Flux[… in IFunction_NodalSgs()
H A Dchannel.h34 const CeedScalar k = gas.k; in Exact_Channel() local
42 const CeedScalar Pr = mu / (cp * k); in Exact_Channel()
H A Dnewtonian_types.h23 CeedScalar k; member
H A Dutils_eigensolver_jacobi.h254 for (CeedInt k = 0; k < N; k++) SwapScalar(&evec[i * N + k], &evec[i_max * N + k]); in SortRows() local
/honee/doc/
H A Dreferences.bib54 author = {Papanastasiou, TC and Malamataris, N and Ellwood K},
65 …rry M and Wilhelmson, Robert B and Wicker, Louis J and Anderson, John R and Droegemeier, Kelvin K},
148 …thor = {Shur, Michael L. and Spalart, Philippe R. and Strelets, Michael K. and Travin, Andrey K.},
171 author = {Shakib, Farzin and Hughes, Thomas JR and Johan, Zden{\v{e}}k},
H A Druntime_options.md451 - `J/(kg K)`
456 - `J/(kg K)`
473 * - `-k`
476 - `W/(m K)`
H A Dexamples.md156 - `K`
161 - `K`
212 - `K`
323 - `K`
333 - `K`
H A Dtheory.md15 … \nabla \cdot \left( \frac{(E + P)\bm{U}}{\rho} -\bm{u} \cdot \bm{\sigma} - k \nabla T \right) - \…
20 … matrix, $\bm{b}$ is a body force vector (e.g., gravity vector $\bm{g}$), $k$ the thermal conduct…
63 - \bm{u} \cdot \bm{\sigma} - k \nabla T
79 \bm{q}_N (\bm{x},t)^{(e)} = \sum_{k=1}^{P}\psi_k (\bm{x})\bm{q}_k^{(e)}
595 kn[k^n];
641 | $\{\kappa^n\}_{n=1}^N$ | k^n | No | Yes |
/honee/tests/
H A Dchannel_divdiff_verify.yaml38 k: 5.07e-5
/honee/examples/
H A Dtaylor_green_vortex.yaml18 k: 0.8836267606 #thermal conductivity
H A Dvortexshedding.yaml40 k: 14.34 # thermal conductivity, Pr = 0.71 typical of air
/honee/problems/
H A Ddensitycurrent.c30 CeedScalar theta0 = 300.; // K in NS_DENSITY_CURRENT()
31 CeedScalar thetaC = -15.; // K in NS_DENSITY_CURRENT()
H A Dnewtonian.c395 .k = 0.02638, in NS_NEWTONIAN_IG()
422 …PetscCall(PetscOptionsScalar("-k", "Thermal conductivity", NULL, newtonian_ig_ctx->gas.k, &newtoni… in NS_NEWTONIAN_IG()
480 newtonian_ig_ctx->gas.k *= units->W_per_m_K; in NS_NEWTONIAN_IG()
735 const CeedScalar kg = units->kilogram, m = units->meter, sec = units->second, K = units->Kelvin; in UnitTests_Newtonian() local
739 const CeedScalar T = 200 * K; in UnitTests_Newtonian()
H A Dblasius.c150 CeedScalar T_inf = 288.; // K in NS_BLASIUS()
151 CeedScalar T_wall = 288.; // K in NS_BLASIUS()
/honee/
H A D.gitlab-ci.yml137 …- NPROC_TEST=1 make -k -j$((NPROC_CPU / NPROC_POOL)) CEED_BACKENDS="/cpu/self/memcheck" JUNIT_BATC…
195 …- NPROC_TEST=1 make -k -j$((NPROC_CPU / NPROC_POOL / 1)) CEED_BACKENDS="/cpu/self" JUNIT_BATCH="cp…
198 …- NPROC_TEST=1 make -k -j$((NPROC_CPU / NPROC_POOL / 1)) CEED_BACKENDS="/cpu/self" JUNIT_BATCH="cp…
255 …- NPROC_TEST=1 make -k -j$((NPROC_CPU / NPROC_POOL / 1)) CEED_BACKENDS="/cpu/self" JUNIT_BATCH="cp…
258 …- NPROC_TEST=1 make -k -j$((NPROC_CPU / NPROC_POOL / 1)) CEED_BACKENDS="/cpu/self" JUNIT_BATCH="cp…
317 # - NPROC_TEST=1 make -k -j$((NPROC_GPU / NPROC_POOL / 1)) CEED_BACKENDS="/gpu/hip/shared" JUNI…
377 …- NPROC_TEST=1 make -k -j$((NPROC_GPU / NPROC_POOL / 1)) CEED_BACKENDS="/gpu/cuda/shared" JUNIT_BA…
/honee/src/
H A Ddm-utils.c347 for (PetscInt k = 0; k < num_points; k++) { in PetscDTUniformTensorQuadrature() local
348 coords[((i * num_points + j) * num_points + k) * dim + 0] = coords_1d[i]; in PetscDTUniformTensorQuadrature()
349 coords[((i * num_points + j) * num_points + k) * dim + 1] = coords_1d[j]; in PetscDTUniformTensorQuadrature()
350 coords[((i * num_points + j) * num_points + k) * dim + 2] = coords_1d[k]; in PetscDTUniformTensorQuadrature()
351 …for (PetscInt c = 0; c < num_comp; c++) weights[((i * num_points + j) * num_points + k) * num_comp… in PetscDTUniformTensorQuadrature()

12