.. include:: /images.rst
MMS04: No-slip wall BC, single-phase, 3D, curl-based functions
--------------------------------------------------------------
.. _description-6:
Description
~~~~~~~~~~~
The no-slip wall boundary condition in MFIX is verified using the techniques described in :cite:`Ani2014`; the manufactured solution is selected such that it satisfies both the divergence-free constraint and the no-slip wall boundary condition. Specifically, the no-slip wall boundary condition requires that the velocity at the (stationary) no-slip wall is zero. The manufactured solution is generated using the curl-based derivation to ensure divergence-free velocity fields :cite:`Ani2014` along with the technique given in :cite:`Bond2007` to ensure that the velocity component functions approach the value of zero at each boundary tested. The manufactured solution for the velocity field used for the verification of no-slip wall is given as :cite:`Ani2016`:
.. math::
\overrightarrow{V} = S^{2}\left( \overrightarrow{\nabla} \times \overrightarrow{H} \right) + 2S\left( \nabla S \times \overrightarrow{H} \right)
:label: mms4eq1
where, :math:`\overrightarrow{V}` is the velocity field vector, :math:`S` is the mathematical equation of the boundary being tested, and :math:`\overrightarrow{H}` is a general vector field consisting of sinusoidal expressions. The manufactured solution for pressure is selected as in :eq:`a1eq1` since there are no constraints on pressure with this boundary condition.
.. _setup-6:
Setup
~~~~~
This case is setup for single-phase flows on a domain with unit dimensions; the boundary tested is the West boundary (i.e., :math:`x = 0`).
.. _table16:
.. csv-table:: MMS-04 Setup, Initial and Boundary Conditions.
:widths: auto
:header: "Computational/Physical model", " ", " "
"3D, Steady-state, incompressible", " ", " "
"Single-phase (no solids)", " ", " "
"No gravity", " ", " "
"Thermal energy equations are not solved", " ", " "
"Turbulence equations are not solved (Laminar)", " ", " "
"Non-uniform mesh", " ", " "
"Central scheme", " ", " "
" ", " ", " "
"**Geometry**", " ", " "
"Coordinate system", "Cartesian", " "
"Domain length, :math:`L` \(x\)", "1.0", "\(m\)"
"Domain height, :math:`H` \(y\)", "1.0", "\(m\)"
"Domain width, :math:`W` \(z\)", "1.0", "\(m\)"
" ", " ", " "
"**Material** :sup:`†`", " ", " "
"Fluid density, :math:`\rho_{g}`", "1.0", "(kg·m\ :sup:`-3`)"
"Fluid viscosity, :math:`\mu_{g}`", "1.0", "(Pa·s)"
" ", " ", " "
"**Initial Conditions**", " ", " "
"Pressure *(gauge)*, :math:`P_{g}`", "MMS", "(Pa)"
"Fluid x-velocity, :math:`u_{g}`", "5.0", "(m·s\ :sup:`-1`)"
"Fluid y-velocity, :math:`v_{g}`", "5.0", "(m·s\ :sup:`-1`)"
"Fluid z-velocity, :math:`w_{g}`", "5.0", "(m·s\ :sup:`-1`)"
" ", " ", " "
"**Boundary Conditions** :sup:`‡`", " ", " "
"West boundary", "No-slip wall", " "
"All other boundaries", "Mass inflow", " "
**†** Material properties selected to ensure comparable contribution from convection and diffusion terms.
:sup:`‡` The manufactured solution is imposed on all boundaries (i.e., Dirichlet specification).
.. _results-6:
Results
~~~~~~~
Numerical solutions were obtained using the Central discretization scheme for 8x8, 16x16, 32x32, 64x64, and 128x128 grid meshes. Iterative convergence was not achieved for this case when pressure was solved. Hence, the pressure variable (:math:`P_{g}`) was fixed by specifying pressure using the manufactured solution in the initial conditions routine and discarding the pressure solution in the main solver routine. The observed order of accuracy matches the formal order as shown in :numref:`fig18` for the velocity variables.
.. _fig18:
.. figure:: /media/image26.png
:align: center
:figwidth: 50%
Observed orders of accuracy for no-slip wall verification (3D, single-phase flows) using :math:`\mathbf{L}_{\mathbf{2}}` and :math:`\mathbf{L}_{\mathbf{\infty}}` norms of the discretization error.