/*?duct.c?*************************************\
*??-?Turbulent?Flow?in?Duct?modeled?by?LES?-???*
*????????two-stage?time?differencing???????????*
*??written?by?Andrei?Chernousov?Mar?16?2001???*
*??E-mail:?????????????????*
*??http://www.geocities.com/andrei_chernousov??*
\**********************************************/

#include?
#include?
#include??

#include?

/*?continuation?flag?（to?be?changed?by?SIGINT?handler）?*/
int?continFlag?=?1;

/*?SIGINT?handler?*/
#pragma?argsused
void?termHandler（int?sigType）
{
??sigType?=?sigType;
??if?（continFlag--）;?/*?1?->?0?*/
??fprintf（stderr?“caught?[Ctrl+C]?terminating...“）;
}

/*---?Combinations?of?the?ratio?of?specific?heats?---*/
#define?K????????1.4
#define?K_1??????0.4
#define?_K_1?????（-0.4）
#define?_05K?????0.7
#define?K_K_1????3.5
#define?K_1_2????0.2
#define?_1_2_K_1?1.25

#define?Cv??????717.75??/*?specific?constant?volume?heat?capacity?J/（kg*K）?*/
#define?R_VOZD??287.10??/*?gas?constant?for?air?*/
#define?KR_VOZD?401.94??/*?R_VOZD?*?K?*/

#define?Fx1?xU1?/*?aliases?for?x-fluxes?*/
#define?Fx2?xU2
#define?Fx3?xU3
#define?Fx4?xU4
#define?Fx5?xU5

#define?Fy1?yU1?/*?aliases?for?y-fluxes?*/
#define?Fy2?yU2
#define?Fy3?yU3
#define?Fy4?yU4
#define?Fy5?yU5

#define?Fz1?zU1?/*?aliases?for?z-fluxes?*/
#define?Fz2?zU2
#define?Fz3?zU3
#define?Fz4?zU4
#define?Fz5?zU5

/*===?Global?variables?===*/
float
??***U1_?***U2_?***U3_?***U4_?***U5_?/*?substitutional?pointers?*/
??/*?Cons.?variables?at?the?cells‘?centroids?*/
??***U1??***U2??***U3??***U4??***U5?
??***U1p?***U2p?***U3p?***U4p?***U5p
??/*?Cons.?variables?at?cell?boundaries?in?x-direction?*/
??***xU1?***xU2?***xU3?***xU4?***xU5?
??***U1x?***U2x?***U3x?***U4x?***U5x
??/*?Cons.?variables?at?cell?boundaries?in?y-direction?*/
??***yU1?***yU2?***yU3?***yU4?***yU5?
??***U1y?***U2y?***U3y?***U4y?***U5y
??/*?Cons.?variables?at?cell?boundaries?in?z-direction?*/
??***zU1?***zU2?***zU3?***zU4?***zU5?
??***U1z?***U2z?***U3z?***U4z?***U5z
???
??R0?/*?initial?density?*/
??R?U?V?W?P?C???????/*?vector?of?primitive?flow?parameters?*/
??u1?u2?u3?u4?u5?????/*?conservative?flow?variables?????????*/
??deltaX?deltaY?deltaZ?/*?cell?spacings???????????????????????*/
??deltaT?????????????????/*?step?in?time????????????????????????*/
??deltaT_X?deltaT_Y?deltaT_Z?/*?convenient?ratios?????????????*/
??_deltaX?_deltaY?_deltaZ
??_4deltaX?_4deltaY?_4deltaZ
??_2deltaX?_2deltaY?_2deltaZ;

char?
??Stage?/*?indicator?of?the?current?stage?*/
??str[80];?/*?buffer?string?*/

unsigned
??LEN?HIG?DEP?/*?numbers?o?cells?in?x-?y-?z-directions?for?3d?box?*/
??LENN?HIGG?DEPP
??Answer??/*?flag?of?continuing?solution?*/
??i??j??k?i__??j__??k__?/*?indices?*/
??l
???
??mem?????/*?total?amount?of?dynamic?memory?required?（in?bytes）?*/
??step????/*?current?time?step?*/
??f_step??/*?frame?taking?step?*/
??numstep;?/*?overall?prescribed?number?of?time?steps

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????66752??2009-08-06?23:20??les_code\example\duct\100

?????文件??????71110??2009-08-06?23:20??les_code\example\duct\1000

?????文件??????71310??2009-08-06?23:20??les_code\example\duct\1100

?????文件??????64865??2009-08-06?23:20??les_code\example\duct\12

?????文件??????71485??2009-08-06?23:20??les_code\example\duct\1200

?????文件??????71601??2009-08-06?23:20??les_code\example\duct\1300

?????文件??????71675??2009-08-06?23:20??les_code\example\duct\1400

?????文件??????71779??2009-08-06?23:20??les_code\example\duct\1500

?????文件??????71788??2009-08-06?23:20??les_code\example\duct\1600

?????文件??????71813??2009-08-06?23:20??les_code\example\duct\1700

?????文件??????71894??2009-08-06?23:20??les_code\example\duct\1800

?????文件??????71937??2009-08-06?23:20??les_code\example\duct\1900

?????文件??????67592??2009-08-06?23:20??les_code\example\duct\200

?????文件??????71979??2009-08-06?23:20??les_code\example\duct\2000

?????文件??????71955??2009-08-06?23:20??les_code\example\duct\2100

?????文件??????71956??2009-08-06?23:20??les_code\example\duct\2200

?????文件??????72009??2009-08-06?23:20??les_code\example\duct\2300

?????文件??????72002??2009-08-06?23:20??les_code\example\duct\2400

?????文件??????72119??2009-08-06?23:20??les_code\example\duct\2500

?????文件??????72182??2009-08-06?23:20??les_code\example\duct\2600

?????文件??????72224??2009-08-06?23:20??les_code\example\duct\2700

?????文件??????72332??2009-08-06?23:20??les_code\example\duct\2800

?????文件??????72185??2009-08-06?23:20??les_code\example\duct\2900

?????文件??????68570??2009-08-06?23:20??les_code\example\duct\300

?????文件??????72277??2009-08-06?23:20??les_code\example\duct\3000

?????文件??????72185??2009-08-06?23:20??les_code\example\duct\3100

?????文件??????72273??2009-08-06?23:20??les_code\example\duct\3200

?????文件??????72302??2009-08-06?23:20??les_code\example\duct\3300

?????文件??????72359??2009-08-06?23:20??les_code\example\duct\3400

?????文件??????72348??2009-08-06?23:20??les_code\example\duct\3500

............此處省略100個文件信息