
[Sponsors] 
May 10, 2011, 07:11 
Streamwise periodic heat trasnfer

#1 
Senior Member
Robert Sawko
Join Date: Mar 2009
Posts: 117
Rep Power: 22 
Dear All,
I am trying to implement heat transfer in interFoam but coupled with periodic boundary conditions in streamwise direction. In this case the temperature field is not periodic in any obvious way. But if I focus on a 2D channel case and impose a constant wall temperature and a bulk inlet temperature a renormalised temperature field \theta = (T  Tw)/(Tbin  Tw) is selfsimilar in the sense that if I rescale the profile at the inlet it should coincide with the profile at the outlet. My question is: how can I do it with interFoam on a cyclic patch? Can I apply anything else than a cyclic boundary condition on this patch? The rescaling would be fairly easy if I had a direct access to values. Is there any easy way of doing it? Something like writing your new BC applicable on cyclic patch would be great. I would be obliged for any comments or hints . So far I've been using this thread as a guidance: http://www.cfdonline.com/Forums/ope...terfoam3.html 

May 10, 2011, 09:34 

#2 
Senior Member
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 15 
I cannot help you with this problem.
But why do you want to do it like this? Normally (in literature) temperature is made periodic using a source term (I performed 3D periodic channel and pipe flows using this source). References are Kawamura, Tiselj, Bergant (channel) and Saad (pipe). Is there a specific reason why to avoid their approach? If you want I can send you one of the references. 

May 10, 2011, 10:50 

#3 
Senior Member
Robert Sawko
Join Date: Mar 2009
Posts: 117
Rep Power: 22 
Thank you for your comment!
I have to confess that ignorance occurs to me as the only reason. To which papers do you refer exactly? I downloaded the papers by Kawamura 1998 and Tiselj 2001. The latter one entitled "Effect of wall boundary condition on scalar transfer in a fully developed turbulent flume" talks about constant wall temperature so I am looking at this one at the moment. I can't see exactly why \theta is to be periodic after inclusion of the source term u^+/2u^+_b in eq (4). In constant heat flux scenario it's correct but for constant wall temperature... shouldn't it be different? 

May 11, 2011, 07:09 

#4 
Senior Member
Robert Sawko
Join Date: Mar 2009
Posts: 117
Rep Power: 22 
I thought I'll post a clarification after my yesterday evening reading.
@stevenvanharen: You were right about the source term. This seems to be an expedient way of doing it. But I still can't see how these references can be applied to isothermal wall condition. But I've found this: S. V. Patankar, C. H. Liu, and E. M. Sparrow. Fully Developed Flow and Heat Transfer in Ducts Having StreamwisePeriodic Variations of CrossSectional Area. ASME J. of Heat Transfer, 99:180186, 1977. which was actually my initial reading on streamwise periodic heat flow. The problem of 'selfsimilar' profiles is welladdressed there. The difference is that the source term for normalised temperature has to be obtained by solving another equation (no 38 in the publication). 

May 18, 2011, 11:20 

#5 
Senior Member
Robert Sawko
Join Date: Mar 2009
Posts: 117
Rep Power: 22 
I had a major breakthrough in my problem and now I can gladly report that it actually works. Not everything is sorted out yet, but it seems to be progressing. Steven was right about the source term and my initial post and idea were probably misleading. It's much easier to incorporate a source term in the equations.
My only difficulty was that I didn't know the form of the source term and I couldn't agree with the form given in papers on constant heat flux. Fortunately Patankar paper nails it down. It's really just a bit of calculus Instaed of solving for tempearture you solve for \theta which is \theta = (T  Tw)/(Tbx  Tw) where Tbx is the bulk temperature in the cross section. Moreover you need to define another field lambda which is \lambda = dTbx/dx/(Tbx  Tw) Now you can solve the equations for \theta and \lambda and these quantites happen to be periodic. u d\theta/dx v d\theta/dy = \kappa (d2\theta/dx2 + d2theta/dy2 ) + \sigma Source term \sigma = [2 \kappa d \theta/dx  u \theta] \lambda + \kappa \theta [d lambda/dx + lambda^2] It might look a bit longwinded but it actually makes sense. The solution requires you to iteratively solve for \theta, correct so that bulk \theta is one in each crosssection and than correct \lambda. Currently I am working on the last bit i.e. correcting lambda. But I attach some results with this post. This is interFoam with heat transfer and periodic boundaries, laminar flow. 

December 1, 2011, 11:50 

#6 
Senior Member
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 16 
Thanks to both Robert Sawko & Steven van Haren for this useful topic.
There is also the Steven van Haren M.Sc. thesis that might be helpful: "Testing DNS capability of OpenFOAM and STARCCM+" That you can download in DELFT UNIVERSITY OF TECHNOLOGY: http://www.lr.tudelft.nl/fileadmin/F..._van_Haren.pdf 

May 3, 2012, 06:00 

#7 
New Member
Gabriele G.
Join Date: May 2012
Location: Italy
Posts: 1
Rep Power: 0 
Dear Robert,
I am really interested in this thread since it has been the subject of my M. Sc. Thesis. At present, I would like to implement in OpenFOAM the model proposed by Patankar for the constant wall temperature case (maybe in the modified form of Stalio and Piller, where only the average lambda over one periodic module is computed by using the integral form of the energy conservation equation). Did you succeed in correcting lambda? Can you give me any advice to perform this task avoiding numerical instability? Thanks in advance, Gabriele 

July 18, 2016, 15:27 

#8  
New Member
Dantong Shi
Join Date: Nov 2015
Posts: 6
Rep Power: 10 
Quote:
Hi, Thank you for your sharing with this problem and your progress. I am also doing a similar problem related to periodic boundary for constant wall temperature case. But I don't know how to add the source term and the iteration for lamda. Can you post your modified solver for this or can you give me any suggestions on that? Thank! 

May 9, 2018, 06:55 

#9  
Member
HK
Join Date: Oct 2015
Location: Madras
Posts: 31
Rep Power: 10 
Quote:
Or Modified method by Wang et al. (PRE,2017) ? 

February 4, 2020, 01:39 

#10  
Member
HK
Join Date: Oct 2015
Location: Madras
Posts: 31
Rep Power: 10 
Quote:
We successfully implemented the method proposed by Wang et al. (2017) in our paper. Those who are interested can go through it. 1. Simulation of Fully Developed Flow and Heat Transfer in Wavy Channels Using OpenFOAM (ResearchGate) 2. Unsteady Flow and Heat Transfer Characteristics of Primary and Secondary Corrugated Channels (ResearchGate) 

March 31, 2022, 13:57 

#11 
New Member
Philipp
Join Date: Mar 2022
Posts: 6
Rep Power: 4 
Hello,
i'm currently facing a similar problem as described in this thread. I am trying to simulate a fully developed channel flow with cyclic boundary conditions on inlet and outlet using rhoPimpleFoam and constant heat flux as a wall BC. Since i am relatively new to Openfoam and i never wrote my own solver i was wondering if it is possible to get rhoPimpleFoam to solve the Energy Equation for a variable like ThetaPlus = T  T_w / T_tau since this solver is not using Temperature but Enthalpy in the Energy Equation and the temperature is then calculated later. I just want to get an idea of how complex it would be to do this and was wondering if anyone has probably done something similar. Thanks in advance! 

Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Conjugate heat transfer with periodic boundaries  Suresh  FLUENT  0  February 23, 2009 09:51 
Streamwise periodic heat transfer  Novice  CFX  0  March 4, 2006 08:14 
Convective Heat Transfer  Heat Exchanger  Mark  CFX  6  November 15, 2004 15:55 
Heat conduction through solid rib in periodic duct  Dieter Fauconnier  FLUENT  0  October 7, 2004 08:22 
Conjugate heat transfer (vs) periodic boundary  mp  FLUENT  1  January 13, 2003 01:38 