Wall Model - IBM

Wall Model - IBM#

To couple the wall model to fluid flow simulation, the immersed boundary method (IBM) operates with a velocity correction calculated from the selected wall function.

From the a Lagrangian mesh some auxiliar points must be defined as illustrated below ( Shi et al.[1], Du et al.[2], Wang et al.[3]):

../../_images/ibm_wall_model.svg

where the point P is the reference point \(\mathrm{P_{N}}\) at normal direction in which the velocity will be linearly interpolated from nearest points. It is usual to fix a distance \(H\) of at least \(3\Delta x\) from the wall, which will be located in a point external to IBM diffusive layer. This velocity profile between P and the wall is modeled and the velocity at Q is calculated. The distance from W to Q is usually fixed \(h=0.1 - 0.25\Delta x\).

The point Q is also distanced from the wall towards normal direction and is an auxiliar point where a calculated velocity will be imposed. The correction will be applied using a interpolation-spread procedure from IBM at Q point.

Important

It is important that the simulation first reaches an stable flow with the standard IBM before activating the wall model. After activating the wall model, the standard IBM ceases its execution and is replaced by the wall model.

Hence the algorithm for the WM-IB-LBM executes the following sequence of operations:

  • Reference velocity \(\mathbf{u}\) is interpolated at P using a linear interpolation between nearest nodes and decomposed in its tangential \(t\) and normal \(n\) components with:

(1)#\[\mathbf{u}_{n} = \left(\mathbf{u}\cdot \mathbf{n}\right)\mathbf{n}\]
(2)#\[\mathbf{u}_{t} = \mathbf{u}-\mathbf{u}_{n}\]

  • The magnitude of \(\mathbf{u}_{t}\) is then assumed as \(u_{\mathrm{t,P_{N}}}\) and replaced in the wall function for the calculation of tangential velocity at Q.

  • The magnitude from \(u_{\mathrm{t,Q}}\) is multiplied by the tensor \(\mathbf{u}_{t}/ |\mathbf{u}_{t}|\) to obtain \(\mathbf{u}_{\mathrm{t,Q}}\)

  • The normal velocity at Q is given by a quadratic relation (Du et al.[2]):

(3)#\[\mathbf{u}_{\mathrm{n,Q}}=\left(\frac{h}{H}\right)^{2}\mathbf{u}_{n}\]

  • Then, the target velocity at Q is given by \(\mathbf{u}_{\mathrm{Q}}=\mathbf{u}_{\mathrm{t,Q}}+\mathbf{u}_{\mathrm{n,Q}}\).

  • The velocity correction is implemented using a traditional interpolation-spread procedure from Q.