CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Phase depending absorption & emission model properties

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By OPFO

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 11, 2020, 11:03
Question Phase depending absorption & emission model properties
  #1
New Member
 
Hektor
Join Date: Feb 2020
Posts: 3
Rep Power: 2
hektorWartgard is on a distinguished road
Hello,

I am trying to model laser melting using OpenFOAM1912, with the icoReactingMultiphaseInterFoam solver. This solver has the radiation model laserDTRM built in.

I have managed to solve simple problems nicely, but the optical properties are significantly different in the gas (almost completely transmissive)/liquid(highly absorbing)/solid phases(highly absorbing and reflective).

Is it possible to set the absorptivity/emissivity/E to different values in different phases? Does any of the built in models support this?

I have tried looking around in the source code, but I do not quite follow the radiation models.


Any help on the topic is greatly appreciated.


EDIT: It seems that localDensityAbsorptionEmission does the trick, and it does more or less what I would expect.



Best regards
Hektor


Attachments:

radiationProperties
Code:
radiation       on;
radiationModel  laserDTRM;
solverFreq 3;

absorptionEmissionModel constantAbsorptionEmission;
    mode    "Gaussian";
    nTheta    20;
    nr     100;
    opd    0.4425;
    focalLaserPosition constant     (.0002 .00020 0.0002);
    laserDirection constant     (0 -1 0);
    focalLaserRadius         0.000025;
    laserPower            200;
    sigma                25e-6;
    //qualityBeamLaser        1.0;
    absorptivity    absorptivity    [0 -1 0 0 0 0 0] 1e5;
    emissivity      emissivity      [0 -1 0 0 0 0 0] 0.04;
    E               E               [1 -1 -3 0 0 0 0] 0;

scatterModel    none;
sootModel       none;
transmissivityModel none;
Included absorptionEmissionModels
Quote:
binaryAbsorptionEmission
constantAbsorptionEmission
greyMeanAbsorptionEmission
greyMeanSolidAbsorptionEmission
localDensityAbsorptionEmission
multiBandAbsorptioEmission
multiBandZoneAbsorptionEmission
none
wideBandAbsorptionEmission

Last edited by hektorWartgard; February 12, 2020 at 09:25. Reason: Found solution
hektorWartgard is offline   Reply With Quote

Old   February 18, 2020, 10:29
Default LaserDTRM
  #2
New Member
 
Bill
Join Date: Jun 2019
Posts: 16
Rep Power: 2
OPFO is on a distinguished road
Hi Hektor,

I am also trying to work with LaserDTRM for laser melting. Unfortunately, I don't know much about your problem, but I wonder if it is connected to a problem I am having. I also have three phases set up (gas, liquid, solid), however, the laser source will not travel through the gas phase and into the solid phase. I have attached quick snapshots from paraview. I use a simple blockmesh and then initialize the phases using setFields.

Have you also encountered this problem?

Cheers
Attached Images
File Type: jpg Phases.JPG (22.1 KB, 3 views)
File Type: jpg Temp.JPG (20.9 KB, 3 views)
OPFO is offline   Reply With Quote

Old   February 18, 2020, 10:37
Default
  #3
New Member
 
Hektor
Join Date: Feb 2020
Posts: 3
Rep Power: 2
hektorWartgard is on a distinguished road
My solution has been to use the "localDensityAbsorptionEmission" absorption/emission model, there you can set the absorptivity/emissivity/E of the different phases.

The value can be set to 0 for the gas, and metals typically have around 10^8 1/m. Find good values for your material at (https://refractiveindex.info/).


The other option is to to skip the gas phase and use a boundary condition to model the laser. This will work if you have a high absorptivity and a simple geometry.
hektorWartgard is offline   Reply With Quote

Old   February 18, 2020, 11:09
Default LaserDTRM
  #4
New Member
 
Bill
Join Date: Jun 2019
Posts: 16
Rep Power: 2
OPFO is on a distinguished road
Hi Hektor,

Thanks for your quick reply. I'm glad you were able to figure your problem out. I'm trying to incorporate your comments, but I don't know how exactly to code in the coefficients. Here's what I have:

localDensityAbsorptionEmissionCoeffs
{
alphaNames alpha.gas;
absorptivity absorptivity [0 -1 0 0 0 0 0] 0;
emissivity emissivity [0 -1 0 0 0 0 0] 0;
E E [1 -1 -3 0 0 0 0] 0;

alphaNames alpha.liquid;
absorptivity absorptivity [0 -1 0 0 0 0 0] 0;
emissivity emissivity [0 -1 0 0 0 0 0] 0;
E E [1 -1 -3 0 0 0 0] 0;

alphaNames alpha.solid;
absorptivity absorptivity [0 -1 0 0 0 0 0] 0;
emissivity emissivity [0 -1 0 0 0 0 0] 0;
E E [1 -1 -3 0 0 0 0] 0;
}

Any insight?
Thanks in advance!
OPFO is offline   Reply With Quote

Old   February 19, 2020, 05:19
Default LaserDTRM
  #5
New Member
 
Bill
Join Date: Jun 2019
Posts: 16
Rep Power: 2
OPFO is on a distinguished road
I think I was able to figure it out after trial and error. For anyone else interested, here's what ended up working for me:

localDensityAbsorptionEmissionCoeffs
{
alphaNames (alpha.gas alpha.liquid alpha.solid);
aCoeff (0 1e3 1e10);
eCoeff (0 0.5 0.5);
ECoeff (0 0.5 0.5);
}

Please let me know if I made any errors or left something out.

Thanks!
HPE likes this.
OPFO is offline   Reply With Quote

Old   Yesterday, 06:50
Default
  #6
New Member
 
Hektor
Join Date: Feb 2020
Posts: 3
Rep Power: 2
hektorWartgard is on a distinguished road
Yes this looks correct.
hektorWartgard is offline   Reply With Quote

Reply

Tags
laser melting, phase change material, radiation model

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Access phase fraction field in a phase change model hend OpenFOAM Programming & Development 0 October 17, 2019 06:03
What viscosity to use for solid phase in VOF model redwanamit034 FLUENT 0 June 14, 2017 18:09
Error in Two phase (condensation) modeling adilsyyed CFX 15 June 24, 2015 20:42
Udf to define new mixing law between properties of different phases in Mixture model saurabhg2014 Fluent UDF and Scheme Programming 0 February 24, 2015 04:23
Eulerian Model - properties ssamton FLUENT 1 March 16, 2012 02:34


All times are GMT -4. The time now is 08:50.