|
[Sponsors] |
August 6, 2007, 12:20 |
Hi Hrvoje,
Sorry to not rep
|
#1 |
Member
Hoang-Lam
Join Date: Mar 2009
Location: Lausanne, Switzerland
Posts: 60
Rep Power: 17 |
Hi Hrvoje,
Sorry to not replying you for the last time (about a_N coefficients. I re-read your thesys, but it stil not ok for me), but in emmergency, I want to know where I'm wrong in the implementation of a Stokes solver: I begin with the simpleFOam.C file: so I just removed the div(phi,U) and changed some "minus" by "plus", according to he discretisation of Navier-Stokes, explained in your thesys. (in my simulation, I put turbulence OFF) **************************** tmp<fvvectormatrix> UEqn (turbulence->divR(U)); UEqn().relax(); solve(UEqn() == -fvc::grad(p)); p.boundaryField().updateCoeffs(); volScalarField AU = UEqn().A(); U = UEqn().H()/AU; UEqn.clear(); phi = fvc::interpolate(U) & mesh.Sf(); adjustPhi(phi, U, p); for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++) { fvScalarMatrix pEqn ( // I put a "minus" here, because now: Up=[H(U)+grad(p)]/ap, so the discretisation of div(U)=0 gives me the equation of pressure below: fvm::laplacian(1.0/AU, p) == -fvc::div(phi)); pEqn.setReference(pRefCell, pRefValue); pEqn.solve(); if (nonOrth == nNonOrthCorr) { // "+=" instead of "-=" because F = S.Uf = S.[H(U)/ap]f + S.[grad(p)/ap]f phi += pEqn.flux(); } } # include "continuityErrs.H" p.relax(); // "+=" instead of "-=" because Up=[H(U)+grad(p)]/ap U += fvc::grad(p)/AU; U.correctBoundaryConditions(); } ******************** the end of the program is the same. It seems so easy, but I can't get by myself. I just need your patience to tell me where I'm wrong. Thanks in advance, Lam |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
anything wrong in this udf | brusly | FLUENT | 0 | November 26, 2007 06:02 |
what's wrong with my UDF??? | lk zhu | FLUENT | 2 | June 16, 2006 02:52 |
what's wrong with 3.2 | F.K. | Siemens | 8 | August 9, 2004 02:44 |
What's wrong with my UDF? | olivia | FLUENT | 1 | June 23, 2001 17:06 |
wrong y+ | jurek | Main CFD Forum | 1 | August 26, 2000 21:01 |