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

segmentation error(again)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 13, 2016, 14:54
Default segmentation error(again)
  #1
Boh
Member
 
Join Date: Sep 2016
Posts: 35
Rep Power: 3
Boh is on a distinguished road
Hi! I have again a problem with the segmentation error....the function works fine until the 3 face loop...but when there is also the third loop the segmentation error appears...does anyone have any idea?



#include "udf.h"
#define Diff 2.82e-05
#define k 0.0264
#define cp 1005
#define Dh 11.85e-3
#define ht 16.8
#define s 1e-3
#define Cads 0.001

DEFINE_PROFILE(Conc_boundary, t, i)
{
real Tin = 0;
real Tout = 0;
real Tm;
real Cinter;
real Cinf;
real xy[ND_ND];
Domain *d;
real Le;
real rho;
real a = 0;
real b = 0;
real hm;
Thread *t_interface;
Thread *t_air_symmetry;
Thread *t_inlet;
Thread *t_outlet;
Thread *t0_air;
cell_t c0;
Thread *t0;
cell_t c0_air;
Thread *t0_inter;
cell_t c0_inter;
face_t f;
t_interface = Lookup_Thread(d, 8);
t_air_symmetry = Lookup_Thread(d, 18);
t_inlet = Lookup_Thread(d, 15);
t_outlet = Lookup_Thread(d, 22);
begin_f_loop(f, t_inlet)
{
c0 = F_C0(f, t_inlet);
t0 = THREAD_T0(t_inlet);
F_CENTROID(xy, c0, t0);
Tin+= F_T(c0,t0);
a++;
}
end_f_loop(f,t_inlet)
begin_f_loop(f, t_outlet)
{

c0 = F_C0(f, t_outlet);
t0 = THREAD_T0(t_outlet);
F_CENTROID(xy, c0, t0);
Tout+= F_T(c0, t0);
b++;
}
end_f_loop(f, t_outlet)
Tin = Tin / a;
Tout = Tout / b;
Tm = 0.5*(Tout + Tin);
rho = pow(10, 5) / (287 * Tm);
Le = k / (rho*cp*Diff);
hm = ht / (pow(Le, 2 / 3)*rho*cp);
begin_f_loop(f, t)
{

c0_air = F_C0(f, t_air_symmetry);
t0_air = THREAD_T0(t_air_symmetry);
F_CENTROID(xy, c0_air, t0_air);
Cinf = F_YI(c0_air, t0_air, i);
c0_inter = F_C0(f, t_interface);
t0_inter = THREAD_T0(t_interface);
F_CENTROID(xy, c0_inter, t0_inter);
Cinter = F_YI(c0_inter, t0_inter,i);
F_PROFILE(f, t, i) = Cinter - s*hm / Diff*(Cinf - Cinter);
}
end_f_loop(f, t)

}
Boh is offline   Reply With Quote

Old   November 14, 2016, 05:06
Default
  #2
Boh
Member
 
Join Date: Sep 2016
Posts: 35
Rep Power: 3
Boh is on a distinguished road
The error appears when I try to hook the function to the boundary
Boh is offline   Reply With Quote

Old   November 14, 2016, 20:15
Default
  #3
`e`
Senior Member
 
Join Date: Mar 2015
Posts: 800
Rep Power: 11
`e` is on a distinguished road
The third face loop is using 't' instead of 't_air_symmetry'.

Code:
begin_f_loop(f, t)
{
    c0_air = F_C0(f, t_air_symmetry);
    ...
`e` 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
Segmentation fault when running dieselFoam or dieselEngineFoam in parallel francesco OpenFOAM Bugs 4 May 2, 2017 21:59
sigSegv error, segmentation error? ollebapur OpenFOAM Running, Solving & CFD 0 August 18, 2014 05:20
Segmentation Fault w/ compiled OF 2.2.0 - motorBike example sudo OpenFOAM Running, Solving & CFD 3 April 2, 2013 17:27
segmentation fault when installing OF-2.1.1 on a cluster Rebecca513 OpenFOAM Installation 9 July 31, 2012 15:06
Segmentation Violation Corentin FLUENT 1 February 13, 2011 02:07


All times are GMT -4. The time now is 10:07.