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

Calculus syntax

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 28, 2006, 12:25
Default I am wondering if some is able
  #1
Member
 
chris book
Join Date: Mar 2009
Posts: 85
Rep Power: 8
chris1980 is on a distinguished road
I am wondering if some is able to "translate" me the code in a more mathematican way.

phi = fvc::interpolate(rho)
*((fvc::interpolate(U) & mesh.Sf()) - fvc::meshPhi(rho, U));

> what is phi in this context? I do not understand how it is computed!

for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
{
fvScalarMatrix pEqn
(
fvm::ddt(psi, p)
+ fvc::div(phi)
- fvm::laplacian(rho*rUA, p)
);

>How to determine this equation? I think this is PISO related stuff.

pEqn.solve();
chris1980 is offline   Reply With Quote

Old   April 29, 2006, 09:26
Default Hi Chris, in foam usually a
  #2
Senior Member
 
Markus Hartinger
Join Date: Mar 2009
Posts: 102
Rep Power: 8
hartinger is on a distinguished road
Hi Chris,

in foam usually all variables are stored at the cell centre. But for calculating the divergence of a field the face values are needed. So fvc::interpolate returns a field containing the face values interpolated from the cell-centre values. The interpolation scheme is defined in system/fvScheme by the entry interpolationScheme, usually linear.

phi is in your example the mass-flux across the faces corrected for mesh-motion by fvc::meshPhi.

the pEqn is part of PISO, probably you look in one of the foam thesis' for it's derivation or ferziger and peric

regards
Markus
hartinger is offline   Reply With Quote

Reply

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
Param. & calculus in Fluent's Journal: Possible ? Kiril FLUENT 0 August 24, 2006 09:12
Syntax r2d2 OpenFOAM 1 December 13, 2005 14:19
Syntax What does mean evgenii OpenFOAM Pre-Processing 2 November 23, 2005 09:42
Syntax changes between versions 1.0.x and 1.1 OpenFOAM discussion board administrator OpenFOAM Announcements from ESI-OpenCFD 0 March 11, 2005 12:48
mesh and calculus servicies ro-design team Main CFD Forum 0 March 7, 2005 06:14


All times are GMT -4. The time now is 12:39.