
[Sponsors] 
Regarding the adaptation of scalarTransport Foam as well as general openFOAM logic. 

LinkBack  Thread Tools  Display Modes 
August 10, 2010, 13:05 
Regarding the adaptation of scalarTransport Foam as well as general openFOAM logic.

#1 
New Member
Hunter
Join Date: Jul 2010
Posts: 2
Rep Power: 0 
Hello openFOAM Users!
I am attempting to implement the scalarTransport Foam to model gas separation using high speed centrifuges. I believe that I am close, real close, but I have a few unresolved issues with the fluid velocity field that I could use some help with. Here is the equation I am trying to solve: solve ( fvc::ddt(rho) fvm::laplacian(DT,rho) +fvm::div(phi,rho) ); Where: rho: fluid density at a given point DT: diffusion coefficient phi: fluid velocity (This I have some questions on as well) In the most simple form the problem can be modeled in one dimension. Therefore, I am declaring the y and z faces to be "empty" (this I learned from the user guide). Under these conditions the fluid velocity is described by the following term: U=A*x (where A is a collection of constants with units of (1/s)) Under more complicated conditions (i.e. with circulation of the fluid in the cylinder) the fluid velocity terms will take the form: U=(A*x+v(x), 0, v(z)) (where v(x) and v(z) are velocity components in the x and z direction, the problem assumes there is always symmetry about the y plane). //==========================// Questions: First, how is U related to phi? I can see there has to be a relationship, perhaps something like phi is generated from U, but I can't find the code that does it (perhaps I don't need/want to?). Second, how do I make the x component of U a function of x? I think I understand how to make U constant by using "uniform". I would like to do something like this: /////////////////////////////////// FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 1 0 0 0 0]; internalField uniform (A*x 0 0); ////////////////////////////////// But this obviously doesn't work. How do I go about generating U as a function of x? Finally, there is another form of my equation which looks more or less something like this: solve ( fvc::ddt(N) fvm::laplacian(DT,N) +fvm::div(phi,N(1N)) ); Where N is the number of particles. However, openFOAM won't compile the this as a solver. In fact it won't compile any solver in which the term, (N) in this case, has any operator on it (i.e 2*N, N+1, NN^2). I am positive that there is some sound mathematical/CFD reasoning for this, and I would like to understand it. If there is a way around it I would really like to know about it. The answers to these questions would really help me get my project off the ground. Please feel free to contact me with any questions, or comments if I have not made something clear. Thank you for your time. 

Tags 
scalartransportfoam, velocity field 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
CGNS converters available  mbeaudoin  OpenFOAM Meshing & Mesh Conversion  123  April 9, 2013 09:53 
BlockMesh FOAM warning  gaottino  OpenFOAM Native Meshers: blockMesh  7  July 19, 2010 14:11 
Import problem  ARC  Open Source Meshers: Gmsh, Netgen, CGNS, ...  0  February 27, 2010 10:56 
Axisymmetrical mesh  Rasmus Gjesing (Gjesing)  OpenFOAM Native Meshers: blockMesh  10  April 2, 2007 14:00 
Import gmsh msh to Foam  adorean  Open Source Meshers: Gmsh, Netgen, CGNS, ...  24  April 27, 2005 08:19 