cs_lagr_coupling.cpp File Reference

#include "base/cs_defs.h"
#include <limits.h>
#include <stdio.h>
#include <stddef.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <ctype.h>
#include <float.h>
#include <assert.h>
#include "base/cs_array.h"
#include "base/cs_base.h"
#include "base/cs_math.h"
#include "mesh/cs_mesh.h"
#include "mesh/cs_mesh_quantities.h"
#include "bft/bft_mem.h"
#include "bft/bft_error.h"
#include "base/cs_physical_constants.h"
#include "base/cs_time_step.h"
#include "turb/cs_turbulence_model.h"
#include "lagr/cs_lagr.h"
#include "lagr/cs_lagr_particle.h"
#include "lagr/cs_lagr_coupling.h"

Include dependency graph for cs_lagr_coupling.cpp:

Functions
void	cs_lagr_coupling (const cs_real_t **taup, const cs_real_t tempct[], cs_real_t tsfext[], const cs_real_3_t *force_p)
	Compute source terms for Lagrangian 2-way coupling.

Function Documentation

cs_lagr_coupling()

void cs_lagr_coupling	(	const cs_real_t **	taup,
		const cs_real_t	tempct[],
		cs_real_t	tsfext[],
		const cs_real_3_t *	force_p )

Compute source terms for Lagrangian 2-way coupling.

Remarks

Source terms are computed for the starting cell of a particle during a given iteration. Even if particle exits the domain, it s necessary to compute a source term matching the exchange between the carrier fluid and the particle at the beginning of the time step. If cs_glob_lagr_time_step->nor == 2 and the particle interacts with a boundary, then the source terms are computed as if nor == 1.

Parameters

[in]	taup	dynamic characteristic time
[in]	tempct	thermal characteristic time
[out]	tsfext	external forces
[in]	force_p	forces per mass unit on particles (m/s^2)