[alexeym]

Hi,

You are trying to install software system-wide while being simple user. You should either use 'sudo wmake' (though in this case I am not quite sure about environment variables being set right), or change $(FOAM_APPBIN) to $(FOAM_USER_APPBIN) in Make/files file to install software in user's folder.

[student666]

Thank you very much for helping; I change Make/file by changing into FOAM_USER_APPBIN.

Regards.

[Scram_1]

Hi,
I'm new to openFOAM and thought of checking if OF gives the right value of the Nu for laminar flow through a pipe.
I used the buoyantBoussinesqSimpleFoam solver by turning off turbulence and g (so that the buoyancy effects are turned off). The velocity and temperature fields seem correct. Now, I want to check the Nu. I tried creating the function NusseltCalc() as detailed in one of your previous posts. However during compilation, I'm getting an error. I tried compiling the wallHeatFlux function as well just to check if it compiles properly (obviously by creating a copy of it first). Even that is giving the same error. In fact even when I tried adding temperature field to icoFoam and tried compiling, its giving me the same error.

Here's the error message I'm getting

csriram91@SRIRAM-PC /opt/OpenFOAM/OpenFOAM-2.3.x/applications/utilities/postProcessing/wall/wallHeatFlux - Copy (2)
$ wclean

csriram91@SRIRAM-PC /opt/OpenFOAM/OpenFOAM-2.3.x/applications/utilities/postProcessing/wall/wallHeatFlux - Copy (2)
$ wmake
Making dependency files...
In file included from /opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/faceZoneMesh.H:33:0,
from /opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/polyMesh.H:56,
from /opt/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/fvMesh.H:50,
from /opt/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/fvCFD.H:7,
from wallHeatFlux.C:34:
/opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/faceZone.H:46:32: fatal error: primitiveFacePatch.H: No such file or directory
#include "primitiveFacePatch.H"
^
compilation terminated.
SOURCE=wallHeatFlux.C ; x86_64-w64-mingw32-g++ -m64 -Dcygwin64 -DWM_DP -DMSWIN -DLITTLE_ENDIAN -Wall -Wextra -Wno-unused-parameter -Wnon-virtual-dtor -O3 -DNDEBUG -DNoRepository -ftemplate-depth-100 -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/turbulenceModels -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/turbulenceModels/compressible/turbulenceModel/lnInclude -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/thermophysicalModels/specie/lnInclude -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/thermophysicalModels/reactionThermo/lnInclude -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/thermophysicalModels/basic/lnInclude -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/thermophysicalModels/solidThermo/lnInclude -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/OSspecific/MSWindows/lnInclude -c $SOURCE -o Make/cygwin64mingw-w64DPOpt/wallHeatFlux.o
In file included from /opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/faceZoneMesh.H:33:0,
from /opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/polyMesh.H:56,
from /opt/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/fvMesh.H:50,
from /opt/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/fvCFD.H:7,
from wallHeatFlux.C:34:
/opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/faceZone.H:46:32: fatal error: primitiveFacePatch.H: No such file or directory
#include "primitiveFacePatch.H"
^
compilation terminated.
wallHeatFlux.dep:3: recipe for target 'Make/cygwin64mingw-w64DPOpt/wallHeatFlux.o' failed
make: *** [Make/cygwin64mingw-w64DPOpt/wallHeatFlux.o] Error 1


I'm running openFOAM 2.3x on cygwin if that's any help.

Thanks!!
Ram

[djamila77]

hi; all

thank you for the answer,

I'll try your suggestions and I'll share the results with you.

Cordially

[djamila77]

Originally Posted by giovanni10

What does (T_hot-T_ini) stand for ?
-- wall_localNusselt.boundaryField()[patchi] = length*patchHeatFlux[patchi]/((T_hot-T_ini)*k); --

I want to calculate the localNusselt number at a wall which is called "right". I have a rectangular cavity which has 4 walls with the following names: bottom, right, left and top. So, the Nusselt number is calculated as the integration of dT/dx for limits from y=0 to 1, at x=0, which is the "right" wall.
How can I alter your code?
I am using an incompressible Magnetohydrodynamics` solver.
Thanks a lot for your help!

Dear Sir
I don't know exactly what your kind of problem research but the Code Fluent is useful in kind of convection in a box headed by the below wall and cladded by the top
you can use it for calculate the Nusselt number and other thing

[acs]

Originally Posted by Scram_1

Hi,
I'm new to openFOAM and thought of checking if OF gives the right value of the Nu for laminar flow through a pipe.
I used the buoyantBoussinesqSimpleFoam solver by turning off turbulence and g (so that the buoyancy effects are turned off). The velocity and temperature fields seem correct. Now, I want to check the Nu. I tried creating the function NusseltCalc() as detailed in one of your previous posts. However during compilation, I'm getting an error. I tried compiling the wallHeatFlux function as well just to check if it compiles properly (obviously by creating a copy of it first). Even that is giving the same error. In fact even when I tried adding temperature field to icoFoam and tried compiling, its giving me the same error.

Here's the error message I'm getting

csriram91@SRIRAM-PC /opt/OpenFOAM/OpenFOAM-2.3.x/applications/utilities/postProcessing/wall/wallHeatFlux - Copy (2)
$ wclean

csriram91@SRIRAM-PC /opt/OpenFOAM/OpenFOAM-2.3.x/applications/utilities/postProcessing/wall/wallHeatFlux - Copy (2)
$ wmake
Making dependency files...
In file included from /opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/faceZoneMesh.H:33:0,
from /opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/polyMesh.H:56,
from /opt/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/fvMesh.H:50,
from /opt/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/fvCFD.H:7,
from wallHeatFlux.C:34:
/opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/faceZone.H:46:32: fatal error: primitiveFacePatch.H: No such file or directory
#include "primitiveFacePatch.H"
^
compilation terminated.
SOURCE=wallHeatFlux.C ; x86_64-w64-mingw32-g++ -m64 -Dcygwin64 -DWM_DP -DMSWIN -DLITTLE_ENDIAN -Wall -Wextra -Wno-unused-parameter -Wnon-virtual-dtor -O3 -DNDEBUG -DNoRepository -ftemplate-depth-100 -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/turbulenceModels -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/turbulenceModels/compressible/turbulenceModel/lnInclude -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/thermophysicalModels/specie/lnInclude -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/thermophysicalModels/reactionThermo/lnInclude -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/thermophysicalModels/basic/lnInclude -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/thermophysicalModels/solidThermo/lnInclude -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude -I/opt/OpenFOAM/OpenFOAM-2.3.x/src/OSspecific/MSWindows/lnInclude -c $SOURCE -o Make/cygwin64mingw-w64DPOpt/wallHeatFlux.o
In file included from /opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/faceZoneMesh.H:33:0,
from /opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/polyMesh.H:56,
from /opt/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/fvMesh.H:50,
from /opt/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/fvCFD.H:7,
from wallHeatFlux.C:34:
/opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/faceZone.H:46:32: fatal error: primitiveFacePatch.H: No such file or directory
#include "primitiveFacePatch.H"
^
compilation terminated.
wallHeatFlux.dep:3: recipe for target 'Make/cygwin64mingw-w64DPOpt/wallHeatFlux.o' failed
make: *** [Make/cygwin64mingw-w64DPOpt/wallHeatFlux.o] Error 1


I'm running openFOAM 2.3x on cygwin if that's any help.

Thanks!!
Ram

Hi, Ram, I've got the same problem. Have you solved it? What did you do?

Thanx!

[djamila77]

Hi there
I'd like asking about the Lenght in Report - Reference Values.
is it the Depth that means Hydraulic diameter for 3D or note

[djamila77]

Hi there
I'd like add this program to calculate Nusselt number, but I don't know how to do that
please help me

[djamila77]

Originally Posted by Goutam

The code was correct. NusseltCalc is working now.

Hi
how can I include the program
please with Details

[pille]

If someone wants to implement the code of Johan Magnusson (shown in post #3) in OpenFOAM 4.1, one only need to change

 
forAll(wall_localNusselt.boundaryField(),patchi
{    
     wall_localNusselt.boundaryField()[patchi] = length*
               patchHeatFlux[patchi]/((T_hot-T_ini)*k);     
}

to

 
forAll(wall_localNusselt.boundaryFieldRef(),patchi
{    
     wall_localNusselt.boundaryFieldRef()[patchi] = length*
               patchHeatFlux[patchi]/((T_hot-T_ini)*k);     
}

and everything works fine.
A detailed description of the implementation is given in his report
("conjugateHeatFoam with explanational tutorial together with a buoyancy driven flow tutorial and a convective conductive tutorial")

[MFN]

Originally Posted by pille

If someone wants to implement the code of Johan Magnusson (shown in post #3) in OpenFOAM 4.1, one only need to change

 
forAll(wall_localNusselt.boundaryField(),patchi
{    
     wall_localNusselt.boundaryField()[patchi] = length*
               patchHeatFlux[patchi]/((T_hot-T_ini)*k);     
}

to

 
forAll(wall_localNusselt.boundaryFieldRef(),patchi
{    
     wall_localNusselt.boundaryFieldRef()[patchi] = length*
               patchHeatFlux[patchi]/((T_hot-T_ini)*k);     
}

and everything works fine.
A detailed description of the implementation is given in his report
("conjugateHeatFoam with explanational tutorial together with a buoyancy driven flow tutorial and a convective conductive tutorial")

does it work in openfoam 4.x ?
I will test and re check that

[Hgholami]

Hi
I modify icoFluid solver of fsiFoam and add this code to it as:

wallHeatFlux_
(
IOobject
(
"wallHeatFlux",
runTime().timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
dimensionedScalar("wallHeatFlux",dimMass/pow3(dimTime), 0.0)
),

forAll(wallHeatFlux_.boundaryField(),patchi)
{
if (mesh.boundary()[patchi].isWall())
{
wallHeatFlux_.boundaryField()[patchi] = (DT_*rho_*Cp_)*fvc::snGrad(T_);
}
}

to calculate heat flux during the simulation. The problem is that the output file "wallheatFlux" in each timeStep is different size
for example: "size 10441 is not equal to the given value of 117"
the code make same size for each boundary surfaces (for inlet, top,bottom etc with different size, the wallHeatFlux is same size 10441). Do you know what is the mistake?

[Daniel_Khazaei]

Originally Posted by Hgholami

Hi
I modify icoFluid solver of fsiFoam and add this code to it as:

to calculate heat flux during the simulation. The problem is that the output file "wallheatFlux" in each timeStep is different size
for example: "size 10441 is not equal to the given value of 117"
the code make same size for each boundary surfaces (for inlet, top,bottom etc with different size, the wallHeatFlux is same size 10441). Do you know what is the mistake?

Assuming (DT_*rho_*Cp_) results in a dimensionedScalar, what you need is surface normal gradient of temperature at boundaries:

wallHeatFlux_.boundaryField()[patchi] = (DT_*rho_*Cp_)*T_.boundaryField()[patchi].snGrad();

[Hgholami]

Originally Posted by Goutam

Dear Members,

This is the code for local Nusselt number:

volScalarField wall_localNusselt
(
IOobject
(
"wall_localNusselt",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
dimensionedScalar("wall_localNusselt ",dimless,0.0)
);

forAll(wall_localNusselt.boundaryField(),patchi)
{
wall_localNusselt.boundaryField()[patchi] = length*patchHeatFlux[patchi]/((T_hot-T_ini)*k);
}

wall_localNusselt.write();
}

I want to write a code for Average Nusselt Number.

Dear Goutam
Do you have any code for calculation Nusselt number with bulk temperature?

I can't use swak4Foam, in order, I should calculate bulk temperature and add it to Nusselt calculation

[chberbelali]

hello everyone,

is there anyone who can help me?!

How can i calculate the local Nusselt number for a vertical wall in cavity using a Matlab implimentation?

[piu58]

The Nußelt number is connected to the thickness of the boundary layer d. If you may measure this thickness, you may calculate the Nußelt number stepwise (but not pointwise). You have to select al length along L the surface for which you plan getting Nu. The you have

Nu = L / d

[chberbelali]

Originally Posted by piu58

The Nußelt number is connected to the thickness of the boundary layer d. If you may measure this thickness, you may calculate the Nußelt number stepwise (but not pointwise). You have to select al length along L the surface for which you plan getting Nu. The you have

Nu = L / d

Details here:

Nusselt number

Matlab: Local nusselt number computing

My problem figure:


Best regards