Programming in OpenFOAM: Adding energy equation Part 1

[@sherlockholmes1832]

What would be my additional term if I had to include radiation mode of heat transfer as well?

[@MrSadra]

its amazing. Thanks a lot for this video

[@1414dmk]

Hi. Thanks for the video. Have you used this solver for simulating mixed convection? something like flow over a heated cylinder? true that ready-made solvers like buoyantBoussinesqPimpleFoam's already available. I have been facing issues with that. please let me know.

[@hyperlyceum60]

Hi Damu
Here, we have focused on solver modification, So, I have never checked in mixed convection as I know that it will not work fine for that case.
You can also share your issues by Email with us.
Peymanhavaej@ hyperlyceum.com

[@antoine1407]

I have a dumb question, adding energy equation can I turn an incompressible solver ( such as InterFoam ) into an compressible one? Thanks for your kidness

[@STWNoman]

Amazing

[@durgeshkumar-ml9nq]

sir, please tell the XE syllabus

[@guijiesang477]

Why the second term on the left side of the energy equation is expressed as fvm::div(phi, T), rather than fvm::div(u, T)?

[@yl96]

I have same question, do you have any idea on it?

[@prakashverma7053]

phi is a scalar field interpolated from vector u field. Basically it is a scalar form of U field.

Invaluable!

[@khaoulalouali1121]

Hello,
I am studying on regenerator which consists of two regions solid and fluid.So, I'm trying to develop new solver which combines two preexisting solvers (icoFoam and chtMultiRegionFoam) for solving a system of equations : Continuity, momentum and energy equations of fluid and solid. if u could help me i will be grateful.
Best Regards

[@francescopignatelli7235]

Hello,
thank you very much for the explanation.
I followed each step of this video and now I am trying to compile the code by typing "wake" but terminal reports me the following error:
Making dependency list for source file thermalPimpleFoam.C
g++ -std=c++11 -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -Wno-attributes -O3 -DNoRepository -ftemplate-depth-100 -I/opt/openfoam7/src/TurbulenceModels/turbulenceModels/lnInclude -I/opt/openfoam7/src/TurbulenceModels/incompressible/lnInclude -I/opt/openfoam7/src/transportModels -I/opt/openfoam7/src/transportModels/incompressible/singlePhaseTransportModel -I/opt/openfoam7/src/finiteVolume/lnInclude -I/opt/openfoam7/src/sampling/lnInclude -I/opt/openfoam7/src/dynamicFvMesh/lnInclude -I/opt/openfoam7/src/dynamicMesh/lnInclude -I/opt/openfoam7/src/meshTools/lnInclude -IlnInclude -I. -I/opt/openfoam7/src/OpenFOAM/lnInclude -I/opt/openfoam7/src/OSspecific/POSIX/lnInclude -fPIC -c thermalPimpleFoam.C -o Make/linux64GccDPInt32Opt/thermalPimpleFoam.o
make: *** No rule to make target 'Make/linux64GccDPInt32Opt///source.dep', needed by 'Make/linux64GccDPInt32Opt///source.o'. Stop.
Can you please help me to fix this issue? I am currently using OpenFOAM 7 if it is of influence.
Thank you very much Peyman.

[@jccasava2035]

Me too!!! By any chance could you solve the issue, Francesco?

[@abderrahmanemokrani2542]

hello sir
please
how i can add a function to the differential equation in openfoam, especially the scale function,
*like here :
fvm::ddt(T) - (1/(1+dfra(T)*Lat/(CP*eT)))*fvm::laplacian(DT, T)
==
fvOptions(T)
*for the continuous functions I managed to define them by this program in createFields.H
volScalarField dfra
(
IOobject
(
"dfra",
runTime.timeName(),
mesh
),
{
0.5-0*T
}
);
*but for the scale function i can't difinde (dfra(T))
*the function what i want add it is this:
int fraction(int n);
using namespace std;
double dfra (double a)
{
double r;
if(a>-0.5 && a

[@cristobalus]

Hi Peyman, thanks for the usseful explanation. Nevetheless, I have problems with wmake command, I recieve "mkdir: cannot create directory ‘/opt/openfoam7/platforms/linux64GccDPInt32Opt/applications’: Permission denied .... ". Now Im using the ubuntu's virtual machine for Windows, Could be that ?

[@hyperlyceum60]

Hi
It could not.
See this page for changing permission in Ubuntu.
help.ubuntu.com/community/FilePermissions
It could be better to make a directory by "mkdir" in your home directory.

[@BoBbY-gj8jz]

It's a great video, learnt alot from it.
please make video on adding reaction in the solver, or modifying combustion solver for other reactions.

[@linabelhadjsenini9085]

I did exactly the same but its not working i always get this message "bash: thermalPimpleFoam : commande introuvable"

[@hyperlyceum60]

Something probably went wrong with the compilation of your new solver.
You can share your own solver through Email. I will compile and test it on my system.
[e-mail]

[@linabelhadjsenini9085]

@@hyperlyceum60 thank you i solved the problem its .c file the should be .C

[@psn5623]

wmake command not working

[@hyperlyceum60]

Please check the following items:
1- Load OpenFOAM package.
2- Are you sure that you in the developed solver directory?
you can check the current directory by "pwd" command and
3- "option" and "file" in the make directory
In case that you still have trouble with the compile of your solver, please contact me by email:
[e-mail]

[@miguelangelampuerosuarez8985]

why do you define twice #include "CorrectPhi.H" and #include "correctPhi.H" and declared in the file (thermalPipeFoam.C)
Here,
*---------------------------------------------------------------------------*/
...
#include "pimpleControl.H"
#include "CorrectPhi.H"
#include "fvOptions.H"
#include "localEulerDdtScheme.H"
#include "fvcSmooth.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
and here,
...
if (correctPhi)
{
// Calculate absolute flux
// from the mapped surface velocity
phi = mesh.Sf() & Uf();
#include "correctPhi.H"
// Make the flux relative to the mesh motion
fvc::makeRelative(phi, U);
}
it has an error,
In file included from thermalPimpleFoam.C:106:0:
correctPhi.H:10:1: error: no matching function for call to ‘CorrectPhi(Foam::volVectorField&, Foam::surfaceScalarField&, Foam::volScalarField&, Foam::dimensionedScalar, Foam::geometricZeroField, Foam::pimpleControl&, bool)’
);
thanks Mrs.Peyman!!

[@hyperlyceum60]

'correctPhi' is a Switch variable while "correctPhi.H" is a file included at beginning of solver.