|
[Sponsors] |
MULESTemplates.C for compressible two phase flow |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
September 10, 2011, 16:23 |
MULESTemplates.C for compressible two phase flow
|
#1 |
Senior Member
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,267
Blog Entries: 1
Rep Power: 25 |
i have a two phase flow with variable density which change with pressure, temperature and chemical conversion!
look following step to derive a formulation for alphal1 1)ddt(rho) + div.(rho,U) = 0 and rho = rho1*alpha1+rho2*(1-alpha1)=(rho1-rho2)*alpha1+rho2 so: 2)ddt(rho_*alpha1) + div.(rho_,U) + ddt(rho2) + div.(rho2,U) = 0 and rho_ = rho1 - rho2 from MULESTemplates.C : fvScalarMatrix psiConvectionDiffusion ( fvm::ddt(rho, psi) + fv::gaussConvectionScheme<scalar>(mesh, phi, UDs).fvmDiv(phi, psi) //- fv::gaussLaplacianScheme<scalar, scalar>(mesh, CDs, snGrads) //.fvmLaplacian(Dpsif, psi) - fvm::Sp(Sp, psi) - Su ); and from incompressibleInterFoam or interPhaseChangeFoam : MULES::implicitSolve(oneField(), alphalpha1, phi, phiAlphalpha1, Sp, Su, 1, 0); can i use rho_ instead of oneField() in above structure? |
|
July 14, 2015, 08:12 |
|
#2 |
Senior Member
Przemek
Join Date: Jun 2011
Posts: 249
Rep Power: 16 |
Yes you have to do it.
__________________
best regards pblasiak |
|
Tags |
compressible, mules, two phase flow |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Two Phase Darcy flow with gravity -> problem | Cyp | OpenFOAM Programming & Development | 0 | November 17, 2010 08:58 |
two phase flow + ice phase and breaking | soeren87 | OpenFOAM Running, Solving & CFD | 0 | September 19, 2010 11:20 |
Validation for laminar Disperse phase flow | shashwat | Main CFD Forum | 0 | April 4, 2008 03:20 |
how 2 freeze 1 phase flow field & start lagrangian | KK | CFX | 5 | February 14, 2008 17:48 |
Two phase Flow | Ashok | Main CFD Forum | 0 | February 14, 2004 07:37 |