Use of Proper Closure Equations in Finite Volume Discretization Schemes

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "Use of Proper Closure Equations in Finite Volume Discretization Schemes"

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 thrmo-fluid 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 prssur-basd primitiv variabl formulation (Acharya, 2007), in which th flow continuity quation is considrd as a constraint on fluid prssur. rssur-basd 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 co-locatd grid arrangmnt. Th discussion of finit volum discrtization schms in this articl is limitd to prssur-basd formulations on co-locatd 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 two-dimnsional grids in an lmnt-basd 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 lvl-1 approximation, Eq. (1) is convrtd to a smi-discrt form comprising of both nodal and intgration point variabls. Lvl-1 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 lvl-2 approximation, Eq. (3) is convrtd to a fully-discrt 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 intr-quation couplings should also b numrically modld. Two-dimnsional, incomprssibl Navir-Stoks 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 x-momntum balanc quation, v in th y-momntum 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 lmnt-basd 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 intr-rlatd. 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 lvl-2 approximation of Navir-Stoks quations in th contxt of co-locatd, prssur-basd 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 lvl-2 approximations ar rviwd and th mthod of propr closur quations (MCE), proposd by th first author, is introducd. To avoid th complxitis associatd with multi-dimnsional problms and to focus on th natur of th approximations in a simpl stting, on dimnsional quations and tst cass ar discussd in dtails. Two-dimnsional rsults ar only brifly prsntd to show th applicability of th MCE in multi-dimnsional problms. 2. Th smi-discrt form of quations in a 1D contxt A on-dimnsional variabl-ara duct is shown in Fig. 2. Th complt st of govrning quations for unstady 1D-inviscid 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 smi-discrt 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 on-dimnsional variabl ara duct. 3. Closur quations for intgration point variabls Aftr obtaining th smi-discrt 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 smi-discrt quations ar convrtd to fully-discrt quations. Ths ip-quations 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 ip-valu in a uniform grid: On non-uniform 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 2nd-ordr accurat, but is unboundd so that non-physical oscillations may appar in rgions with high gradints. To show th problm associatd with this schm, considr th 1-D inviscid incomprssibl flow in a constant ara duct shown in Fig. 3. Th smi-discrt 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 1-D Co-locatd grid in a constant cross sctional ara duct. Using simpl symmtric intrpolation, Eqs. (13) and (14) ar convrtd to th following fully-discrt 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 taylor-sris xpansion It is logical to assum that convctd quantitis in a thrmo-fluid 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 co-locatd 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 mass-carrying ( û ) and convctd or transportd ( u ) vlocity componnts at th intgration points. Th convcting vlocity, usd in th smi-discrt 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 prssur-vlocity 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 co-locatd grids. Thy rtaind a unifid dfinition for th intgration-point vlocity componnts, but argud that th closur for u should not b obtaind solly by purly mathmatical upwind schms. Thy suggstd th non-consrvativ 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 1-D 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 co-locatd 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 1-D tst cas show that this mthod works wll and succssfully supprsss any oscillatory numrical solution vn in tim-dpndd 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 prssur-basd 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 usr-dfind 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 momntum-variabl 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) On-dimnsional 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 1-D inviscid incomprssibl and comprssibl flow problms (Ashrafizadh t al., 2008, 2009) and stady 2-D viscous incomprssibl flow problms on both structurd and unstructurd grids (Alinia, 2011; Bakhtiari, 2008). In this mthod, following th physical-basd approach proposd by Schnidr and Raw (Schnidr t al., 1987) and modifications proposd by Karimian and Darbandi (Darbandi, 1997; Karimian, 1994a), a combination of non-consrvativ form of momntum and continuity quations is usd to obtain an intrpolation formula for th ip-vlocity: 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 CDS-basd 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 (Rzvani-2010) On-dimnsional 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 unstady-comprssibl 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 unstady-comprssibl 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 non-consrvativ 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 non-consrvativ 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 smi-implicit (sgrgatd) solution stratgy. p f

13 Us of ropr Closur Equations in Finit Volum Discrtization Schms Two-dimnsional 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. On-dimnsional 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 prssur-vlocity coupling in MCE. Quasi 1D flow in a convrgnt-divrgnt 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 cll-facs 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 control-volum 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 non-zro valu of p. Not that th notation N-D in ths figurs stands for Non-Dimnsional 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 control-volum 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 control-volum momntum sourc/sink Th quasi-1-d tst cas Idal flow in a convrging-divrging nozzl is anothr tst cas to xamin th prformanc of MCE. Figur 7 compars th non-dimnsional vlocity and prssur distributions along th nozzl with th analytical solutions. Excllnt agrmnt is obsrvd. Fig. 7. Vlocity and prssur fild in a 1D convrging-divrging 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 convrgnt-divrgnt nozzl, stady comprssibl flow with normal shock in convrgnt-divrgnt nozzl and unstady comprssibl flow with normal shock and xpansion wavs in a shock tub Stady subsonic comprssibl flow in a convrgnt-divrgnt nozzl Rsults of solving subsonic comprssibl flow in a symmtric convrgnt-divrgnt 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, non-dimnsional 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 convrgnt-divrgnt 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, non-dimnsional 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 non-dimnsionalizd using th low prssur sid valus. Good agrmnt with th xact solution is obsrvd. Fig. 9. Stady flow with shock wavs in a convrgnt-divrgnt nozzl. 5. Two-dimnsional 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 lid-drivn cavity Flow in a lid-drivn cavity, shown in Fig. 12, is a widly usd tst cas for 2D stady incomprssibl flows govrnd by Navir-Stoks 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. Lid-drivn cavity boundary conditions with an schmatic of flow faturs. Fig. 13. Gnratd unstructurd msh for lid-drivn cavity tst cas. 5.2 Flow ovr a backward facing-stp 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 r-attachmnt 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 r-attachmnt 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 backward-stp 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 facing-stp flow problm.

20 80 Finit Volum Mthod owrful Mans of Enginring Dsign Fig. 15. Calculatd stramlins in th backward-stp 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 prssur-basd finit volum discrtization schms for th solution of fluid flow problms ar xplaind and compard in th contxt of co-locatd 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 co-locatd grids which mploys only on dfinition for th intgration point vlocity and succssfully solvs incomprssibl as wll as comprssibl flow problms. Simpl on-dimnsional 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 Finit-Volum 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 Cll-Fac Vlocitis and Thir Effcts on th rssur-vlocity 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 Smi-Coupld 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 Co-locatd rssur-basd Finit Volum Mthod, 16th Annual Confrnc of th CFD Socity of Canada, Saskatoon, Saskatchwan, Canada Ashrafizadh, A.; Rzvani M.; & Bakhtiari, B. (2009). rssur-vlocity 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). Two-Dimnsional 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 rssur-basd 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 2-D Stady Incomprssibl Drivn Cavity Flow at High Rynolds Numbrs, Int. J. Numr. Mth. Fluids, Vol. 48, pp Erturk, E. (2008). Numrical Solutions of 2-D 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 Control-Volum-Basd Finit-Elmnt Formulation to th Shock Tub roblm, AIAA Journal, Vol.33 No.1, pp Karimian, S. M. H. & Schnidr, G. E. (1994b). rssur-basd 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 rssur-basd 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 Intr-Equation Couplings in All-Spd 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 Finit-Elmnt Mthod For hat Transfr and Fluid Flow Using Colocatd Variabls-1. Computational rocdur. Numrical Hat Transfr, part B, Vol. 11, pp Schnidr, G. E. & Raw, M. J. (1987b). Control Volum Finit-Elmnt Mthod For hat Transfr and Fluid Flow Using Colocatd Variabls-2. 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 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 information

Finite element discretization of Laplace and Poisson equations

Finite 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 2D-Discrtization

More information

3 Finite Element Parametric Geometry

3 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 information

Elements of Statistical Thermodynamics

Elements 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 information

Analysis of Convection-Diffusion Problems at Various Peclet Numbers Using Finite Volume and Finite Difference Schemes Anand Shukla

Analysis of Convection-Diffusion Problems at Various Peclet Numbers Using Finite Volume and Finite Difference Schemes Anand Shukla Mathmatical Thory and Modling.iist.org ISSN 4-5804 (apr) ISSN 5-05 (Onlin) Vol., No.6, 01-Slctd from Intrnational Confrnc on Rcnt Trnds in Applid Scincs ith nginring Applications Analysis of Convction-iffusion

More information

Nusselt number correlations for simultaneously developing laminar duct flows of liquids with temperature dependent properties

Nusselt 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 information

u 3 = u 3 (x 1, x 2, x 3 )

u 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 information

Sliding Mode Flow Rate Observer Design

Sliding 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 information

The Matrix Exponential

The 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 information

4.2 Design of Sections for Flexure

4.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 information

16. Electromagnetics and vector elements (draft, under construction)

16. 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 information

Element connectivity parameterization method for the stress based topology optimization for geometrically nonlinear structure

Element connectivity parameterization method for the stress based topology optimization for geometrically nonlinear structure 0 th World Congrss on Structural and Multidisciplinary Optimization May 9-4, 03, Orlando, Florida, USA Elmnt connctivity paramtrization mthod for th strss basd topology optimization for gomtrically nonlinar

More information

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

The 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 information

CE 530 Molecular Simulation

CE 530 Molecular Simulation CE 53 Molcular Simulation Lctur 8 Fr-nrgy calculations David A. Kofk Dpartmnt of Chmical Enginring SUNY Buffalo kofk@ng.buffalo.du 2 Fr-Enrgy Calculations Uss of fr nrgy Phas quilibria Raction quilibria

More information

CHAPTER 1. Introductory Concepts Elements of Vector Analysis Newton s Laws Units The basis of Newtonian Mechanics D Alembert s Principle

CHAPTER 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 information

5.80 Small-Molecule Spectroscopy and Dynamics

5.80 Small-Molecule Spectroscopy and Dynamics MIT OpnCoursWar http://ocw.mit.du 5.80 Small-Molcul 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 information

MCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17)

MCB137: 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 information

On the Hamiltonian of a Multi-Electron Atom

On the Hamiltonian of a Multi-Electron Atom On th Hamiltonian of a Multi-Elctron 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 information

1 General boundary conditions in diffusion

1 General boundary conditions in diffusion Gnral boundary conditions in diffusion Πρόβλημα 4.8 : Δίνεται μονοδιάτατη πλάκα πάχους, που το ένα άκρο της κρατιέται ε θερμοκραία T t και το άλλο ε θερμοκραία T 2 t. Αν η αρχική θερμοκραία της πλάκας

More information

Observer Bias and Reliability By Xunchi Pu

Observer 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 information

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator.

Exam 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 information

10. The Discrete-Time Fourier Transform (DTFT)

10. The Discrete-Time Fourier Transform (DTFT) Th Discrt-Tim Fourir Transform (DTFT Dfinition of th discrt-tim Fourir transform Th Fourir rprsntation of signals plays an important rol in both continuous and discrt signal procssing In this sction w

More information

MEASURING HEAT FLUX FROM A COMPONENT ON A PCB

MEASURING 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 information

A dynamic boundary model for implementation of boundary conditions in lattice-boltzmann method

A dynamic boundary model for implementation of boundary conditions in lattice-boltzmann method Journal of Mchanical Scinc and Tchnology Journal of Mchanical Scinc and Tchnology 22 (28) 92~2.springrlink.com/contnt/738-494x A dynamic boundary modl for implmntation of boundary conditions in lattic-boltzmann

More information

Available online at ScienceDirect. Procedia Engineering 126 (2015 )

Available 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 information

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design

MAE4700/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 information

Types of Transfer Functions. Types of Transfer Functions. Types of Transfer Functions. Ideal Filters. Ideal Filters

Types 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 tim-domain classification of an LTI digital transfr function is basd on th lngth of its impuls rspons h[n]:

More information

Introduction to Condensed Matter Physics

Introduction to Condensed Matter Physics Introduction to Condnsd Mattr Physics pcific hat M.P. Vaughan Ovrviw Ovrviw of spcific hat Hat capacity Dulong-Ptit Law Einstin modl Dby modl h Hat Capacity Hat capacity h hat capacity of a systm hld at

More information

is an appropriate single phase forced convection heat transfer coefficient (e.g. Weisman), and h

is 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 Jns-Lotts Corrlation b) T clad surfac tmpratur assuming t Tom Corrlation c) T clad surfac tmpratur assuming

More information

Differentiation of Exponential Functions

Differentiation 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 information

A Prey-Predator Model with an Alternative Food for the Predator, Harvesting of Both the Species and with A Gestation Period for Interaction

A Prey-Predator 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 Pry-Prdator 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 information

Search sequence databases 3 10/25/2016

Search 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 information

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches. Subjct Chmistry Papr No and Titl Modul No and Titl Modul Tag 8/ Physical Spctroscopy / Brakdown of th Born-Oppnhimr approximation. Slction ruls for rotational-vibrational transitions. P, R branchs. CHE_P8_M

More information

A nonequilibrium molecular dynamics simulation of evaporation

A 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 information

Modeling of the Plasma Flow in High-Power TAL

Modeling of the Plasma Flow in High-Power TAL 39th AIAA/ASME/SAE/ASEE Joint Propulsion Confrnc and Exhibit 20-23 July 2003, Huntsvill, Alabama AIAA 2003-4701 AIAA-2003-4701 Modling of th Plasma Flow in High-Powr TAL Michal Kidar, Iain D. Boyd and

More information

2.3 Matrix Formulation

2.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 information

CO-ORDINATION OF FAST NUMERICAL RELAYS AND CURRENT TRANSFORMERS OVERDIMENSIONING FACTORS AND INFLUENCING PARAMETERS

CO-ORDINATION OF FAST NUMERICAL RELAYS AND CURRENT TRANSFORMERS OVERDIMENSIONING FACTORS AND INFLUENCING PARAMETERS CO-ORDINATION 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 information

What are those βs anyway? Understanding Design Matrix & Odds ratios

What 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 information

Lecture # 12: Shock Waves and De Laval Nozzle

Lecture # 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 information

BINOMIAL COEFFICIENTS INVOLVING INFINITE POWERS OF PRIMES

BINOMIAL 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 information

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.

SECTION 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 information

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002

ME 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 information

Calculus concepts derivatives

Calculus 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 information

A 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 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 information

Sec 2.3 Modeling with First Order Equations

Sec 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 information

2. Laser physics - basics

2. 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 information

Collisions between electrons and ions

Collisions 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 Fokkr-Planck collision

More information

The pn junction: 2 Current vs Voltage (IV) characteristics

The 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 p-typ Dpltion rgion n-typ Elctron movmnt across th junction: 1. n

More information

Differential Equations

Differential 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 information

Ultimate strength analysis & design of residential slabs on reactive soil

Ultimate 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 information

A SPLITTING METHOD USING DISCONTINUOUS GALERKIN FOR THE TRANSIENT INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

A SPLITTING METHOD USING DISCONTINUOUS GALERKIN FOR THE TRANSIENT INCOMPRESSIBLE NAVIER-STOKES 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 NAVIER-STOKES

More information

PROOF OF FIRST STANDARD FORM OF NONELEMENTARY FUNCTIONS

PROOF 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 ISSN-319-8354(E) PROOF OF FIRST STANDARD FORM OF NONELEMENTARY FUNCTIONS 1 Dharmndra

More information

EFFECT OF BALL PROPERTIES ON THE BALL-BAT COEFFICIENT OF RESTITUTION

EFFECT OF BALL PROPERTIES ON THE BALL-BAT COEFFICIENT OF RESTITUTION EFFECT OF BALL PROPERTIES ON THE BALL-BAT COEFFICIENT OF RESTITUTION A. M. NATHAN 1 AND L. V. SMITH 2 1 Univrsity of Illinois, 1110 W. Grn Strt, Urbana, IL 61801, USA, E-mail: a-nathan@illinois.du 2 Washington

More information

Proximity Templates for Modeling of Skin and Proximity Effects on Packages and High Frequency Interconnect

Proximity 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 Sangiovanni-incntlli Univ. of California, Brkly

More information

VALIDATION OF FINITE ELEMENT PROGRAM FOR JOURNAL BEARINGS -- STATIC AND DYNAMIC PROPERTIES

VALIDATION 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 information

INC 693, 481 Dynamics System and Modelling: The Language of Bound Graphs Dr.-Ing. Sudchai Boonto Assistant Professor

INC 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 information

Definition1: The ratio of the radiation intensity in a given direction from the antenna to the radiation intensity averaged over all directions.

Definition1: 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 information

PHA 5127 Answers Homework 2 Fall 2001

PHA 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 information

Pipe flow friction, small vs. big pipes

Pipe 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 information

Waves in cavities such as vehicle compartments, rooms or ducts

Waves 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 information

Engineering 323 Beautiful HW #13 Page 1 of 6 Brown Problem 5-12

Engineering 323 Beautiful HW #13 Page 1 of 6 Brown Problem 5-12 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 information

Determination of Vibrational and Electronic Parameters From an Electronic Spectrum of I 2 and a Birge-Sponer Plot

Determination of Vibrational and Electronic Parameters From an Electronic Spectrum of I 2 and a Birge-Sponer Plot 5 J. Phys. Chm G Dtrmination of Vibrational and Elctronic Paramtrs From an Elctronic Spctrum of I 2 and a Birg-Sponr Plot 1 15 2 25 3 35 4 45 Dpartmnt of Chmistry, Gustavus Adolphus Collg. 8 Wst Collg

More information

VII. Quantum Entanglement

VII. 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 information

Bifurcation Theory. , a stationary point, depends on the value of α. At certain values

Bifurcation 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 information

Gradebook & Midterm & Office Hours

Gradebook & Midterm & Office Hours Your commnts So what do w do whn on of th r's is 0 in th quation GmM(1/r-1/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 information

REGISTER!!! The Farmer and the Seeds (a parable of scientific reasoning) Class Updates. The Farmer and the Seeds. The Farmer and the Seeds

REGISTER!!! 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 information

u 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

u 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 information

Physics 312 First Pledged Problem Set

Physics 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 information

Notes on Fermi-Dirac Integrals 2 nd Edition. Raseong Kim and Mark Lundstrom Network for Computational Nanotechnology Purdue University

Notes on Fermi-Dirac Integrals 2 nd Edition. Raseong Kim and Mark Lundstrom Network for Computational Nanotechnology Purdue University Nots on rmi-dirac Intgrals nd Edition Rasong Kim and Mark Lundstrom Ntwork for Computational Nanotchnology Purdu Univrsity Sptmbr 3, 8 1. Introduction rmi-dirac intgrals appar frquntly in smiconductor

More information

Design Guidelines for Quartz Crystal Oscillators. R 1 Motional Resistance L 1 Motional Inductance C 1 Motional Capacitance C 0 Shunt Capacitance

Design 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 information

Sara Godoy del Olmo Calculation of contaminated soil volumes : Geostatistics applied to a hydrocarbons spill Lac Megantic Case

Sara Godoy del Olmo Calculation of contaminated soil volumes : Geostatistics applied to a hydrocarbons spill Lac Megantic Case wwwnvisol-canadaca 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 information

SUMMER 17 EXAMINATION

SUMMER 17 EXAMINATION (ISO/IEC - 7-5 Crtifid) SUMMER 7 EXAMINATION Modl wr jct Cod: Important Instructions to aminrs: ) Th answrs should b amind by ky words and not as word-to-word as givn in th modl answr schm. ) Th modl answr

More information

Maxwellian Collisions

Maxwellian Collisions Maxwllian Collisions Maxwll ralizd arly on that th particular typ of collision in which th cross-sction varis at Q rs 1/g offrs drastic siplifications. Intrstingly, this bhavior is physically corrct for

More information

ECE602 Exam 1 April 5, You must show ALL of your work for full credit.

ECE602 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 closd-book. 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 information

International Journal of Scientific & Engineering Research, Volume 6, Issue 7, July ISSN

International Journal of Scientific & Engineering Research, Volume 6, Issue 7, July ISSN Intrnational Journal of Scintific & Enginring Rsarch, Volum 6, Issu 7, July-25 64 ISSN 2229-558 HARATERISTIS OF EDGE UTSET MATRIX OF PETERSON GRAPH WITH ALGEBRAI GRAPH THEORY Dr. G. Nirmala M. Murugan

More information

Dealing with quantitative data and problem solving life is a story problem! Attacking Quantitative Problems

Dealing 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 information

Probability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J.

Probability 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 information

The failure of the classical mechanics

The 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 blac-body radiation. - Atomic and molcular spctra. - h particl-li charactr

More information

Iterative approximate solutions of kinetic equations for reversible enzyme reactions

Iterative approximate solutions of kinetic equations for reversible enzyme reactions Vol.5, No.6, 740-755 (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 information

BINOMIAL COEFFICIENTS INVOLVING INFINITE POWERS OF PRIMES. 1. Statement of results

BINOMIAL 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 information

Kernel Integral Images: A Framework for Fast Non-Uniform Filtering

Kernel Integral Images: A Framework for Fast Non-Uniform Filtering Krnl Intgral Imags: A Framwork for Fast Non-Uniform 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 information

Utilization of Universal Controls Analysis Tool and CFD in Parallel for a Cold Nozzle Flow

Utilization of Universal Controls Analysis Tool and CFD in Parallel for a Cold Nozzle Flow 46th AIAA/ASME/SAE/ASEE Joint Propulsion Confrnc & Exhibit 5-8 July 010, Nashvill, TN AIAA 010-6558 Utilization of Univrsal Controls Analysis Tool and CFD in Paralll for a Cold Nozzl Flow Jol M. Faur 1,

More information

Broadband All-Angle Negative Refraction by Phononic Crystals

Broadband All-Angle Negative Refraction by Phononic Crystals Supplmntar Information Broadband All-Angl Ngativ Rfraction b Phononic Crstals Yang Fan Li, Fi Mng, Shiwi Zhou, Ming-Hui Lu and Xiaodong Huang 1 Optimization algorithm and procss Bfor th optimization procss,

More information

Unit 6: Solving Exponential Equations and More

Unit 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 information

7.4 Potential Difference and Electric Potential

7.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 information

INTRODUCTION TO AUTOMATIC CONTROLS INDEX LAPLACE TRANSFORMS

INTRODUCTION 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 information

Estimation of odds ratios in Logistic Regression models under different parameterizations and Design matrices

Estimation of odds ratios in Logistic Regression models under different parameterizations and Design matrices Advancs in Computational Intllignc, Man-Machin Systms and Cybrntics Estimation of odds ratios in Logistic Rgrssion modls undr diffrnt paramtrizations and Dsign matrics SURENDRA PRASAD SINHA*, LUIS NAVA

More information

Thus, 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.

Thus, 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 information

Category Theory Approach to Fusion of Wavelet-Based Features

Category Theory Approach to Fusion of Wavelet-Based Features Catgory Thory Approach to Fusion of Wavlt-Basd Faturs Scott A. DLoach Air Forc Institut of Tchnology Dpartmnt of Elctrical and Computr Enginring Wright-Pattrson AFB, Ohio 45433 Scott.DLoach@afit.af.mil

More information

Another view for a posteriori error estimates for variational inequalities of the second kind

Another 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 information

Section 11.6: Directional Derivatives and the Gradient Vector

Section 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 # 4-5 7 9 9 35 37 odd Th Dirctional Drivativ Rcall that a b Slop o th tangnt lin to th

More information

EAcos θ, where θ is the angle between the electric field and

EAcos θ, 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 information

Multi-objective Optimization of web profile of railway wheel using Bi-directional Evolutionary Structural Optimization

Multi-objective Optimization of web profile of railway wheel using Bi-directional Evolutionary Structural Optimization JCAMECH Vol. XX, No. XX, XX 20XX, pp 307-318 DOI: 10.22059/jcamch.2017.237353.160 Multi-objctiv Optimization of wb profil of railway whl using Bi-dirctional Evolutionary Structural Optimization Ali Asghar

More information

Chapter 7b Electron Spin and Spin- Orbit Coupling

Chapter 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 information

A=P=E M-A=N Alpha particle Beta Particle. Periodic table

A=P=E M-A=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 information

ADAPTIVE LIMIT ANALYSIS USING DEVIATORIC FIELDS

ADAPTIVE 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 information

UNIFIED 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 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 information

Fixed-Point Harmonic-Balanced Method for Nonlinear Eddy Current Problems

Fixed-Point Harmonic-Balanced Method for Nonlinear Eddy Current Problems Intrnational Journal of Enrgy and Powr Enginring 206; 5(-): 37-4 Publishd onlin Octobr 4, 205 (http://www.scincpublishinggroup.com/j/ijp) doi: 0.648/j.ijp.s.2060500.5 ISSN: 2326-957X (Print); ISSN: 2326-960X

More information

CS 6353 Compiler Construction, Homework #1. 1. Write regular expressions for the following informally described languages:

CS 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 information

Analysis of Gas Expansion Required for Hydrate Formation

Analysis 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.22-26 Analysis

More information