Lines Matching refs:CeedScalar
23 CeedScalar theta0; // !< Reference temperature
24 CeedScalar P0; // !< Reference Pressure
25 CeedScalar umax; // !< Centerline velocity
26 CeedScalar center; // !< Y Coordinate for center of channel
27 CeedScalar H; // !< Channel half-height
28 CeedScalar B; // !< Body-force driving the flow
32 CEED_QFUNCTION_HELPER State Exact_Channel(CeedInt dim, CeedScalar time, const CeedScalar X[], CeedI… in Exact_Channel()
34 const CeedScalar theta0 = context->theta0; in Exact_Channel()
35 const CeedScalar P0 = context->P0; in Exact_Channel()
36 const CeedScalar umax = context->umax; in Exact_Channel()
37 const CeedScalar center = context->center; in Exact_Channel()
38 const CeedScalar H = context->H; in Exact_Channel()
40 const CeedScalar cp = gas->cp; in Exact_Channel()
41 const CeedScalar mu = gas->mu; in Exact_Channel()
42 const CeedScalar k = gas->k; in Exact_Channel()
48 const CeedScalar x[3] = {0, X[1], X[2]}; in Exact_Channel()
50 const CeedScalar Pr = mu / (cp * k); in Exact_Channel()
51 const CeedScalar Ec = (umax * umax) / (cp * theta0); in Exact_Channel()
52 const CeedScalar theta = theta0 * (1 + (Pr * Ec / 3) * (1 - Square(Square((x[1] - center) / H)))); in Exact_Channel()
53 CeedScalar Y[5] = {0.}; in Exact_Channel()
66 CEED_QFUNCTION(ICsChannel)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *ou… in ICsChannel()
67 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsChannel()
68 CeedScalar(*q0)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ICsChannel()
74 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsChannel()
76 CeedScalar q[5] = {0}; in ICsChannel()
86 CEED_QFUNCTION(Channel_Inflow)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const… in Channel_Inflow()
87 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in Channel_Inflow()
88 const CeedScalar(*q_data_sur) = in[2]; in Channel_Inflow()
89 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[3]; in Channel_Inflow()
90 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Channel_Inflow()
95 const CeedScalar gamma = HeatCapacityRatio(&context->newtonian_ctx); in Channel_Inflow()
98 CeedScalar wdetJb, norm[3]; in Channel_Inflow()
107 const CeedScalar x[3] = {0, X[1][i], X[2][i]}; in Channel_Inflow()
111 CeedScalar q_exact[5] = {0.}; in Channel_Inflow()
115 CeedScalar q_inside[5] = {0}; in Channel_Inflow()
118 const CeedScalar P = s_inside.Y.pressure; in Channel_Inflow()
121 const CeedScalar e_internal = gas->cv * s_exact.Y.temperature; in Channel_Inflow()
122 const CeedScalar rho_in = P / ((gamma - 1) * e_internal); in Channel_Inflow()
123 const CeedScalar E_kinetic = .5 * rho_in * Dot3(s_exact.Y.velocity, s_exact.Y.velocity); in Channel_Inflow()
124 const CeedScalar E = rho_in * e_internal + E_kinetic; in Channel_Inflow()
130 const CeedScalar u_normal = Dot3(norm, s_exact.Y.velocity); in Channel_Inflow()
148 CEED_QFUNCTION(Channel_Outflow)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *cons… in Channel_Outflow()
149 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in Channel_Outflow()
150 const CeedScalar(*q_data_sur) = in[2]; in Channel_Outflow()
151 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Channel_Outflow()
157 CeedScalar wdetJb, norm[3]; in Channel_Outflow()
161 const CeedScalar rho = q[0][i]; in Channel_Outflow()
162 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Channel_Outflow()
163 const CeedScalar E = q[4][i]; in Channel_Outflow()
170 const CeedScalar P = context->P0; // pressure in Channel_Outflow()
171 const CeedScalar u_normal = Dot3(norm, u); // Normal velocity in Channel_Outflow()