CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT > Fluent UDF and Scheme Programming

Gauss plane heat source in the interface of two-phase flow

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 2, 2021, 21:01
Default Gauss plane heat source in the interface of two-phase flow
  #1
New Member
 
Xu Zhao
Join Date: Jun 2021
Posts: 5
Rep Power: 5
zx145833 is on a distinguished road
Hello,
Everyone, I have a trouble about the loading of plane heat source of welding arc in the interface of the two-phase flow (metal-air system). The partial UDF is listed as below. The delta function is introduced to transform the plane heat source to volume heat source. While it is noted that the forum needs the gradient of volume fraction of metal phase and I'm not sure the transformation is reasonable. In fact, the function is referenced from the level set method. When I test the result, the metal surface temperature doesn't increase and keeps room temperature. Anyone who can give me some suggestions and thanks a lot.

#define PI 3.14159
#define reff 0.00284 /*Arc effective radius*/
#define current 280. /*Welding current*/
#define voltage 29.5 /*Welding voltage*/
#define eta_arc 0.8 /* Arc heat efficiency*/

DEFINE_SOURCE(arc_source,c,t,dS,eqn)
{
real xc[ND_ND],x,y,z,time,Qarc,r,source,var1,delta,grad;
Thread *pri_th, *sec_th;
pri_th=THREAD_SUB_THREAD(t,0);
sec_th=THREAD_SUB_THREAD(t,1);
time=RP_Get_Real("flow-time");
C_CENTROID(xc, c, t);
x=xc[0]-xx2-v*time; /*xx2 is the initial position of arc in x
direction*/
y=xc[1];
z=xc[2];
r=sqrt(x*x+y*y);
Qarc=current*voltage;
grad=sqrt(C_UDMI(c,t,0)*C_UDMI(c,t,0)+C_UDMI(c,t,1 )*C_UDMI(c,t,1)+C_U
DMI(c,t,2)*C_UDMI(c,t,2)+0.00000001);
delta=6*C_VOF(c,sec_th)*(1-C_VOF(c,sec_th))*grad;
if(r<=reff&&C_VOF(c,sec_th)>0.&&C_VOF(c,sec_th)<1. )
{
source=eta_arc*Qarc/2./PI/(reff*reff)*exp(-r*r/(2.*reff*reff));
dS[eqn]=0;
}
else
{
source=0;
dS[eqn]=0;
}
return source*delta;
}
zx145833 is offline   Reply With Quote

Old   July 2, 2021, 23:42
Default
  #2
Member
 
Join Date: Jul 2020
Location: India
Posts: 65
Rep Power: 5
Cooper24 is on a distinguished road
I am also working on a similar problem and had written about it in the forum yesterday. I have sent you a PM. Kindly revert back.
Cooper24 is offline   Reply With Quote

Old   July 7, 2021, 01:32
Default
  #3
Senior Member
 
Alexander
Join Date: Apr 2013
Posts: 2,363
Rep Power: 34
AlexanderZ will become famous soon enoughAlexanderZ will become famous soon enough
compile your code first, read log files, fix errors
__________________
best regards


******************************
press LIKE if this message was helpful
AlexanderZ is offline   Reply With Quote

Reply

Tags
heat source, two-phase flow

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
Table bounds warnings at: END OF TIME STEP CFXer CFX 4 July 16, 2020 23:44
[foam-extend.org] Problems installing foam-extend-4.0 on openSUSE 42.2 and Ubuntu 16.04 ordinary OpenFOAM Installation 19 September 3, 2019 18:13
My radial inflow turbine Abo Anas CFX 27 May 11, 2018 01:44
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 06:28
friction forces icoFoam ofslcm OpenFOAM 3 April 7, 2012 10:57


All times are GMT -4. The time now is 11:16.