[Msure]

Hello everyone,

I wrote an OpenFOAM solver using two steps projection method, not using original PISO or PIMPLE in OpenFOAM. The second-order Adams-bashforth for time scheme is used.

        fvVectorMatrix UEqn
        (

          fvm::ddt(U)
          + 1.5*fvc::div(phi, U) - 0.5*fvc::div(phi.oldTime(), U.oldTime())
          == 1.5*fvm::laplacian(nu, U) - 0.5*fvc::laplacian(nu, U)
          + 0.5*fvc::grad(p)
         
        );

        fvScalarMatrix pEqn
        (
            fvm::laplacian(p) == 2.0/3.0*fvc::div(phi)/dt
        );

I got two questions, thanks for everyone who has some ideas about this:
1. How to add turbulence model?
Is it right that using nuEff() to replace the nu directly?
2. How to add VOF model?
I have no idea about this.

Thanks a lot!

Regards,
Shuo

[fedez91]

Hello Msure

1- You should check out the implementation of turbulence in other models, such as pisoFoam. It is done by creating a pointer in createFields.H that points to the selected turbulent closure. Then the function turbulence->divDevReff(U) is used to obtain the viscous and turbulent stresses. You can replicate this procedure in your solver. Check that other functions are used throughout the solver to, for example, solve the turbulent closure equations. You can take a look a this presentation for more information.

2- To add VOF take a look at the implementation of interFoam. You can check this thesis to see details of the implementation.

I hope this helps,
Federico

[Msure]

thanks Federico， these documents are useful.