|
[Sponsors] |
|
November 22, 2013, 10:44 |
conjugate heat transfer in OpenFOAM
|
#1 |
Member
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12 |
Dear FOAMers
I'm new in OpenFoam and would like to model steady state conjugate heat transfer(solid and fluid region) in a heat exchanger which consists of a few vertical metal plates attached to a horizontal plate with fluid flowing in gaps between plates. My current experience with OpenFoam is restricted by “Lid-driven cavity flow tutorial” and “Stress analysis of a hole” tutorials. I know that there are “chtMultiRegionFoam” and “chtMultiRegionSimpleFoam” cases in /tutorials/heatTransfer/ folder. However, i didn't findt tutorials which would explain these cases and it's hard to understand what is going on in these cases without explanation. I wonder if there is a step-by-step tutorial similar to “Lid-driven cavity” tutorial explaining conjugate heat transfer. I also wonder if there are tutorials explaining - how to create a case containing regions with different properties(fluid and solid), assign regions and boundaries; - how to use conjugate heat transfer solvers; My next step will be to model a steady state conjugate heat transfer in a porous material consisting of orthogonally interlaced cylindrical wires and fluid between them. I have a geometry for this case created in Abaqus. I wonder if there is any way to translate my geometry from Abaqus to Open Foam preserving surfaces and regions definitions. Or what would be appropriate tool consistent with Open FOAM where I can create such geometry and define regions and surfaces? Thank you! |
|
November 22, 2013, 16:30 |
|
#2 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Greetings skuznet and welcome to the forum!
Have a look into this (unofficial) tutorial: http://openfoamwiki.net/index.php/Ge..._-_planeWall2D - and feel free to add information to it, when you've understood better how things work! As for Abaqus, I have no idea if there is a direct way. You can try looking here in the forum. You can even use Google: Code:
site:www.cfd-online.com/Forums/ Abaqus OpenFOAM Best regards, Bruno
__________________
|
|
November 22, 2013, 19:12 |
|
#3 |
Member
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12 |
Bruno,
thank you for your answer. Actually I've run already the planeWall2D case. It didn't work. There was an error message in log.chtMultiRegionSipleFOAM Code:
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.2.2 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.2.2-9739c53ec43f Exec : chtMultiRegionSimpleFoam Date : Nov 22 2013 Time : 19:54:37 Host : "scyld.localdomain" PID : 65269 Case : /home/skuznetsov/OpenFOAM/RUN/planeWall2D nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time --> FOAM FATAL IO ERROR: keyword regions is undefined in dictionary "/home/skuznetsov/OpenFOAM/RUN/planeWall2D/constant/regionProperties" file: /home/skuznetsov/OpenFOAM/RUN/planeWall2D/constant/regionProperties from line 18 to line 20. From function dictionary::lookupEntry(const word&, bool, bool) const in file db/dictionary/dictionary.C at line 402. FOAM exiting and I also had error message when was trying to open results with paraview Code:
ERROR: In ..\..\..\src\VTK\IO\vtkOpenFOAMReader.cxx, line 7200 vtkOpenFOAMReaderPrivate (058864C8): Error reading line 58 of E:\Dropbox\planeWall2D\constant/bottomAir/polyMesh/boundary: Expected number, string or (, found ; ERROR: In ..\..\..\src\VTK\IO\vtkOpenFOAMReader.cxx, line 7200 vtkOpenFOAMReaderPrivate (058865E0): Error reading line 58 of E:\Dropbox\planeWall2D\constant/topAir/polyMesh/boundary: Expected number, string or (, found ; ERROR: In ..\..\..\src\VTK\IO\vtkOpenFOAMReader.cxx, line 7200 vtkOpenFOAMReaderPrivate (058866F8): Error reading line 47 of E:\Dropbox\planeWall2D\constant/wall/polyMesh/boundary: Expected number, string or (, found ; ERROR: In ..\..\..\src\VTK\Filtering\vtkExecutive.cxx, line 756 vtkCompositeDataPipeline (05BAEC40): Algorithm vtkPOpenFOAMReader(05AD8AA8) returned failure for request: vtkInformation (05B304F0) Debug: Off Modified Time: 95148 Reference Count: 1 Registered Events: (none) Request: REQUEST_INFORMATION FORWARD_DIRECTION: 0 ALGORITHM_AFTER_FORWARD: 1 Last edited by skuznet; November 25, 2013 at 15:46. |
|
November 23, 2013, 07:29 |
|
#4 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Hi Sergey,
Have a look into the wiki page again: http://openfoamwiki.net/index.php/Ge..._-_planeWall2D - I've added an adapted case for OpenFOAM 2.2. Best regards, Bruno
__________________
|
|
November 23, 2013, 12:14 |
|
#5 |
Member
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12 |
Awesome! Thank you Bruno!
|
|
November 23, 2013, 18:38 |
|
#6 |
Member
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12 |
Bruno,
the case now seems to run well. But when I'm trying to open it paraview I still can't see anything. I downloaded all case files on my pc and trying to open it with paraview (i tried versions 3.81, 3.12, 3.98.0). -I choose .OpenFOAM file to open, -choose OpenFOAM reader, -check Mesh Regions box -and click Apply and nothing appears. I see message Code:
ERROR: In ..\..\..\src\VTK\IO\vtkOpenFOAMReader.cxx, line 4498 vtkOpenFOAMReaderPrivate (05944570): Error opening E:\Dropbox\planeWall2D\constant/polyMesh/faces.gz: Can't open. If you are trying to read a parallel decomposed case, set Case Type to Decomposed Case. ERROR: In ..\..\..\src\VTK\Filtering\vtkExecutive.cxx, line 756 vtkCompositeDataPipeline (05CB2EA0): Algorithm vtkPOpenFOAMReader(05B9E960) returned failure for request: vtkInformation (05BE93C8) Debug: Off Modified Time: 98373 Reference Count: 1 Registered Events: (none) Request: REQUEST_DATA FROM_OUTPUT_PORT: 0 FORWARD_DIRECTION: 0 ALGORITHM_AFTER_FORWARD: 1 but still don't see anything. The time panel is not active, I can't change time and in Color panel I can't choose anything except Solid Color. If I open the cavity case, then I can see geometry and pressure and velocity distributions. Last edited by skuznet; November 25, 2013 at 15:44. |
|
November 24, 2013, 07:58 |
|
#7 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Hi Sergey,
Please, when posting code and error outputs, please use the "[CODE]" markers, as explained in the second link of my signature. And you can edit your previous posts as well, so please add this marker to them, to make it easier to read the posts. Regarding your question: the "*.gz" files might not be opened with success with ParaView's internal ".foam" reader. First unzip the files; on Linux, you can run this command: Code:
find . -name "*.gz" | xargs gunzip Best regards, Bruno
__________________
|
|
December 4, 2013, 23:49 |
|
#8 |
Member
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12 |
Bruno, thank you for the help!
I installed ubuntu, openfoan and paraview on my homecomputer and now it works well. I can run case planeWall2d and multiRegionHeatre and can open results. I think my previous issues with opening results were because I was opening them in Windows and files were not transferred or opened correctly in paraview under windows. |
|
December 10, 2013, 12:36 |
|
#9 |
Member
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12 |
Dear Bruno:
Thank you for adapting the planeWall2D. I was able to run this case and I got some basic understanding how it works. Based on planeWall2D tutorial I was able to set up and run a simple 3d fin case - there is a fin which is heated from the bottom and the fluid flow removes the heat from the fin. I was able to run it with air and it looks ok. However, when I try to put water instead of air by replacing themoPhysicalProperties with file from HTML Code:
/chtMultiRegionFoam/multiRegionLiquidHeater/constant/bottomWater/thermophysicalProperties Code:
\*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object thermophysicalProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // thermoType { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } mixture { specie { nMoles 1; molWeight 18; } equationOfState { rho 1000; } thermodynamics { Cp 4181; Hf 0; } transport { mu 959e-6; Pr 6.62; } } // ************************************************************************* // Something is going wrong. Open FOAM starts solving it but during iterations my temperature (minimal temperature in solid and in fluid domains) is reducing, which should not occur, since I have a heater. At the end temperature in fluid domain becomes even negative, and OpenFOAM crashes. Can you please tell me if I have to change anything else in my case when I replaced my thermophysicalProperties file with water? I guess the problem might be because thermophysicalProperties for water has different equation of state - rhoConst instead of perfectGas. Thank you! |
|
December 10, 2013, 16:55 |
|
#10 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Hi Sergey,
A few questions:
Bruno
__________________
|
|
December 10, 2013, 18:53 |
|
#11 |
Member
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12 |
1. I'm using chtMultiRegionSimpleFoam solver, but the thermophysicalProperties file was taken from the example for chtMultiRegionFoam.
2. I have only one fluid - water, the inlet velocity is 0.1 m/s, while the size of the fin 1x1x3 m. 3. I didn't check the solid. How can I check? The very first fluid iteration gives also very high temperature 2000, while initial temperature was 300 and temperature of the warmed surface of solid is 350: Code:
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.2.x | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.2.x-37940e0cd1ae Exec : /home/sergkuznet/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam Date : Dec 11 2013 Time : 16:43:49 Host : "ubuntu" PID : 8592 Case : /home/sergkuznet/OpenFOAM/run/fin2R_initial nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create fluid mesh for region fluidDomain for time = 0 Create solid mesh for region solidDomain for time = 0 *** Reading fluid mesh thermophysical properties for region fluidDomain Adding to thermoFluid Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } Adding to rhoFluid Adding to UFluid Adding to phiFluid Adding to gFluid Adding to turbulence Selecting turbulence model type laminar Adding to ghFluid Adding to ghfFluid Radiation model not active: radiationProperties not found Selecting radiationModel none Adding fvOptions No finite volume options present *** Reading solid mesh thermophysical properties for region solidDomain Adding to thermos Selecting thermodynamics package { type heSolidThermo; mixture pureMixture; transport constIso; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } Adding to radiations Radiation model not active: radiationProperties not found Selecting radiationModel none Adding fvOptions No finite volume options present Time = 1 Solving for fluid region fluidDomain DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.041964, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.0001767454, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.0003058668, No Iterations 2 DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.02843725, No Iterations 1 Min/max T:299.9999 2035.089 GAMG: Solving for p_rgh, Initial residual = 0.8476081, Final residual = 0.003588793, No Iterations 11 time step continuity errors : sum local = 0.7781759, global = -0.07076251, cumulative = -0.07076251 Min/max rho:2 2 Solving for solid region solidDomain DICPCG: Solving for h, Initial residual = 1, Final residual = 0.05509895, No Iterations 1 Min/max T:min(T) [0 0 0 1 0 0 0] 300 max(T) [0 0 0 1 0 0 0] 350 ExecutionTime = 1.39 s ClockTime = 1 s Time = 2 Solving for fluid region fluidDomain DILUPBiCG: Solving for Ux, Initial residual = 0.3551669, Final residual = 0.003849144, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.05855359, Final residual = 0.001012355, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.02422332, Final residual = 0.0003484752, No Iterations 1 DILUPBiCG: Solving for h, Initial residual = 0.5433119, Final residual = 0.007134343, No Iterations 1 Min/max T:264.6221 2035.089 GAMG: Solving for p_rgh, Initial residual = 0.04481994, Final residual = 0.0004256873, No Iterations 18 time step continuity errors : sum local = 3.010856, global = 1.085556, cumulative = 1.014794 Min/max rho:2 2 Solving for solid region solidDomain DICPCG: Solving for h, Initial residual = 0.189656, Final residual = 0.014255, No Iterations 1 Min/max T:min(T) [0 0 0 1 0 0 0] 299.7434 max(T) [0 0 0 1 0 0 0] 350 ExecutionTime = 1.92 s ClockTime = 2 s Time = 3 Solving for fluid region fluidDomain DILUPBiCG: Solving for Ux, Initial residual = 0.1454075, Final residual = 0.001948315, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.09988127, Final residual = 0.001433494, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.0186658, Final residual = 0.0001720673, No Iterations 1 DILUPBiCG: Solving for h, Initial residual = 0.3185567, Final residual = 0.01199445, No Iterations 1 Min/max T:244.3615 1008.353 GAMG: Solving for p_rgh, Initial residual = 0.2012721, Final residual = 0.001699132, No Iterations 18 time step continuity errors : sum local = 0.8901295, global = 0.2734903, cumulative = 1.288284 Min/max rho:2 2 Solving for solid region solidDomain DICPCG: Solving for h, Initial residual = 0.09596577, Final residual = 0.00727127, No Iterations 1 Min/max T:min(T) [0 0 0 1 0 0 0] 299.396 max(T) [0 0 0 1 0 0 0] 350 ExecutionTime = 2.44 s ClockTime = 2 s Last edited by skuznet; December 11, 2013 at 17:31. |
|
December 10, 2013, 23:34 |
|
#12 |
Member
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12 |
Bruno,
I can see some inconsistency in what i'm trying to do: the chtMultiRegionSimpleFoam in compressible solver, and in thermophysicalProperties which I used the density is constant. Can it be the reason for temperature decrease and OF crash? |
|
December 15, 2013, 14:28 |
|
#13 | |||
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Hi Sergey,
Quote:
Check the file "constant/wall/thermophysicalProperties". Quote:
Quote:
Best regards, Bruno
__________________
|
||||
December 16, 2013, 15:40 |
|
#14 |
Member
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12 |
Bruno,
I took the file "constant/solidDomain/thermophysicalProperties" from multiRegionHeater tutorial, and it looks OK to me. I was reffering to the file "constant/fluidDomain/thermophysicalProperties", which describes my fluid properties. I can't find out what is wrong with this case. Could you please take a look at it? Thank you! |
|
September 18, 2016, 11:23 |
|
#15 |
Member
Join Date: Sep 2016
Posts: 63
Rep Power: 9 |
Dear Bruno,
thank you very much first of all for your planeWall2D. I downloaded your case files and run it. But the resultant planeWall2D.OpenFOAM file is empty. The Paraview exit automatically when I press on "apply" button to load the data. My OpenFOAM is V4.0. I would appreciate greatly if you can help me on this, although it is a 3 year old post. Best regards, Hai |
|
September 18, 2016, 15:50 |
|
#16 | |||
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Quick answers: First of all, sorry to everyone, but:
Quote:
Now, to answer in detail to a few questions: Quote:
It's as I indicated at the start of this post that you're reading. Quote:
And it looks like you did not properly ready the wiki page. This is explained here: http://openfoamwiki.net/index.php/Ge...ost-Processing
__________________
|
||||
September 19, 2016, 14:42 |
|
#17 |
Member
Join Date: Sep 2016
Posts: 63
Rep Power: 9 |
Dear Bruno,
thank you very much for your prompt reply! I am a beginner for OpenFOAM and can easily make stupid mistakes. I wonder whether I have done anything wrong. After I executed Allrun, I didn't get planeWall2D{bottomAir}.OpenFOAM planeWall2D{topAir}.OpenFOAM planeWall2D{wall}.OpenFOAM but only planeWallD.OpenFOAM and planeWallD.blockMesh. The other 7 files begin with log.. I use Ubundu 14.04, and OpenFOAM 4.0. I would appreciate greatly if you can give me a guide. Thank you very much in advance! Have a nice evening! Sitajeje |
|
September 20, 2016, 07:36 |
|
#18 |
Senior Member
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 22 |
You dont need separate files for separate regions to view results, the file planeWallD.OpenFOAM should be enough.
|
|
September 21, 2016, 11:13 |
|
#19 |
Senior Member
Alex
Join Date: Oct 2013
Posts: 337
Rep Power: 21 |
I think you are wrong Paritosh, actually you need those files in order to be able to see the results in each region, otherwise you can't watch the results!
In order to create them, the command to be executed is Code:
paraFoam -touchAll Alex
__________________
Web site where I present my Master's Thesis: foamingtime.wordpress.com The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. Some basic tutorials are also resolved step by step in the web. If you are interested in these matters, you are invited to come in! |
|
September 21, 2016, 13:56 |
|
#20 |
Member
Join Date: Sep 2016
Posts: 63
Rep Power: 9 |
Dear Alex and Paritosh,
thank you very much for your help! I wonder whether you have the planeWall2D case file which you have run through? I would appreciate if you would share it with me. I still cannot get the planeWall2D{bottomAir}.OpenFOAM, planeWall2D{topAir}.OpenFOAM and planeWall2D{wall}.OpenFOAM files. All the best sitajeje |
|
Tags |
cht, solid-fluid interface |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Conjugate Heat Transfer: Wall Heat Flux at Coupled Walls? | MaxHeat | FLUENT | 4 | September 14, 2017 10:44 |
openfoam for heat transfer | kirankarki | OpenFOAM Running, Solving & CFD | 29 | February 12, 2015 18:46 |
Conjugate Heat transfer in CFX | ksp1717 | CFX | 11 | December 10, 2010 22:07 |
Conjugate Heat Transfer of Motorized EGR | enr_venkat | CFX | 1 | October 12, 2010 18:17 |
best mesh generator for conjugate heat transfer? | phsieh2005 | Main CFD Forum | 1 | June 1, 2007 17:35 |