|
[Sponsors] |
March 4, 2022, 09:18 |
Melting of Steel via 3D Goldak Heat Source
|
#1 |
New Member
Joe Thermal
Join Date: Feb 2022
Posts: 4
Rep Power: 4 |
Hello together,
I'm trying to model a 3D volumetric Goldak heat source that is moving over a simple steel plate. The problem I’m facing is that the simulation stops after a short amount of iterations (about 0.2s simulated time). I used an expression for thermal expansion and phase density difference buoyancy forces, the Kozeny-carman permeability model with min/max clipping. The setup is shown in the attached pictures. The used CEL expressions are : LIBRARY: CEL: EXPRESSIONS: BetaLiq = 4.2333e-05[K^-1] BuoySrcRho = ((LiqFrac*rhoL+(1-LiqFrac)*rhoS)-rhoRef)*g BuoySrcT = -BetaLiq*rhoL*(T-Tref)*LiqFrac*g Kperm = Kstar*(1-LiqFrac)^2/max(LiqFrac^3, 1e-10) Kperm2 = max(min(KpermMax,Kperm), 1e-15 [m^-2]) KpermMax = 1.0E+9[kg s^-1 m^-3] / Visc Kstar = 1e10[m^-2] L = 246800[J/kg] LiqFrac = Liquid Steel.Mass Fraction PL = 2400[W] Q = if(qf1*qf+qr1*qr<=0[kg m^-1 s^-3],1[kg m^-1 s^-3],qf1*qf+qr1*qr) Tcold = 1700[K] Thot = 1800[K] Tmelt = 1800[K] Tref = Tcold Visc = Liquid Steel.Dynamic Viscosity a = 0.004[m] b = 0.005[m] cf = 0.0025[m] cr = 0.0075[m] ff = 0.5 fr = 1.5 kf = (6*sqrt(3)*ff*PL)/(a*b*cf*pi*sqrt(pi)) kr = (6*sqrt(3)*fr*PL)/(a*b*cr*pi*sqrt(pi)) qf = kf*exp(-3*((x-vl*t)^2/cf^2+y^2/a^2+z^2/b^2)) qf1 = if(x>vl*t,1,0) qr = kr*exp(-3*((x-vl*t)^2/cr^2+y^2/a^2+z^2/b^2)) qr1 = if(x<vl*t,1,0) rhoL = Liquid Steel.Density rhoRef = rhoS rhoS = 7854[kg/m^3] vl = 0.005[m s^-1] END END Something is wrong with this model, but I don’t know what it is. Maybe someone has a clue what is wrong with the simulation. Best regards Joe T. |
|
March 4, 2022, 09:20 |
Additional Infos
|
#2 |
New Member
Joe Thermal
Join Date: Feb 2022
Posts: 4
Rep Power: 4 |
More setup details
|
|
March 4, 2022, 14:04 |
More Setup Details II, more setup screenshots
|
#3 |
New Member
Joe Thermal
Join Date: Feb 2022
Posts: 4
Rep Power: 4 |
More details on the run...
|
|
March 4, 2022, 14:07 |
Finally the last screenshots are about to arrive
|
#4 |
New Member
Joe Thermal
Join Date: Feb 2022
Posts: 4
Rep Power: 4 |
Finally the last screenshots are about to arrive
|
|
March 4, 2022, 17:53 |
|
#5 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,862
Rep Power: 144 |
CFX does not have a melting model, so of course it cannot do this.
What are you trying to do here? Why are you doing this? Modelling welding is very complex and difficult and CFX is not ideal for it. But depending on what you are trying to achieve you might be able to do it in CFX.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
March 8, 2022, 14:07 |
|
#6 |
Senior Member
Erik
Join Date: Feb 2011
Location: Earth (Land portion)
Posts: 1,188
Rep Power: 23 |
Just one comment: You seem to have the same reference enthalpy and point for each phase. Meaning your latent heat of melting is zero, unless you defined it elsewhere?
|
|
Tags |
cfx, laser, momentum sources, thermal |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[swak4Foam] funkyDoCalc with OF2.3 massflow | NiFl | OpenFOAM Community Contributions | 14 | November 25, 2020 04:30 |
[OpenFOAM.org] Error creating ParaView-4.1.0 OpenFOAM 2.3.0 | tlcoons | OpenFOAM Installation | 13 | April 20, 2016 18:34 |
Trouble compiling utilities using source-built OpenFOAM | Artur | OpenFOAM Programming & Development | 14 | October 29, 2013 11:59 |
centOS 5.6 : paraFoam not working | yossi | OpenFOAM Installation | 2 | October 9, 2013 02:41 |
[swak4Foam] funkySetFields compilation error | tayo | OpenFOAM Community Contributions | 39 | December 3, 2012 06:18 |