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

maxwell boundary condition

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 3, 2013, 02:31
Default maxwell boundary condition
  #1
New Member
 
morteza
Join Date: Mar 2013
Posts: 2
Rep Power: 0
gashti is on a distinguished road
Hi, I write this UDF for slip flow for fluent 6.3 but i receive this error15Error: FLUENT received fatal signal (ACCESS_VIOLATION)1. Note exact events leading to error.2. Save case/data under new name.3. Exit program and restart to continue.4. Report error to your distributor.Error Object: ()15I write my UDF by your suggestion, please guide me for this errorthanksmy UDF:

#include "udf.h" /* must be at the beginning of every UDF */
#include "mem.h" /* must be at the beginning of this UDF */
#include "math.h" /* must be at the beginning of this UDF */
#include "metric.h" /* must be at the beginning of this UDF */
#define k_B 1.3805e-23 /* Boltzman constant (j/K)*/
#define sigma 419e-12 /* molecular diameter*/
#define sigma_v 1 /* Tangential momentum accommodation coefficient*/
#define sigma_T 1 /*Thermal accommodation coefficient*/
#define pi 3.14159 /* pi number*/
#define sqrt_2 1.41421 /* sqrt(2) number */
DEFINE_PROFILE(wall_slip_velocity, t , index )
{
double T , P, density , MU , K , landa ,a_m , du_da , dT_dx ,u_f ;
double A[ND_ND] ;
face_t f;
Thread *t0;
cell_t c0;

begin_f_loop(f,t)
{
t0 = THREAD_T0(t); /* adjacent cell thread to f */
c0 = F_C0(f, t); /* adjacent cell to f */
T=C_T(c0,t0); /* temperature of the cell */
P=C_P(c0,t0); /* peressure of the cell */
density=C_R(c0,t0); /* density of the cell*/
MU = C_MU_L(c0,t0); /* laminar viscosity of the cell */
K = C_K_L(c0,t0); /* thermal conductivity of the cell */
landa=(k_B*T)/(sqrt_2*pi*sigma*sigma*P); /* mean free path line */

F_AREA(A,f,t);
a_m = NV_MAG(A); /* magnitude of the cell area vector */
du_da=NV_DOT(A,A)/a_m; /* n- component of the cell x-velocity reconstruction gradient(RG) vector */
dT_dx=1 ; /* x-component of the cell temperature reconstruction gradient(RG) vector */

F_PROFILE(f,t,index) = ((2-sigma_v)/sigma_v)*landa * du_da + 0.75*(MU/(density*T))*dT_dx ;
}

end_f_loop(f, t)
}
gashti is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Radiation interface hinca CFX 15 January 26, 2014 18:11
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 17:44
Setting outlet Pressure boundary condition using CAFFA code Mukund Pondkule Main CFD Forum 0 March 16, 2011 04:23
Domain Imbalance HMR CFX 3 March 6, 2011 21:10
How exactly the "pressure outlet" bdry condition compute properties on the boundary? yating9901 FLUENT 3 June 28, 2010 12:26


All times are GMT -4. The time now is 21:33.