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

gradientEnthalpy and wallHeatTransfer

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 9, 2009, 07:31
Default gradientEnthalpy and wallHeatTransfer
  #1
Senior Member
 
Markus Rehm
Join Date: Mar 2009
Location: Erlangen (Germany)
Posts: 178
Rep Power: 8
markusrehm is on a distinguished road
Hello,

I am having problems to use the gradientEnthalpy BC for buoyantSimpleFoam in the hotRoom case at the walls. For 0/T I applied it as follows:

Code:
...
fixedWalls
    {
        type            gradientEnthalpy;
        gradient        uniform -1000.0;
    }
...
It produces a memory error in the hEqn. On the other hand for the wallHeatTransfer BC everything works well.

It doesn't seem to be a problem of the thermos - so much I found out.

In gradientEnthalpyFvPatchScalarField.C:111 Tw.evaluate(); is calling its own updateCoeffs() function again and again.

Any solution to this?

I use the 1.5.x git version. See the output below.

Regards, Markus.

Code:
rem@IEC126:~/OpenFOAM/rem-1.5.x/run/buoyantSimpleFoam/hotRoom> buoyantSimpleFoam
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.5.x                                 |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Exec   : buoyantSimpleFoam
Date   : Jul 09 2009
Time   : 13:28:26
Host   : IEC126
PID    : 529
Case   : /home/rem/OpenFOAM/rem-1.5.x/run/buoyantSimpleFoam/hotRoom
nProcs : 1

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading environmentalProperties
Reading thermophysical properties

Selecting thermodynamics package hThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0.85;
    alphah          1;
    alphak          1;
    alphaEps        0.76923;
}

Calculating field g.h

Creating field pd


Starting time loop

Time = 1

DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 8.16329e-07, No Iterations 6
DILUPBiCG:  Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 8.16329e-07, No Iterations 6
Speicherzugriffsfehler
markusrehm is offline   Reply With Quote

Old   July 12, 2009, 14:05
Default
  #2
Senior Member
 
Henrik Rusche
Join Date: Mar 2009
Location: Braunschweig, Niedersachsen, Germany
Posts: 275
Rep Power: 9
henrik is on a distinguished road
Dear Markus,

The *Enthalpy-boundary conditions are not meant to be used on temperature-fields.

You can understand the internal workings by having a look into $FOAM_SRC/thermophysicalModels/basic/lnInclude/basicThermo.C

Try putting a fixedGradient boundary condition on T.

Henrik
henrik is offline   Reply With Quote

Old   July 13, 2009, 04:06
Default
  #3
Senior Member
 
Markus Rehm
Join Date: Mar 2009
Location: Erlangen (Germany)
Posts: 178
Rep Power: 8
markusrehm is on a distinguished road
Thanks Henrik,

I will have a deeper look. Markus
markusrehm is offline   Reply With Quote

Reply

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



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