Use of Proper Closure Equations in Finite Volume Discretization Schemes


 Homer Campbell
 1 years ago
 Views:
Transcription
1 Us of ropr Closur Equations in Finit Volum Discrtization Schms 3 A. Ashrafizadh, M. Rzvani and B. Alinia K. N. Toosi Univrsity of Tchnology Iran 1. Introduction Finit Volum Mthod (FVM) is a popular fild discrtization approach for th numrical simulation of physical procsss dscribd by consrvation laws. This articl prsnts som advancs in finit volum discrtization of hat and fluid flow problms. In th cas of thrmofluid problms, dcision rgarding th choic of unknown variabls nds to b takn at th continuous lvl of formulation. A popular choic, among many availabl options, is th prssurbasd primitiv variabl formulation (Acharya, 2007), in which th flow continuity quation is considrd as a constraint on fluid prssur. rssurbasd finit volum mthods, in turn, can b catgorizd basd on th computational grid arrangmnt usd for th discrtization of solution domain. In a staggrd grid arrangmnt, vlocity componnts and othr unknown scalars ar dfind at diffrnt nodal points (atankar, 1980). On th contrary, all unknown variabls ar dfind at th sam nodal locations in a colocatd grid arrangmnt. Th discussion of finit volum discrtization schms in this articl is limitd to prssurbasd formulations on colocatd grids. At th continuous lvl of formulation in a finit volum mthod, an intnsiv variabl is constraind by an intgral balanc quation for a control volum CV boundd by th control surfac CS as follows: In Eq. (1) J t dv J nda S dv (1) CV CS CV rprsnts th flux of across th control surfac, S stands for th rat of gnration of within th control volum and n is th outward unit normal vctor to th control surfac. In most of th nginring problms of intrst, thr ar two mchanisms for th transport of, i.. advction and diffusion, and th flux function is mathmatically dscribd as follows: J V In Eq. (2) V is th advctiv mass flux and is th diffusion cofficint associatd with variabl. By proprly dfining, J and S, Eq. (1) can b usd as a gnric quation to (2)
2 62 Finit Volum Mthod owrful Mans of Enginring Dsign mathmatically dscrib th transport of mass, momntum and nrgy in fluid flow problms (atankar, 1980). A typical two dimnsional discrt domain covrd by contiguous control volums (clls) is shown in Fig. 1. Each cll rprsnts a nodal point and is boundd by a numbr of panls which compris th control surfac. Typically, ach panl includs on intgration point. Thr isn't any limitation on th shap and th numbr of panls of a cll. Hr w us a background finit lmnt msh to dfin th control volums. Figur 1 shows structurd and unstructurd twodimnsional grids in an lmntbasd finit volum mthod (EBFVM) contxt (Ashrafizadh t al., 2011). Th objctiv of any finit volum discrtization schm is to us th CV balanc quations, similar to Eq. (1), to obtain a st of algbraic quations which constrains th unknown nodal valus throughout th domain. This objctiv is achivd by carrying out th discrtization in two stags. At th first stag of discrtization, hr calld lvl1 approximation, Eq. (1) is convrtd to a smidiscrt form comprising of both nodal and intgration point variabls. Lvl1 approximat form of Eq. (1) for an intrnal finit volum in Fig. 1 can b writtn as follows: ( ) t np N N ip ip V ( V nda ) ( A ) ip n S V (3) ip ip1 ip1 In Eq. (3), subscripts np and ip stand for nodal and intgration point variabls rspctivly. Closur quations ar now ndd to rlat ip variabls to np variabls. At th scond stag of discrtization, hr calld lvl2 approximation, Eq. (3) is convrtd to a fullydiscrt form, i.. an algbraic quation which constrains th variabl at nodal point. Th computational molcul obtaind at this stag of discrtization is writtn in th following convnint form: N nb a a S (4) nb nb nb1 Hr, subscript nb stands for th nighbor nods and suprscript N nb rfrs to th numbr of influntial nighbor nods. Influnc cofficints a and a nb and th sourc trm S contain trms and paramtrs which modl all physical ffcts rlvant to th valu of at th nodal point. If Eq. (1) is a nonlinar quation, linarization also bcoms ncssary and should b carrid out during th discrtization. In principl, th quation can b linarizd bfor or aftr ach of th approximation lvls and this might affct th proprtis of th final discrt modl. In most nginring problms, including hat and fluid flow ons, a coupld st of nonlinar partial diffrntial quations nds to b discrtizd and solvd. In such cass, th intrquation couplings should also b numrically modld. Twodimnsional, incomprssibl NavirStoks quations, for xampl, can b writtn in a form similar to Eq. (1) in which is qual to 1 in th mass balanc quation, u in th xmomntum balanc quation, v in th ymomntum balanc quation and (spcific nrgy) in th nrgy
3 Us of ropr Closur Equations in Finit Volum Discrtization Schms 63 Fig. 1. Structurd and unstructurd lmntbasd finit volum grids. balanc quation. In gnral, th computational molcul associatd with ach quation of a coupld st includs mor than on variabl to rflct th fact that th quations and th variabls ar intrrlatd. Thrfor, in th cas of a st of quations with two coupld variabls and, th following modifid form of Eq. (4) is usd to constrain at th nodal point : Nnb Nnb nbnb nbnb nb1 nb1 (5) a a a a S In this articl, w focus on lvl2 approximation of NavirStoks quations in th contxt of colocatd, prssurbasd finit volum mthod. This lvl of approximation is particularly important bcaus it dals with th modling of physical influncs such as upwind and
4 64 Finit Volum Mthod owrful Mans of Enginring Dsign downwind ffcts and th couplings btwn variabls. A numbr of famous lvl2 approximations ar rviwd and th mthod of propr closur quations (MCE), proposd by th first author, is introducd. To avoid th complxitis associatd with multidimnsional problms and to focus on th natur of th approximations in a simpl stting, on dimnsional quations and tst cass ar discussd in dtails. Twodimnsional rsults ar only brifly prsntd to show th applicability of th MCE in multidimnsional problms. 2. Th smidiscrt form of quations in a 1D contxt A ondimnsional variablara duct is shown in Fig. 2. Th complt st of govrning quations for unstady 1Dinviscid comprssibl flow in this duct, using th momntum variabl formulation, is as follows (Chtrntal, 1999): 1 da m f S f adx 1 da u f fu p fu S x a dx f pu 1 da pu adx Whr is tim, f u is th momntum variabl, u is vlocity, p is prssur and a is th cross sctional ara. Not that no artificial nrgy sourc trm is usd in th nrgy quation. In ordr to clos th abov st of quations, an auxiliary quation is ndd. Hr, th idal gas quation of stat is usd: (6) p Rt (7) Whr t is tmpratur and R is th gas constant. Nglcting th potntial nrgy trm, th total nrgy pr unit mass ( ) of an idal gas can b writtn as follows: Whr cv is th spcific hat at constant volum. 1 2 cvt u (8) 2 By intgrating Eq. (6) ovr th control volum in Fig. 2 and aftr using divrgnc thorm, th smidiscrt forms of th govrning quations ar obtaind as follows: m V SV fa fwaw V p F F fu p w ( w a) a fu p w a a a S V V EE f pu a fu a w w 0 (9) Not that capital lttrs (lik F and ) ar usd for th nodal valus and small lttrs (lik f and p ) ar usd for intgration point variabls. Hr, V stands for th volum of cll. Closur quations for th intgration point variabls ar discussd nxt.
5 Us of ropr Closur Equations in Finit Volum Discrtization Schms 65 Fig. 2. Grid arrangmnt in a ondimnsional variabl ara duct. 3. Closur quations for intgration point variabls Aftr obtaining th smidiscrt form of govrning quations and linarizing thm, closur quations ar ndd to approximat intgration point (ip) quantitis in trms of th nodal valus. By mploying ths closur quations, th smidiscrt quations ar convrtd to fullydiscrt quations. Ths ipquations play a critical rol in th robustnss and accuracy of a collocatd schm. Thrfor, diffrnt solution stratgis for dfining ipquation ar brifly dscribd in this sction. Hr, w only discuss th closur quations for ast intgration point of a typical cll around nod. Th closur quations at th wst ip ar obtaind similarly. 3.1 Mthods basd on gomtrical intrpolation Simpl symmtric intrpolation is th most trivial choic for an ipvalu in a uniform grid: On nonuniform grids, wightd intrpolation is usd as follows: E (10) 2 (1 ) (11) E Whr is th wight factor for ast ip, calculatd by th following formula: x X E X X (12)
6 66 Finit Volum Mthod owrful Mans of Enginring Dsign This schm is 2ndordr accurat, but is unboundd so that nonphysical oscillations may appar in rgions with high gradints. To show th problm associatd with this schm, considr th 1D inviscid incomprssibl flow in a constant ara duct shown in Fig. 3. Th smidiscrt form of continuity and momntum quations in this cas ar as follows: u 0 Au (13) 0 w mu u Ap p (14) w w Whr A is th cross sctional ara of th duct and m is th mass flow rat basd on th most rcnt availabl valu of vlocity. Fig. 3. A 1D Colocatd grid in a constant cross sctional ara duct. Using simpl symmtric intrpolation, Eqs. (13) and (14) ar convrtd to th following fullydiscrt forms: A U E U W 0 (15) 2 m A 2 2 U U 0 E W E W Combining th abov quations, on obtains th following constraint on th prssur at nodal point : E W (16) (17) Not that prssur at nod is absnt in this quation and this constraint cannot distinguish btwn a uniform prssur fild and th wiggly prssur fild shown in Fig Mthods basd on taylorsris xpansion It is logical to assum that convctd quantitis in a thrmofluid problm ar influncd by th upstram condition. Thrfor, assuming th flow dirction from to E, th following truncatd Taylor sris xpansions ar both valid approximations for : (18) φ x x X (19)
7 Us of ropr Closur Equations in Finit Volum Discrtization Schms 67 Fig. 4. A typical wiggly prssur fild. Equation (18) is calld th simpl or Fist Ordr Upwinding (FOU) and provids highly stabl and unconditionally boundd numrical schms. Howvr, it is only first ordr accurat. Equation (19) provids a Scond Ordr Upwinding (SOU) stratgy. On can show that a numrical schm which mploys Eq. (10) for intgration point prssurs and Eqs. (18) or (19) for intgration point vlocitis in Eqs. (13) and (14), would not b abl to dtct th chckrboard prssur filds. Such numrical schms ar pron to unphysical prssurvlocity dcoupling symptom. 3.3 Mthods basd on momntum intrpolation Rhi and Chow (Rhi, 1983) mployd a colocatd grid to solv th flow fild around an airfoil. To avoid nonphysical, wiggly numrical solutions, thy cam to th conclusion that diffrnt closur quations had to b usd for th convcting or masscarrying ( û ) and convctd or transportd ( u ) vlocity componnts at th intgration points. Th convcting vlocity, usd in th smidiscrt continuity quation, is linkd to th prssur through a momntum intrpolation procdur. Hnc, this discrtization schm is calld th Momntum Intrpolation Mthod (MIM). As compard to th classical schms on staggrd grids, on may assum that nodal momntum balancs ar concptually staggrd to obtain closur quations for convcting vlocity at intgration points. Th discrt momntum balanc for nod can b writtn as follows: bnbunb V au anbunb V 0 U (20) x a a x Th convcting vlocity at is obtaind by staggring th stncil of th computational molcul to th intgration point: Which rsults in: nb nb ˆ u a U V a a x (21) U ˆ UE d u x x x E (22)
8 68 Finit Volum Mthod owrful Mans of Enginring Dsign Whr 1 V VE d 2 a ae (23) Th closur for th convctd vlocity, and any othr transportd variabl, coms from an upwind schm. Th simplst option at th intgration point would b th FOU: u U if m 0 (24) Th intrpolation formula, Eq. (22), is not uniqu as xplaind in (Acharya, 2007) and on might succssfully implmnt th ida using modifid forms of this approach. Th ky point, howvr, is to us diffrnt closur quations for th convctd and convcting vlocity componnts. Th convcting vlocity should b proprly rlatd to th prssur fild. Th closur quation for p has not bn a mattr of concrn, at last for incomprssibl flows. Taking into account th lliptic bhavior of th prssur, th following closur quation is usd in th MIM: E (25) 2 This tchniqu guarants th rquird physical prssurvlocity coupling in th numrical modl and prvnts any wiggly solutions (Ashrafizadh t al., 2009). 3.4 Mthods basd on physical influncs Th plasur of working with all variabls on a singl nodal position, as opposd to th pain of working with staggrd grids limitd to simpl gomtris, has mad th MIM vry popular. Thrfor, attmpts hav bn mad to mploy and improv th MIM for th solution of flows at all spds. Schnidr and Raw (Schnidr & Raw, 1987a, 1987b) proposd th physical influnc schm to unify th intgration point vlocity intrpolation formulas in colocatd grids. Thy rtaind a unifid dfinition for th intgrationpoint vlocity componnts, but argud that th closur for u should not b obtaind solly by purly mathmatical upwind schms. Thy suggstd th nonconsrvativ form of th momntum quation to provid a closur for u : u u u u 0 x x x (26) Th mthod, calld th physical influnc schm, provids th following closur quation in th contxt of th 1D tst cas: x ˆ u u Up 2u x E (27)
9 Us of ropr Closur Equations in Finit Volum Discrtization Schms 69 Whr th ovr bar rfrs to th most rcnt availabl valu in th nonlinar itrations. Equation (25) is usd as th closur for p in this schm. Whil succssfully implmntd in a numbr of multidimnsional viscous problms, th physical influnc schm fails to prvnt prssur chckrboard problm in inviscid incomprssibl flows (Bakhtiari, 2008). This failur rinforcs th blif that two dfinitions for convctd and convcting vlocitis at intgration points ar ncssary whn th calculations ar carrid out on a colocatd grid. 3.5 Mthods basd on modifid physical influncs To fix th problm associatd with th physical influnc schm, Karimian and Schnidr (Karimian, 1994a) accpt th notion of dual vlocitis at th intgration points and mploy a combination of momntum and continuity quations to provid a closur for th convcting vlocity ( û ): Which rsults in: U ˆ UE x u Momntum Error ucontinuity Error 2 2u (28) uˆ U UE E (29) 2 2u Thy also us Eq. (27) for convctd vlocity ( u ) and Eq. (25) as th closur for p. Numrical xprimnts with th 1D tst cas show that this mthod works wll and succssfully supprsss any oscillatory numrical solution vn in timdpndd all spd flows (Karimian, 1994b). This mthod has also bn xtndd to 2D flows on structurd (Alisadghi t al., 2011a) and unstructurd (Alisadghi t al., 2011b) grids. Darbandi (Darbandi t al., 1997) proposd a prssurbasd all spd mthod similar to th Karimian's mthod, in which th momntum componnts ( f u ) ar usd as th flow dpndnt variabls instad of th vlocity componnts ( u ). Thy usd th following form of th momntum quation as th physical constraint on th intgration point convctd vlocity: f f u u f 0 x x x (30) Whr f u is th momntum variabl. Appropriat discrtization of Eq.(30) lads to: 1 f F E (31) 4u In this mthod, a combination of momntum and continuity quations is usd to constrain th intgration point convcting vlocity: f f u f u f u 0 x x x x (32)
10 70 Finit Volum Mthod owrful Mans of Enginring Dsign Whr is a usrdfind wight factor. This quation can b discrtizd to obtain th following formula: ˆ 1 f F E (33) 2 u In this mthod, lik all othr mthods discussd so far, Eq. (25) is usd as th closur for p. Numrical rsults hav shown that this st of closur quations rmovs th possibility of a chckrboard solution and provids strong coupling btwn prssur and vlocity filds. Anothr analysis has shown that th us of momntumvariabl formulation, instad of primitiv variabl formulation, improvs th stability and accuracy of th solution spcially in th solution of comprssibl flows with shock wavs (Darbandi, 2004). 3.6 Mthod of propr closur quations (MCE) Ondimnsional incomprssibl flow Ashrafizadh t al. (Ashrafizadh t al., 2009) hav shown that it is possibl to dvlop a colocatd numrical schm without rsorting to th convcting and convctd vlocity concpts, originally proposd by Rhi and Chow. Th proposd mthod, calld th mthod of propr closur quations (MCE), mploys a propr st of physically rlvant quations to constrain th vlocity and prssur at intgration points. It has alrady bn succssfully implmntd and usd in th solution of stady 1D inviscid incomprssibl and comprssibl flow problms (Ashrafizadh t al., 2008, 2009) and stady 2D viscous incomprssibl flow problms on both structurd and unstructurd grids (Alinia, 2011; Bakhtiari, 2008). In this mthod, following th physicalbasd approach proposd by Schnidr and Raw (Schnidr t al., 1987) and modifications proposd by Karimian and Darbandi (Darbandi, 1997; Karimian, 1994a), a combination of nonconsrvativ form of momntum and continuity quations is usd to obtain an intrpolation formula for th ipvlocity: which rsults in: momntum Eq. u continuity Eq. 0 u (34) u u 2 1 da u u 1 da m u u u S uu u S 0 x x Adx x x Adx UDS 1 CDS CDS 2 CDS (35) u m Hr S and S ar artificial momntum and mass sourc trms, and a and u ar th ipcross sctional ara and a schm control paramtr rspctivly. Thr ar options availabl to th usr at this stag. Using UDS or CDS Schms for trms 1 and 2 andu 0, 1, 2, as discussd in (Rzvani, 2008), diffrnt formulas for u can b obtaind. Th most common formula is as follows: E x u m u U ( S us ) (36) 2u 2u
11 Us of ropr Closur Equations in Finit Volum Discrtization Schms 71 As mntiond arlir, th classical closur choic for th intgration point prssur has bn th linar approximation, i.. Eq.(25). Howvr, th MCE mploys rlvant govrning quations to obtain th rquird closurs for all intgration point quantitis, including th prssur. Th momntum quation is a natural choic for obtaining a propr closur quation for p. Howvr, only th prssur gradint appars in th momntum quation and CDSbasd discrtization, physically propr for th prssur, rsults in wiggly solutions. As discussd in (Rzvani, 2010), by taking th divrgnc of th momntum quation, a prssur poisson quation for th prssur is obtaind. Th prssur oisson quation is an appropriat constraint quation for th prssur discrtization. Thrfor, as compard to Eq. (34), th following closur quation is proposd for th intgration point prssur: Divrgnc of 4u Continuity Eq. 0 Momntum Eq. x (37) Which rsults in: U U A A S 2 2 E u u x u u m E E a 4 x x (38) Whr 0 is a schm control paramtr. Numrical tst rsults hav shown that th mthod works wll with any nonzro positiv valu for (Rzvani2010) Ondimnsional comprssibl flow Extnsion of th MCE to comprssibl flow problms rquirs two additional tasks. First, th nrgy quation should also b solvd in ordr to find th tmpratur fild. Scond, dpnding on th chosn unknown variabls, propr linarization is crucially important. Following approachs prsntd in (Darbandi t al., 2007; Karimian, 1994a), th Nwton Raphson linarization stratgy is mployd to linariz th nonlinar trms in th balanc quations. For xampl, trm pu in th nrgy quation is linarizd as follows: pu Rt u u Rt R ft R ft tf ft (39) Closur quation for f can b asily obtaind using th unstadycomprssibl form of Eq. (34). First ordr upwind diffrnc schm is usd for th discrtization of momntum variabl gradint trms. Us of u 1 in th gnral form of th closur for th momntum variabl at th ast fac of th cll, rsults in th following formula for f : 2C C f C x C x f F S u S 12C u 12C 12C 12C u 12C E u m (40) Whr, C is th Courant numbr dfind as C u x and suprscript stands for th valus from prvious tim stp.
12 72 Finit Volum Mthod owrful Mans of Enginring Dsign Similarly, Closur quation for Eq. (37) as follow: p can b obtaind using th unstadycomprssibl form of p E p fu p u F F A A a E E 2 x u u m p 2 p S ux 4 x x 2 Th nonconsrvativ form of th continuity quation in Eq. (6), in th absnc of its sourc trm, is proposd as a propr closur quation for th dnsity (Darbandi t al., 1997; Karimian, 1994a). Hr, in ordr to control th stability and accuracy of th cod in both subsonic and suprsonic rgims, a smart blnding factor ( ) is usd which maks th quation hyprbolic in suprsonic flow rgions [Karimian, 1994a]: (41) x F f A x 1 2u x a x x 2 (42) m m M M 1 whr M is local Mach numbr at th ast intgration point and m=2 is suggstd for th suprscript m. Th final rquird closur quation is an intrpolation formula for th intgration point tmpratur. Hr, th nonconsrvativ form of th nrgy quation in th absnc of th sourc trm is a natural candidat (Darbandi t al., 1997) and on obtains th following xprssion: 2C Cp F FE C p A AE t t T 12C fc v1 2C E 12C cva 12C Th final discrt computational molculs for solving vlocity, prssur and tmpratur filds ar obtaind aftr closur quations ar substitutd in Eq. (9): pp pp pf pf pt pt p nb nb nb nb nb nb a a a F a F a T a T b (43) (44) (45) (46) ff ff fp fp ft ft f nb nb nb nb nb nb a F a F a a a T a T b (47) tt tt tf tf tp tp t nb nb nb nb nb nb at a T a F a F a a b ff fp In this linar algbraic st, a, a,... ar th influnc cofficints and b, b,... ar right hand sid constant trms. Subscript nb stands for th immdiat nighbors of nod (i.. W and E ). Th first suprscript rfrs to th rlvant quation and th scond suprscript points to th rlvant physical variabl. Th assmbld linar quations can b solvd using any dirct or itrativ solvr. Hr, th Gauss limination mthod is usd to solv th st of simultanous quations. Th coupld st of th algbraic quations might also b solvd using a smiimplicit (sgrgatd) solution stratgy. p f
13 Us of ropr Closur Equations in Finit Volum Discrtization Schms Twodimnsional incomprssibl flow Th MCE has also bn usd to solv 2D, incomprssibl laminar flows on both structurd and unstructurd grids (Bakhtiari, 2008; Alinia, 2011). Th implmntation follows th basic philosophy alrady dscribd in a 1D contxt. Howvr, mor complx intrpolation formulas ar rquird to tak into account th ffct of cll shap and orintation on balanc quations. Implmntation dtails ar not discussd hr for th sak of brvity. 4. Ondimnsional tst cass 4.1 Incomprssibl flow rformanc of th MCE for solving incomprssibl and comprssibl flows is xamind hr through two diffrnt tst cass. Stady incomprssibl flow in a constant ara duct is th first tst cas which is usd to xamin th prssurvlocity coupling in MCE. Quasi 1D flow in a convrgntdivrgnt nozzl is anothr tst cas usd to chck th accuracy of th proposd mthod in inviscid incomprssibl flow problms Stady incomprssibl flow in a constant ara duct Th 1D incomprssibl flow in a constant ara duct is usd to s if oscillatory solutions appar whn th MCE quations govrn th discrt flow fild. For this purpos, artificial lmnt and volum sourcs/sinks ar usd in continuity and momntum quations. In a domain with 30 nods, th artificial lmnt sourc trms ar locatd at th 10th and 20th cllfacs and control volum sourc trms ar activatd at th 10th and 20th nodal positions. Numrical rsults show that th ffcts of sharp variations in vlocity and prssur du to th artificial sourc trms ar quickly dampd out and MCE works wll vn with p 0. Figurs 5 and 6 show th ffcts of controlvolum artificial sourcs on th numrical solution. Ths rsults show that numrical oscillations ar dvlopd in both vlocity and prssur filds for p 0. Howvr, oscillations ar strongly dampd for any nonzro valu of p. Not that th notation ND in ths figurs stands for NonDimnsional and vlocity is normalizd with th inlt vlocity and prssur is normalizd using th xit prssur. Fig. 5. Vlocity and prssur filds in a constant ara duct with controlvolum mass sourc/sink.
14 74 Finit Volum Mthod owrful Mans of Enginring Dsign Fig. 6. Vlocity and prssur filds in a constant ara duct with controlvolum momntum sourc/sink Th quasi1d tst cas Idal flow in a convrgingdivrging nozzl is anothr tst cas to xamin th prformanc of MCE. Figur 7 compars th nondimnsional vlocity and prssur distributions along th nozzl with th analytical solutions. Excllnt agrmnt is obsrvd. Fig. 7. Vlocity and prssur fild in a 1D convrgingdivrging nozzl ( p 1 ). 4.2 Comprssibl flow For comprssibl flow, prformanc of MCE is xamind in thr diffrnt tst problms. Stady subsonic comprssibl flow in convrgntdivrgnt nozzl, stady comprssibl flow with normal shock in convrgntdivrgnt nozzl and unstady comprssibl flow with normal shock and xpansion wavs in a shock tub Stady subsonic comprssibl flow in a convrgntdivrgnt nozzl Rsults of solving subsonic comprssibl flow in a symmtric convrgntdivrgnt nozzl with aspct ratio 2 using 51 nods is shown in this sction. For th boundary condition implmntation, mass flow rat and th static tmpratur ar spcifid at th inlt and th
15 Us of ropr Closur Equations in Finit Volum Discrtization Schms 75 static prssur is providd at th outlt. In ordr to assss th accuracy of th proposd mthod, numrical solutions corrsponding to inlt Mach numbrs of 0.05, 0.1, 0.15, 0.2, and 0.25 ar prsntd in Fig. 8. Flow with th inlt Mach numbr of , for which th nozzl is chokd, is also considrd. In this tst cas, p 1 is usd and prssur, tmpratur, and dnsity ar nondimnsionalizd with th outlt prssur, inlt tmpratur, and inlt dnsity, rspctivly Stady flow with shock wavs In th divrgnt sction of th prvious tst cas, normal shock wavs appar for crtain inlt boundary conditions whil th position and th strngth of th shock can b controlld by rgulating th back prssur. Figur 9 prsnts th Mach numbr, nondimnsional prssur, tmpratur and dnsity distributions along th nozzl for th prssur ratios ( outlt inlt ) of 0.65, 0.75 and 0.85 using 81 nods. It is sn that th normal shock is wll capturd using a fw nods in all cass and no numrical oscillation is dvlopd along th nozzl. Not that this rsults ar computd using m 2 in Eq.(43) and p taks valus lowr than 1 in th suprsonic parts of th solution domain and valus highr than 1 in th subsonic rgions. Fig. 8. Stady subsonic comprssibl flow in a convrgntdivrgnt nozzl.
16 76 Finit Volum Mthod owrful Mans of Enginring Dsign Unstady comprssibl flow in a shock tub As th final unstady 1D tst cas with both comprssion and xpansion wavs, shock tub problm is discussd. Th gomtry shown in Fig.10 is a constant ara duct dividd into two rgions sparatd by a diaphragm at th middl of th duct. At 0 th gas tmpratur is 298 K, th prssur at th right sid of th mmbran is 100 Ka and th prssur at th lft sid of th mmbran is 1000 Ka. Th tim stp for th marching numrical solution is 6 qual to 4 10 sconds. Figur 11 shows th numrically calculatd distributions of Mach numbr, nondimnsional prssur, tmpratur, and dnsity aftr sconds. Hr, th Courant numbr is which is computd using th maximum magnitud of vlocity in th domain. Th prssur, tmpratur and dnsity ar nondimnsionalizd using th low prssur sid valus. Good agrmnt with th xact solution is obsrvd. Fig. 9. Stady flow with shock wavs in a convrgntdivrgnt nozzl. 5. Twodimnsional incomprssibl flow tst cass Th MCE has also bn usd to solv many 2D bnchmark tst problms. Hr, only th rsults for two standard tst cass ar prsntd.
17 Us of ropr Closur Equations in Finit Volum Discrtization Schms 77 Fig. 10. Gomtry and th flow pattrn for th shock tub problm. Fig. 11. Rsults of unstady comprssibl flow in shock tub
18 78 Finit Volum Mthod owrful Mans of Enginring Dsign 5.1 Flow in a 2D liddrivn cavity Flow in a liddrivn cavity, shown in Fig. 12, is a widly usd tst cas for 2D stady incomprssibl flows govrnd by NavirStoks quations. Numrical solution via MCE on an unstructurd grid with triangular lmnts, shown in Fig. 13, is compard to th numrical rsults of Erturk (Erturk, 2005). Flow pattrns at diffrnt Rynolds numbrs ar shown in Fig. 14. Excllnt agrmnt btwn th calculatd vlocitis at cavity cntrlins and th bnchmark rsults in Fig. 14 is also obsrvd. Fig. 12. Liddrivn cavity boundary conditions with an schmatic of flow faturs. Fig. 13. Gnratd unstructurd msh for liddrivn cavity tst cas. 5.2 Flow ovr a backward facingstp Th backward stp flow is anothr commonly usd tst cas for numrical solution algorithms. Th gomtry of this tst cas is shown in Fig. 15. Forty nods ar usd across th channl, 300 nods ar usd from th stp to 30h and 50 nods ar usd from 30h to 50h along th channl. Nods ar distributd uniformly. Figur 16 shows th calculatd stramlins at R=800, in which a rcirculation zon appars at th uppr wall. Normalizd locations of th rattachmnt point at th lowr wall (X1), th sparation point at th uppr
19 Us of ropr Closur Equations in Finit Volum Discrtization Schms 79 wall (X2), and th rattachmnt point at th uppr wall (X3) ar compard with th rsults obtaind by Gartling (Gartling, 1990) for flow at R = 800. Comparison rsults ar shown in Tabl 1 and xcllnt agrmnt is obvious (Erturk, 2008). A comparison has also bn mad btwn numrical rsults via MCE and th bnchmark rsults in (Gartling, 1990) and shown in Fig. 18. Again, xcllnt agrmnt is obsrvd. Similar calculations and comparisons hav also bn mad for a backwardstp flow modl with an inlt channl. Excllnt rsults, not shown hr, ar also obtaind in this tst cas. (a) Stramlins at R = (b) Vlocity rofils at R = (c) Stramlins at R = (d) Vlocity rofils at R = Fig. 14. Schmatic viw of th backward facingstp flow problm.
20 80 Finit Volum Mthod owrful Mans of Enginring Dsign Fig. 15. Calculatd stramlins in th backwardstp flow problm at R = 800. MCE Erturk X X X Tabl 1. Comparison of normalizd X1, X2 and X3 locations and lngth of th first rcirculating rgion on th uppr wall (X3 X2) for R = 800 X3 X Fig. 16. Horizontal vlocity profils across th channl at X/h = 14 and 30 for R = Conclusion Diffrnt prssurbasd finit volum discrtization schms for th solution of fluid flow problms ar xplaind and compard in th contxt of colocatd grid arrangmnts and a nw schm, calld th Mthod of ropr Closur Equations (MCE), is proposd. Th proposd mthod is th only discrtization schm on colocatd grids which mploys only on dfinition for th intgration point vlocity and succssfully solvs incomprssibl as wll as comprssibl flow problms. Simpl ondimnsional tst cass ar prsntd to discuss diffrnt availabl schms and th building blocks of MCE. Implmntation of th MCE on 2D structurd and unstructurd grids has also bn carrid out but th dtails ar not rportd hr for th sak of brvity. Numrical solutions via MCE, howvr, ar prsntd for both 1D and 2D problms which clarly show th satisfactory prformanc of th MCE in th numrical solution of fluid flow problms.
21 Us of ropr Closur Equations in Finit Volum Discrtization Schms Rfrncs Acharya, S; Baliga, B. R; Karki, K.; Murthy, J. Y.; rakash, C. & Vanka, S.. (2007). rssur Basd FinitVolum Mthods in Computational Fluid Dynamics. Journal of Hat Transfr, Vol. 129, pp Alinia, B. (2011). Implmntation of th Mthod of ropr Closur Equations on Two Dimnsional Unstructurd Grids for Stady, Incomprssibl Viscous roblms, M.Sc. Thsis, Dpartmnt of Mchanical Enginring, K.N.Toosi Univrsity of Tchnology, Thran, Iran, (In rsian) Alisadghi, H. & Karimian, S. M. H. (2011a). Diffrnt Modlings of CllFac Vlocitis and Thir Effcts on th rssurvlocity Coupling, Accuracy and Convrgnc of Solution, Intrnational Journal for Numrical Mthods in Fluids, Vol. 65, pp Alisadghi, H. & Karimian, S. M. H. (2011b). Comparison of Diffrnt Solution Algorithms for Collocatd Mthod of MCIM to Calculat Stady and Unstady Incomprssibl Flows on Unstructurd Grids, Computrs and Fluids, Vol. 46, Ashrafizadh, A.; Okhovat, S.; ourbagian, M. & Raithby, G. D. (2011). A SmiCoupld Solution Algorithm in Arodynamic Invrs Shap Dsign, Invrs roblms in Scinc and Enginring, Vol. 19, No. 4, pp , ISSN Ashrafizadh, A.; Rzvani M.; & Bakhtiari, B. (2008). A Nw Colocatd rssurbasd Finit Volum Mthod, 16th Annual Confrnc of th CFD Socity of Canada, Saskatoon, Saskatchwan, Canada Ashrafizadh, A.; Rzvani M.; & Bakhtiari, B. (2009). rssurvlocity Coupling On Co Locatd Grids Using th Mthod of ropr Closur Equations, Numrical Hat Transfr; art B, Vol. 56, pp Bakhtiari, B. (2008). Dcoupling of rssur and Vlocity in th Numrical Simulation of Incomprssibl Flows: A rviw. M.Sc. Thsis, K.N.Toosi Univrsity of Tchnology, Thran, Iran, (In rsian) Chtrntal, I. (1999). TwoDimnsional Eulr Equations Solvr, M.Sc. Thsis, Graduat Dpartmnt of Arospac Scinc and Enginring, Univrsity of Toronto, Toronto, Canada Darbandi, M. & Schnidr, G. E. (1997). Momntum Variabl rocdur for Solving Comprssibl and Incomprssibl Flows. AIAA Journal, Vol. 35, pp Darbandi, M. & Mokarizadh, V. (2004). A Modifid rssurbasd Algorithm to Solv Flow Filds With Shock and Expansion Wavs, Numrical Hat Transfr, art B, Vol. 46, pp Darbandi, M.; Roohi, E. & Mokarizadh, V. (2007). Concptual Linarization of Eulr Govrning Equations to Solv High Spd Comprssibl Flow Using a rssur Basd Mthod, Numrical Mthods for artial Diffrntial Equations, DOI , pp Erturk, E.; Cork, T. C. & Gökçöl, C. (2005). Numrical Solutions of 2D Stady Incomprssibl Drivn Cavity Flow at High Rynolds Numbrs, Int. J. Numr. Mth. Fluids, Vol. 48, pp Erturk, E. (2008). Numrical Solutions of 2D Stady Incomprssibl Flow Ovr a Backward Facing Stp, art I: High Rynolds Numbr Solutions, Int. J. Computrs & Fluids, Vol. 37, pp Gartling, D. K. (1990). A Tst roblm for Outflow Boundary Conditions Flow Ovr a Backward Facing Stp, Int. J. Numr. Mth. Fluids, Vol. 11, pp
22 82 Finit Volum Mthod owrful Mans of Enginring Dsign Karimian, S. M. H. & Schnidr, G. E. (1994a). Application of a ControlVolumBasd FinitElmnt Formulation to th Shock Tub roblm, AIAA Journal, Vol.33 No.1, pp Karimian, S. M. H. & Schnidr, G. E. (1994b). rssurbasd Computational Mthod for Comprssibl and Incomprssibl Flows, Journal of Thrmophysics and Hat Transfr, Vol.8, No.2, pp atankar, S. V. (1980). Numrical Hat Transfr and Fluid Flow, Hmisphr ublishing Corporation, Washington, D. C. Rzvani, M. (2008). A Nw rssurbasd Computational Mthod for Solving Fluid flow at all spds, M.Sc. Thsis, Dpartmnt of Mchanical Enginring, K.N.Toosi Univrsity of Tchnology, Thran, Iran, (In rsian) Rzvani, M. & Ashrafizadh, A. (2010). Numrical Simulation of th IntrEquation Couplings in AllSpd Flows via th Mthod of ropr Closur Equations, Numrical Hat Transfr; art A, Vol. 58, pp Rhi, C.M. & Chow, W. L. (1983). A Numrical Study of th Turbulnt Flow ast an Isolatd airfoil With Trailing Edg Sparation. AIAA Journal, Vol. 21, pp Schnidr, G. E. & Raw, M. J. (1987a). Control Volum FinitElmnt Mthod For hat Transfr and Fluid Flow Using Colocatd Variabls1. Computational rocdur. Numrical Hat Transfr, part B, Vol. 11, pp Schnidr, G. E. & Raw, M. J. (1987b). Control Volum FinitElmnt Mthod For hat Transfr and Fluid Flow Using Colocatd Variabls2. Application and Validation. Numrical Hat Transfr, part B, Vol. 11, pp
23 Finit Volum Mthod  owrful Mans of Enginring Dsign Editd by hd. Radostina trova ISBN Hard covr, 370 pags ublishr InTch ublishd onlin 28, March, 2012 ublishd in print dition March, 2012 W hop that among ths chaptrs you will find a topic which will rais your intrst and ngag you to furthr invstigat a problm and build on th prsntd work. This book could srv ithr as a txtbook or as a practical guid. It includs a wid varity of concpts in FVM, rsult of th fforts of scintists from all ovr th world. Howvr, just to hlp you, all book chaptrs ar systmizd in thr gnral groups: Nw tchniqus and algorithms in FVM; Solution of particular problms through FVM and Application of FVM in mdicin and nginring. This book is for vryon who wants to grow, to improv and to invstigat. How to rfrnc In ordr to corrctly rfrnc this scholarly work, fl fr to copy and past th following: A. Ashrafizadh, M. Rzvani and B. Alinia (2012). Us of ropr Closur Equations in Finit Volum Discrtization Schms, Finit Volum Mthod  owrful Mans of Enginring Dsign, hd. Radostina trova (Ed.), ISBN: , InTch, Availabl from: InTch Europ Univrsity Campus ST Ri Slavka Krautzka 83/A Rijka, Croatia hon: +385 (51) Fax: +385 (51) InTch China Unit 405, Offic Block, Hotl Equatorial Shanghai No.65, Yan An Road (Wst), Shanghai, , China hon: Fax:
ME469A Numerical Methods for Fluid Mechanics
ME469A Numrical Mthods for Fluid Mchanics Handout #5 Gianluca Iaccarino Finit Volum Mthods Last tim w introducd th FV mthod as a discrtization tchniqu applid to th intgral form of th govrning quations
More informationFinite element discretization of Laplace and Poisson equations
Finit lmnt discrtization of Laplac and Poisson quations Yashwanth Tummala Tutor: Prof S.Mittal 1 Outlin Finit Elmnt Mthod for 1D Introduction to Poisson s and Laplac s Equations Finit Elmnt Mthod for 2DDiscrtization
More information3 Finite Element Parametric Geometry
3 Finit Elmnt Paramtric Gomtry 3. Introduction Th intgral of a matrix is th matrix containing th intgral of ach and vry on of its original componnts. Practical finit lmnt analysis rquirs intgrating matrics,
More informationElements of Statistical Thermodynamics
24 Elmnts of Statistical Thrmodynamics Statistical thrmodynamics is a branch of knowldg that has its own postulats and tchniqus. W do not attmpt to giv hr vn an introduction to th fild. In this chaptr,
More informationAnalysis of ConvectionDiffusion Problems at Various Peclet Numbers Using Finite Volume and Finite Difference Schemes Anand Shukla
Mathmatical Thory and Modling.iist.org ISSN 45804 (apr) ISSN 505 (Onlin) Vol., No.6, 01Slctd from Intrnational Confrnc on Rcnt Trnds in Applid Scincs ith nginring Applications Analysis of Convctioniffusion
More informationNusselt number correlations for simultaneously developing laminar duct flows of liquids with temperature dependent properties
Journal of Physics: Confrnc Sris OPEN ACCESS Nusslt numbr corrlations for simultanously dvloping laminar duct flows of liquids with tmpratur dpndnt proprtis To cit this articl: Stfano Dl Giudic t al 2014
More informationu 3 = u 3 (x 1, x 2, x 3 )
Lctur 23: Curvilinar Coordinats (RHB 8.0 It is oftn convnint to work with variabls othr than th Cartsian coordinats x i ( = x, y, z. For xampl in Lctur 5 w mt sphrical polar and cylindrical polar coordinats.
More informationSliding Mode Flow Rate Observer Design
Sliding Mod Flow Rat Obsrvr Dsign Song Liu and Bin Yao School of Mchanical Enginring, Purdu Univrsity, Wst Lafaytt, IN797, USA liu(byao)@purdudu Abstract Dynamic flow rat information is ndd in a lot of
More informationThe Matrix Exponential
Th Matrix Exponntial (with xrciss) by D. Klain Vrsion 207.0.05 Corrctions and commnts ar wlcom. Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix A A k I + A + k!
More information4.2 Design of Sections for Flexure
4. Dsign of Sctions for Flxur This sction covrs th following topics Prliminary Dsign Final Dsign for Typ 1 Mmbrs Spcial Cas Calculation of Momnt Dmand For simply supportd prstrssd bams, th maximum momnt
More information16. Electromagnetics and vector elements (draft, under construction)
16. Elctromagntics (draft)... 1 16.1 Introduction... 1 16.2 Paramtric coordinats... 2 16.3 Edg Basd (Vctor) Finit Elmnts... 4 16.4 Whitny vctor lmnts... 5 16.5 Wak Form... 8 16.6 Vctor lmnt matrics...
More informationElement connectivity parameterization method for the stress based topology optimization for geometrically nonlinear structure
0 th World Congrss on Structural and Multidisciplinary Optimization May 94, 03, Orlando, Florida, USA Elmnt connctivity paramtrization mthod for th strss basd topology optimization for gomtrically nonlinar
More informationThe van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012
Th van dr Waals intraction D. E. Sopr 2 Univrsity of Orgon 20 pril 202 Th van dr Waals intraction is discussd in Chaptr 5 of J. J. Sakurai, Modrn Quantum Mchanics. Hr I tak a look at it in a littl mor
More informationCE 530 Molecular Simulation
CE 53 Molcular Simulation Lctur 8 Frnrgy calculations David A. Kofk Dpartmnt of Chmical Enginring SUNY Buffalo kofk@ng.buffalo.du 2 FrEnrgy Calculations Uss of fr nrgy Phas quilibria Raction quilibria
More informationCHAPTER 1. Introductory Concepts Elements of Vector Analysis Newton s Laws Units The basis of Newtonian Mechanics D Alembert s Principle
CHPTER 1 Introductory Concpts Elmnts of Vctor nalysis Nwton s Laws Units Th basis of Nwtonian Mchanics D lmbrt s Principl 1 Scinc of Mchanics: It is concrnd with th motion of matrial bodis. odis hav diffrnt
More information5.80 SmallMolecule Spectroscopy and Dynamics
MIT OpnCoursWar http://ocw.mit.du 5.80 SmallMolcul Spctroscopy and Dynamics Fall 008 For information about citing ths matrials or our Trms of Us, visit: http://ocw.mit.du/trms. Lctur # 3 Supplmnt Contnts
More informationMCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17)
MCB37: Physical Biology of th Cll Spring 207 Homwork 6: Ligand binding and th MWC modl of allostry (Du 3/23/7) Hrnan G. Garcia March 2, 207 Simpl rprssion In class, w drivd a mathmatical modl of how simpl
More informationOn the Hamiltonian of a MultiElectron Atom
On th Hamiltonian of a MultiElctron Atom Austn Gronr Drxl Univrsity Philadlphia, PA Octobr 29, 2010 1 Introduction In this papr, w will xhibit th procss of achiving th Hamiltonian for an lctron gas. Making
More information1 General boundary conditions in diffusion
Gnral boundary conditions in diffusion Πρόβλημα 4.8 : Δίνεται μονοδιάτατη πλάκα πάχους, που το ένα άκρο της κρατιέται ε θερμοκραία T t και το άλλο ε θερμοκραία T 2 t. Αν η αρχική θερμοκραία της πλάκας
More informationObserver Bias and Reliability By Xunchi Pu
Obsrvr Bias and Rliability By Xunchi Pu Introduction Clarly all masurmnts or obsrvations nd to b mad as accuratly as possibl and invstigators nd to pay carful attntion to chcking th rliability of thir
More informationExam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator.
Exam N a m : _ S O L U T I O N P U I D : I n s t r u c t i o n s : It is important that you clarly show your work and mark th final answr clarly, closd book, closd nots, no calculator. T i m : h o u r
More information10. The DiscreteTime Fourier Transform (DTFT)
Th DiscrtTim Fourir Transform (DTFT Dfinition of th discrttim Fourir transform Th Fourir rprsntation of signals plays an important rol in both continuous and discrt signal procssing In this sction w
More informationMEASURING HEAT FLUX FROM A COMPONENT ON A PCB
MEASURING HEAT FLUX FROM A COMPONENT ON A PCB INTRODUCTION Elctronic circuit boards consist of componnts which gnrats substantial amounts of hat during thir opration. A clar knowldg of th lvl of hat dissipation
More informationA dynamic boundary model for implementation of boundary conditions in latticeboltzmann method
Journal of Mchanical Scinc and Tchnology Journal of Mchanical Scinc and Tchnology 22 (28) 92~2.springrlink.com/contnt/738494x A dynamic boundary modl for implmntation of boundary conditions in latticboltzmann
More informationAvailable online at ScienceDirect. Procedia Engineering 126 (2015 )
Availabl onlin at www.scincdirct.com ScincDirct Procdia Enginring 26 (25 ) 628 632 7t Intrnational Confrnc on Fluid Mcanics, ICFM7 Applications of ig ordr ybrid DG/FV scms for twodimnsional RAS simulations
More informationMAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design
MAE4700/5700 Finit Elmnt Analysis for Mchanical and Arospac Dsign Cornll Univrsity, Fall 2009 Nicholas Zabaras Matrials Procss Dsign and Control Laboratory Sibly School of Mchanical and Arospac Enginring
More informationTypes of Transfer Functions. Types of Transfer Functions. Types of Transfer Functions. Ideal Filters. Ideal Filters
Typs of Transfr Typs of Transfr x[n] X( LTI h[n] H( y[n] Y( y [ n] h[ k] x[ n k] k Y ( H ( X ( Th timdomain classification of an LTI digital transfr function is basd on th lngth of its impuls rspons h[n]:
More informationIntroduction to Condensed Matter Physics
Introduction to Condnsd Mattr Physics pcific hat M.P. Vaughan Ovrviw Ovrviw of spcific hat Hat capacity DulongPtit Law Einstin modl Dby modl h Hat Capacity Hat capacity h hat capacity of a systm hld at
More informationis an appropriate single phase forced convection heat transfer coefficient (e.g. Weisman), and h
For t BWR oprating paramtrs givn blow, comput and plot: a) T clad surfac tmpratur assuming t JnsLotts Corrlation b) T clad surfac tmpratur assuming t Tom Corrlation c) T clad surfac tmpratur assuming
More informationDifferentiation of Exponential Functions
Calculus Modul C Diffrntiation of Eponntial Functions Copyright This publication Th Northrn Albrta Institut of Tchnology 007. All Rights Rsrvd. LAST REVISED March, 009 Introduction to Diffrntiation of
More informationA PreyPredator Model with an Alternative Food for the Predator, Harvesting of Both the Species and with A Gestation Period for Interaction
Int. J. Opn Problms Compt. Math., Vol., o., Jun 008 A PryPrdator Modl with an Altrnativ Food for th Prdator, Harvsting of Both th Spcis and with A Gstation Priod for Intraction K. L. arayan and. CH. P.
More informationSearch sequence databases 3 10/25/2016
Sarch squnc databass 3 10/25/2016 Etrm valu distribution Ø Suppos X is a random variabl with probability dnsity function p(, w sampl a larg numbr S of indpndnt valus of X from this distribution for an
More information22/ Breakdown of the BornOppenheimer approximation. Selection rules for rotationalvibrational transitions. P, R branches.
Subjct Chmistry Papr No and Titl Modul No and Titl Modul Tag 8/ Physical Spctroscopy / Brakdown of th BornOppnhimr approximation. Slction ruls for rotationalvibrational transitions. P, R branchs. CHE_P8_M
More informationA nonequilibrium molecular dynamics simulation of evaporation
Intrnational Confrnc Passiv and Low Enrgy Cooling 543 A nonquilibrium molcular dynamics simulation of vaporation Z.J. Wang, M. Chn and Z.Y. Guo Dpartmnt of Enginring Mchanics, Tsinghua Univrsity, Bijing
More informationModeling of the Plasma Flow in HighPower TAL
39th AIAA/ASME/SAE/ASEE Joint Propulsion Confrnc and Exhibit 2023 July 2003, Huntsvill, Alabama AIAA 20034701 AIAA20034701 Modling of th Plasma Flow in HighPowr TAL Michal Kidar, Iain D. Boyd and
More information2.3 Matrix Formulation
23 Matrix Formulation 43 A mor complicatd xampl ariss for a nonlinar systm of diffrntial quations Considr th following xampl Exampl 23 x y + x( x 2 y 2 y x + y( x 2 y 2 (233 Transforming to polar coordinats,
More informationCOORDINATION OF FAST NUMERICAL RELAYS AND CURRENT TRANSFORMERS OVERDIMENSIONING FACTORS AND INFLUENCING PARAMETERS
COORDINATION OF FAST NUMERICAL RELAYS AND CURRENT TRANSFORMERS OVERDIMENSIONING FACTORS AND INFLUENCING PARAMETERS Stig Holst ABB Automation Products Swdn Bapuji S Palki ABB Utilitis India This papr rports
More informationWhat are those βs anyway? Understanding Design Matrix & Odds ratios
Ral paramtr stimat WILD 750  Wildlif Population Analysis of 6 What ar thos βs anyway? Undrsting Dsign Matrix & Odds ratios Rfrncs Hosmr D.W.. Lmshow. 000. Applid logistic rgrssion. John Wily & ons Inc.
More informationLecture # 12: Shock Waves and De Laval Nozzle
ArE 311L & ArE343L Lctur Nots Lctur # 1: Shock Wavs and D Laval Nozzl Dr. Hui H Hu Dpartmnt of Arospac Enginring Iowa Stat Univrsity Ams, Iowa 50011, U.S.A ArE311L Lab#3: rssur Masurmnts in a d Laval Nozzl
More informationBINOMIAL COEFFICIENTS INVOLVING INFINITE POWERS OF PRIMES
BINOMIAL COEFFICIENTS INVOLVING INFINITE POWERS OF PRIMES DONALD M. DAVIS Abstract. If p is a prim (implicit in notation and n a positiv intgr, lt ν(n dnot th xponnt of p in n, and U(n n/p ν(n, th unit
More informationSECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.
SETION 6. 57 6. Evaluation of Dfinit Intgrals Exampl 6.6 W hav usd dfinit intgrals to valuat contour intgrals. It may com as a surpris to larn that contour intgrals and rsidus can b usd to valuat crtain
More informationME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002
3.4 Forc Analysis of Linkas An undrstandin of forc analysis of linkas is rquird to: Dtrmin th raction forcs on pins, tc. as a consqunc of a spcifid motion (don t undrstimat th sinificanc of dynamic or
More informationCalculus concepts derivatives
All rasonabl fforts hav bn mad to mak sur th nots ar accurat. Th author cannot b hld rsponsibl for any damags arising from th us of ths nots in any fashion. Calculus concpts drivativs Concpts involving
More informationA POSTERIORI ERROR ESTIMATION FOR AN INTERIOR PENALTY TYPE METHOD EMPLOYING H(DIV) ELEMENTS FOR THE STOKES EQUATIONS
A POSTERIORI ERROR ESTIMATION FOR AN INTERIOR PENALTY TYPE METHOD EMPLOYING H(DIV) ELEMENTS FOR THE STOKES EQUATIONS JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This papr stablishs a postriori rror
More informationSec 2.3 Modeling with First Order Equations
Sc.3 Modling with First Ordr Equations Mathmatical modls charactriz physical systms, oftn using diffrntial quations. Modl Construction: Translating physical situation into mathmatical trms. Clarly stat
More information2. Laser physics  basics
. Lasr physics  basics Spontanous and stimulatd procsss Einstin A and B cofficints Rat quation analysis Gain saturation What is a lasr? LASER: Light Amplification by Stimulatd Emission of Radiation "light"
More informationCollisions between electrons and ions
DRAFT 1 Collisions btwn lctrons and ions Flix I. Parra Rudolf Pirls Cntr for Thortical Physics, Unirsity of Oxford, Oxford OX1 NP, UK This rsion is of 8 May 217 1. Introduction Th FokkrPlanck collision
More informationThe pn junction: 2 Current vs Voltage (IV) characteristics
Th pn junction: Currnt vs Voltag (V) charactristics Considr a pn junction in quilibrium with no applid xtrnal voltag: o th V E F E F V ptyp Dpltion rgion ntyp Elctron movmnt across th junction: 1. n
More informationDifferential Equations
UNIT I Diffrntial Equations.0 INTRODUCTION W li in a world of intrrlatd changing ntitis. Th locit of a falling bod changs with distanc, th position of th arth changs with tim, th ara of a circl changs
More informationUltimate strength analysis & design of residential slabs on reactive soil
Ultimat strngth analysis & dsign of rsidntial slabs on ractiv soil This documnt prsnts an ovrviw of thory undrlying ultimat strngth analysis and dsign of stiffnd raft and waffl raft slabs, as commonly
More informationA SPLITTING METHOD USING DISCONTINUOUS GALERKIN FOR THE TRANSIENT INCOMPRESSIBLE NAVIERSTOKES EQUATIONS
ESAIM: MAN Vol. 39, N o 6, 005, pp. 5 47 DOI: 0.05/man:005048 ESAIM: Mathmatical Modlling and Numrical Analysis A SPLITTING METHOD USING DISCONTINUOUS GALERKIN FOR THE TRANSIENT INCOMPRESSIBLE NAVIERSTOKES
More informationPROOF OF FIRST STANDARD FORM OF NONELEMENTARY FUNCTIONS
Intrnational Journal Of Advanc Rsarch In Scinc And Enginring http://www.ijars.com IJARSE, Vol. No., Issu No., Fbruary, 013 ISSN3198354(E) PROOF OF FIRST STANDARD FORM OF NONELEMENTARY FUNCTIONS 1 Dharmndra
More informationEFFECT OF BALL PROPERTIES ON THE BALLBAT COEFFICIENT OF RESTITUTION
EFFECT OF BALL PROPERTIES ON THE BALLBAT COEFFICIENT OF RESTITUTION A. M. NATHAN 1 AND L. V. SMITH 2 1 Univrsity of Illinois, 1110 W. Grn Strt, Urbana, IL 61801, USA, Email: anathan@illinois.du 2 Washington
More informationProximity Templates for Modeling of Skin and Proximity Effects on Packages and High Frequency Interconnect
Proximity Tmplats for Modling of Skin and Proximity Effcts on Packags and High Frquncy Intrconnct Luca Danil Massachustts Instit. of Tch. dluca@mit.du Albrto Sangiovanniincntlli Univ. of California, Brkly
More informationVALIDATION OF FINITE ELEMENT PROGRAM FOR JOURNAL BEARINGS  STATIC AND DYNAMIC PROPERTIES
Univrsity of Kntucky UKnowldg Univrsity of Kntucky Mastr's Thss Graduat School 2004 VALIDATION OF FINITE ELEMENT PROGRAM FOR JOURNAL BEARINGS  STATIC AND DYNAMIC PROPERTIES Raja Shkar Balupari Univrsity
More informationINC 693, 481 Dynamics System and Modelling: The Language of Bound Graphs Dr.Ing. Sudchai Boonto Assistant Professor
INC 693, 48 Dynamics Systm and Modlling: Th Languag o Bound Graphs Dr.Ing. Sudchai Boonto Assistant Prossor Dpartmnt o Control Systm and Instrumntation Enginring King Mongkut s Unnivrsity o Tchnology
More informationDefinition1: The ratio of the radiation intensity in a given direction from the antenna to the radiation intensity averaged over all directions.
Dirctivity or Dirctiv Gain. 1 Dfinition1: Dirctivity Th ratio of th radiation intnsity in a givn dirction from th antnna to th radiation intnsity avragd ovr all dirctions. Dfinition2: Th avg U is obtaind
More informationPHA 5127 Answers Homework 2 Fall 2001
PH 5127 nswrs Homwork 2 Fall 2001 OK, bfor you rad th answrs, many of you spnt a lot of tim on this homwork. Plas, nxt tim if you hav qustions plas com talk/ask us. Thr is no nd to suffr (wll a littl suffring
More informationPipe flow friction, small vs. big pipes
Friction actor (t/0 t o pip) Friction small vs larg pips J. Chaurtt May 016 It is an intrsting act that riction is highr in small pips than largr pips or th sam vlocity o low and th sam lngth. Friction
More informationWaves in cavities such as vehicle compartments, rooms or ducts
7.1 Wavs in cavitis such as vhicl compartmnts, rooms or ducts Sound propagation from sourcs into th fr fild is rlativly simpl. At a rciving position in a distanc to th sourc th sound will arriv dlayd by
More informationEngineering 323 Beautiful HW #13 Page 1 of 6 Brown Problem 512
Enginring Bautiful HW #1 Pag 1 of 6 5.1 Two componnts of a minicomputr hav th following joint pdf for thir usful liftims X and Y: = x(1+ x and y othrwis a. What is th probability that th liftim X of th
More informationDetermination of Vibrational and Electronic Parameters From an Electronic Spectrum of I 2 and a BirgeSponer Plot
5 J. Phys. Chm G Dtrmination of Vibrational and Elctronic Paramtrs From an Elctronic Spctrum of I 2 and a BirgSponr Plot 1 15 2 25 3 35 4 45 Dpartmnt of Chmistry, Gustavus Adolphus Collg. 8 Wst Collg
More informationVII. Quantum Entanglement
VII. Quantum Entanglmnt Quantum ntanglmnt is a uniqu stat of quantum suprposition. It has bn studid mainly from a scintific intrst as an vidnc of quantum mchanics. Rcntly, it is also bing studid as a basic
More informationBifurcation Theory. , a stationary point, depends on the value of α. At certain values
Dnamic Macroconomic Thor Prof. Thomas Lux Bifurcation Thor Bifurcation: qualitativ chang in th natur of th solution occurs if a paramtr passs through a critical point bifurcation or branch valu. Local
More informationGradebook & Midterm & Office Hours
Your commnts So what do w do whn on of th r's is 0 in th quation GmM(1/r1/r)? Do w nd to driv all of ths potntial nrgy formulas? I don't undrstand springs This was th first lctur I actually larnd somthing
More informationREGISTER!!! The Farmer and the Seeds (a parable of scientific reasoning) Class Updates. The Farmer and the Seeds. The Farmer and the Seeds
How dos light intract with mattr? And what dos (this say about) mattr? REGISTER!!! If Schrödingr s Cat walks into a forst, and no on is around to obsrv it, is h rally in th forst? sourc unknown Phys 1010
More informationu r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C
Tchniqus of Intgration c Donald Kridr and Dwight Lahr In this sction w ar going to introduc th first approachs to valuating an indfinit intgral whos intgrand dos not hav an immdiat antidrivativ. W bgin
More informationPhysics 312 First Pledged Problem Set
Physics 31 First Pldgd Problm St 1. Th ground stat of hydrogn is dscribd by th wavfunction whr a is th Bohr radius. (a) Comput th charg dnsity Ã (r) = 1 p ¼ µ 1 a 3 r=a ; ½ (r) = jã (r)j : and plot 4¼r
More informationNotes on FermiDirac Integrals 2 nd Edition. Raseong Kim and Mark Lundstrom Network for Computational Nanotechnology Purdue University
Nots on rmidirac Intgrals nd Edition Rasong Kim and Mark Lundstrom Ntwork for Computational Nanotchnology Purdu Univrsity Sptmbr 3, 8 1. Introduction rmidirac intgrals appar frquntly in smiconductor
More informationDesign Guidelines for Quartz Crystal Oscillators. R 1 Motional Resistance L 1 Motional Inductance C 1 Motional Capacitance C 0 Shunt Capacitance
TECHNICAL NTE 30 Dsign Guidlins for Quartz Crystal scillators Introduction A CMS Pirc oscillator circuit is wll known and is widly usd for its xcllnt frquncy stability and th wid rang of frquncis ovr which
More informationSara Godoy del Olmo Calculation of contaminated soil volumes : Geostatistics applied to a hydrocarbons spill Lac Megantic Case
wwwnvisolcanadaca Sara Godoy dl Olmo Calculation of contaminatd soil volums : Gostatistics applid to a hydrocarbons spill Lac Mgantic Cas Gostatistics: study of a PH contamination CONTEXT OF THE STUDY
More informationSUMMER 17 EXAMINATION
(ISO/IEC  75 Crtifid) SUMMER 7 EXAMINATION Modl wr jct Cod: Important Instructions to aminrs: ) Th answrs should b amind by ky words and not as wordtoword as givn in th modl answr schm. ) Th modl answr
More informationMaxwellian Collisions
Maxwllian Collisions Maxwll ralizd arly on that th particular typ of collision in which th crosssction varis at Q rs 1/g offrs drastic siplifications. Intrstingly, this bhavior is physically corrct for
More informationECE602 Exam 1 April 5, You must show ALL of your work for full credit.
ECE62 Exam April 5, 27 Nam: Solution Scor: / This xam is closdbook. You must show ALL of your work for full crdit. Plas rad th qustions carfully. Plas chck your answrs carfully. Calculators may NOT b
More informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 7, July ISSN
Intrnational Journal of Scintific & Enginring Rsarch, Volum 6, Issu 7, July25 64 ISSN 2229558 HARATERISTIS OF EDGE UTSET MATRIX OF PETERSON GRAPH WITH ALGEBRAI GRAPH THEORY Dr. G. Nirmala M. Murugan
More informationDealing with quantitative data and problem solving life is a story problem! Attacking Quantitative Problems
Daling with quantitati data and problm soling lif is a story problm! A larg portion of scinc inols quantitati data that has both alu and units. Units can sa your butt! Nd handl on mtric prfixs Dimnsional
More informationProbability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J.
Probability and Stochastic Procsss: A Frindly Introduction for Elctrical and Computr Enginrs Roy D. Yats and David J. Goodman Problm Solutions : Yats and Goodman,4.3. 4.3.4 4.3. 4.4. 4.4.4 4.4.6 4.. 4..7
More informationThe failure of the classical mechanics
h failur of th classical mchanics W rviw som xprimntal vidncs showing that svral concpts of classical mchanics cannot b applid.  h blacbody radiation.  Atomic and molcular spctra.  h particlli charactr
More informationIterative approximate solutions of kinetic equations for reversible enzyme reactions
Vol.5, No.6, 740755 (0) doi:0.46/ns.0.5609 Natural Scinc Itrativ approximat solutions of kintic quations for rvrsibl nzym ractions Sarbaz H. A. Khoshnaw Dpartmnt of Mathmatics, Univrsity of Licstr, Licstr,
More informationBINOMIAL COEFFICIENTS INVOLVING INFINITE POWERS OF PRIMES. 1. Statement of results
BINOMIAL COEFFICIENTS INVOLVING INFINITE POWERS OF PRIMES DONALD M. DAVIS Abstract. If p is a prim and n a positiv intgr, lt ν p (n dnot th xponnt of p in n, and u p (n n/p νp(n th unit part of n. If α
More informationKernel Integral Images: A Framework for Fast NonUniform Filtering
Krnl Intgral Imags: A Framwork for Fast NonUniform Filtring Mohamd Hussin Dpt. of Computr Scinc Univrsity of Maryland mhussin@cs.umd.du Fatih Porikli Mitsubishi Elctric Rsarch Labs Cambridg, MA 0239 fatih@mrl.com
More informationUtilization of Universal Controls Analysis Tool and CFD in Parallel for a Cold Nozzle Flow
46th AIAA/ASME/SAE/ASEE Joint Propulsion Confrnc & Exhibit 58 July 010, Nashvill, TN AIAA 0106558 Utilization of Univrsal Controls Analysis Tool and CFD in Paralll for a Cold Nozzl Flow Jol M. Faur 1,
More informationBroadband AllAngle Negative Refraction by Phononic Crystals
Supplmntar Information Broadband AllAngl Ngativ Rfraction b Phononic Crstals Yang Fan Li, Fi Mng, Shiwi Zhou, MingHui Lu and Xiaodong Huang 1 Optimization algorithm and procss Bfor th optimization procss,
More informationUnit 6: Solving Exponential Equations and More
Habrman MTH 111 Sction II: Eonntial and Logarithmic Functions Unit 6: Solving Eonntial Equations and Mor EXAMPLE: Solv th quation 10 100 for. Obtain an act solution. This quation is so asy to solv that
More information7.4 Potential Difference and Electric Potential
7.4 Potntial Diffrnc and Elctric Potntial In th prvious sction, you larnd how two paralll chargd surfacs produc a uniform lctric fild. From th dfinition of an lctric fild as a forc acting on a charg, it
More informationINTRODUCTION TO AUTOMATIC CONTROLS INDEX LAPLACE TRANSFORMS
adjoint...6 block diagram...4 clod loop ytm... 5, 0 E()...6 (t)...6 rror tady tat tracking...6 tracking...6...6 gloary... 0 impul function...3 input...5 invr Laplac tranform, INTRODUCTION TO AUTOMATIC
More informationEstimation of odds ratios in Logistic Regression models under different parameterizations and Design matrices
Advancs in Computational Intllignc, ManMachin Systms and Cybrntics Estimation of odds ratios in Logistic Rgrssion modls undr diffrnt paramtrizations and Dsign matrics SURENDRA PRASAD SINHA*, LUIS NAVA
More informationThus, because if either [G : H] or [H : K] is infinite, then [G : K] is infinite, then [G : K] = [G : H][H : K] for all infinite cases.
Homwork 5 M 373K Solutions Mark Lindbrg and Travis Schdlr 1. Prov that th ring Z/mZ (for m 0) is a fild if and only if m is prim. ( ) Proof by Contrapositiv: Hr, thr ar thr cass for m not prim. m 0: Whn
More informationCategory Theory Approach to Fusion of WaveletBased Features
Catgory Thory Approach to Fusion of WavltBasd Faturs Scott A. DLoach Air Forc Institut of Tchnology Dpartmnt of Elctrical and Computr Enginring WrightPattrson AFB, Ohio 45433 Scott.DLoach@afit.af.mil
More informationAnother view for a posteriori error estimates for variational inequalities of the second kind
Accptd by Applid Numrical Mathmatics in July 2013 Anothr viw for a postriori rror stimats for variational inqualitis of th scond kind Fi Wang 1 and Wimin Han 2 Abstract. In this papr, w giv anothr viw
More informationSection 11.6: Directional Derivatives and the Gradient Vector
Sction.6: Dirctional Drivativs and th Gradint Vctor Practic HW rom Stwart Ttbook not to hand in p. 778 # 4 p. 799 # 45 7 9 9 35 37 odd Th Dirctional Drivativ Rcall that a b Slop o th tangnt lin to th
More informationEAcos θ, where θ is the angle between the electric field and
8.4. Modl: Th lctric flux flows out of a closd surfac around a rgion of spac containing a nt positiv charg and into a closd surfac surrounding a nt ngativ charg. Visualiz: Plas rfr to Figur EX8.4. Lt A
More informationMultiobjective Optimization of web profile of railway wheel using Bidirectional Evolutionary Structural Optimization
JCAMECH Vol. XX, No. XX, XX 20XX, pp 307318 DOI: 10.22059/jcamch.2017.237353.160 Multiobjctiv Optimization of wb profil of railway whl using Bidirctional Evolutionary Structural Optimization Ali Asghar
More informationChapter 7b Electron Spin and Spin Orbit Coupling
Wintr 3 Chm 356: Introductory Quantum Mchanics Chaptr 7b Elctron Spin and Spin Orbit Coupling... 96 H atom in a Magntic Fild: Elctron Spin... 96 Total Angular Momntum... 3 Chaptr 7b Elctron Spin and
More informationA=P=E MA=N Alpha particle Beta Particle. Periodic table
Nam Pr. Atomic Structur/Nuclar Chmistry (Ch. 3 & 21) OTHS Acadmic Chmistry Objctivs: Undrstand th xprimntal dsign and conclusions usd in th dvlopmnt of modrn atomic thory, including Dalton's Postulats,
More informationADAPTIVE LIMIT ANALYSIS USING DEVIATORIC FIELDS
VI Intrnational Confrnc on Adaptiv Modling and Simulation ADMOS 2013 J. P. Moitinho d Almida, P. Díz, C. Tiago and N. Parés (Eds) ADAPTIVE LIMIT ANALYSIS USING DEVIATORIC FIELDS ANDREI V. LYAMIN *, KRISTIAN
More informationUNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS
UNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This papr is concrnd with rsidual typ a postriori rror stimators
More informationFixedPoint HarmonicBalanced Method for Nonlinear Eddy Current Problems
Intrnational Journal of Enrgy and Powr Enginring 206; 5(): 374 Publishd onlin Octobr 4, 205 (http://www.scincpublishinggroup.com/j/ijp) doi: 0.648/j.ijp.s.2060500.5 ISSN: 2326957X (Print); ISSN: 2326960X
More informationCS 6353 Compiler Construction, Homework #1. 1. Write regular expressions for the following informally described languages:
CS 6353 Compilr Construction, Homwork #1 1. Writ rgular xprssions for th following informally dscribd languags: a. All strings of 0 s and 1 s with th substring 01*1. Answr: (0 1)*01*1(0 1)* b. All strings
More informationAnalysis of Gas Expansion Required for Hydrate Formation
J. Rajnauth: Analysis of Gas Expansion Rquird for Hydrat Formation 22 ISSN 1 7924 Th Journal of th Association of Profssional Enginrs of Trinidad and Tobago Vol.42, No.1, April/May 214, pp.2226 Analysis
More information