# Need Confirmation - Eulerian solver for tiny particules

 Register Blogs Members List Search Today's Posts Mark Forums Read

 August 29, 2012, 10:26 Need Confirmation - Eulerian solver for tiny particules #1 Senior Member   HECKMANN Frédéric Join Date: Jul 2010 Posts: 236 Rep Power: 8 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()" 3) solve the continuity equation using the following code: Code: fvScalarMatrix cEqn ( fvm::ddt(alpha) + fvc::div(alpha_phi) ); cEqn.solve() 4) update the flux phi_alpha with: Code: alpha_phi -= cEqn.flux(); 5) compute the relative Reynolds number, the function and 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() My questions are: 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 #2 Senior Member   Gerhard Holzinger Join Date: Feb 2012 Location: Austria Posts: 173 Rep 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

 Tags equations, euler, particles

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Tuhin Rakshit Main CFD Forum 0 June 16, 2005 10:03 Suteh CFX 54 November 7, 2001 21:12 Roued CFX 1 October 2, 2001 16:49 tokai CFX 10 July 17, 2001 16:25 cfd guy CFX 4 May 8, 2001 06:04

All times are GMT -4. The time now is 19:56.