|
[Sponsors] |
June 22, 2013, 18:47 |
Doubts in SimpleFoam.....
|
#1 |
New Member
Abhiroop
Join Date: Jun 2013
Posts: 26
Rep Power: 12 |
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); } |
|
June 25, 2013, 06:58 |
Doubts in scalarTransportFoam
|
#2 |
New Member
Abhiroop
Join Date: Jun 2013
Posts: 26
Rep Power: 12 |
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; } |
|
July 2, 2013, 14:22 |
|
#3 |
New Member
Chen Yang
Join Date: Oct 2012
Posts: 3
Rep Power: 13 |
I also have the same questions. Does anyboody know the answers?
|
|
July 4, 2013, 06:47 |
|
#4 |
New Member
Abhiroop
Join Date: Jun 2013
Posts: 26
Rep Power: 12 |
Hey cyang,, i dont have the answers yet....but hv a look here ,,might help
http://www.cfd-online.com/Forums/ope...sportfoam.html |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
interFoam vs. simpleFoam channel flow comparison | DanM | OpenFOAM Running, Solving & CFD | 12 | January 31, 2020 15:26 |
Laminar simpleFoam and inviscid simpleFoam | herenger | OpenFOAM Running, Solving & CFD | 7 | July 11, 2013 06:27 |
Trying to run a benchmark case with simpleFoam | spsb | OpenFOAM | 3 | February 24, 2012 09: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 |