CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

Melting of Steel via 3D Goldak Heat Source

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 4, 2022, 09:18
Default Melting of Steel via 3D Goldak Heat Source
  #1
thj
New Member
 
Joe Thermal
Join Date: Feb 2022
Posts: 4
Rep Power: 2
thj is on a distinguished road
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.
Attached Images
File Type: jpg 01.jpg (51.7 KB, 10 views)
File Type: jpg 02.jpg (50.9 KB, 9 views)
File Type: jpg 03.jpg (72.9 KB, 8 views)
File Type: jpg 04.jpg (63.6 KB, 9 views)
File Type: png 05.png (17.3 KB, 3 views)
thj is offline   Reply With Quote

Old   March 4, 2022, 09:20
Default Additional Infos
  #2
thj
New Member
 
Joe Thermal
Join Date: Feb 2022
Posts: 4
Rep Power: 2
thj is on a distinguished road
More setup details
Attached Images
File Type: png 06.png (18.6 KB, 5 views)
File Type: png 07.png (16.3 KB, 4 views)
File Type: png 08.png (36.1 KB, 4 views)
File Type: png 09.png (16.0 KB, 1 views)
File Type: png 10.png (33.8 KB, 1 views)
thj is offline   Reply With Quote

Old   March 4, 2022, 14:04
Default More Setup Details II, more setup screenshots
  #3
thj
New Member
 
Joe Thermal
Join Date: Feb 2022
Posts: 4
Rep Power: 2
thj is on a distinguished road
More details on the run...
Attached Images
File Type: png 11.png (18.2 KB, 2 views)
File Type: png 12.png (26.5 KB, 3 views)
File Type: jpg 13.jpg (120.1 KB, 0 views)
File Type: png 14.png (12.4 KB, 1 views)
File Type: jpg 15.jpg (136.4 KB, 1 views)
thj is offline   Reply With Quote

Old   March 4, 2022, 14:07
Default Finally the last screenshots are about to arrive
  #4
thj
New Member
 
Joe Thermal
Join Date: Feb 2022
Posts: 4
Rep Power: 2
thj is on a distinguished road
Finally the last screenshots are about to arrive
Attached Images
File Type: jpg 16.jpg (101.3 KB, 3 views)
File Type: jpg 17.jpg (97.3 KB, 2 views)
File Type: jpg 18.jpg (132.8 KB, 1 views)
thj is offline   Reply With Quote

Old   March 4, 2022, 17:53
Default
  #5
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,231
Rep Power: 135
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
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.
ghorrocks is offline   Reply With Quote

Old   March 8, 2022, 14:07
Default
  #6
Senior Member
 
Erik
Join Date: Feb 2011
Location: Earth (Land portion)
Posts: 1,116
Rep Power: 21
evcelica is on a distinguished road
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?
evcelica is offline   Reply With Quote

Reply

Tags
cfx, laser, momentum sources, thermal

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
[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


All times are GMT -4. The time now is 17:09.