Monday, 07/24: Writing a pm-code - AIP

PM: main code blocks 21
• Triangular Shaped Cloud (TSC):
S(x) = 1
∆x
{
1 − |x|/∆x, |x| < ∆x
0, otherwise
The fraction of particle’s mass assigned to a cell ijk is the shape function
averaged over this cell:
W (x p − x ijk ) =
x ijk +∆x/2
∫
dx ′ S(x p − x ′ );
x ijk −∆x/2
W (r p − r ijk ) = W (x p − x ijk )W (y p − y ijk )W (z p − z ijk );
The density in a cell ijk is then
ρ ijk =
N p
∑
p=1
m p W (r p − r ijk )