|August 29, 2012, 10:26||
Need Confirmation - Eulerian solver for tiny particules
Join Date: Jul 2010
Posts: 236Rep 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).
- 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 :
The algorithm I have in mind :
1) setup a "phi_alpha" flux equal to
2) initialize this flux in the createfield with:
"linearInterpolate(alpha)*linearInterpolate(u) & mesh.Sf()"
fvScalarMatrix cEqn ( fvm::ddt(alpha) + fvc::div(alpha_phi) ); cEqn.solve()
alpha_phi -= cEqn.flux();
6) solve the momentum equation using:
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 10:28. Reason: Tiny corrections
|September 13, 2012, 02:24||
Join Date: Feb 2012
Posts: 165Rep Power: 14
have a look into twoPhaseEulerFoam or bubbleFoam and the thesis of H. Rusche. They all deal with a Eulerian formulation of two-phase 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 02:27. Reason: added clarifiaction
|equations, euler, particles|
|Thread||Thread Starter||Forum||Replies||Last Post|
|interfacing a fluid solver with abaqus||Tuhin Rakshit||Main CFD Forum||0||June 16, 2005 10: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 16:49|
|Setting a B.C using UserFortran in 4.3||tokai||CFX||10||July 17, 2001 16:25|
|Error during Solver||cfd guy||CFX||4||May 8, 2001 06:04|