problems about surfaceTension implement in reactingmultiphaseFoam

yfClark · July 24, 2019, 22:30

Hello.I am studying the reactingmultiphaseEulerFoam,
the code for surfaceTension implement is:

Code:

      surfaceScalarField ghSnGradRho
       (
          "ghSnGradRho",
           ghf*fvc::snGrad(rho)*mesh.magSf()
        ); //g*h*\Delta \codt \rho 
        phigFs.set
        (
            phasei,
            (
                alpharAUfs[phasei]
               *(
                   ghSnGradRho
                 - (fvc::interpolate(phase.rho() - rho))*(g & mesh.Sf())
                 - fluid.surfaceTension(phase)*mesh.magSf()
                )
            ).ptr()
        );

finnaly,the velocity is updated by:

Code:

            mSfGradp = pEqnIncomp.flux()/rAUf;

            forAll(phases, phasei)
            {
                phaseModel& phase = phases[phasei];

                phase.U() =
                    HbyAs[phasei]
                  + fvc::reconstruct
                    (
                        alpharAUfs[phasei]*mSfGradp
                      - phigFs[phasei]
                    );
                phase.U().correctBoundaryConditions();
                fvOptions.correct(phase.U());
            }

My derivation mathematical formula is:
$\mathbf{U}_{i}=\mathbf{H} \mathbf{b y A}_{i}+\frac{-\alpha_{i} \nabla p_{r g h}}{A_{P, i}}-\frac{\left[\left(\alpha_{i} \mathbf{g} \cdot \mathbf{h} \nabla \rho-\alpha_{i}\left(\rho_{i}-\rho\right) \mathbf{g}\right)-\sigma \kappa \nabla \alpha_{i}\right]}{A_{p, i}}$
the code seems to be:
$\mathbf{U}_{i}=\mathbf{H} \mathbf{b y A}_{i}+\frac{-\alpha_{i} \nabla p_{r g h}}{A_{P, i}}-\frac{\left[\left(\alpha_{i} \mathbf{g} \cdot \mathbf{h} \nabla \rho-\alpha_{i}\left(\rho_{i}-\rho\right) \mathbf{g}\right)-\alpha_{i}\sigma \kappa \nabla \alpha_{i}\right]}{A_{p, i}}$
an addition alpha_{i} is add to surfaceTension,any one can show me what's the problem?

July 24, 2019, 22:30	problems about surfaceTension implement in reactingmultiphaseFoam	#1
yfClark New Member Yang Fan Join Date: Jul 2019 Posts: 1 Rep Power: 0	Hello.I am studying the reactingmultiphaseEulerFoam, the code for surfaceTension implement is: Code: surfaceScalarField ghSnGradRho ( "ghSnGradRho", ghffvc::snGrad(rho)mesh.magSf() ); //gh\Delta \codt \rho phigFs.set ( phasei, ( alpharAUfs[phasei] ( ghSnGradRho - (fvc::interpolate(phase.rho() - rho))(g & mesh.Sf()) - fluid.surfaceTension(phase)mesh.magSf() ) ).ptr() ); finnaly,the velocity is updated by: Code: mSfGradp = pEqnIncomp.flux()/rAUf; forAll(phases, phasei) { phaseModel& phase = phases[phasei]; phase.U() = HbyAs[phasei] + fvc::reconstruct ( alpharAUfs[phasei]mSfGradp - phigFs[phasei] ); phase.U().correctBoundaryConditions(); fvOptions.correct(phase.U()); } My derivation mathematical formula is: $\mathbf{U}_{i}=\mathbf{H} \mathbf{b y A}_{i}+\frac{-\alpha_{i} \nabla p_{r g h}}{A_{P, i}}-\frac{\left[\left(\alpha_{i} \mathbf{g} \cdot \mathbf{h} \nabla \rho-\alpha_{i}\left(\rho_{i}-\rho\right) \mathbf{g}\right)-\sigma \kappa \nabla \alpha_{i}\right]}{A_{p, i}}$ the code seems to be: $\mathbf{U}_{i}=\mathbf{H} \mathbf{b y A}_{i}+\frac{-\alpha_{i} \nabla p_{r g h}}{A_{P, i}}-\frac{\left[\left(\alpha_{i} \mathbf{g} \cdot \mathbf{h} \nabla \rho-\alpha_{i}\left(\rho_{i}-\rho\right) \mathbf{g}\right)-\alpha_{i}\sigma \kappa \nabla \alpha_{i}\right]}{A_{p, i}}$ an addition alpha_{i} is add to surfaceTension,any one can show me what's the problem?

Thread Tools	Search this Thread
Show Printable Version Email this Page	Search this Thread: Advanced Search
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[ICEM] Problems with coedge curves and surfaces	tommymoose	ANSYS Meshing & Geometry	6	December 1, 2020 12:12
Needed Benchmark Problems for FSI	Mechstud	Main CFD Forum	4	July 26, 2011 13:13
modified SurfaceTension (with filmthickness)	Martin D!	Fluent UDF and Scheme Programming	0	March 31, 2010 12:13
Help required to solve Hydraulic related problems	aero	CFX	0	October 30, 2006 12:00