
[Sponsors] 
Multiphase flow and conjugate heat transfer simulation 

LinkBack  Thread Tools  Search this Thread  Display Modes 
August 27, 2009, 04:39 
Multiphase flow and conjugate heat transfer simulation

#1 
Member
Jitao Liu
Join Date: Mar 2009
Location: Jinan , China
Posts: 64
Rep Power: 17 
Hi all,
I am trying to create a solver based on conjugateHeatFoam called conjugateInterHeatFoam. It is a transient solver for incompressible, laminar, two phase flow of nonisothermal, nonnewtonian fluids with conjugate heat transfer. All fields are created in createFields.H and createSolidFields.H. The solving sequence is: (1) gammaEqnSubCycle.H // solving the volume fraction equation (2) solveFluid.H // solving the flow equations (3)solveEnergy.H // solving the coupled energy equation // my solveEnergy.H { // Decoupled patches # include "attachPatches.H" // Solid side # include "readSolidControls.H" volScalarField mu = twoPhaseProperties.mu(); for (int nonOrth = 0; nonOrth <= nNonOrthCorr; nonOrth++) { coupledFvScalarMatrix TEqns(2); // Add fluid equation TEqns.set ( 0, new fvScalarMatrix ( rho * Cp * fvm::ddt(T) + rho * Cp * fvm::div(phi, T)  fvm::laplacian(kT, T)  mu*magSqr(fvc::grad(U)) //viscous disipation ) ); // Add solid equation TEqns.set ( 1, new fvScalarMatrix ( rho * Cp * fvm::ddt(Tsolid)  fvm::laplacian(kTsolid, Tsolid) ) ); TEqns.solve(); } } My questions are: (1) Do I need to couple rho, cp and kT ? Or just kT ? (2) How can I couple the conductivity kT ? In twophase flow, kT is a weighted average of the conductivity for both phase based on gamma (the volume fraction ): gamma*kT1 + (scalar(1)  gamma)*kT2 How to couple this weighted kT with the solid conductivity kT ? // kT definition in createFields.H // Thermal conductivity volScalarField kT ( IOobject ( "kT", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), gamma*kT1 + (scalar(1)  gamma)*kT2, gamma.boundaryField().types() ); //kT definition in createSolidFields.H Info<< "Reading solid conductivity kT\n" << endl; volScalarField kTsolid ( IOobject ( "kT", runTime.timeName(), solidMesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), solidMesh ); // On creation of both conductivities, couple the sides kT.correctBoundaryConditions(); kTsolid.correctBoundaryConditions(); Is it right? (3)How to treat with the boundary conditions at the twophase interface? Any suggestions would be appreciated. Thanks in advance:). Best Regards Jitao P.S. I could compile through the new solver with wmake. But the calculation of a 2D test stopped with a "segmentation fault" error. Debugging with gdb: Program received signal SIGSEGV, Segmentation fault. 0x00000000004628a5 in Foam::fvMatrix<double>::operator*=(Foam::Dimension edField<double, Foam::volMesh> const&) () Current language: auto; currently asm Last edited by awacs; August 31, 2009 at 03:13. 

August 27, 2009, 05:31 

#2 
Member
Jitao Liu
Join Date: Mar 2009
Location: Jinan , China
Posts: 64
Rep Power: 17 
The modification is done under OpenFOAM1.5dev. I want to simulate the process of viscous polymer melt filling into the heated mold cavity. The mold cavity can be heated both sides or just one side.The coupled temperature of the plastic part and the mold should be simulated.(see attachment)
Regards, Jitao 

August 28, 2009, 09:43 

#3 
Member
Jitao Liu
Join Date: Mar 2009
Location: Jinan , China
Posts: 64
Rep Power: 17 
I could moldfiy the nonNewtonianIcoFoam to conjugateNonNewtonianIcoHeatFoam by adding the tenperature equation . It is a transient solver for incompressible, laminar, single phase flow of nonisothermal, nonnewtonian fluids with conjugate heat transfer.
But I am stilled confused with twophase flows, especially the coupling of the twophase weighted kT with the solid conductivity kT ? Can any one help me ? Thanks a lot. Regards, Jitao ps Results of the modified conjugateNonNewtonianIcoHeatFoam. The effect of viscous dissipation is neglected. 2sTT.jpg 2sUT.jpg 2sUT2.jpg 2sTnu.jpg Last edited by awacs; August 29, 2009 at 23:41. 

September 1, 2009, 06:46 

#4 
Member
Jitao Liu
Join Date: Mar 2009
Location: Jinan , China
Posts: 64
Rep Power: 17 
I have settled some of the issues. The modified conjugateInterHeatFoam solver runs now.
A 2D test: 02TSolidAndPFluid.jpg 012TAndGamma.jpg 012TSolidAndTFluid.jpg 012TSolidAndUFluid.jpg 012TSolidAndViscosityFluid.jpg But I am still confused with the coupling of the twophase weighted kT with the solid conductivity kT.The following code gives a uniform kT. fluidTop { type regionCoupling; value uniform 0.03; remoteField kT; } But the twophase weighted kT is varied as the twofluid interface moving forward: // kT definition in createFields.H // Thermal conductivity volScalarField kT ( IOobject ( "kT", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), gamma*kT1 + (scalar(1)  gamma)*kT2, gamma.boundaryField().types() ); How to initialize this weighted kT in 0/kT ? Regards， Jitao Last edited by awacs; September 1, 2009 at 22:59. 

October 22, 2009, 14:12 

#5 
Member
Edin Berberovic
Join Date: Mar 2009
Posts: 31
Rep Power: 17 
Hello
You can initialize it normally using setFields (or funkySetFields), depending on your initial phase fraction distribution. But you will probably need to update the conductivity before solving energy. Something like kT == twoPhaseProperties.k(), if you have such a member. It will also probably be necessary to apply the coupling again, by correcting boundary conditions. Regards, Edin. 

November 9, 2009, 06:49 

#6  
Member
Jitao Liu
Join Date: Mar 2009
Location: Jinan , China
Posts: 64
Rep Power: 17 
Quote:
Thanks a lot for your reply. I will try it out. Regards, Jitao 

March 10, 2010, 20:56 

#7 
Senior Member
PeiYing Hsieh
Join Date: Mar 2009
Posts: 334
Rep Power: 18 

March 10, 2010, 21:56 

#8 
Senior Member
Sandy Lee
Join Date: Mar 2009
Posts: 213
Rep Power: 18 
I am using OF 1.6.x, are there some new develpments about this kind of case? Who can give a note?


March 1, 2013, 05:25 

#9 
Member
nadine moussa
Join Date: Mar 2012
Posts: 30
Rep Power: 14 
hello,
Are there any news about your implementation? I am also working with a multiphase flow in a mold so my goal is to couple conjugate heat transfer to multiphase flow! Any suggestion! Nadine Moussa 

Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
conjugate heat transfer of a turbine blade  StevenGY  FLUENT  1  February 23, 2013 14:31 
[ICEM] Conjugate CFD  Heat Transfer mesh  Subhadeep  ANSYS Meshing & Geometry  3  May 31, 2011 12:14 
Conjugate heat transfer with periodic boundaries  Suresh  FLUENT  0  February 23, 2009 09:51 
2D conjugate heat transfer  Barbaros  CFX  7  August 12, 2006 01:45 
Conjugate heat transfer for filmcooled flat plate  Michele  FLUENT  0  July 3, 2006 08:42 