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

maxwell boundary condition

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 3, 2013, 02:31
Default maxwell boundary condition
  #1
New Member
 
morteza hajati
Join Date: Mar 2013
Posts: 6
Rep Power: 13
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


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
Domain Imbalance HMR CFX 5 October 10, 2016 06:57
Radiation interface hinca CFX 15 January 26, 2014 18:11
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 18:44
Setting outlet Pressure boundary condition using CAFFA code Mukund Pondkule Main CFD Forum 0 March 16, 2011 04:23
How exactly the "pressure outlet" bdry condition compute properties on the boundary? yating9901 FLUENT 3 June 28, 2010 13:26


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