30 f_y, f_z, c_Xh, gs_Xh, bc_prs_surface,bc_sym_surface, Ax, bd, dt, mu, rho)
31 type(
field_t),
intent(inout) :: p, u, v, w
32 type(
field_t),
intent(inout) :: u_e, v_e, w_e
33 type(
field_t),
intent(inout) :: p_res
34 type(
field_t),
intent(inout) :: f_x, f_y, f_z
35 type(
coef_t),
intent(inout) :: c_Xh
36 type(
gs_t),
intent(inout) :: gs_Xh
39 class(
ax_t),
intent(inout) :: Ax
40 real(kind=
rp),
intent(inout) :: bd
41 real(kind=
rp),
intent(in) :: dt
42 real(kind=
rp),
intent(in) :: mu
43 real(kind=
rp),
intent(in) :: rho
47 type(
field_t),
pointer :: ta1, ta2, ta3, wa1, wa2, wa3, work1, work2
48 integer :: temp_indices(8)
62 c_xh%h1(i,1,1,1) = 1.0_rp / rho
63 c_xh%h2(i,1,1,1) = 0.0_rp
67 call curl(ta1, ta2, ta3, u_e, v_e, w_e, work1, work2, c_xh)
68 call curl(wa1, wa2, wa3, ta1, ta2, ta3, work1, work2, c_xh)
71 wa1%x(i,1,1,1) = (wa1%x(i,1,1,1) * mu / rho) * c_xh%B(i,1,1,1)
72 wa2%x(i,1,1,1) = (wa2%x(i,1,1,1) * mu / rho) * c_xh%B(i,1,1,1)
73 wa3%x(i,1,1,1) = (wa3%x(i,1,1,1) * mu / rho) * c_xh%B(i,1,1,1)
77 ta1%x(i,1,1,1) = f_x%x(i,1,1,1) / rho - wa1%x(i,1,1,1)
78 ta2%x(i,1,1,1) = f_y%x(i,1,1,1) / rho - wa2%x(i,1,1,1)
79 ta3%x(i,1,1,1) = f_z%x(i,1,1,1) / rho - wa3%x(i,1,1,1)
82 call gs_xh%op(ta1, gs_op_add)
83 call gs_xh%op(ta2, gs_op_add)
84 call gs_xh%op(ta3, gs_op_add)
87 ta1%x(i,1,1,1) = ta1%x(i,1,1,1) * c_xh%Binv(i,1,1,1)
88 ta2%x(i,1,1,1) = ta2%x(i,1,1,1) * c_xh%Binv(i,1,1,1)
89 ta3%x(i,1,1,1) = ta3%x(i,1,1,1) * c_xh%Binv(i,1,1,1)
92 call ax%compute(p_res%x,p%x,c_xh,p%msh,p%Xh)
94 call cdtp(wa1%x, ta1%x, c_xh%drdx, c_xh%dsdx, c_xh%dtdx, c_xh)
95 call cdtp(wa2%x, ta2%x, c_xh%drdy, c_xh%dsdy, c_xh%dtdy, c_xh)
96 call cdtp(wa3%x, ta3%x, c_xh%drdz, c_xh%dsdz, c_xh%dtdz, c_xh)
98 p_res%x(i,1,1,1) = (-p_res%x(i,1,1,1)) &
99 + wa1%x(i,1,1,1) + wa2%x(i,1,1,1) + wa3%x(i,1,1,1)
106 wa1%x(i,1,1,1) = 0.0_rp
107 wa2%x(i,1,1,1) = 0.0_rp
108 wa3%x(i,1,1,1) = 0.0_rp
111 call bc_sym_surface%apply_surfvec(wa1%x,wa2%x,wa3%x,ta1%x, ta2%x, ta3%x, n)
115 ta1%x(i,1,1,1) = 0.0_rp
116 ta2%x(i,1,1,1) = 0.0_rp
117 ta3%x(i,1,1,1) = 0.0_rp
120 call bc_prs_surface%apply_surfvec(ta1%x, ta2%x, ta3%x, u%x, v%x, w%x, n)
123 p_res%x(i,1,1,1) = p_res%x(i,1,1,1) &
124 - (dtbd * (ta1%x(i,1,1,1) + ta2%x(i,1,1,1) + ta3%x(i,1,1,1)))&
125 - (wa1%x(i,1,1,1) + wa2%x(i,1,1,1) + wa3%x(i,1,1,1))
132 p, f_x, f_y, f_z, c_Xh, msh, Xh, mu, rho, bd, dt, n)
133 class(
ax_t),
intent(in) :: Ax
134 type(
mesh_t),
intent(inout) :: msh
135 type(
space_t),
intent(inout) :: Xh
136 type(
field_t),
intent(inout) :: p, u, v, w
137 type(
field_t),
intent(inout) :: u_res, v_res, w_res
138 type(
field_t),
intent(inout) :: f_x, f_y, f_z
139 type(
coef_t),
intent(inout) :: c_Xh
140 real(kind=
rp),
intent(in) :: mu
141 real(kind=
rp),
intent(in) :: rho
142 real(kind=
rp),
intent(in) :: bd
143 real(kind=
rp),
intent(in) :: dt
144 integer,
intent(in) :: n
145 integer :: temp_indices(3)
146 type(
field_t),
pointer :: ta1, ta2, ta3
150 c_xh%h1(i,1,1,1) = mu
151 c_xh%h2(i,1,1,1) = rho * (bd / dt)
155 call ax%compute(u_res%x, u%x, c_xh, msh, xh)
156 call ax%compute(v_res%x, v%x, c_xh, msh, xh)
157 call ax%compute(w_res%x, w%x, c_xh, msh, xh)
163 call opgrad(ta1%x, ta2%x, ta3%x, p%x, c_xh)
166 u_res%x(i,1,1,1) = (-u_res%x(i,1,1,1)) - ta1%x(i,1,1,1) + f_x%x(i,1,1,1)
167 v_res%x(i,1,1,1) = (-v_res%x(i,1,1,1)) - ta2%x(i,1,1,1) + f_y%x(i,1,1,1)
168 w_res%x(i,1,1,1) = (-w_res%x(i,1,1,1)) - ta3%x(i,1,1,1) + f_z%x(i,1,1,1)
Defines a Matrix-vector product.
Dirichlet condition applied in the facet normal direction.
integer, parameter, public rp
Global precision used in computations.
subroutine, public cdtp(dtx, x, dr, ds, dt, coef)
Apply D^T to a scalar field, where D is the derivative matrix.
subroutine, public opgrad(ux, uy, uz, u, coef, es, ee)
Compute the gradient of a scalar field.
subroutine, public curl(w1, w2, w3, u1, u2, u3, work1, work2, coef)
Residuals in the Pn-Pn formulation (SX version)
subroutine pnpn_vel_res_sx_compute(Ax, u, v, w, u_res, v_res, w_res, p, f_x, f_y, f_z, c_Xh, msh, Xh, mu, rho, bd, dt, n)
subroutine pnpn_prs_res_sx_compute(p, p_res, u, v, w, u_e, v_e, w_e, f_x, f_y, f_z, c_Xh, gs_Xh, bc_prs_surface, bc_sym_surface, Ax, bd, dt, mu, rho)
Defines Pressure and velocity residuals in the Pn-Pn formulation.
Defines a registry for storing and requesting temporary fields This can be used when you have a funct...
type(scratch_registry_t), target, public neko_scratch_registry
Global scratch registry.
Defines a function space.
Base type for a matrix-vector product providing .
Coefficients defined on a given (mesh, ) tuple. Arrays use indices (i,j,k,e): element e,...
Dirichlet condition in facet normal direction.
Abstract type to compute pressure residual.
Abstract type to compute velocity residual.
The function space for the SEM solution fields.