# MULESTemplates.C for compressible two phase flow

 Register Blogs Members List Search Today's Posts Mark Forums Read

 September 10, 2011, 15:23 MULESTemplates.C for compressible two phase flow #1 Senior Member   Nima Samkhaniani Join Date: Sep 2009 Location: Tehran, Iran Posts: 1,233 Blog Entries: 1 Rep Power: 18 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(mesh, phi, UDs).fvmDiv(phi, psi) //- fv::gaussLaplacianScheme(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, 07:12 #2 Senior Member   Przemek Join Date: Jun 2011 Posts: 162 Rep Power: 8 Yes you have to do it. __________________ best regards pblasiak

 Tags compressible, mules, two phase flow

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Cyp OpenFOAM Programming & Development 0 November 17, 2010 08:58 soeren87 OpenFOAM Running, Solving & CFD 0 September 19, 2010 10:20 shashwat Main CFD Forum 0 April 4, 2008 02:20 KK CFX 5 February 14, 2008 17:48 Ashok Main CFD Forum 0 February 14, 2004 07:37

All times are GMT -4. The time now is 13:39.