Spectral Particle Methods
2018 SPHeric meeting
Galway, IE
Daniel Duque & Javier Calderon-Sanchez
CEHINAV, UPM, Spain
daniel.duque@upm.es
Why Fourier
dtdu=−∇p+ν∇2u+g
Pressure:
∇⋅u=0
dtduq=−iqpq−νq2uq+gq
q⋅uq=0
... but, of course
du/dt is a total derivative (hence, use particles)
uq should be computed on a mesh
(hence, use a mesh)
A compromise: use both, Particle-in-Cell style (also, pFEM v2.0)
Quick check: TG vortices
Kolmogorov flow
External periodic shear force:
g=f0cos(λ2πy)e^x
Kolmogorov force
Kolmogorov force
Onset of 2D turbulence
Original idea: P. Branson et al. CC BY-NC 4.0.
Ending: J.L. Cercos-Pita
2D turbulence: velocity
2D turbulence: acceleration
2D turbulence: power spectrum
Membranes: 2D fluids
Flow in membranes
Overdamped regime
![lo Re]()
Stokes flow
Flow driven only by forces. Stokes equation:
η∇2u−∇p=−f
∇⋅u=0
Solution:
uq,a=ηq21b∑[δab−q2qaqb]fq,a=:b∑Tabfq,b
Tab:=ηq21[δab−q2qaqb]
Phase segregation in membranes
![cover]()
Modeling phase segregation: CH
dtdϕ=M∇2μ
μ=δϕδ∫[−2aϕ2+4bϕ4+2c∣∇ϕ∣2]dV
Landau double-well
Equilibrium profile: μ=0
Closure: connection with 2D Stokes flow
Force is due to concentration gradients:
f=−ϕ∇μ
Also, the surrounding liquid modifies the Oseen tensor:
Tab:=ηm(q2+q/LSD)1[δab−q2qaqb]LSD=2ηfηm
If interested
Spectral Particle Methods
2018 SPHeric meeting
Galway, IE
Daniel Duque & Javier Calderon-Sanchez
CEHINAV, UPM, Spain
daniel.duque@upm.es