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

Doubts in SimpleFoam.....

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 22, 2013, 18:47
Default Doubts in SimpleFoam.....
  #1
New Member
 
Abhiroop
Join Date: Jun 2013
Posts: 26
Rep Power: 4
AB08 is on a distinguished road
Hi,, am very new to OpenFoam....... can sm1 please tell me

Q1) What is the exact function of fvOptions(U) in the UEqn.H file of simpleFoam? (I assumed that it gave the guessed value of the velocity(U) field (am I right) and if so from where does the initial guess come) ??



tmp<fvVectorMatrix> UEqn
(
fvm::div(phi, U)
+ turbulence->divDevReff(U)
==
fvOptions(U)
);

UEqn().relax();

fvOptions.constrain(UEqn());

solve(UEqn() == -fvc::grad(p));

fvOptions.correct(U);


Q2) What is the function of adjustPhi(phiHbyA, U, p); in the pEqn.H of simpleFoam?



{
volScalarField rAU(1.0/UEqn().A());
volVectorField HbyA("HbyA", U);
HbyA = rAU*UEqn().H();
UEqn.clear();

surfaceScalarField phiHbyA("phiHbyA", fvc::interpolate(HbyA) & mesh.Sf());
adjustPhi(phiHbyA, U, p);


fvOptions.relativeFlux(phiHbyA);

// Non-orthogonal pressure corrector loop
while (simple.correctNonOrthogonal())
{
fvScalarMatrix pEqn
(
fvm::laplacian(rAU, p) == fvc::div(phiHbyA)
);

pEqn.setReference(pRefCell, pRefValue);

pEqn.solve();

if (simple.finalNonOrthogonalIter())
{
phi = phiHbyA - pEqn.flux();
}
}

#include "continuityErrs.H"

// Explicitly relax pressure for momentum corrector
p.relax();

// Momentum corrector
U = HbyA - rAU*fvc::grad(p);
U.correctBoundaryConditions();
fvOptions.correct(U);
}
AB08 is offline   Reply With Quote

Old   June 25, 2013, 06:58
Default Doubts in scalarTransportFoam
  #2
New Member
 
Abhiroop
Join Date: Jun 2013
Posts: 26
Rep Power: 4
AB08 is on a distinguished road
Hi,,, i m new to OpenFoam nd need sm help .......

What is the function of fvOptions in the following????

#include "fvCFD.H"
#include "fvIOoptionList.H"
#include "simpleControl.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
{
#include "setRootCase.H"
#include "createTime.H"
#include "createMesh.H"
#include "createFields.H"
#include "createFvOptions.H"

simpleControl simple(mesh);

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Info<< "\nCalculating scalar transport\n" << endl;

#include "CourantNo.H"

while (simple.loop())
{
Info<< "Time = " << runTime.timeName() << nl << endl;

while (simple.correctNonOrthogonal())
{
solve
(
fvm::ddt(T)
+ fvm::div(phi, T)
- fvm::laplacian(DT, T)
==
fvOptions(T)
);
}

runTime.write();
}

Info<< "End\n" << endl;

return 0;
}
AB08 is offline   Reply With Quote

Old   July 2, 2013, 14:22
Default
  #3
New Member
 
Chen Yang
Join Date: Oct 2012
Posts: 3
Rep Power: 4
cyang is on a distinguished road
I also have the same questions. Does anyboody know the answers?
cyang is offline   Reply With Quote

Old   July 4, 2013, 06:47
Default
  #4
New Member
 
Abhiroop
Join Date: Jun 2013
Posts: 26
Rep Power: 4
AB08 is on a distinguished road
Hey cyang,, i dont have the answers yet....but hv a look here ,,might help

Doubt in scalarTransportFoam...
AB08 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
Laminar simpleFoam and inviscid simpleFoam herenger OpenFOAM Running, Solving & CFD 7 July 11, 2013 06:27
interFoam vs. simpleFoam channel flow comparison DanM OpenFOAM Running, Solving & CFD 11 January 5, 2013 07:21
Trying to run a benchmark case with simpleFoam spsb OpenFOAM 3 February 24, 2012 10:07
Naca0012 k-e mpirun gives fpe whereas simpleFoam not Pierpaolo OpenFOAM 1 May 8, 2010 03:08
Error running simpleFoam in parallel skabilan OpenFOAM Running, Solving & CFD 2 August 29, 2008 09:42


All times are GMT -4. The time now is 23:29.