CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   potentialFoam solver (https://www.cfd-online.com/Forums/openfoam-programming-development/96410-potentialfoam-solver.html)

dshawul January 22, 2012 09:44

potentialFoam solver
 
Hello,
I am confused about the basic potential foam equation in OF
Code:

fvm::laplacian
            (
                dimensionedScalar
                (
                    "1",
                    dimTime/p.dimensions()*dimensionSet(0, 2, -2, 0, 0),
                    1
                ),
                p
            )
                ==
            fvc::div(phi)

a) why is (1/ap) set to 1 unlike in other solvers such as simplefoam, pisofoam etc . I gather from Hrvoje's thesis (1/ap) should be used due to rhie-chow interpolation. Isn't the rhie-chow interpolation not used here?

b) Is p really pressure or velocity potential ? If it is velocity potential , then it would not require the rhs (i.e laplace(psi) = 0). Are we trying to solve stead-state euler equations Or the laplace equation for the velocity potential ?

I am really confused about this. please reply.
thanks

dshawul January 27, 2012 22:14

I have finally understood the algorithm after weeks of searching! But this board is the most unhelpful discussion board I have ever come across by far.
Just in case some one has similar difficulties understanding the algorithm , I don't want you to go through same diffculty so here it is.

* It is first of all solving the laplace equation laplace(psi) = 0
* There is a "hack" in createFields.H that zeros both p and U in the internal field so that only the U at the boundary will be used as a source term for the above equation.

The odd placement of that code in createFields.H was the main problem for me. Here is a link to the document that explained to me how the solver works. Thanks to the writer.

Antony_Bosco November 23, 2015 10:36

Hello everyone ,

Please let me know how to calculate the total pressure field (= static p + dynamic p) in case of potentialFoam.
i used potentialFoam -writep.

i have a doubt if p is actually p or something else.

Thanks in Advance,
antony


All times are GMT -4. The time now is 05:35.