Hello everyone,<\/p>\n
I am trying to simulate Selective Laser Melting in Ansys Fluent. Its a two phase problem. So I am using VOF model.<\/p>\n
my domain dimension is :<\/p>\n
L = 0.8 mm (X axis)
W = 0.2 mm (Y axis)
H = 0.28 mm (Z axis)<\/p>\n
In this domain , I patched the domain ( x = 0.8 mm, y = 0.2 mm , z= 0.2 mm) as solid and the rest is gas.<\/p>\n
I am trying to model the heat source as surface heat source and place it at x = 0 mm , y = 0.1 mm and z is undefined as I will use if statement to find the interface ( where 0<VOF<1)<\/p>\n
My UDF code is given below. After Compiling and initialising , when I try to view the user defined volumetric heat source under the contour tab , FLUENT Crashes instantly. I am not sure what went wrong. Yes I also created 5 user defined memory locations. Can anyone help?<\/p>\n
#include “udf.h”
#include “sg_mphase.h”
#include “mem.h”
#include “sg_mem.h”
#include “math.h”
#include “flow.h”
#include “unsteady.h”
#include “metric.h”<\/p>\n
\/\/ Constants
#define A 0.4 \/\/ Absorption coefficient
#define P 200 \/\/ Laser power (W)
#define R 40e-6 \/\/ Spot radius (m)
#define v 0.5 \/\/ Scan speed of laser (m\/s)
#define Pi 3.1415926535 \/\/ Pi constant
#define x0 0.0 \/\/ Initial x position of the laser (m)
#define y0 0.1e-3 \/\/ Initial y position of the laser (m)<\/p>\n
\/\/ UDF for adjusting the gradient heat
DEFINE_ADJUST(adjust_gradient_heat, domain)
{
Thread *t;
Thread **pt;
cell_t c;
int phase_domain_index = 3.0; \/\/ thjats my metal domain
Domain *pDomain = DOMAIN_SUB_DOMAIN(domain, phase_domain_index);<\/p>\n
Alloc_Storage_Vars(pDomain, SV_VOF_RG, SV_VOF_G, SV_NULL);
Scalar_Reconstruction(pDomain, SV_VOF, -1, SV_VOF_RG, NULL);
Scalar_Derivatives(pDomain, SV_VOF, -1, SV_VOF_G, SV_VOF_RG, Vof_Deriv_Accumulate);<\/p>\n
mp_thread_loop_c(t, domain, pt)
if (FLUID_THREAD_P(t))
{
Thread *ppt = pt[phase_domain_index];<\/p>\n
begin_c_loop(c, t)
{
C_UDMI(c, t, 0) = C_VOF_G(c, ppt)[0];
C_UDMI(c, t, 1) = C_VOF_G(c, ppt)[1];
C_UDMI(c, t, 2) = C_VOF_G(c, ppt)[2];
C_UDMI(c, t, 3) = sqrt(C_UDMI(c, t, 0) * C_UDMI(c, t, 0) +
C_UDMI(c, t, 1) * C_UDMI(c, t, 1) +
C_UDMI(c, t, 2) * C_UDMI(c, t, 2)); \/\/ Magnitude of gradient of volume fraction<\/p>\n
}
end_c_loop(c, t)
}<\/p>\n
Free_Storage_Vars(pDomain, SV_VOF_RG, SV_VOF_G, SV_NULL);
}<\/p>\n
\/\/ UDF for defining the heat source
DEFINE_SOURCE(heat_source, c, t, dS, eqn)
{
Thread *pri_th; \/\/ Gas phase
Thread *sec_th; \/\/ solid phase
real source;
real x[ND_ND], time;
time = CURRENT_TIME; \/\/ Acquire time from Fluent solver
C_CENTROID(x, c, t); \/\/ Acquire the cell centroid location
real T = C_T(c, t);<\/p>\n
pri_th = THREAD_SUB_THREAD(t, 0);
sec_th = THREAD_SUB_THREAD(t, 1);<\/p>\n
real rho = C_R(c,t);
real Cp = C_CP(c,t);<\/p>\n
real rhom = C_R(c,sec_th);
real Cpm = C_CP(c,sec_th);<\/p>\n
real rhog = C_R(c,pri_th);
real Cpg = C_CP(c,pri_th);
real factor = (2 * rho * Cp) \/ (rhom * Cpm + rhog * Cpg);<\/p>\n
real r = sqrt(pow(x[0] – x0 – v * time, 2.0) + pow(x[1] – y0, 2.0));<\/p>\n
if (C_VOF(c, t) > 0.05 && C_VOF(c, t) < 1)
{
source = ((2 * A * P) \/ (Pi * R * R)) * exp((-2 * (r * r)) \/ (R * R))*factor * C_UDMI(c,t,3);
dS[eqn] = 0.0;
}
else
{
source = 0.0;
dS[eqn] = 0.0;
}<\/p>\n
return source;
}<\/p>\n","protected":false},"template":"","class_list":["post-377965","topic","type-topic","status-publish","hentry"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_bbp_subscription":["440817","199"],"_bbp_author_ip":["118.179.2.146"]," _bbp_last_reply_id":["0"]," _bbp_likes_count":["0"],"_btv_view_count":["165"],"_bbp_topic_status":["unanswered"],"_bbp_topic_id":["377965"],"_bbp_forum_id":["27792"],"_bbp_engagement":["199","440817"],"_bbp_voice_count":["2"],"_bbp_reply_count":["15"],"_bbp_last_reply_id":["379982"],"_bbp_last_active_id":["379982"],"_bbp_last_active_time":["2024-08-23 10:08:13"]},"test":"rakibul-buet19gmail-com"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/377965","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics"}],"about":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/types\/topic"}],"version-history":[{"count":0,"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/377965\/revisions"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/media?parent=377965"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}