Add source term in icoReactingMultiphaseInterFoam VOF equation

CHUIKOV · July 14, 2022, 07:02

Dear Foamers

I'm looking into icoReactingMultiphaseInterFoam for a laser melting case with added mass (blown powder or wire).

I had a look at the solver structure but cant find the alphaEqn.H file. I would like add a source term for the liquid phase in the corresponding VOF equation.

Any idea how to proceed?

many thanks in advance,

Chuikov

überschwupper · July 15, 2022, 02:24

Peace,

the alpha equation is integrated inside of the >>fluid.solve()<< routine of the multiphaseSystem class. The solve() function then calls the function >>solveAlphas()<< function.

My ideas would be:
1.) calculate your source term and add it to the equation. Ensure that this term is zero for the other phases
2.) Conditional statement and just add your source term if the condition is met

Copy that library in you rown library folder

recompile that library into FOAM_USER_LIBBIN (best under a different name)

change inside the Make/options file inside the icoReactingMultiphaseInterFoam directory the source directory of your library and the library name

CHUIKOV · July 19, 2022, 02:12

many thanks for your kind assistance!

looks like quite a complex function, I will try my luck with a coded source via fvOptions

cheers!

July 14, 2022, 07:02	Add source term in icoReactingMultiphaseInterFoam VOF equation	#1
CHUIKOV New Member Jon Join Date: Mar 2013 Posts: 14 Rep Power: 13	Dear Foamers I'm looking into icoReactingMultiphaseInterFoam for a laser melting case with added mass (blown powder or wire). I had a look at the solver structure but cant find the alphaEqn.H file. I would like add a source term for the liquid phase in the corresponding VOF equation. Any idea how to proceed? many thanks in advance, Chuikov

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[Other] Tabulated thermophysicalProperties library	chriss85	OpenFOAM Community Contributions	62	October 2, 2022 03:50
Source Term due to evaporation in energy transport equation	styleworker	OpenFOAM Programming & Development	3	September 7, 2022 03:09
[swak4Foam] swak4foam for OpenFOAM 4.0	mnikku	OpenFOAM Community Contributions	80	May 17, 2022 08:06
Using PengRobinsonGas EoS with sprayFoam	Jabo	OpenFOAM Running, Solving & CFD	35	April 29, 2022 15:35
friction forces icoFoam	ofslcm	OpenFOAM	3	April 7, 2012 10:57

July 15, 2022, 02:24		#2
überschwupper Member Join Date: Jan 2022 Location: Germany Posts: 72 Rep Power: 4	Peace, the alpha equation is integrated inside of the >>fluid.solve()<< routine of the multiphaseSystem class. The solve() function then calls the function >>solveAlphas()<< function. My ideas would be: 1.) calculate your source term and add it to the equation. Ensure that this term is zero for the other phases 2.) Conditional statement and just add your source term if the condition is met Copy that library in you rown library folder recompile that library into FOAM_USER_LIBBIN (best under a different name) change inside the Make/options file inside the icoReactingMultiphaseInterFoam directory the source directory of your library and the library name

July 19, 2022, 02:12		#3
CHUIKOV New Member Jon Join Date: Mar 2013 Posts: 14 Rep Power: 13	many thanks for your kind assistance! looks like quite a complex function, I will try my luck with a coded source via fvOptions cheers!