TALK=T;RUN(1,1) DISPLAY This case invloves steady, turbulent flow over a backward-facing step. The flow is 2d and incompressible. The edge of the step provides a fixed point of flow separation. This example uses the low Reynolds form of the Spalart-Allmaras turbulence model. The case is based on the NASA validation case https://turbmodels.larc.nasa.gov/backstep_val.html. The step is of height h, in a two-dimensional channel of exit width 9h, with expansion ratio of 1.125. The inlet is at a distance of 110h upstream of the step, with a fixed-pressure boundary at a distance 50h downstream of the step. The Reynolds number based on step height h is 36,000, with uniform profiles of u and ENTI specified at the inlet. The large distance upstream of the step allows the boundary layer to fully develop along the wall before reaching the step. The model predicts a reattachment length XR/h=6.4, compared with experimental measurement XR/h=7.1. //////////////////////// wall /////////////////////// ----------------------------------------------------- Pressure Inlet -------> ---------> fixed at zero ________________ ////////////////| Exit wall /| Recirculation /| <---- -----> /|____________________________________ ////////////////////////////////////// wall ENDDIS ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(Low Reynolds S-A 2D Backward Facing Step) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 400 XULAST =2. XFRAC(1)=2.704782E-05 ;XFRAC(2)=5.680043E-05 XFRAC(3)=8.952829E-05 ;XFRAC(4)=1.255289E-04 XFRAC(5)=1.651297E-04 ;XFRAC(6)=2.086905E-04 XFRAC(7)=2.566073E-04 ;XFRAC(8)=3.093159E-04 XFRAC(9)=3.672953E-04 ;XFRAC(10)=4.310726E-04 XFRAC(11)=5.012277E-04 ;XFRAC(12)=5.783984E-04 XFRAC(13)=6.632861E-04 ;XFRAC(14)=7.566625E-04 XFRAC(15)=8.593765E-04 ;XFRAC(16)=9.72362E-04 XFRAC(17)=1.096646E-03 ;XFRAC(18)=1.233358E-03 XFRAC(19)=1.383742E-03 ;XFRAC(20)=1.549164E-03 XFRAC(21)=1.731129E-03 ;XFRAC(22)=1.931289E-03 XFRAC(23)=2.151466E-03 ;XFRAC(24)=2.393661E-03 XFRAC(25)=2.660075E-03 ;XFRAC(26)=2.95313E-03 XFRAC(27)=3.275491E-03 ;XFRAC(28)=3.630088E-03 XFRAC(29)=4.020144E-03 ;XFRAC(30)=4.449207E-03 XFRAC(31)=4.921176E-03 ;XFRAC(32)=5.440341E-03 XFRAC(33)=6.011424E-03 ;XFRAC(34)=6.639614E-03 XFRAC(35)=7.330623E-03 ;XFRAC(36)=8.090734E-03 XFRAC(37)=8.926854E-03 ;XFRAC(38)=9.846589E-03 XFRAC(39)=0.010858 ;XFRAC(40)=0.011971 XFRAC(41)=0.013195 ;XFRAC(42)=0.014542 XFRAC(43)=0.016023 ;XFRAC(44)=0.017653 XFRAC(45)=0.019445 ;XFRAC(46)=0.021416 XFRAC(47)=0.023585 ;XFRAC(48)=0.025971 XFRAC(49)=0.028595 ;XFRAC(50)=0.031481 XFRAC(51)=0.034656 ;XFRAC(52)=0.038149 XFRAC(53)=0.041991 ;XFRAC(54)=0.046217 XFRAC(55)=0.050866 ;XFRAC(56)=0.05598 XFRAC(57)=0.061605 ;XFRAC(58)=0.067792 XFRAC(59)=0.074598 ;XFRAC(60)=0.082085 XFRAC(61)=0.090321 ;XFRAC(62)=0.09938 XFRAC(63)=0.109345 ;XFRAC(64)=0.120307 XFRAC(65)=0.132364 ;XFRAC(66)=0.145628 XFRAC(67)=0.160218 ;XFRAC(68)=0.176266 XFRAC(69)=0.19392 ;XFRAC(70)=0.213339 XFRAC(71)=0.2347 ;XFRAC(72)=0.258197 XFRAC(73)=0.284044 ;XFRAC(74)=0.312475 XFRAC(75)=0.34375 ;XFRAC(76)=0.375025 XFRAC(77)=0.403456 ;XFRAC(78)=0.429303 XFRAC(79)=0.4528 ;XFRAC(80)=0.474161 XFRAC(81)=0.49358 ;XFRAC(82)=0.511234 XFRAC(83)=0.527282 ;XFRAC(84)=0.541872 XFRAC(85)=0.555136 ;XFRAC(86)=0.567194 XFRAC(87)=0.578155 ;XFRAC(88)=0.58812 XFRAC(89)=0.597179 ;XFRAC(90)=0.605415 XFRAC(91)=0.612902 ;XFRAC(92)=0.619708 XFRAC(93)=0.625896 ;XFRAC(94)=0.63152 XFRAC(95)=0.636634 ;XFRAC(96)=0.641283 XFRAC(97)=0.645509 ;XFRAC(98)=0.649351 XFRAC(99)=0.652844 ;XFRAC(100)=0.656019 XFRAC(101)=0.658905 ;XFRAC(102)=0.66153 XFRAC(103)=0.663915 ;XFRAC(104)=0.666084 XFRAC(105)=0.668055 ;XFRAC(106)=0.669848 XFRAC(107)=0.671477 ;XFRAC(108)=0.672958 XFRAC(109)=0.674305 ;XFRAC(110)=0.675529 XFRAC(111)=0.676642 ;XFRAC(112)=0.677653 XFRAC(113)=0.678573 ;XFRAC(114)=0.679409 XFRAC(115)=0.680169 ;XFRAC(116)=0.68086 XFRAC(117)=0.681489 ;XFRAC(118)=0.68206 XFRAC(119)=0.682579 ;XFRAC(120)=0.683051 XFRAC(121)=0.68348 ;XFRAC(122)=0.68387 XFRAC(123)=0.684225 ;XFRAC(124)=0.684547 XFRAC(125)=0.68484 ;XFRAC(126)=0.685106 XFRAC(127)=0.685349 ;XFRAC(128)=0.685569 XFRAC(129)=0.685769 ;XFRAC(130)=0.685951 XFRAC(131)=0.686116 ;XFRAC(132)=0.686267 XFRAC(133)=0.686403 ;XFRAC(134)=0.686528 XFRAC(135)=0.686641 ;XFRAC(136)=0.686743 XFRAC(137)=0.686837 ;XFRAC(138)=0.686922 XFRAC(139)=0.686999 ;XFRAC(140)=0.687069 XFRAC(141)=0.687133 ;XFRAC(142)=0.687191 XFRAC(143)=0.687243 ;XFRAC(144)=0.687291 XFRAC(145)=0.687335 ;XFRAC(146)=0.687375 XFRAC(147)=0.687411 ;XFRAC(148)=0.687443 XFRAC(149)=0.687473 ;XFRAC(150)=0.6875 XFRAC(151)=0.687545 ;XFRAC(152)=0.68759 XFRAC(153)=0.687636 ;XFRAC(154)=0.687684 XFRAC(155)=0.687732 ;XFRAC(156)=0.687781 XFRAC(157)=0.687831 ;XFRAC(158)=0.687882 XFRAC(159)=0.687935 ;XFRAC(160)=0.687988 XFRAC(161)=0.688042 ;XFRAC(162)=0.688098 XFRAC(163)=0.688154 ;XFRAC(164)=0.688212 XFRAC(165)=0.688271 ;XFRAC(166)=0.688331 XFRAC(167)=0.688392 ;XFRAC(168)=0.688454 XFRAC(169)=0.688518 ;XFRAC(170)=0.688583 XFRAC(171)=0.688649 ;XFRAC(172)=0.688716 XFRAC(173)=0.688785 ;XFRAC(174)=0.688856 XFRAC(175)=0.688927 ;XFRAC(176)=0.689 XFRAC(177)=0.689075 ;XFRAC(178)=0.689151 XFRAC(179)=0.689228 ;XFRAC(180)=0.689308 XFRAC(181)=0.689388 ;XFRAC(182)=0.689471 XFRAC(183)=0.689555 ;XFRAC(184)=0.68964 XFRAC(185)=0.689728 ;XFRAC(186)=0.689817 XFRAC(187)=0.689908 ;XFRAC(188)=0.69 XFRAC(189)=0.690095 ;XFRAC(190)=0.690191 XFRAC(191)=0.69029 ;XFRAC(192)=0.69039 XFRAC(193)=0.690492 ;XFRAC(194)=0.690597 XFRAC(195)=0.690703 ;XFRAC(196)=0.690812 XFRAC(197)=0.690923 ;XFRAC(198)=0.691036 XFRAC(199)=0.691151 ;XFRAC(200)=0.691269 XFRAC(201)=0.691389 ;XFRAC(202)=0.691511 XFRAC(203)=0.691636 ;XFRAC(204)=0.691763 XFRAC(205)=0.691893 ;XFRAC(206)=0.692025 XFRAC(207)=0.69216 ;XFRAC(208)=0.692298 XFRAC(209)=0.692439 ;XFRAC(210)=0.692582 XFRAC(211)=0.692728 ;XFRAC(212)=0.692877 XFRAC(213)=0.693029 ;XFRAC(214)=0.693184 XFRAC(215)=0.693343 ;XFRAC(216)=0.693504 XFRAC(217)=0.693669 ;XFRAC(218)=0.693837 XFRAC(219)=0.694008 ;XFRAC(220)=0.694183 XFRAC(221)=0.694361 ;XFRAC(222)=0.694543 XFRAC(223)=0.694728 ;XFRAC(224)=0.694917 XFRAC(225)=0.69511 ;XFRAC(226)=0.695307 XFRAC(227)=0.695508 ;XFRAC(228)=0.695712 XFRAC(229)=0.695921 ;XFRAC(230)=0.696134 XFRAC(231)=0.696351 ;XFRAC(232)=0.696573 XFRAC(233)=0.696799 ;XFRAC(234)=0.697029 XFRAC(235)=0.697265 ;XFRAC(236)=0.697504 XFRAC(237)=0.697749 ;XFRAC(238)=0.697999 XFRAC(239)=0.698253 ;XFRAC(240)=0.698513 XFRAC(241)=0.698778 ;XFRAC(242)=0.699048 XFRAC(243)=0.699323 ;XFRAC(244)=0.699604 XFRAC(245)=0.699891 ;XFRAC(246)=0.700183 XFRAC(247)=0.700481 ;XFRAC(248)=0.700786 XFRAC(249)=0.701096 ;XFRAC(250)=0.701412 XFRAC(251)=0.701735 ;XFRAC(252)=0.702064 XFRAC(253)=0.7024 ;XFRAC(254)=0.702743 XFRAC(255)=0.703092 ;XFRAC(256)=0.703449 XFRAC(257)=0.703812 ;XFRAC(258)=0.704183 XFRAC(259)=0.704561 ;XFRAC(260)=0.704947 XFRAC(261)=0.705341 ;XFRAC(262)=0.705742 XFRAC(263)=0.706151 ;XFRAC(264)=0.706569 XFRAC(265)=0.706995 ;XFRAC(266)=0.707429 XFRAC(267)=0.707872 ;XFRAC(268)=0.708324 XFRAC(269)=0.708785 ;XFRAC(270)=0.709256 XFRAC(271)=0.709735 ;XFRAC(272)=0.710225 XFRAC(273)=0.710724 ;XFRAC(274)=0.711233 XFRAC(275)=0.711752 ;XFRAC(276)=0.712282 XFRAC(277)=0.712822 ;XFRAC(278)=0.713373 XFRAC(279)=0.713935 ;XFRAC(280)=0.714508 XFRAC(281)=0.715093 ;XFRAC(282)=0.715689 XFRAC(283)=0.716297 ;XFRAC(284)=0.716918 XFRAC(285)=0.717551 ;XFRAC(286)=0.718196 XFRAC(287)=0.718855 ;XFRAC(288)=0.719527 XFRAC(289)=0.720212 ;XFRAC(290)=0.72091 XFRAC(291)=0.721623 ;XFRAC(292)=0.72235 XFRAC(293)=0.723092 ;XFRAC(294)=0.723848 XFRAC(295)=0.72462 ;XFRAC(296)=0.725407 XFRAC(297)=0.726209 ;XFRAC(298)=0.727028 XFRAC(299)=0.727863 ;XFRAC(300)=0.728715 XFRAC(301)=0.729584 ;XFRAC(302)=0.73047 XFRAC(303)=0.731374 ;XFRAC(304)=0.732296 XFRAC(305)=0.733237 ;XFRAC(306)=0.734196 XFRAC(307)=0.735174 ;XFRAC(308)=0.736172 XFRAC(309)=0.73719 ;XFRAC(310)=0.738229 XFRAC(311)=0.739288 ;XFRAC(312)=0.740368 XFRAC(313)=0.74147 ;XFRAC(314)=0.742594 XFRAC(315)=0.74374 ;XFRAC(316)=0.74491 XFRAC(317)=0.746103 ;XFRAC(318)=0.747319 XFRAC(319)=0.74856 ;XFRAC(320)=0.749826 XFRAC(321)=0.751117 ;XFRAC(322)=0.752434 XFRAC(323)=0.753777 ;XFRAC(324)=0.755147 XFRAC(325)=0.756545 ;XFRAC(326)=0.75797 XFRAC(327)=0.759424 ;XFRAC(328)=0.760907 XFRAC(329)=0.76242 ;XFRAC(330)=0.763963 XFRAC(331)=0.765537 ;XFRAC(332)=0.767142 XFRAC(333)=0.768779 ;XFRAC(334)=0.77045 XFRAC(335)=0.772153 ;XFRAC(336)=0.773891 XFRAC(337)=0.775663 ;XFRAC(338)=0.777471 XFRAC(339)=0.779315 ;XFRAC(340)=0.781196 XFRAC(341)=0.783114 ;XFRAC(342)=0.785071 XFRAC(343)=0.787067 ;XFRAC(344)=0.789103 XFRAC(345)=0.791179 ;XFRAC(346)=0.793298 XFRAC(347)=0.795458 ;XFRAC(348)=0.797662 XFRAC(349)=0.79991 ;XFRAC(350)=0.802202 XFRAC(351)=0.804541 ;XFRAC(352)=0.806926 XFRAC(353)=0.809359 ;XFRAC(354)=0.811841 XFRAC(355)=0.814373 ;XFRAC(356)=0.816955 XFRAC(357)=0.819588 ;XFRAC(358)=0.822275 XFRAC(359)=0.825015 ;XFRAC(360)=0.827809 XFRAC(361)=0.83066 ;XFRAC(362)=0.833568 XFRAC(363)=0.836534 ;XFRAC(364)=0.839559 XFRAC(365)=0.842645 ;XFRAC(366)=0.845792 XFRAC(367)=0.849003 ;XFRAC(368)=0.852277 XFRAC(369)=0.855617 ;XFRAC(370)=0.859024 XFRAC(371)=0.862499 ;XFRAC(372)=0.866044 XFRAC(373)=0.869659 ;XFRAC(374)=0.873347 XFRAC(375)=0.877109 ;XFRAC(376)=0.880945 XFRAC(377)=0.884859 ;XFRAC(378)=0.88885 XFRAC(379)=0.892922 ;XFRAC(380)=0.897075 XFRAC(381)=0.901311 ;XFRAC(382)=0.905632 XFRAC(383)=0.910039 ;XFRAC(384)=0.914534 XFRAC(385)=0.91912 ;XFRAC(386)=0.923797 XFRAC(387)=0.928567 ;XFRAC(388)=0.933433 XFRAC(389)=0.938396 ;XFRAC(390)=0.943459 XFRAC(391)=0.948622 ;XFRAC(392)=0.953889 XFRAC(393)=0.959262 ;XFRAC(394)=0.964742 XFRAC(395)=0.970331 ;XFRAC(396)=0.976032 XFRAC(397)=0.981847 ;XFRAC(398)=0.987779 XFRAC(399)=0.993829 ;XFRAC(400)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 265 YVLAST =0.1125 YFRAC(1)=7.7278E-05 ;YFRAC(2)=1.622838E-04 YFRAC(3)=2.557902E-04 ;YFRAC(4)=3.586472E-04 YFRAC(5)=4.717899E-04 ;YFRAC(6)=5.96247E-04 YFRAC(7)=7.331496E-04 ;YFRAC(8)=8.837427E-04 YFRAC(9)=1.049395E-03 ;YFRAC(10)=1.231613E-03 YFRAC(11)=1.432052E-03 ;YFRAC(12)=1.652535E-03 YFRAC(13)=1.895066E-03 ;YFRAC(14)=2.161851E-03 YFRAC(15)=2.455314E-03 ;YFRAC(16)=2.778124E-03 YFRAC(17)=3.133214E-03 ;YFRAC(18)=3.523814E-03 YFRAC(19)=3.953473E-03 ;YFRAC(20)=4.426098E-03 YFRAC(21)=4.945986E-03 ;YFRAC(22)=5.517862E-03 YFRAC(23)=6.146927E-03 ;YFRAC(24)=6.838899E-03 YFRAC(25)=7.600067E-03 ;YFRAC(26)=8.437351E-03 YFRAC(27)=9.358364E-03 ;YFRAC(28)=0.010371 YFRAC(29)=0.011486 ;YFRAC(30)=0.012712 YFRAC(31)=0.01406 ;YFRAC(32)=0.015544 YFRAC(33)=0.017175 ;YFRAC(34)=0.01897 YFRAC(35)=0.020944 ;YFRAC(36)=0.023116 YFRAC(37)=0.025505 ;YFRAC(38)=0.028133 YFRAC(39)=0.031023 ;YFRAC(40)=0.034203 YFRAC(41)=0.0377 ;YFRAC(42)=0.041548 YFRAC(43)=0.04578 ;YFRAC(44)=0.050435 YFRAC(45)=0.055556 ;YFRAC(46)=0.060676 YFRAC(47)=0.065332 ;YFRAC(48)=0.069564 YFRAC(49)=0.073411 ;YFRAC(50)=0.076908 YFRAC(51)=0.080088 ;YFRAC(52)=0.082979 YFRAC(53)=0.085606 ;YFRAC(54)=0.087995 YFRAC(55)=0.090167 ;YFRAC(56)=0.092141 YFRAC(57)=0.093936 ;YFRAC(58)=0.095568 YFRAC(59)=0.097051 ;YFRAC(60)=0.098399 YFRAC(61)=0.099625 ;YFRAC(62)=0.10074 YFRAC(63)=0.101753 ;YFRAC(64)=0.102674 YFRAC(65)=0.103511 ;YFRAC(66)=0.104272 YFRAC(67)=0.104964 ;YFRAC(68)=0.105593 YFRAC(69)=0.106165 ;YFRAC(70)=0.106685 YFRAC(71)=0.107158 ;YFRAC(72)=0.107587 YFRAC(73)=0.107978 ;YFRAC(74)=0.108333 YFRAC(75)=0.108656 ;YFRAC(76)=0.108949 YFRAC(77)=0.109216 ;YFRAC(78)=0.109459 YFRAC(79)=0.109679 ;YFRAC(80)=0.10988 YFRAC(81)=0.110062 ;YFRAC(82)=0.110227 YFRAC(83)=0.110378 ;YFRAC(84)=0.110515 YFRAC(85)=0.110639 ;YFRAC(86)=0.110752 YFRAC(87)=0.110855 ;YFRAC(88)=0.110949 YFRAC(89)=0.111034 ;YFRAC(90)=0.111111 YFRAC(91)=0.111122 ;YFRAC(92)=0.111133 YFRAC(93)=0.111146 ;YFRAC(94)=0.11116 YFRAC(95)=0.111176 ;YFRAC(96)=0.111193 YFRAC(97)=0.111212 ;YFRAC(98)=0.111232 YFRAC(99)=0.111255 ;YFRAC(100)=0.11128 YFRAC(101)=0.111308 ;YFRAC(102)=0.111338 YFRAC(103)=0.111371 ;YFRAC(104)=0.111408 YFRAC(105)=0.111448 ;YFRAC(106)=0.111492 YFRAC(107)=0.111541 ;YFRAC(108)=0.111595 YFRAC(109)=0.111654 ;YFRAC(110)=0.111719 YFRAC(111)=0.11179 ;YFRAC(112)=0.111868 YFRAC(113)=0.111955 ;YFRAC(114)=0.11205 YFRAC(115)=0.112154 ;YFRAC(116)=0.112269 YFRAC(117)=0.112396 ;YFRAC(118)=0.112535 YFRAC(119)=0.112688 ;YFRAC(120)=0.112856 YFRAC(121)=0.113041 ;YFRAC(122)=0.113244 YFRAC(123)=0.113468 ;YFRAC(124)=0.113715 YFRAC(125)=0.113986 ;YFRAC(126)=0.114284 YFRAC(127)=0.114612 ;YFRAC(128)=0.114972 YFRAC(129)=0.115369 ;YFRAC(130)=0.115805 YFRAC(131)=0.116285 ;YFRAC(132)=0.116814 YFRAC(133)=0.117394 ;YFRAC(134)=0.118033 YFRAC(135)=0.118736 ;YFRAC(136)=0.119509 YFRAC(137)=0.12036 ;YFRAC(138)=0.121295 YFRAC(139)=0.122324 ;YFRAC(140)=0.123456 YFRAC(141)=0.124701 ;YFRAC(142)=0.126071 YFRAC(143)=0.127577 ;YFRAC(144)=0.129234 YFRAC(145)=0.131057 ;YFRAC(146)=0.133063 YFRAC(147)=0.135268 ;YFRAC(148)=0.137695 YFRAC(149)=0.140364 ;YFRAC(150)=0.1433 YFRAC(151)=0.146529 ;YFRAC(152)=0.150081 YFRAC(153)=0.153989 ;YFRAC(154)=0.158288 YFRAC(155)=0.163016 ;YFRAC(156)=0.168217 YFRAC(157)=0.173938 ;YFRAC(158)=0.180231 YFRAC(159)=0.187154 ;YFRAC(160)=0.194769 YFRAC(161)=0.203145 ;YFRAC(162)=0.212359 YFRAC(163)=0.222495 ;YFRAC(164)=0.233644 YFRAC(165)=0.245907 ;YFRAC(166)=0.259398 YFRAC(167)=0.274237 ;YFRAC(168)=0.29056 YFRAC(169)=0.308516 ;YFRAC(170)=0.328267 YFRAC(171)=0.349993 ;YFRAC(172)=0.373892 YFRAC(173)=0.40018 ;YFRAC(174)=0.429098 YFRAC(175)=0.460907 ;YFRAC(176)=0.495897 YFRAC(177)=0.534387 ;YFRAC(178)=0.576725 YFRAC(179)=0.615214 ;YFRAC(180)=0.650204 YFRAC(181)=0.682013 ;YFRAC(182)=0.710931 YFRAC(183)=0.73722 ;YFRAC(184)=0.761118 YFRAC(185)=0.782844 ;YFRAC(186)=0.802595 YFRAC(187)=0.820551 ;YFRAC(188)=0.836874 YFRAC(189)=0.851713 ;YFRAC(190)=0.865204 YFRAC(191)=0.877468 ;YFRAC(192)=0.888616 YFRAC(193)=0.898752 ;YFRAC(194)=0.907966 YFRAC(195)=0.916342 ;YFRAC(196)=0.923957 YFRAC(197)=0.93088 ;YFRAC(198)=0.937173 YFRAC(199)=0.942894 ;YFRAC(200)=0.948095 YFRAC(201)=0.952824 ;YFRAC(202)=0.957122 YFRAC(203)=0.96103 ;YFRAC(204)=0.964582 YFRAC(205)=0.967811 ;YFRAC(206)=0.970747 YFRAC(207)=0.973416 ;YFRAC(208)=0.975843 YFRAC(209)=0.978048 ;YFRAC(210)=0.980054 YFRAC(211)=0.981877 ;YFRAC(212)=0.983534 YFRAC(213)=0.98504 ;YFRAC(214)=0.98641 YFRAC(215)=0.987655 ;YFRAC(216)=0.988787 YFRAC(217)=0.989816 ;YFRAC(218)=0.990752 YFRAC(219)=0.991602 ;YFRAC(220)=0.992375 YFRAC(221)=0.993078 ;YFRAC(222)=0.993717 YFRAC(223)=0.994298 ;YFRAC(224)=0.994826 YFRAC(225)=0.995306 ;YFRAC(226)=0.995742 YFRAC(227)=0.996139 ;YFRAC(228)=0.996499 YFRAC(229)=0.996827 ;YFRAC(230)=0.997125 YFRAC(231)=0.997396 ;YFRAC(232)=0.997643 YFRAC(233)=0.997867 ;YFRAC(234)=0.99807 YFRAC(235)=0.998255 ;YFRAC(236)=0.998423 YFRAC(237)=0.998576 ;YFRAC(238)=0.998715 YFRAC(239)=0.998842 ;YFRAC(240)=0.998957 YFRAC(241)=0.999061 ;YFRAC(242)=0.999156 YFRAC(243)=0.999243 ;YFRAC(244)=0.999321 YFRAC(245)=0.999393 ;YFRAC(246)=0.999457 YFRAC(247)=0.999516 ;YFRAC(248)=0.99957 YFRAC(249)=0.999619 ;YFRAC(250)=0.999663 YFRAC(251)=0.999703 ;YFRAC(252)=0.99974 YFRAC(253)=0.999773 ;YFRAC(254)=0.999803 YFRAC(255)=0.999831 ;YFRAC(256)=0.999856 YFRAC(257)=0.999879 ;YFRAC(258)=0.999899 YFRAC(259)=0.999918 ;YFRAC(260)=0.999935 YFRAC(261)=0.999951 ;YFRAC(262)=0.999965 YFRAC(263)=0.999978 ;YFRAC(264)=0.999989 YFRAC(265)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 1 ZWLAST =0.1 ZFRAC(1)=1. ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(1)=P1 ;NAME(3)=U1 NAME(5)=V1 ;NAME(137)=PRPS NAME(138)=STRS ;NAME(139)=DVDX NAME(140)=DUDY ;NAME(141)=CP NAME(142)=WDIS ;NAME(143)=VOR1 NAME(144)=CF ;NAME(145)=YPLS NAME(146)=LTLS ;NAME(147)=ENTI NAME(148)=DEN1 ;NAME(149)=EL1 NAME(150)=ENUT * Y in SOLUTN argument list denotes: * 1-stored 2-solved 3-whole-field * 4-point-by-point 5-explicit 6-harmonic averaging SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(U1,Y,Y,Y,N,N,Y) SOLUTN(V1,Y,Y,Y,N,N,Y) SOLUTN(PRPS,Y,N,N,N,N,Y) SOLUTN(STRS,Y,N,N,N,N,Y) SOLUTN(DVDX,Y,N,N,N,N,Y) SOLUTN(DUDY,Y,N,N,N,N,Y) SOLUTN(CP,Y,N,N,N,N,Y) SOLUTN(WDIS,Y,N,N,N,N,Y) SOLUTN(VOR1,Y,N,N,N,N,Y) SOLUTN(CF,Y,N,N,N,N,Y) SOLUTN(YPLS,Y,N,N,N,N,Y) SOLUTN(LTLS,Y,Y,Y,N,N,Y) SOLUTN(ENTI,Y,Y,Y,N,N,Y) SOLUTN(DEN1,Y,N,N,N,N,Y) SOLUTN(EL1,Y,N,N,N,N,Y) SOLUTN(ENUT,Y,N,N,N,N,Y) DEN1 = 148 VIST = 150 LEN1 = 149 PRPS = 137 ************************************************************ Group 8. Terms & Devices * Y in TERMS argument list denotes: * 1-built-in source 2-convection 3-diffusion 4-transient * 5-first phase variable 6-interphase transport TERMS(P1,Y,Y,Y,N,Y,Y) TERMS(U1,Y,Y,Y,Y,Y,Y) TERMS(V1,Y,Y,Y,Y,Y,Y) TERMS(LTLS,N,N,Y,N,Y,Y) TERMS(ENTI,N,Y,Y,Y,N,Y) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = F NEWENT = F ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1 ************************************************************ Group 9. Properties used if PRPS is not stored, and where PRPS = -1.0 if it is! RHO1 =1.189 ;TMP1 = GRND1 ;EL1 =0. TSURR =0. ;TEMP0 =0. ;PRESS0 =1.01325E+05 DVO1DT =3.41E-03 ;DRH1DP =0. TMP1A =300. ;TMP1B =0. ;TMP1C =0. TMP2A =0. ;TMP2B =0. ;TMP2C =0. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. ENUL =1.544E-05 ;ENUT = GRND2 ENUTA =0. ;ENUTB =0. ;ENUTC =0. IENUTA = 22 PRNDTL(U1)=1. ;PRNDTL(V1)=1. PRNDTL(LTLS)=1. ;PRNDTL(ENTI)=0.666667 PRT(U1)=1. ;PRT(V1)=1. PRT(LTLS)=1. ;PRT(ENTI)=0.666667 CP1 =1005. ;CP2 =1. ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1)=0. ;FIINIT(U1)=44.400002 FIINIT(V1)=0.1 ;FIINIT(PRPS)=-1. FIINIT(STRS)=1.0E-10 ;FIINIT(DVDX)=1.0E-10 FIINIT(DUDY)=1.0E-10 ;FIINIT(CP)=1.0E-10 FIINIT(WDIS)=1.0E-02 ;FIINIT(VOR1)=1.0E-10 FIINIT(CF)=1.0E-10 ;FIINIT(YPLS)=1.0E-10 FIINIT(LTLS)=1.0E-10 ;FIINIT(ENTI)=1.0E-03 FIINIT(DEN1)=1.0E-10 ;FIINIT(EL1)=1.0E-10 FIINIT(ENUT)=1.0E-10 PATCH(^OB1 ,INIVAL, 1, 0, 0, 0, 0, 0, 1, 1) INIT(^OB1 ,PRPS,0. ,198. ) INIADD = F FSWEEP = 1 NAMFI =CHAM ************************************************************ Group 12. Patchwise adjustment of terms Patches for this group are printed with those for Group 13. Their names begin either with GP12 or & ************************************************************ Group 13. Boundary & Special Sources PATCH(SASOPRO ,PHASEM, 1, 400, 1, 265, 1, 1, 1, 1) COVAL(SASOPRO ,ENTI, FIXFLU , GRND4 ) PATCH(SASODIS ,PHASEM, 1, 400, 1, 265, 1, 1, 1, 1) COVAL(SASODIS ,ENTI, GRND4 ,0. ) PATCH(SASODIF ,PHASEM, 1, 400, 1, 265, 1, 1, 1, 1) COVAL(SASODIF ,ENTI, FIXFLU , GRND4 ) PATCH(^OB2 ,WEST , 2, 0, 0, 0, 0, 0, 1, 1) COVAL(^OB2 ,P1 , FIXFLU ,52.791603 ) COVAL(^OB2 ,U1 ,0. ,44.400002 ) COVAL(^OB2 ,V1 ,0. ,0. ) COVAL(^OB2 ,ENTI,0. ,7.720001E-05 ) PATCH(^OB3 ,EAST , 3, 0, 0, 0, 0, 0, 1, 1) COVAL(^OB3 ,P1 ,1000. ,0. ) COVAL(^OB3 ,U1 ,0. ,0. ) COVAL(^OB3 ,V1 ,0. ,0. ) COVAL(^OB3 ,ENTI,0. ,7.72E-05 ) PATCH(^OB4 ,CELL , 4, 7, 0, 0, 0, 0, 1, 1) COVAL(^OB4 ,U1 , GRND2 ,0. ) COVAL(^OB4 ,ENTI,1. ,0. ) PATCH(^OB5 ,CELL , 5, 8, 0, 0, 0, 0, 1, 1) COVAL(^OB5 ,U1 , GRND2 ,0. ) COVAL(^OB5 ,ENTI,1. ,0. ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 3000 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC =1.0E-04 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(U1)=10 LITER(V1)=10 ;LITER(LTLS)=20 LITER(ENTI)=20 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03 ENDIT(V1)=1.0E-03 ;ENDIT(LTLS)=1.0E-03 ENDIT(ENTI)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(U1,FALSDT,1.0E-03) RELAX(V1,FALSDT,1.0E-03) RELAX(PRPS,LINRLX,1.) RELAX(STRS,LINRLX,1.) RELAX(DVDX,LINRLX,1.) RELAX(DUDY,LINRLX,1.) RELAX(CP,LINRLX,1.) RELAX(WDIS,LINRLX,1.) RELAX(VOR1,LINRLX,1.) RELAX(CF,LINRLX,1.) RELAX(YPLS,LINRLX,1.) RELAX(LTLS,LINRLX,1.) RELAX(ENTI,LINRLX,0.3) RELAX(DEN1,LINRLX,1.) RELAX(EL1,LINRLX,1.) RELAX(ENUT,LINRLX,1.) OVRRLX =0. EXPERT = F ;NNORSL = F ************************************************************ Group 18. Limits VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10 VARMAX(STRS)=1.0E+10 ;VARMIN(STRS)=-1.0E+10 VARMAX(DVDX)=1.0E+10 ;VARMIN(DVDX)=-1.0E+10 VARMAX(DUDY)=1.0E+10 ;VARMIN(DUDY)=-1.0E+10 VARMAX(CP)=1.0E+10 ;VARMIN(CP)=-1.0E+10 VARMAX(WDIS)=1.0E+10 ;VARMIN(WDIS)=-1.0E+10 VARMAX(VOR1)=1.0E+10 ;VARMIN(VOR1)=-1.0E+10 VARMAX(CF)=1.0E+10 ;VARMIN(CF)=-1.0E+10 VARMAX(YPLS)=1.0E+10 ;VARMIN(YPLS)=-1.0E+10 VARMAX(LTLS)=1.0E+10 ;VARMIN(LTLS)=-1.0E+10 VARMAX(ENTI)=1.0E+10 ;VARMIN(ENTI)=1.0E-10 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10 VARMAX(EL1)=1.0E+10 ;VARMIN(EL1)=-1.0E+10 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10 ************************************************************ Group 19. Data transmitted to GROUND GENK = T YPLS = T PARSOL = F IENUTA = 22 ISG62 = 1 SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,0) SPEDAT(SET,MAKE,UINF,C,=0.0!IF(ISWEEP.EQ.1)) SPEDAT(SET,STORED,UINF,C,=ABS(U1{1.325&0.05625&0.05})!SWPFIN) SPEDAT(SET,PRINT,UINF,C,=UINF) SPEDAT(SET,STORED,CF,C,=2*STRS*(ABS(U1)/U1)/(UINF^2)!SWPFIN) SPEDAT(SET,STORED,CP,C,=P1/(0.5*DEN1*UINF^2)!SWPFIN) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,OBJNAM,^OB1,C,STEP) SPEDAT(SET,OBJTYP,^OB1,C,BLOCKAGE) SPEDAT(SET,STEP,MATERIAL,R,198.) SPEDAT(SET,OBJNAM,^OB2,C,INLET) SPEDAT(SET,OBJTYP,^OB2,C,INLET) SPEDAT(SET,OBJNAM,^OB3,C,OUT) SPEDAT(SET,OBJTYP,^OB3,C,OUTLET) SPEDAT(SET,ARATIO,^OB3,R,1.) SPEDAT(SET,OBJNAM,^OB4,C,WALL_UPPER) SPEDAT(SET,OBJTYP,^OB4,C,PLATE) SPEDAT(SET,OBJNAM,^OB5,C,WALL_LOWER) SPEDAT(SET,OBJTYP,^OB5,C,PLATE) SPEDAT(SET,FACETDAT,NUMOBJ,I,5) SPEDAT(SET,MATERIAL,198,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=107.699997 ;EX(U1)=22.99 EX(V1)=0.4821 ;EX(PRPS)=0.8726 EX(STRS)=0.03026 ;EX(DVDX)=68.970001 EX(DUDY)=3.582E+04 ;EX(CP)=0.08241 EX(WDIS)=7.038E-03 ;EX(VOR1)=3.587E+04 EX(CF)=2.749E-05 ;EX(YPLS)=1.673E-03 EX(LTLS)=2.953E-04 ;EX(ENTI)=1.327E-03 EX(DEN1)=1.038 ;EX(EL1)=0. EX(ENUT)=1.313E-03 ************************************************************ Group 21. Print-out of Variables INIFLD = F ;SUBWGR = F * Y in OUTPUT argument list denotes: * 1-field 2-correction-eq. monitor 3-selective dumping * 4-whole-field residual 5-spot-value table 6-residual table OUTPUT(P1,N,N,Y,Y,Y,Y) OUTPUT(U1,N,N,Y,Y,Y,Y) OUTPUT(V1,N,N,Y,Y,Y,Y) OUTPUT(PRPS,N,N,Y,N,N,N) OUTPUT(STRS,N,N,Y,N,N,N) OUTPUT(DVDX,N,N,Y,N,N,N) OUTPUT(DUDY,N,N,Y,N,N,N) OUTPUT(CP,N,N,Y,N,N,N) OUTPUT(WDIS,N,N,Y,N,N,N) OUTPUT(VOR1,N,N,Y,N,N,N) OUTPUT(CF,N,N,Y,N,N,N) OUTPUT(YPLS,N,N,Y,N,N,N) OUTPUT(LTLS,N,N,Y,Y,Y,Y) OUTPUT(ENTI,N,N,Y,Y,Y,Y) OUTPUT(DEN1,N,N,Y,N,N,N) OUTPUT(EL1,N,N,Y,N,N,N) OUTPUT(ENUT,N,N,Y,N,N,N) WALPRN = T ************************************************************ Group 22. Monitor Print-Out IXMON = 379 ;IYMON = 177 ;IZMON = 1 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1 UWATCH = T ;USTEER = T HIGHLO = F ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 ;NUMCLS = 5 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 3 ;IPROF = 1 ABSIZ =0.5 ;ORSIZ =0.4 NTZPRF = 1 ;NCOLPF = 50 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20 No PATCHes yet used for this Group ************************************************************ Group 24. Dumps For Restarts SAVE = T ;NOWIPE = F NSAVE =CHAM STOP