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

Conjugate Heat Transfer with Immersed Boundary Method (fe-3.2)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 30, 2016, 09:30
Default Conjugate Heat Transfer with Immersed Boundary Method (fe-3.2)
  #1
New Member
 
Join Date: Oct 2012
Posts: 8
Rep Power: 5
SeMar is on a distinguished road
Dear foamers,

I'm trying to modify the icoIbFoam in the latest foam-extend release 3.2. The idea is, adding a temperature diffusion/convection-equation to the code in order to simulate a conjugate-heat-transfer problem. So far I've implemented a modified "DT" Field
Code:
DT = cellIbMask * lambdaF / (rhoF * CpF) + (1-cellIbMask) * lambdaS / (rhoS * CpS);
All the values xF, xS (fluid, solid) are lookup values from the transportPropertiesDict.

In the solverloop I added
Code:
fvScalarMatrix TEqn
        (
            fvm::ddt(T)
            + fvm::div(phi, T)
            - fvm::laplacian(DT, T)
        );
TEqn.solve();
T.correctBoundaryConditions();
While temperature is being calculated in the fluid cells, the value of the Ib-cells does not change. I tried several combinations of 'setDeadCellValue' and 'fixesValue' in the 0/T BC file and also erased the 'T.correctBoundaryConditions();' without success. Setting a temperature manually in the code however works. So I guess there must be something preventing the immersed boundary cells from being solved. Does anybody have a clue, where to find the corresponding lines of code?

Thanks in advance

Sebastian
SeMar is offline   Reply With Quote

Old   June 22, 2016, 06:12
Talking Have you had success?
  #2
New Member
 
XuDongNa
Join Date: Feb 2012
Location: Harbin Heilongjiang Province,China
Posts: 7
Rep Power: 6
supersoldier is on a distinguished road
Send a message via ICQ to supersoldier
Recently , I also want to do some work like yours, so have you had success?
supersoldier is offline   Reply With Quote

Old   June 22, 2016, 10:09
Default
  #3
New Member
 
Join Date: Oct 2012
Posts: 8
Rep Power: 5
SeMar is on a distinguished road
No, unfortunately not. But I'm still very interested in a solution.
SeMar is offline   Reply With Quote

Old   June 24, 2016, 01:05
Default
  #4
New Member
 
XuDongNa
Join Date: Feb 2012
Location: Harbin Heilongjiang Province,China
Posts: 7
Rep Power: 6
supersoldier is on a distinguished road
Send a message via ICQ to supersoldier
Quote:
Originally Posted by SeMar View Post
Dear foamers,

I'm trying to modify the icoIbFoam in the latest foam-extend release 3.2. The idea is, adding a temperature diffusion/convection-equation to the code in order to simulate a conjugate-heat-transfer problem. So far I've implemented a modified "DT" Field
Code:
DT = cellIbMask * lambdaF / (rhoF * CpF) + (1-cellIbMask) * lambdaS / (rhoS * CpS);
All the values xF, xS (fluid, solid) are lookup values from the transportPropertiesDict.

In the solverloop I added
Code:
fvScalarMatrix TEqn
        (
            fvm::ddt(T)
            + fvm::div(phi, T)
            - fvm::laplacian(DT, T)
        );
TEqn.solve();
T.correctBoundaryConditions();
While temperature is being calculated in the fluid cells, the value of the Ib-cells does not change. I tried several combinations of 'setDeadCellValue' and 'fixesValue' in the 0/T BC file and also erased the 'T.correctBoundaryConditions();' without success. Setting a temperature manually in the code however works. So I guess there must be something preventing the immersed boundary cells from being solved. Does anybody have a clue, where to find the corresponding lines of code?

Thanks in advance

Sebastian
what do you mean by "Setting a temperature manually in the code however works" ? I posted a thread at Immersed Boundary Method in OpenFOAM-3.1-ext . #67.

Can you give me some advance?
supersoldier is offline   Reply With Quote

Old   June 24, 2016, 03:22
Default
  #5
New Member
 
Join Date: Oct 2012
Posts: 8
Rep Power: 5
SeMar is on a distinguished road
By manually I mean assigning a value to the field.
Code:
T = 300.0;
This would result in every cell of field T including immersedboundary cells having a value of 300.0.
SeMar is offline   Reply With Quote

Reply

Tags
foam extend 3.2, heat exchange, immersed boundary, programming

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
Setting the height of the stream in the free channel kevinmccartin CFX 10 July 9, 2015 21:36
Difficulty In Setting Boundary Conditions Moinul Haque CFX 4 November 25, 2014 18:30
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 07:00
what's the type of conjugate heat transfer method in cfx? loose or direct coumpling? skinga281 CFX 8 November 29, 2012 17:07
Domain Imbalance HMR CFX 3 March 6, 2011 21:10


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