CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Multiphase flow and conjugate heat transfer simulation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 27, 2009, 04:39
Default Multiphase flow and conjugate heat transfer simulation
  #1
Member
 
Jitao Liu
Join Date: Mar 2009
Location: Jinan , China
Posts: 64
Rep Power: 8
awacs is on a distinguished road
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 non-isothermal, non-newtonian 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.
awacs is offline   Reply With Quote

Old   August 27, 2009, 05:31
Default
  #2
Member
 
Jitao Liu
Join Date: Mar 2009
Location: Jinan , China
Posts: 64
Rep Power: 8
awacs is on a distinguished road
The modification is done under OpenFOAM-1.5-dev. 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
Attached Images
File Type: jpg RHCM-2D.jpg (49.7 KB, 86 views)
awacs is offline   Reply With Quote

Old   August 28, 2009, 09:43
Default
  #3
Member
 
Jitao Liu
Join Date: Mar 2009
Location: Jinan , China
Posts: 64
Rep Power: 8
awacs is on a distinguished road
I could moldfiy the nonNewtonianIcoFoam to conjugateNonNewtonianIcoHeatFoam by adding the tenperature equation . It is a transient solver for incompressible, laminar, single phase flow of non-isothermal, non-newtonian fluids with conjugate heat transfer.

But I am stilled confused with two-phase flows, especially the coupling of the two-phase 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.
2s-T-T.jpg 2s-U-T.jpg

2s-U-T2.jpg 2s-T-nu.jpg

Last edited by awacs; August 29, 2009 at 23:41.
awacs is offline   Reply With Quote

Old   September 1, 2009, 06:46
Default
  #4
Member
 
Jitao Liu
Join Date: Mar 2009
Location: Jinan , China
Posts: 64
Rep Power: 8
awacs is on a distinguished road
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 two-phase 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 two-phase weighted kT is varied as the two-fluid 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.
awacs is offline   Reply With Quote

Old   October 22, 2009, 14:12
Default
  #5
New Member
 
Edin Berberovic
Join Date: Mar 2009
Posts: 29
Rep Power: 8
eberberovic is on a distinguished road
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.
eberberovic is offline   Reply With Quote

Old   November 9, 2009, 07:49
Default
  #6
Member
 
Jitao Liu
Join Date: Mar 2009
Location: Jinan , China
Posts: 64
Rep Power: 8
awacs is on a distinguished road
Quote:
Originally Posted by eberberovic View Post
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.
Hi Edin,

Thanks a lot for your reply. I will try it out.

Regards,
Jitao
awacs is offline   Reply With Quote

Old   March 10, 2010, 21:56
Default
  #7
Senior Member
 
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 272
Rep Power: 9
phsieh2005 is on a distinguished road
Hi, Jitao,

I am wondering if you are successful in developing the new solver? I am also interested in such solver.

Pei

Quote:
Originally Posted by awacs View Post
Hi Edin,

Thanks a lot for your reply. I will try it out.

Regards,
Jitao
phsieh2005 is offline   Reply With Quote

Old   March 10, 2010, 22:56
Default
  #8
Senior Member
 
Sandy Lee
Join Date: Mar 2009
Posts: 207
Rep Power: 9
sandy is on a distinguished road
I am using OF 1.6.x, are there some new develpments about this kind of case? Who can give a note?
sandy is offline   Reply With Quote

Old   March 1, 2013, 06:25
Default
  #9
Member
 
nadine moussa
Join Date: Mar 2012
Posts: 30
Rep Power: 5
nadine is on a distinguished road
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
nadine is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
conjugate heat transfer of a turbine blade Steven-GY FLUENT 1 February 23, 2013 15: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 10:51
2-D conjugate heat transfer Barbaros CFX 7 August 12, 2006 01:45
Conjugate heat transfer for film-cooled flat plate Michele FLUENT 0 July 3, 2006 08:42


All times are GMT -4. The time now is 04:34.