Programming in OpenFOAM: Adding energy equation Part 1

In this video, you will learn how to define a new variable and implement a new PDE in OpenFOAM. If you need to develop a new CFD solver in OpenFOAM, I would definitely recommend seeing this video.
In the first part of this product, we learn the implementation of the energy equation, while in the second part we going to use and evaluate the developed solver in a planar Couette flow simulation.
This product includes a video file, the OpenFOAM solver, a planar Couette flow test case, and a presentation file.
Please visit our website for more downloading the model files:
www.Hyperlyceum.com
Or the following link:
hyperlyceum.co...
For more information please send me an Email:
[e-mail]

[@MrSadra]

its amazing. Thanks a lot for this video

[@STWNoman]

Amazing

[@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.

[@sherlockholmes1832]

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

Invaluable!

[@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.

[@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

[@durgeshkumar-ml9nq]

sir, please tell the XE syllabus

[@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.

[@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

[@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.