
[Sponsors] 
August 29, 2012, 11:26 
Need Confirmation  Eulerian solver for tiny particules

#1 
Senior Member
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 235
Rep Power: 7 
Hi everyone, it would be nice if somebody can give me some confirmations on the general aspect of my algorithm and the "code" that follows (especially the translation of the equations). I want to setup an unsteady Eulerian solver to compute the behavior of tiny particles in a fluid (only the fluid impact the particules motion through the drag).
Nomenclature :  volume fraction of the particles  velocity of the particles  function of the relative Reynolds number  time respond of the particle  velocity of the fluid Continuity equation : Momentum equation The algorithm I have in mind : 1) setup a "phi_alpha" flux equal to 2) initialize this flux in the createfield with: Code:
"linearInterpolate(alpha)*linearInterpolate(u) & mesh.Sf()" Code:
fvScalarMatrix cEqn ( fvm::ddt(alpha) + fvc::div(alpha_phi) ); cEqn.solve() Code:
alpha_phi = cEqn.flux(); 6) solve the momentum equation using: Code:
fvVectorMatrix uEqn ( fvm::ddt(alpha, u) + fvm::div(alpha_phi, u) == fRe / Tau * alpa * (V – u) ); uEqn.solve() i) I'm not sure of the step 4... I kind of reproduce the simpleFoam code. Is there any documentation somewhere about the command ".flux" ? ii) I'm not sure of the translation of the momentum equation step 6. Do I need to use a specific formulation for the source term ? Does the solver identify the last "u" by itself ? Thank you in advance for reading this. Last edited by fredo490; August 29, 2012 at 11:28. Reason: Tiny corrections 

September 13, 2012, 03:24 

#2 
Senior Member
Gerhard Holzinger
Join Date: Feb 2012
Location: Austria
Posts: 146
Rep Power: 13 
have a look into twoPhaseEulerFoam or bubbleFoam and the thesis of H. Rusche. They all deal with a Eulerian formulation of twophase fluid flows.
I would recommend to get into bubbleFoam or twoPhaseEulerFoam and remove everything you don't need. Particles with constant diameter and only drag as momentum exchange should not be a big problem. Last edited by GerhardHolzinger; September 13, 2012 at 03:27. Reason: added clarifiaction 

Tags 
equations, euler, particles 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
interfacing a fluid solver with abaqus  Tuhin Rakshit  Main CFD Forum  0  June 16, 2005 11:03 
Could you comare StarCD with CFX 5?Help, please...  Suteh  CFX  54  November 7, 2001 21:12 
CFX 5.5  Roued  CFX  1  October 2, 2001 17:49 
Setting a B.C using UserFortran in 4.3  tokai  CFX  10  July 17, 2001 17:25 
Error during Solver  cfd guy  CFX  4  May 8, 2001 07:04 