Programming in OpenFOAM: writing a Custom Boundary Condition

In this video, we will show how to implement a new boundary condition step by step in OpenFOAM. If you need to develop a new boundary condition in OpenFOAM, I would definitely recommend seeing this video.
The fluid flow in a rectangular channel (Hagen-Poiseuille flow) is investigated in two different conditions:
1- A constant velocity profile at the inlet
2- A fully developed velocity profile imposed at the inlet of the channel
For the second test case, we have written the parabolic velocity boundary condition employing “codedFixedvalue” in OpenFOAM.

Please visit our website for more downloading test cases and source codes:

Or the following link:

For more information please send me an Email:

[TheOnlyRaceEngineer]

The most perfect explanation ever. Extraordinary. You define each term in very simple terms and that's the highest skill for an engineer. Thank you Hyper Lyceum for your assistance. Videos like this must get above million views and likes. Great job. SUBSCRIBED no doubt.

[luismir]

Great video!
It was really helpful for any OpenFOAM user.
Best regards from Brazil!

[asifraj]

Very nice tutorials. it works perfectly. Thank you.

[ulrichbenash]

hello I am a young enthusiast of fluid mechanics...I am doing a simulation regarding atmospheric flows and I would like to know how to write the C++ code to have the normal speed at the level of a patch and not throughout the domain...I tried everything but nothing I need help please

Nice video, thank you! In version > 8, "redirectType" needs to be replaced by the keyword "name".

[samueltwain]

helpful video, thanks for your effort. after openfoamv1912, there is a new utility called setExprBoundaryFields, I am not sure whether It could reach the same results.

[ramakantgadhewal]

Thank you for your video. I want to apply codedfixedValue boundary conditions in the steady-state solver. Will it work on the steady-state for oscillating inlet?

[DiegoAndrade]

Bravo and thank you, how can you add a 2d function ?

[jacquesenboit]

What's the difference btwn assigning values to the fvPatchField directly, as here, versus initializing a new surface field of the same size & then putting it ob the RHS of 'operator==' as is shown in the OF UG?

[ramkumars]

hi, it is an excellent video on programming.. may i know if similar stuff can be done on interpolation functions by any chance?.. i need one by the way.

[lorenleskovar]

how would you do an open channel flow, what would be the bounary condition for the open surface?

[HM-orso]

I want to apply temperature only at discreetly located coordinates in a 2D plane on heated wall. I don't have any function for temperature profile as it is not continuous but discreet profile. I can not find a way to do that in internet. Can anyone please help me here.

[davidcarnero]

hello everyone, I'm trying to create a time varying boundary condition with the next function: U(x_i, t)=U_0( 1+ k cos(tw + phi)), where w is a oscillating characteristic of the flow, but i´m not sure of how to implement it, I hope someone could help me. Tnks

[hamanisofiane]

Hello, Thanks for your vedeo
I am a beginner on  the use of openFoam. I want to simulate a cyclic open channel flow using interFoam solver.
So I start by creating an other solver called "myinterfoam" in order to modify "fvOpition" term in the UEqn :

fvVectorMatrix UEqn
    (
        fvm::ddt(rho, U) + fvm::div(rhoPhi, U)
      + MRF.DDt(rho, U)
      + turbulence->divDevTau(rho, U)
      ==  fvOptions(rho, U)

I want introduce a condition to the term "fvOptions(rho, U)" by multiplying it by the volume fraction "alpha1" like this:
if (alpha1<1)== 0*fvOptions(rho, U)
else== fvOptions(rho, U)
endif
like this the source term when using Ubar in the file fvOption will be available only when the phase1 is present in the computational domain, means when alpha1=1.
Thanks in advance