Fixed ghost particles mirroring process. More...

#include "resources/Scripts/types/types.h"
#include "resources/Scripts/KernelFunctions/Kernel.h"

Include dependency graph for PreInteractions.cl:

Functions
__kernel void	entry (const __global uint iset, const __global int imove, const __global uint associations, const __global vec r, const __global vec normal, const __global vec dudt, const __global vec lap_u, const __global float gp_p, const __global vec gp_u, __global float rho, __global float p, __global vec u, __constant float *visc_dyn, uint N, vec g)
	Pressure and velocity extensions for the computation of the differential operators, except the Laplacian of the velocity already computed.

Detailed Description

Fixed ghost particles mirroring process.

Function Documentation

◆ entry()

__kernel void entry	(	const __global uint *	iset,
		const __global int *	imove,
		const __global uint *	associations,
		const __global vec *	r,
		const __global vec *	normal,
		const __global vec *	dudt,
		const __global vec *	lap_u,
		const __global float *	gp_p,
		const __global vec *	gp_u,
		__global float *	rho,
		__global float *	p,
		__global vec *	u,
		__constant float *	visc_dyn,
		uint	N,
		vec	g
	)

Pressure and velocity extensions for the computation of the differential operators, except the Laplacian of the velocity already computed.

For the Divergence of the velocity, the normmal velocity is mirrored (ASM), while the mirrored interpolated value is used for the tangential velocity (SSM).

For the pressure the following expression is used:

\( p(\mathbf{r}) = p(\hat{\mathbf{r}}) + \frac{\partial{p}}{\partial{\mathbf{n_b}}} \cdot \left( \hat{\mathbf{r}} - \mathbf{r}\right) \)

where the derivative of the pressure is computed as:

\( \frac{\partial{p}}{\partial{\mathbf{n_b}}} = \rho \left[ \mathbf{g} - \frac{\mathrm{d}\mathbf{u_b}}{\mathrm{d}t} + \frac{\mu}{\rho} \Delta \mathbf{u} \right] \cdot \mathbf{n_b} \)

Parameters

iset	Set of particles index.
imove	Moving flags. imove > 0 for regular fluid particles. imove = 0 for sensors. imove < 0 for boundary elements/particles.
associations	Mirroring particles.
r	Position \( \mathbf{r} \).
normal	Normal \( \mathbf{n} \).
dudt	Velocity rate of change \( \left. \frac{d \mathbf{u}}{d t} \right\vert_{n+1} \).
lap_u	Velocity laplacian \( \frac{\Delta \mathbf{u}}{rho} \).
gp_p	Interpolated pressure in the mirrored position \( p \).
gp_u	Interpolated velocity in the mirrored position \( \mathbf{u} \).
rho	Density \( \rho \).
p	Pressure \( p \).
u	Velocity \( \mathbf{u} \).
visc_dyn	Dynamic viscosity \( \mu \).
N	Number of particles.
g	Gravity acceleration \( \mathbf{g} \).