CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

conjugate heat transfer in OpenFOAM

Register Blogs Community New Posts Updated Threads Search

Like Tree28Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 22, 2013, 10:44
Default conjugate heat transfer in OpenFOAM
  #1
Member
 
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12
skuznet is on a distinguished road
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!
awacs, ivd, ordinary and 2 others like this.
skuznet is offline   Reply With Quote

Old   November 22, 2013, 16:30
Default
  #2
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
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
If there is no such direct conversion, have a look into this wiki page and the links therein: http://openfoamwiki.net/index.php/SnappyHexMesh

Best regards,
Bruno
tmik likes this.
__________________
wyldckat is offline   Reply With Quote

Old   November 22, 2013, 19:12
Default
  #3
Member
 
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12
skuznet is on a distinguished road
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.
skuznet is offline   Reply With Quote

Old   November 23, 2013, 07:29
Default
  #4
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
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
__________________
wyldckat is offline   Reply With Quote

Old   November 23, 2013, 12:14
Default
  #5
Member
 
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12
skuznet is on a distinguished road
Awesome! Thank you Bruno!
skuznet is offline   Reply With Quote

Old   November 23, 2013, 18:38
Default
  #6
Member
 
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12
skuznet is on a distinguished road
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
I can get rid off this message by choosing Case Type: Decomposed Case
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.
skuznet is offline   Reply With Quote

Old   November 24, 2013, 07:58
Default
  #7
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
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
Then you can open the case with ParaView's internal reader.

Best regards,
Bruno
__________________
wyldckat is offline   Reply With Quote

Old   December 4, 2013, 23:49
Default
  #8
Member
 
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12
skuznet is on a distinguished road
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.
wyldckat likes this.
skuznet is offline   Reply With Quote

Old   December 10, 2013, 12:36
Default
  #9
Member
 
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12
skuznet is on a distinguished road
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!
skuznet is offline   Reply With Quote

Old   December 10, 2013, 16:55
Default
  #10
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Hi Sergey,

A few questions:
  1. Which specific multi-region solver are you using?
  2. What is the flow speed for each fluid?
  3. Have you checked the solid's characteristics? It might be bleeding energy due to very low or high heat transfer factors.
Best regards,
Bruno
__________________
wyldckat is offline   Reply With Quote

Old   December 10, 2013, 18:53
Default
  #11
Member
 
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12
skuznet is on a distinguished road
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
Kummi likes this.

Last edited by skuznet; December 11, 2013 at 17:31.
skuznet is offline   Reply With Quote

Old   December 10, 2013, 23:34
Default
  #12
Member
 
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12
skuznet is on a distinguished road
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?
skuznet is offline   Reply With Quote

Old   December 15, 2013, 14:28
Default
  #13
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Hi Sergey,

Quote:
Originally Posted by skuznet View Post
2. I have only one fluid - water, the inlet velocity is 0.1 m/s, while the size of the fin 1x1x3 m.
The fin is quite big. It's easy for it to either cool down too much or warm up too much, depending on the thermodynamic properties defined for it.

Quote:
Originally Posted by skuznet View Post
3. I didn't check the solid. How can I check?
Check the file "constant/wall/thermophysicalProperties".

Quote:
Originally Posted by skuznet View Post
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
Then there is something very wrong with the boundary conditions.

Quote:
Originally Posted by skuznet View Post
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?
Which specific file are you referring to? Because the file "constant/thermophysicalProperties" is not used by these multi-region solvers!

Best regards,
Bruno
__________________
wyldckat is offline   Reply With Quote

Old   December 16, 2013, 15:40
Default
  #14
Member
 
Sergey
Join Date: Nov 2013
Posts: 87
Rep Power: 12
skuznet is on a distinguished road
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!
Attached Files
File Type: gz fin2.tar.gz (7.2 KB, 66 views)
skuznet is offline   Reply With Quote

Old   September 18, 2016, 11:23
Default
  #15
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 9
sitajeje is on a distinguished road
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
sitajeje is offline   Reply With Quote

Old   September 18, 2016, 15:50
Default
  #16
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Quick answers: First of all, sorry to everyone, but:
  1. I haven't had enough free time to look into your questions and I won't be able to answer to them all on this post right now.
  2. I can't find it right now, but I know I wrote a post... last year, I think... in which I stated that the tutorial was not working properly as of OpenFOAM 2.3. This situation still stands, as I haven't managed to find enough time to study things properly.
The best I can suggest right now is suggest what I suggested to another forum member a few days ago:
Quote:
Originally Posted by wyldckat View Post
so the best I can do is suggest that you study the course presentation and cases that are provided here: http://www.openfoamworkshop.org/courses - look for the presentation "Learning how to use free surface flows in OpenFOAM 3.0" and click on the download link to the right of the entry "Training Material".
The free surface flows are also using the same kind of pressure management as in the chtMultiRegion*Foam solvers.



Now, to answer in detail to a few questions:

Quote:
Originally Posted by jaydeep View Post
I guess I could narrow down the errors in ParaView to initial condition.
Yes, the internal reader in ParaView does not support all boundary conditions that OpenFOAM supports at the moment. Therefore, ParaView will complain when it doesn't understand a boundary condition


Quote:
Originally Posted by jaydeep View Post
P.S. Still the results look unreasonable
It's as I indicated at the start of this post that you're reading.

Quote:
Originally Posted by sitajeje View Post
But the resultant planeWall2D.OpenFOAM file is empty. The Paraview exit automatically when I press on "apply" button to load the data.
Yes, the file is meant to be empty, it's just sort-of a tag for ParaView to know which folder to use.
And it looks like you did not properly ready the wiki page. This is explained here: http://openfoamwiki.net/index.php/Ge...ost-Processing
__________________
wyldckat is offline   Reply With Quote

Old   September 19, 2016, 14:42
Default
  #17
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 9
sitajeje is on a distinguished road
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
sitajeje is offline   Reply With Quote

Old   September 20, 2016, 07:36
Default
  #18
Senior Member
 
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 22
vasava will become famous soon enough
Quote:
Originally Posted by sitajeje View Post
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..
You dont need separate files for separate regions to view results, the file planeWallD.OpenFOAM should be enough.
vasava is offline   Reply With Quote

Old   September 21, 2016, 11:13
Default
  #19
Senior Member
 
Alex
Join Date: Oct 2013
Posts: 337
Rep Power: 21
zfaraday will become famous soon enough
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
Regards,

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!
zfaraday is offline   Reply With Quote

Old   September 21, 2016, 13:56
Default
  #20
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 9
sitajeje is on a distinguished road
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
sitajeje is offline   Reply With Quote

Reply

Tags
cht, solid-fluid interface


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 Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 00:38.