CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Periodic BC for Heated 2D Rib Channel

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

LinkBack Thread Tools Display Modes
Old   November 16, 2015, 09:24
Default Periodic BC for Heated 2D Rib Channel
Join Date: May 2015
Posts: 68
Rep Power: 4
hcl734 is on a distinguished road
I try to simulated a 2D channel with a rib and periodic boundary condtions. But I need to keep

I am using the cyclic boundary condition to couple inflow and outflow patch.

But since pressure drops and temperature rises from rib to rib I have to modify the inlet conditions to fit.

I would like to avoid modifying the solver, because I am not really familiar with the source code.

I found an option to keep my velocity (=mass flow, incompressible, same area, dm/dt = u*A*rho) constant using fvOptions.

fieldNames (U);
Ubar ( 1 0 0 );
gradPini gradPini [0 2 -2 0 0] 0; // initial pressure gradient
flowDir (1 0 0); // flow direction
But I also need to add a constant term (dq/dt)/((dm/dt)*cp*L) to compensate for the energy added by the heated walls.
How can I do this in OpenFOAM?

My idea so far is to set <<scalarExplicitSource>> for energy equation to subtract the term mentioned above, but I am not sure how to do it.

Something like:

    type            scalarExplicitSource; 
    active          true; 
    timeStart       0.2; // What are those entries for? 
    duration        2.0; 
    selectionMode   points; 
//[all points???] 

        volumeMode      absolute; 
            T         -(dq/dt)/((dm/dt)*cp*L); 

The idea of subtracting the source term I got from

Fully Developed Flow and Heat Transfer in Ducts Having Streamwise-Periodic Variations of Cross-Sectional Area. S. V. Patankar, C. H. Liu and E. M. Sparrow

Every advice is appreciated!

hcl734 is offline   Reply With Quote

Old   January 26, 2016, 07:23
New Member
Chris Coutinho
Join Date: Jan 2015
Location: Netherlands
Posts: 23
Rep Power: 4
cbcoutinho is on a distinguished road
Hello hcl734,

I am also working on streamwise-periodic flow with mass/heat transfer. I'm interested in implementing this in openFOAM, and it looks like you are working on this as well. Have you made any progress since you asked this question?

cbcoutinho is online now   Reply With Quote

Old   January 26, 2016, 08:13
Join Date: May 2015
Posts: 68
Rep Power: 4
hcl734 is on a distinguished road
Hello Chris,

pressureGradientExplicitSourceCoeffs works fine to set a constant mass flux.

And I implemented the heat source the way patankar suggests in the mentioned paper following this thread

It works for me on one grid (mass flow and temperature become constant) but fails at another which is strange and I don't understand yet.
I can share my calculations with you privately but can't publish them here, since they are part of my thesis, pn me.

hcl734 is offline   Reply With Quote


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
periodic channel flow roadRunner STAR-CCM+ 0 June 25, 2015 07:32
Periodic channel flow with time dependent mass flow rate QBeast FLUENT 3 May 10, 2013 13:14
[ANSYS Meshing] Meshing a Periodic channel with high curvature 2D airfoil ergometer ANSYS Meshing & Geometry 4 February 22, 2012 10:13
[ICEM] how to define periodic in icem ? mingersai ANSYS Meshing & Geometry 1 February 3, 2012 18:46
Heat conduction through solid rib in periodic duct Dieter Fauconnier FLUENT 0 October 7, 2004 08:22

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