October 14, 2022 at 4:07 amDubey92Subscriber
I am modelling a Laser melting problem using Solidification/melting. I am applying the laser at a face using UDF. There is no issue in compilation. But the simulation does not starts and throws multiple warnings and errors. I have attache the screenshot. I have checked the properties also. They are fine. Any help is much appreciated. Thanks in advance.
October 14, 2022 at 1:34 pmFederico Alzamora PrevitaliSubscriber
floating point exception typically results from the solver trying to divide by zero. I see the warning states that Cp is zero at two given temperatures. Can you tell us how was Copper defined in your material properties? For example, is it using constant Cp? or other?
October 15, 2022 at 10:17 amDubey92Subscriber
Hi Federico. I am defining copper properties using UDF. Also when I tried to apply a constant flux on the surface instead of the Laser, it works. I am applying the Laser using UDF and that is what creating the problem. I am using DEFINE_PROFILE for applying the Laser. Here is my UDF:
DEFINE_PROFILE(Laser, t, i) // The name of the UDF is Laser
real x[ND_ND], time; // Define face centroid vector, time
time = RP_Get_Real("flow-time"); // Acquire time from Fluent solver
face_t f; // face thread
real r = sqrt(pow(x-x0-v*time,2.0) + pow(x-y0,2.0));
real T = F_T(f,t);
F_PROFILE(f,t,i) = (((2*eta*P)/(Pi*R*R))*exp((-2*(r*r))/(R*R))) - ha*(T-Ta) - (s*e*(pow(T,4) - pow(Ta,4)));
F_UDMI(f,t,0) = F_PROFILE(f,t,i);
F_PROFILE(f,t,i) = - ha*(T-Ta) - (s*e*(pow(T,4) - pow(Ta,4)));
October 17, 2022 at 3:41 pmRobAnsys Employee
If you're adjusting the cp using a UDF you might want to read https://ansyshelp.ansys.com/account/secured?returnurl=/Views/Secured/corp/v222/en/flu_udf/flu_udf_ModelSpecificDEFINE.html%23flu_udf_sec_define_property as opposed to DEFINE_PROFILE Otherwise, check the exact curve you get with the cp values: it must be over 0 at all times. The 1K warning is generally when the energy equation is going horribly wrong, unless you're melting copper at VERY low temperatures....
October 18, 2022 at 4:08 amDubey92Subscriber
Thanks Rob. I figured out the issue. In my DEFINE_PROFILE for Laser, the terms consist of Laser, Convective heat loss and Radiative heat loss. I ran the simulation with only Laser and Convective Loss and it worked. The issue of Invalid Cp, Floating point exception and Divergence in the solver comes when I add the radiative loss also. I don't know why this is happening.
October 18, 2022 at 9:54 amRobAnsys Employee
Plot the radiative loss over the full range of temperatures (say 50K below the minimum your're expecting to 50K above) using the equation in the UDF. I wonder if you've sucked too much heat out of the domain somewhere.
- You must be logged in to reply to this topic.
Simulation World 2022
Check out more than 70 different sessions now available on demand. Get inspired as you hear from visionary companies, leading researchers and educators from around the globe on a variety of topics from life-saving improvements in healthcare, to bold new realities of space travel. Take a leap of certainty and check out a session today here.
Earth Rescue – An Ansys Online Series
The climate crisis is here. But so is the human ingenuity to fight it. Earth Rescue reveals what visionary companies are doing today to engineer radical new ideas in the fight against climate change. Click here to watch the first episode.
Subscribe to the Ansys Blog to get great new content about the power of simulation delivered right to your email on a weekly basis. With content from Ansys experts, partners and customers you will learn about product development advances, thought leadership and trends and tips to better use Ansys tools. Sign up here.
- Suppress Fluent to open with GUI while performing in journal file
- Heat transfer coefficient
- What are the differences between CFX and Fluent?
- Floating point exception in Fluent
- The solver failed with a non-zero exit code of : 2
- Getting graph and tabular data from result in workbench mechanical
- Difference between K-epsilon and K-omega Turbulence Model
- Time Step Size and Courant Number
- Mesh Interfaces in ANSYS FLUENT
- error in cfd post
© 2022 Copyright ANSYS, Inc. All rights reserved.