Lines Matching refs:CeedScalar
25 CeedScalar delta0; // !< Boundary layer height at inflow
27 …CeedScalar T_wall; // !< Temperature at the w…
28 …CeedScalar x_inflow; // !< Location of inflow i…
29 …CeedScalar n_cheb; // !< Number of Chebyshev …
30 …CeedScalar *X; // !< Chebyshev polynomial…
31 …CeedScalar eta_max; // !< Maximum eta in the d…
32 …CeedScalar Tf_cheb[BLASIUS_MAX_N_CHEBYSHEV]; // !< Chebyshev coefficien…
33 …CeedScalar Th_cheb[BLASIUS_MAX_N_CHEBYSHEV - 1]; // !< Chebyshev coefficien…
74 …iusSolution)(const BlasiusContext blasius, const CeedScalar x[3], const CeedScalar x0, const CeedS… in BlasiusSolution()
75 const CeedScalar rho_infty, CeedScalar *t12) { in BlasiusSolution()
77 CeedScalar mu = blasius->newtonian_ctx.mu; in BlasiusSolution()
79 CeedScalar nu = mu / rho_infty; in BlasiusSolution()
80 CeedScalar U_infty = sqrt(Dot3(S_infty.Y.velocity, S_infty.Y.velocity)); in BlasiusSolution()
81 CeedScalar eta = x[1] * sqrt(U_infty / (nu * (x0 + x[0] - x_inflow))); in BlasiusSolution()
82 CeedScalar X = 2 * (eta / blasius->eta_max) - 1.; in BlasiusSolution()
83 CeedScalar Rd = GasConstant(&blasius->newtonian_ctx); in BlasiusSolution()
85 CeedScalar f[4], h[4]; in BlasiusSolution()
91 CeedScalar Y[5]; in BlasiusSolution()
103 CEED_QFUNCTION(ICsBlasius)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *ou… in ICsBlasius()
104 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsBlasius()
105 CeedScalar(*q0)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ICsBlasius()
109 const CeedScalar mu = context->newtonian_ctx.mu; in ICsBlasius()
110 const CeedScalar delta0 = context->delta0; in ICsBlasius()
111 const CeedScalar x_inflow = context->x_inflow; in ICsBlasius()
112 CeedScalar t12; in ICsBlasius()
115 const CeedScalar U_infty = sqrt(Dot3(S_infty.Y.velocity, S_infty.Y.velocity)); in ICsBlasius()
117 const CeedScalar x0 = U_infty * S_infty.U.density / (mu * 25 / Square(delta0)); in ICsBlasius()
120 const CeedScalar x[3] = {X[0][i], X[1][i], X[2][i]}; in ICsBlasius()
122 CeedScalar q[5] = {0}; in ICsBlasius()
131 CEED_QFUNCTION(Blasius_Inflow)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const… in Blasius_Inflow()
133 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in Blasius_Inflow()
134 const CeedScalar(*q_data_sur) = in[2]; in Blasius_Inflow()
135 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[3]; in Blasius_Inflow()
136 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Blasius_Inflow()
137 CeedScalar(*jac_data_sur) = context->newtonian_ctx.is_implicit ? out[1] : NULL; in Blasius_Inflow()
142 const CeedScalar rho_0 = S_infty.U.density; in Blasius_Inflow()
143 const CeedScalar U_infty = sqrt(Dot3(S_infty.Y.velocity, S_infty.Y.velocity)); in Blasius_Inflow()
144 …const CeedScalar x0 = U_infty * rho_0 / (gas->mu * 25 / Square(context->del… in Blasius_Inflow()
145 const CeedScalar zeros[11] = {0.}; in Blasius_Inflow()
148 CeedScalar wdetJb, norm[3]; in Blasius_Inflow()
153 const CeedScalar x[3] = {X[0][i], X[1][i], 0.}; in Blasius_Inflow()
154 CeedScalar t12; in Blasius_Inflow()
156 CeedScalar qi[5]; in Blasius_Inflow()
172 const CeedScalar stress[3][3] = { in Blasius_Inflow()
177 const CeedScalar Fe[3] = {0}; // TODO: viscous energy flux needs grad temperature in Blasius_Inflow()
178 CeedScalar Flux[5]; in Blasius_Inflow()
187 …NCTION(Blasius_Inflow_Jacobian)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in Blasius_Inflow_Jacobian()
188 const CeedScalar(*dq)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in Blasius_Inflow_Jacobian()
189 const CeedScalar(*q_data_sur) = in[2]; in Blasius_Inflow_Jacobian()
190 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[3]; in Blasius_Inflow_Jacobian()
191 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Blasius_Inflow_Jacobian()
196 const CeedScalar Rd = GasConstant(gas); in Blasius_Inflow_Jacobian()
197 const CeedScalar gamma = HeatCapacityRatio(gas); in Blasius_Inflow_Jacobian()
199 const CeedScalar rho_0 = S_infty.U.density; in Blasius_Inflow_Jacobian()
200 const CeedScalar U_infty = sqrt(Dot3(S_infty.Y.velocity, S_infty.Y.velocity)); in Blasius_Inflow_Jacobian()
201 …const CeedScalar x0 = U_infty * rho_0 / (gas->mu * 25 / Square(context->del… in Blasius_Inflow_Jacobian()
204 CeedScalar wdetJb, norm[3]; in Blasius_Inflow_Jacobian()
209 const CeedScalar x[3] = {X[0][i], X[1][i], X[2][i]}; in Blasius_Inflow_Jacobian()
210 CeedScalar t12; in Blasius_Inflow_Jacobian()
214 CeedScalar drho, dE, dP; in Blasius_Inflow_Jacobian()
218 CeedScalar dE_internal = drho * gas->cv * S_infty.Y.temperature; 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()