[TeresaT]

Hey Dhar

yes, it is included. In the thermophysicalProperties.x files you can find it as Hf under thermodynamics.

Regards,
Teresa

[LucaL]

Hello all,

since some time I work within the icoReactingMultiphaseInterFoam solver with the laserDTRM model.
In doing so, I have come up with a few questions.
These are related to the code "DTRMParticle.C" in lines 222 to 226.

Question 1:
Is there a literature reference to this equation? How was the equation for the DTRM model transformed to the present equation?

Question 2:
The part: (ds*e*sigma*pow4(T)/mathematical:i)/ (1 + ds*a) reminds strongly of the radiation of a gray body.
Does the multiplication (ds*e)/(1 + ds*a) therefore establish the relation to a black body? If so, with certain choice of the coefficients e and a, the body could radiate more strongly than a black body.

Question 3:
What is the significance of the coefficient E? In the code "localDensityAbsorptionEmission.H" in line 72 and 73 this is titled as emission contributions.
Through several simulations (with significantly increased ECoeff values, up to 1e12), I could not find any change in the temperature profile within a solid metal.

Hopefully the questions were understandable.
Looking forward to answers and new thoughts.

Greetings!

[VINIT]

Hello Everyone

I am trying to simulate melt pool dynamics in OpenFOAM using IcoReactingMultiphaseInterfoam but I am not getting the correct results. The temperature is not what it should be and the laser is not applied to top surface. I am attaching the results.

Could anyone help me what is wrong?

[Krapf]

Originally Posted by VINIT

[...] the laser is not applied to top surface.

That's how it is supposed to work. In the radiationProperties you define the absorption coefficient [1/m]. If you choose a value large enough, you will see that the laser beam does not penetrate the material completely. This will also change the temperature, because more energy is absorbed in fewer cells.

[VINIT]

Originally Posted by Krapf

That's how it is supposed to work. In the radiationProperties you define the absorption coefficient [1/m]. If you choose a value large enough, you will see that the laser beam does not penetrate the material completely. This will also change the temperature, because more energy is absorbed in fewer cells.

Hello Krapf

Thanks for the response. I tried your suggestion, and I am getting the melt pool but still the temperature is still low. It should be around 3000 K but I am getting 440 K. What am I doing wrong?

I have attached my radiation file and the result.

Thanks

[Krapf]

Do you get the expected values when you plot the absorbed power in ParaView?

[VINIT]

How do you check the absorbed power?

[Krapf]

IIRC the Q field (power/volume) is written. If so you can plot it in ParaView and use the "Integrate Variables" filter to get the power (you'll get one value).


EDIT: If you use the "Mesh Quality" filter and select the volume as measure, you can use the calculator to multiply Q by the volume of each cell. Then you can look at the distribution of the power.

[VINIT]

Originally Posted by Krapf

IIRC the Q field (power/volume) is written. If so you can plot it in ParaView and use the "Integrate Variables" filter to get the power (you'll get one value).


EDIT: If you use the "Mesh Quality" filter and select the volume as measure, you can use the calculator to multiply Q by the volume of each cell. Then you can look at the distribution of the power.

Thanks for the reply. I was able to get correct temperature. I had to set the correct sigma value for laser. Now I am trying to incorporate convection, vaporization, radiation losses. Do you have any idea about any of these?

[nixiaonan]

Originally Posted by VINIT

Thanks for the reply. I was able to get correct temperature. I had to set the correct sigma value for laser. Now I am trying to incorporate convection, vaporization, radiation losses. Do you have any idea about any of these?

What is your sigma value？thank you。

[Krapf]

Originally Posted by VINIT

Thanks for the reply. I was able to get correct temperature. I had to set the correct sigma value for laser. Now I am trying to incorporate convection, vaporization, radiation losses. Do you have any idea about any of these?

For mass transfer you can use Lee, kineticGasTheory (https://www.openfoam.com/news/main-n...er-and-physics) or interfaceHeatResistance (https://www.openfoam.com/news/main-n...rs-and-physics).

[VINIT]

Hello Everyone

I am trying to simulate LBPF and I have developed the custom solver and geometry. I want to add powder layer but I don't have much knowledge about it. Do you have any tutorial that you can share to create powder layer using CFD-DEM. I have coupled liggghts openfoam but I am not able to find any tutorial for the same.

Thanks in advance