 |
 |
 |
|
|
|
[Sponsors] | ||||
problems about surfaceTension implement in reactingmultiphaseFoam |
 |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 24, 2019, 21:30
|
problems about surfaceTension implement in reactingmultiphaseFoam
|
#1 |
|
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",
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()
);
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());
}
![\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}}](/Forums/vbLatex/img/01797a06abdeab4ec42b15fe6bff50d1-1.gif "\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}}](/Forums/vbLatex/img/056728ea5daf97bf53a4847a6fe6dbaa-1.gif "\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? |
|
|
|
|
|
| Tags |
| reactingmultiphaseeuler, surfacetension |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| [ICEM] Problems with coedge curves and surfaces | tommymoose | ANSYS Meshing & Geometry | 6 | December 1, 2020 11:12 |
| Needed Benchmark Problems for FSI | Mechstud | Main CFD Forum | 4 | July 26, 2011 12:13 |
| modified SurfaceTension (with filmthickness) | Martin D! | Fluent UDF and Scheme Programming | 0 | March 31, 2010 11:13 |
| Help required to solve Hydraulic related problems | aero | CFX | 0 | October 30, 2006 11:00 |
Linear Mode
