# Doubts in SimpleFoam.....

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

 June 22, 2013, 18:47 Doubts in SimpleFoam..... #1 New Member   Abhiroop Join Date: Jun 2013 Posts: 26 Rep Power: 5 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 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: 5 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: 5 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: 5 Hey cyang,, i dont have the answers yet....but hv a look here ,,might help Doubt in scalarTransportFoam...

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post herenger OpenFOAM Running, Solving & CFD 7 July 11, 2013 06:27 DanM OpenFOAM Running, Solving & CFD 11 January 5, 2013 07:21 spsb OpenFOAM 3 February 24, 2012 10:07 Pierpaolo OpenFOAM 1 May 8, 2010 03:08 skabilan OpenFOAM Running, Solving & CFD 2 August 29, 2008 09:42

All times are GMT -4. The time now is 19:38.