Rates.cl File Reference

Velocity and density variation rates computation. More...

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

Include dependency graph for Rates.cl:

Functions
__kernel void	entry (const __global uint *iset, const __global int *imove, const __global float *rho, const __global vec *grad_p, const __global vec *lap_u, const __global float *div_u, __global vec *dudt, __global float *drhodt, __constant float *visc_dyn, usize N, vec g)
	Velocity and density variation rates computation.

Detailed Description

Velocity and density variation rates computation.

Function Documentation

entry()

__kernel void entry	(	const __global uint *	iset,
		const __global int *	imove,
		const __global float *	rho,
		const __global vec *	grad_p,
		const __global vec *	lap_u,
		const __global float *	div_u,
		__global vec *	dudt,
		__global float *	drhodt,
		__constant float *	visc_dyn,
		usize	N,
		vec	g )

Velocity and density variation rates computation.

The mass conservation and momentum equations are applied from the already computed differential operators:

• \( \frac{\mathrm{d} \mathbf{u}}{\mathrm{d} t} = - \frac{\nabla p}{rho} + \frac{\mu}{rho} \Delta \mathbf{u} + \mathbf{g}\)
• \( \frac{\mathrm{d} \rho}{\mathrm{d} t} = - \rho \nabla \cdot \mathbf{u} + \delta \Delta t \frac{\rho_a}{\rho_0} \Delta p\)

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.
rho	Density \( \rho_{n+1} \).
grad_p	Pressure gradient \( \frac{\nabla p}{rho} \).
lap_u	Velocity laplacian \( \frac{\Delta \mathbf{u}}{rho} \).
div_u	Velocity divergence \( \rho \nabla \cdot \mathbf{u} \).
dudt	Velocity rate of change \( \left. \frac{d \mathbf{u}}{d t} \right\vert_{n+1} \).
drhodt	Density rate of change \( \left. \frac{d \rho}{d t} \right\vert_{n+1} \).
visc_dyn	Dynamic viscosity \( \mu \).
N	Number of particles.
g	Gravity acceleration \( \mathbf{g} \).