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

question about gravity

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 13, 2009, 05:18
Default question about gravity
  #1
New Member
 
vivien
Join Date: Mar 2009
Posts: 18
Rep Power: 8
vivien is on a distinguished road
Hi, everybody,

I would like to simulate the contrast agent transport in blood flow in human vessel. Currently contrast agent is assumed to be a passive scalar transported by blood. So I use icoFoam coupled with caFoam in this way:
#include "fvCFD.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
int main(int argc, char *argv[])
{
# include "setRootCase.H"
# include "createTime.H"
# include "createMesh.H"
# include "createFields.H"
# include "initContinuityErrs.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Info<< "\nStarting time loop\n" << endl;
for (runTime++; !runTime.end(); runTime++)
{
Info<< "Time = " << runTime.timeName() << nl << endl;
# include "readPISOControls.H"
# include "CourantNo.H"
fvVectorMatrix UEqn
(
fvm::ddt(U)
+ fvm::div(phi, U)
- fvm::laplacian(nu, U)
);
solve(UEqn == -fvc::grad(p));
// --- PISO loop
for (int corr=0; corr<nCorr; corr++)
{
volScalarField rUA = 1.0/UEqn.A();
U = rUA*UEqn.H();
phi = (fvc::interpolate(U) & mesh.Sf())
+ fvc::ddtPhiCorr(rUA, U, phi);
adjustPhi(phi, U, p);
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
{
fvScalarMatrix pEqn
(
fvm::laplacian(rUA, p) == fvc::div(phi)
);
pEqn.setReference(pRefCell, pRefValue);
pEqn.solve();
if (nonOrth == nNonOrthCorr)
{
phi -= pEqn.flux();
}
}
# include "continuityErrs.H"
U -= rUA*fvc::grad(p);
U.correctBoundaryConditions();
}
// part added : copied from scalarTransportFoam.C
// beginning of the added part #####################################
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
{
solve
(
fvm::ddt(T)
+ fvm::div(phi, T)
- fvm::laplacian(DT, T)
);
}
// end of the added part ###########################################
runTime.write();
Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
<< nl << endl;
}
Info<< "End\n" << endl;
return(0);
}

However, we had observed some settings of CA at the bottom of vessel wall due to gravity. Does anybody have ideas how can I include gravity term into the solver for such case? should I solve for 2-phase flow(as I think it may due to the slightly different density of blood and CA which reult in the settings of CA)?

Thanks!!!

vivien
vivien is offline   Reply With Quote

Old   July 13, 2009, 08:50
Default
  #2
New Member
 
vivien
Join Date: Mar 2009
Posts: 18
Rep Power: 8
vivien is on a distinguished road
Can anybody give me suggestions?

Thanks in advance!

Vivien
vivien is offline   Reply With Quote

Old   July 13, 2009, 09:14
Default
  #3
Senior Member
 
Henrik Rusche
Join Date: Mar 2009
Location: Braunschweig, Niedersachsen, Germany
Posts: 275
Rep Power: 9
henrik is on a distinguished road
Dear Vivien,

settlingFoam might be a good starting point for you since the separation occurs due to a difference in density between the CA and the background fluid.

Henrik

Last edited by henrik; July 13, 2009 at 09:30.
henrik is offline   Reply With Quote

Old   July 13, 2009, 10:00
Default
  #4
New Member
 
vivien
Join Date: Mar 2009
Posts: 18
Rep Power: 8
vivien is on a distinguished road
Hi, Henrik,

Thanks for your reply. Maybe I can start to solve my problem step by step...because CA is always injected through a catheter into the vessel and then mixed with blood, do you think if I can still use icoFoam+passivescalartransport (as I posted at first)to solve the problem if CA here is assumed again to have the same density as blood(the geometry will be similar as a Y, 2 inlet2, 1 outlet, the CA and the blood will be mixed at bifurcation). My aim is to first simulate uniform mixing of 2 flow (same property, but different inlets, one inlet with CA). My idea is to give time varing velocity (because blood flow is pulsatile )for vessel inlet, and give fixed velocity to the catheter+ passivescalar transport. Do you think it is the correct way to implement?

Thanks!!

Vivien
vivien 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
how to consider gravity in CFX shrimp CFX 4 September 8, 2008 20:41
help!! Question about activating gravity mahe Phoenics 0 June 7, 2008 05:03
Gravity question zhou FLUENT 0 February 25, 2004 14:31
Help: gravity in CFX Dejun Jing CFX 2 July 22, 2002 08:58
On gravity modelling... Drona CFX 7 November 22, 2001 11:28


All times are GMT -4. The time now is 22:29.