Began at --> Wed Oct 6 16:40:30 EDT 1999 Proc 0 reporting for duty *************************************************** aaaa dddd ppppp aaaa ccccc a a d d p p a a c aaaaa d d ppppp aaaaa c a a d d p a a c a a dddd p a a ccccc version ----> ADPAC V1.0 adpac 3-d multiple-block euler/navier-stokes solver (a)dvanced (d)ucted (p)ropfan (a)nalysis (c)odes Serial/Parallel Production Code ADPAC V1.0 release date Feb., 1998 *************************************************** version V1.0 notes: ******************* Available turbulence models: Mixing-Length (algebraic) Baldwin-Lomax (algebraic) Spalart-Allmaras (1-Equation) Goldberg k-R (2-Equation) (wall functions are available) wall functions now ON in separated flow thermodynamic properties are for air with variable or constant gamma 2-d or 3-d mesh blocks supported serial/parallel from same source implicit algorithm is available FIXED B.C. is OK in parallel PINT B.C. is OK with multigrid/parallel 2-D and 3-D yield identical results (corrects previous diffs due to turbulence models) Other bug corections up to 2/02/98 EJH Searching Input File for ADPAC Case Name... ******************************************* data item identified --->CASENAME = box data item identified --->FCART = 1.000000 data item identified --->FMULTI = 1.000000 case name defined as : box multigrid level set as : 1 grid type set as : Cartesian Solution type is explicit: IIMPLIC = 0 1-Equation turbulence model is disabled: I1EQ = 0 2-Equation turbulence model is disabled: I2EQ = 0 ************************************************** begin grid input ************************************************** mesh file name: box.mesh *************** number of grids: (nblks) = 1 **************** grid block dimensions : *********************** block imx jmx kmx ----- --- --- --- 1 25 13 13 block/processor assignment information ************************************** Block/Processor Assignment file box.blkproc does not exist Default block/processor assignment used IPR( Block # ) Processor ----------------- --------- ipr( 1) = 0 initialize storage arrays and pointer locations *********************************************** block 3-d ijk 2-d ij 1-d iii length total length 1 1 1 1 5096 5096 26 array sizes required for this mesh (No multi-grid) ********************************** nra3d (3-d array size) >= 5096 nra1d (1-d array size) >= 26 nbl2d (2-d block size) >= 676 program array dimensions: ************************* nra3d (3-d array size) = 600001 nra1d (1-d array size) = 3201 nbl2d (2-d block size) = 100001 begin reading grid block coordinates ************************************ grid block 1 size 26 14 14 grid block 1 coordinates successfully read ************************************************** begin flow input ************************************************** input file read on unit 5 *********************** data item identified --->CASENAME = box data item identified --->DIAM = 1.000000 data item identified --->RMACH = 0.5000000 data item identified --->CFL = -5.000000 data item identified --->FINVVI = 0.0000000E+00 data item identified --->FCART = 1.000000 data item identified --->FREST = 0.0000000E+00 data item identified --->FFULMG = 0.0000000E+00 data item identified --->FMULTI = 1.000000 data item identified --->FNCMAX = 10.00000 inviscid flow calculation selected ********************************** 4-stage runge_kutta algorithm selected The following debug triggers have been enabled ********************************************** input file data *************** rm ( reference mach number ) = 0.5000 gamma ( specific heat ratio ) = 1.4000 cfl ( cfl number (+ time accurate ) = 5.0000 cflmax ( max cfl number time-step scheme) = 2.5000 zetarat( exponent for diss/res.sm. ratio) = 0.6000 fvtsfac( viscous time step factor ) = 2.5000 vt ( =1, time acc., =0, local time ) = 0.0000 vis2 ( 2nd order dissipation factor ) = 0.5000 vis4 ( 4th order dissipation factor ) = 0.0156 fncmax ( maximum iterations ) = 10.0000 frest ( 1 = restart solution ) = 0.0000 save ( 1 = save solution ) = 1.0000 p3dprt ( 0 = no plot3d output ) = 1.0000 fitchk ( checkpoint iteration interval ) = 100.0000 viscg2 ( coarse grid 2nd order diss coef) = 0.1250 fistep ( recalculate time step ) = 1.0000 fnprnt ( 1 = print flowfield ) = 1.0000 fiprnt ( 1 = print partial flowfield ) = 1.0000 finvvi ( 0 = inviscid flow solution ) = 0.0000 ( 1 = viscous flow solution ) = epsx ( i implicit smoothing factor ) = 1.0000 epsy ( j implicit smoothing factor ) = 1.0000 epsz ( k implicit smoothing factor ) = 1.0000 ftotsm ( 1 = post multigrid smoothing ) = 0.0000 epstot ( post multigrid smoothing coef. ) = 0.2000 iunint ( # it. for unsteady plot3d file ) = 0.0000 itimei ( #it between time step update ) = 1.0000 iturbi ( #it between turb. model update ) = 1.0000 iturbb ( #it before starting turb model ) = 10.0000 iturbf ( #it before freezing turb model ) = 0.0000 ibconf ( #it before freezing bound. cond) = 0.0000 iwallf ( wall function trigger (on=1) ) = 1.0000 imassav( mxng pln avrgng (area=1,mass=2)) = 1.0000 itimerm( 1= quit if too little time left) = 0.0000 fdeltat( implicit time step (sec.) ) = 0.0000 fsolve ( solution trigger ) = 1.0000 iimplicit (implicit algorithm trigger ) = 0 fturbcht( 1) (C/O Turb. Model Trigger ) = 0.0000 igrafix( trigger for graphics display ) = 0.0000 igrafint(#it between graphics update ) = 0.0000 iimgsav( trigger for image capture ) = 0.0000 iimgint( #it between image captures ) = 0.0000 fcart ( 0 = Cylindrical, 1 = Cartesian ) = 1.0000 diam ( blade diameter (ft) ) = 1.0000 tref ( total temperature (deg. r) ) = 518.6700 pref ( total pressure (lbf/ft**2) ) = 2116.2200 rgas ( gas constant ((ft/s)**2/deg r) ) = 1716.2600 pr ( prandtl number ) = 0.7000 prt ( turbulent prandtl number ) = 0.9000 non-dimensional initial values calculated as: ******************************************** rho0 ( initial density ) = 0.8852 u0 ( initial axial velocity ) = 0.5774 v0 ( initial radial velocity ) = 0.0000 w0 ( initial theta velocity ) = 0.0000 ei0 ( initial internal energy ) = 2.5476 h0 ( initial enthalpy ) = 3.5000 p0 ( initial pressure ) = 0.8430 t0 ( initial temperature ) = 0.9524 dmu0 ( initial viscosity ) = 0.0000 dimensional reference values calculated as: ******************************************* rho0 ( density (slugs/ft**3) ) = 0.210433E-02 u0 ( axial velocity (ft/s) ) = 544.724 v0 ( radial velocity (ft/s) ) = 0.000000E+00 w0 ( theta velocity (ft/s) ) = 0.000000E+00 ei0 ( static int.energy (ft-lbf/lbm) ) = 0.226782E+07 h0 ( static enthalpy (ft-lbf/lbm) ) = 0.311560E+07 p0 ( static pressure (lbf/ft**2) ) = 1784.01 t0 ( static temperature(deg. r) ) = 493.971 dmu0 ( viscosity (lbf-s/ft**2)) = 0.359781E-06 values used for nondimensionalization: ************************************** rho ( density (slugs/ft**3) ) = 0.237731E-02 u ( axial velocity (ft/s) ) = 943.490 ei ( total int. energy (ft-lbf/lbm) ) = 890173. p ( total pressure (lbf/ft**2) ) = 2116.22 t ( total temperature (deg. r) ) = 518.670 dmu ( viscosity (lbf-s/ft**2 ) = 2.24297 omega ( rotational speed (rad/s) ) = 943.490 length ( length (ft) ) = 1.00000 time ( time (seconds) ) = 0.105990E-02 multigrid solution parameters specified as: ******************************************* nmulti ( number of multi-grid grid level) = 1 nsubit ( multigrid coarse grid subiterat) = 1 nfulmg ( full multigrid solution trigger) = 0 ncoag1 ( full multigrid initial level ) = 1 ncoag2 ( full multigrid final level ) = 2 nitfmg ( full multigrid coarse mesh iter) = 0 nrdmul ( multigrid bc read trigger (1-y)) = 0 Block Solution Information Table ******************************** Block RPM Dimensional Non-Dimen. Cyl/Cart. C-Grid Processor Num Revolutions Omega Omega Trigger Turb Trig Number Per Minute (Rad/s) 1-Cart 1-On0-Off MYID ----- ------------- ------------- ------------- --------- --------- --------- 1 0.000000E+00 0.000000E+00 0.000000E+00 1 0 0 input data successfully entered ******************************* ************************************************** begin boundary data input ************************************************** Opening casename.boundata, file is box.boundata INLETG 1 1 I I P P J K 1 1 1 13 1 13 1 13 1 13 b.c. data item identified --->INLETG PTOT TTOT 1.0 1.0 1.000000 1.000000 0.0000000E+00 1.0000000E-03 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 Stored emass = 0.0000000E+00 Stored relax = 0.0000000E+00 # j=1 is y=1.0 # outward normal is +y SSIN 1 1 J J P P I K 1 1 1 25 1 13 1 25 1 13 b.c. data item identified --->SSIN # if this were a 'free' bc, you'd need this PTOT TTOT EMINF ALPHA 1.0 1.0 0.5 0.0 # j=13 is y=0.0 # outward normal is -y FREE 1 1 J J M M I K 13 13 1 25 1 13 1 25 1 13 b.c. data item identified --->FREE PTOT TTOT EMINF ALPHA 1.000000 1.000000 0.5000000 0.0000000E+00 0.0000000E+00 1.0000000E-03 # k=1 is z=0.0 # outward normal is -z SSIN 1 1 K K P P I J 1 1 1 25 1 13 1 25 1 13 b.c. data item identified --->SSIN # if this were a 'free' bc, you'd need this PTOT TTOT EMINF ALPHA 1.0 1.0 0.5 0.0 # k=13 is z=1.0 # outward normal is +z FREE 1 1 K K M M I J 13 13 1 25 1 13 1 25 1 13 b.c. data item identified --->FREE PTOT TTOT EMINF ALPHA 1.000000 1.000000 0.5000000 0.0000000E+00 0.0000000E+00 1.0000000E-03 EXITG 1 1 I I M M J K 25 25 1 13 1 13 1 13 1 13 b.c. data item identified --->EXITG PEXIT EMDOT PRELAX 0.843019 40.000000 0.000000 b.c. data item identified --->ENDDATA boundary data successfully entered ********************************** N B L L L L L L L L L L M M N N M M N N U C B B F F D D S S 1 2 1 1 1 1 2 2 2 2 M T L L A A I I P P L L L L L L L L L L B Y O O C C R R E E I I I I I I I I I I E P C C E E 1 2 C C M M M M M M M M M M R E K K 1 2 1 2 1 2 1 2 1 2 1 2 1 2 ------ -- -- -- -- -- -- -- -- --- --- --- --- --- --- --- --- --- --- 1 INLETG 1 1 I I P P J K 1 1 1 13 1 13 1 13 1 13 2 SSIN 1 1 J J P P I K 1 1 1 25 1 13 1 25 1 13 3 FREE 1 1 J J M M I K 13 13 1 25 1 13 1 25 1 13 4 SSIN 1 1 K K P P I J 1 1 1 25 1 13 1 25 1 13 5 FREE 1 1 K K M M I J 13 13 1 25 1 13 1 25 1 13 6 EXITG 1 1 I I M M J K 25 25 1 13 1 13 1 13 1 13 ***************************************** Examining all Block Outer Boundary Conditions for Completeness ***************************************** set up block data: ****************** cell area terms calculated for grid block # 1 calculating cell volumes for grid block # 1 Mesh Cell Expansion Ratio Table ******************************** Block I Index I J K J Index I J K K Index I J K | Numbr Expans. Max Max Max Expans. Max Max Max Expans. Max Max Max | N Ratio Ratio Ratio V ----- -------- ---- ---- ---- -------- ---- ---- ---- -------- ---- ---- ---- -- 1 1.00000 15 2 5 1.00000 9 2 2 1.00000 4 2 8 * - WARNING!! ********* Solution Accuracy and Convergence May Be Adversely Affected By Large Mesh Cell Expansion Ratios --- Mesh Blocks indicated by a * in the final column Possess Coordinate Mesh Cell Expansion Ratios Which Exceed the Recommended Maximum Value of 1.300000 starting imploop 1 calling impupd convergence history ---> mesh level 1 **************************************** cycle log10 <--- max ----> log10 #ss mass mass number max error rms pts. flow flow error bl# i j k error in out ------ --------- --- --- --- --- --------- ----- -------------- ------------ 1 -6.9550 1 2 13 13 -9.2393 0 36.880405 36.880394 2 -6.8251 1 24 12 12 -9.1769 0 36.880405 36.880383 3 -6.8189 1 2 10 13 -9.1783 0 36.880394 36.880390 4 -6.6871 1 25 12 11 -9.1450 0 36.880390 36.880383 5 -6.8182 1 7 12 10 -9.1347 0 36.880394 36.880383 6 -6.6993 1 2 5 8 -9.1566 0 36.880398 36.880379 7 -6.8091 1 25 9 10 -9.1531 0 36.880398 36.880379 8 -6.8251 1 15 5 12 -9.1446 0 36.880398 36.880375 9 -6.6988 1 2 13 5 -9.1315 0 36.880398 36.880383 10 -6.8172 1 23 6 6 -9.1218 0 36.880402 36.880379 new restart file name box.restart.new number of blocks: (nblks) = 1 *** plot3d absolute output on unit 9 *** *** file box.p3dabs *** plot3d relative output on unit 8 *** *** file box.p3drel Approximate CPU Time Breakdown (Parallel Runs Only) --------------------------------------------------- TIMING INIT tinit 0.0000000000000000E+00 Initialization Time TIMING GRID tgrid 0.0000000000000000E+00 Grid Input Time TIMING COMM tcomm 0.0000000000000000E+00 Boundary Update Time TIMING SOLV tsolv 0.0000000000000000E+00 Flow Solver Time TIMING REST trest 0.0000000000000000E+00 Restart File Time TIMING OUT tout 0.0000000000000000E+00 Output file Time TIMING TOTAL tottim 0.0000000000000000E+00 processing completed execution normally terminated 3.7u 0.1s 0:06 62% 0+0k 91+25io 70pf+0w Finished at --> Wed Oct 6 16:40:41 EDT 1999