CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

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

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 10, 2010, 13:05
Default 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
HFonten1 is on a distinguished road
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(1-N))
);

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, N-N^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.
HFonten1 is offline   Reply With Quote

Reply

Tags
scalartransportfoam, velocity field

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 16:06.