CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   How to obtain force value (http://www.cfd-online.com/Forums/openfoam-programming-development/90128-how-obtain-force-value.html)

ynsudo July 1, 2011 05:28

How to obtain force value
 
Dear all,
I am struggling with the following.
I am using interFoam (version 1.71) and want to obtain the each force-term appered in the code UEqn.H for a certain cell to know which term is large.

Those force terms are

Code:

        fvm::ddt(rho, U)
      + fvm::div(rhoPhi, U)
      - fvm::laplacian(muEff, U)
      - (fvc::grad(U) & fvc::grad(muEff))

                (
                    fvc::interpolate(interface.sigmaK())*fvc::snGrad(alpha1)
                  - ghf*fvc::snGrad(rho)
                  - fvc::snGrad(p_rgh)
                ) * mesh.magSf()

What should I add in the file UEqn.H to get those force values?
Any hints are appreciated.
Best regards,
Yu


All times are GMT -4. The time now is 09:48.