Comparative Study of Netonian Sinusoidal Blood Flow through Normal and Stenosed Carotid Artery

Size: px
Start display at page:

Download "Comparative Study of Netonian Sinusoidal Blood Flow through Normal and Stenosed Carotid Artery"

Transcription

1 IOSR Jornal of Mahemacs (IOSR-JM) e-issn: , p-issn: X. Volme 13, Isse 6 Ver. III (Nov. - Dec. 2017), PP Comparave Sdy of Neonan Snsodal Blood Flow hrogh Normal and Senosed Carod Arery Mohammad Mar Rahman 1,a),Md. Anwar Hossan 1,b), Kharzzaman Mamn 1,c), Mos. Nasrn Aher 1,d). 1 (Deparmen of Mahemacs, Dhaa Unversy of Engneerng and Technology, Gazpr-1700, Bangladesh) a)correspondng ahor: marrahman1975@gmal.com b) anwarmah23@gmal.com,c) presdenmamn@gmal.com,d) nasrn @yahoo.com Absrac: A nmercal smlaon s performed o nvesgae Newonan snsodal flows behavor on hree dmensonal dealzed carod arery () and sngle senosed (75% by area) carod arery(s). The wall vessel s se as rgd drng smlaon. Bfrcaed blood vessel are smlaed by sng hree-dmensonal flow analyss. Snsodal and parabolc velocy profles are se o o f he condons of nle bondares of arery. The nvesgaon has a Reynolds nmber range of 90 o Low Reynolds nmber model has been sed as governng eqaon. The nvesgaon has been carred o o characerze he flow behavor of blood n wo geomery, namely, () Normal carod arery () and () Senosed carod arery (S). The Newonan model has been sed o sdy he physcs of fld. The fndngs of he wo models are horoghly compared n order o observe here behavoral seqence of flows. The nmercal resls were presened n erms of velocy, pressre, wall shear sress dsrbons and cross seconal veloces as well as he sreamlnes conor. Senoss dsrbs he normal paern of blood flow hrogh he arery as redced area. A senoss regon velocy and pea Reynolds nmber rapdly ncrease and Reynolds nmber reach ransonal and rblen regon. These flow flcaon and rblence have bad effec o he blood vessel whch maes o accelerae he progress of senoss. Keywords : Aheroscleross, Carod arery, Snsodal blood flow, Senoss, Vscosy Dae of Sbmsson: Dae of accepance: I. Inrodcon Sroe s a prevalen and severe medcal condon ha mposes an mmense loss of healh and heavy economc brden worldwde. Globally, sroe s he second leadng case of deah and he leadng case of long-erm dsably [1, 2]. The adverse effecs are boh eensve and long-lasng, as one hrd of new sroes resl n deah and mos sroe srvvors are afflced by a nerologcal dsably ha reqres rehablaon or long-erm care [3]. Moreover, once a sroe has occrred, he rs of sroe recrrence ncreases sbsanally [4, 5]. By far, he mos common form of sroe s de o schema a resrcon or shorage n cerebral blood spply whch accons for 78% of sroes [6]. Ischema s mos ofen cased by a hrombs or clo ha has emblazed and ravelled o he bran casng areral occlson of a narrow cerebral vessel [7, 10] The hrombs ypcally orgnaes from a promal locaon n he body, sch as he hear (cardogenc embolsm) or from 2 aheroscleroc plaqe n a large promal arery (arery-o-arery embolsm), commonly he aora or nernal carod arery (I). Aheroscleross localzed n he carod areres s no only prevalen, b also represens a sgnfcan rs for clncal cerebrovasclar evens. The carod arery s he mos freqenly nvolved se n large arery aherohromboc sroe. Carod arery aherohromboss s esmaed o be responsble for 10 30% of all sroes [4, 6-8]. he geomery of he flow vessel s a ey deermnan of s flow characerscs. I has been well esablshed ha senoses, bfrcaons, and srface roghness are assocaed wh flow dsrbances and comple flow [9-11].The nflence of vessel geomery on flow s eemplfed n he areral vasclare. Even n he normal, healhy vasclare, vessel bends and bfrcaons, as well as he plsale flow of he cardac cycle, nrodce secondary flows and oscllang shear sresses. Aheroscleross adms ye more complcaons, parclarly where areral dsease s severe enogh o case a vessel senoss [15]. Whereas flow n normal areres s comple b ypcally lamnar, senoss are assocaed wh shear layers, eddy formaon, dsal flow dsrbances [15], and velocy flcaons [16]. Toros, rreglar vessel srfaces are lely o frher conrbe o flow dsrbances [9, 10]. A he senosed carod arery, flow frs enconers an asymmerc bfrcang ncon, hen mmedaely goes hrogh a senoc flow consrcon, followed downsream by a sdden epanson no he carod blb. The bfrcaon mposes a flow dvson ha affecs he bondary layers along boh he nsde and osde bends, whch s dependen on he eac branch geomery. DOI: Page

2 The comparson of senosed flow behavor wh he normal one can provde he proper ndersandng of nderlyng mechansm behnd he developmen of aheroscleross. The flow rns o be abnormal n he redced cross seconal area of he arery saed by Kader and Shenory [7]. The flow behavor n he senosed arery s qe dfferen n comparson o he normal one. Sress and ressance o flow s mch hgher n senosed arery. Cha and Shread [9] fond ha he flow hrogh he consrced be s characerzed by hgh velocy e generaed a consrced regon. Kader and Shenory [7] fond he resls from nmercal smlaon whch demonsraed ha velocy and senoc e lengh ncreases n ncreasng he severy of senoss. Ther resls also demonsrae ha he 3D senoc CFD model s capable of predcng he changes n flow behavor for ncreased severy of senoss. Yong e al. [17] sded he wall shear sress and pressre graden n he senoss and evalaed he case of plaqe rpre. The ahors sded plsale blood flow hrogh he senoss wh elasc wall o observe he lmen movemen. Accordng o her sdy he pea WSS occrred s before mnmm lmen poson. Pno e al. [18] performed nmercal smlaon assmng a physologcal plsale flow hrogh dfferen models of senoss. In case of sbec-specfc anaomcally realsc senosed carod bfrcaon sbeced o plsale nle condon, he smlaon resls demonsraed he rapd flcaon of velocy and pressre n pos-senoc regon by S. Lee and S. Lee [16]. Ahmed and Gddens [19] sded boh seady and plsale flow hrogh 25%, 50% and 75% consrcon of a rgd be where Reynolds nmber ranged from 500 o In anoher sdy, he senoss havng dfferen geomerc profles le rapezm, sem-ellpse and rangle were analyzed consderng he Non-Newonan behavor by Lorenzn and Casalena [20]. Mos of he researcher consdered senosed carod arery or normal one. B n hs sdy we se normal and senosed carod arery ogeher o ndersand he flow behavor correcly. The nvesgaon has a Reynolds nmber range of 90 o The nmercal resls are presened n erms of velocy, pressre, wall shear sress dsrbons.ths paper s organzed as follows. In secon Model descrpon we dscss he geomery and blood properes. In secon Governng eqaon and Nmercal mehod we dscss he governng eqaon, bondary condon, grd ndependence chec, valdaon and nmercal scheme. In secon Resls and dscsson we dscss and smmares he wss, velocy and pressre dsrbons hrogh he carod arery. II. Model descrpon 2.1 Geomery Three dmensonal carod arery () and 75% senosed(by area) carod arery (S) are sed as geomery for hs sdy shown n Fg. 1 (a). Bfrcaon angles I and E 18.5 ; he rads of he blood vessel C(4mm), mddle of senoss (2mm), end of I(2.5mm), end of E(2mm); lengh of he arery from VI() are a poson () 72mm, () 80mm, (l) 88mm, (m) and (n) 103mm, (o)and (p) 118mm and (q) and (r) 142mm; percen of senoss (75% by area). Fg. 1(b) shows he mesh n cross-seconal plane of a carod arery. (a) (b) Fgre 1: (a) Model of a carod arery ; (b) mesh n cross seconal plane of a carod arery. 2.2 Blood properes The densy of he blood s 1050gm 3. In a Newonan model for he blood vscosy, he vale of μ s 3 reaed as a consan sally se o pa. s. DOI: Page

3 DOI: Page III. Governng Eqaon And Nmercal Mehod 3.1 Governng eqaon De o consrcon, blood passes hrogh he hroa and pos senoc regon wh hgh velocy. Flow velocy a pos senoc regon ncreases b pressre of ha regon decreases. So, moderae and severe areral senoss resrcs he flow severely wh creang rblen flow characerscs. For hs regon, rblen flow s epeced a pos senoc regon especally for 75% senoc arery. Neher lamnar flow modelng nor sandard wo-eqaon models are sable for hs nd of blood flow. So, Wlco low-re rblence model s more accepable for flow analyss fond by Varghese and Franel [22]. Therefore, low Re rblen model s aen for calclaon. Now, he Naver-Soes eqaon can be gven by- 0 (1) p 2 1 (2) Snce each erm of hs eqaon s me averaged, he eqaon s referred o as a Reynolds averaged Naver-Soes (RANS) eqaon. Drng hs procedre, several addonal nnown parameers appear whch reqre addonal eqaons o be nrodced as rblence models. The se of RANS eqaons are- 0 (3) p (4) In hs eqaon s an addonal erm nown as he Reynolds s sress ensor, whch can be appromaed by sng Bossnesq s hypohess- 3 2 (5) Eddy vscosy can be modeled as a fncon of he rblence nec energy and specfc dsspaon rae ; herefore s referred o as he wo-eqaon rblen model. The rblence nec energy and specfc dsspaon rae of sandard model are deermned by followng wo eqaons: The eqaon: S Y G (6) The eqaon: S Y G (7) In hese eqaons, G represens he generaon of rblence nec energy de o mean velocy gradens. G represens he generaon of ω. and represen he effecve dffsvy of and ω, respecvely. Y and Y represen he dsspaon of and ω de o rblence. S and S are ser-defned sorce erms. A low Reynolds nmber correcon facor conrols he nflence on he overall srcre of he flow feld, dependng pon local condons, and s gven as-

4 rr Re R 0 1 Re R (8) Where, Re, R 8, 0, 0.072, 1. 4 Closre Coeffcen for he Transonal Model are- 1, 0.52, , 0.09, 0.072, 8and R 8 Dealed descrpons of seen harmonc coeffcens are shown n TABLE 1. Here, Reynolds nmber vares from 90 o Snce cardac plse cycle s 0.82sec, s plsaon, fond from he calclaon rad 0.82 sec R. Fg. 2(a) and 2(b) show oscllaory physologcal waveform and parabolc nle velocy profle respecvely. In fgre 2(a), a, b, and c represen he posons of early sysole (0.0615sec), pea sysole ( sec), and dasole (0.533 sec) respecvely. Fgre 2: (a) oscllaory snsodal waveform, (b) parabolc nle velocy profle, and (c) velocy dsrbon n 75% senoc arery from dfferen mesh szes. 3.2 Grd ndependence chec An eensve es s carred o wh dfferen szes of mesh sch as mesh 1 ( elemen), mesh2 ( elemen) and mesh3 ( elemen) respecvely. Fg 2(c) shows he velocy dsrbons for 75% senosed carod arery wh menoned mesh szes. In all cases, he velocy dsrbons are same. I mples ha he solon s grd ndependence. 3.3 Valdaon Before sarng of presen nvesgaon he nmercal smlaon s needed o be valdaed. Valdaon of he presen nmercal compaon s done by plong he seady velocy profle a downsream from he senoss hroa and comparng wh he velocy profle of Varghese and Franel [22]. For hs case, a parabolc velocy profle s assmed as nle bondary condon. The mean nle velocy corresponds o Reynolds nmber 500 and he flow s assmed o be seady. The resls are shown n he Fgre 3, where a good agreemen can be fond wh Varghese and Franel [22] Low Re -ω (presen) Varghese and Franel Fgre 3: Valdaon wh he velocy profle of Varghese and Franel [22]. DOI: Page

5 3.4 Nmercal scheme The nmercal smlaons are performed by well-nown sofware ANSYS Flen A pressre based algorhm s chosen as he solver ype. Ths solver s generally seleced for an ncompressble fld. As here s no hea ransfer n he blood flow process, energy eqaon s no solved. Snce rblen s epeced n 75% senoc carod arery a pos senoc and bfrcaon regon, a low Reynolds nmber rblen model s sed hrogho he wor. In solon mehods, he SIMPLE algorhm s seleced for pressre-velocy coplng. Frs Order Upwnd scheme s employed as a nmercal scheme for dscrezaon of he momenm eqaon. The me sep s se o sec wh 2000 nmber of oal me seps. Mamm 10 eraons are performed per each me sep. IV. Resls And Dscsson Of Snsodal Blood Flow The compaonal resls of snsodal blood flow are condced o sdy he nflence of senoss on he flow behavor. The flow parameers le velocy, pressre, WSS are observed from ransverse and longdnal conors a specfc nsans of plse cycle for comparng he flow varaon. The dscsson s caegorzed wh he observaons of flow varaon sarng from early sysole, pea sysole and dasole, respecvely Wall Shear Sress(WSS) dsrbon WSS s an mporan facor o deermne he effec of severy of areral senoss. I depends on he vscosy of he fld and velocy graden and, s defned as d, where s he vscosy and d s he velocy dr dr graden. Fgres, 4(a), 4(b), 4(c) show he conors of WSS n normal carod arery and 75% (by area) senosed carod arery and Fgs.5(a), 5(b), 5(c) show he comparsons of WSS n boh senoc arery a earlysysole, pea sysole and dasole respecvely. In boh cases, WSS s nformly dsrbed a C regons. The resls reveal ha he flow paern changes abrply n he regon of senoss hroa and bfrcaon onwards n he downsream whereas no sgnfcan change s observed n he psream and s sgnfcanly observed ha he changes of WSS s de o senoss and bfrcaon. Srprsngly a sharp rase of WSS and mmedaely declne are observed n he senoc hroa. From he comparsons, or fndng s ha mamm WSS s observed a pea sysole de o senoss. Afer bfrcaon WSS becomes hgh de o snoss a all cases and hen WSS becomes more low a pea sysole and dasole. Fgre 4:Comparson beween Carod arery and Senosed carod arery WSS conor a (a) early sysole (b) pea sysole and (c) dasole. DOI: Page

6 Fgre 5:Comparson beween Carod arery and Senosed carod arery WSS dsrbon a (a) early sysole (b) pea sysole (c) dasole. 4.2.Velocy dsrbon In fgre 6, nle velocy n boh geomeres begns wh same magnde for all case. Veloces a he hroa and pos senoc regon for all cases of S are relavely hgher han ha of accordng arery. Recrclaon regon s observed n S a pos senoc regon for all cases. In he oher hand, velocy s lamnar n C regon for all cases. The hghes velocy s observed n S a pea sysole and he lowes velocy s observed n a dasole. Fgre 6: Comparson beween Carod arery and Senosed carod arery a dfferen plse cycle. Comparson beween Carod arery and Senosed carod arery velocy conors a cross-secon (a) early sysole (b) pea sysole and (c) dasole a dsance from velocy-nle () 0.0mm (VI),() 0.072mm, () 0.08mm (mddle of senoss), (l) 0.088mm (bfrcaon pon) and (m) 0.103mm are shown n fgre In all cases we see no change n velocy dsrbon before senoc regon. Fgre 4.4 shows de o bfrcaon velocy decreases n all cases and velocy s hgher n nner sde han he oer sde of he arery. The presence of senoss ncreases he phenomena,so he possbly of formaon plaqe n oer wall s hgh. Carod Arery Senosed Carod Arery DOI: Page

7 pressre(pa) pressre(pa) pressre(pa) Fgre 7:Comparson beween Carod arery and Senosed carod arery velocy conor a cross-secon (a) early sysole (b) pea sysole and (c) dasole a dsance from velocy-nle () 0.0mm (VI),() 0.072mm, () 0.08mm (mddle of senoss), (l) 0.088mm (bfrcaon pon) and (m) 0.103mm Pressre dsrbons Pressre conor n boh areres a early-sysole, pea-sysole, and dasole are shown n Fgre 8 and cross-seconal pressre dsrbon along dameer are shown n fgre 9. Fgre 8(a) shows he conor of pressre n and S a early-sysole. In he pressre s nform and s varyng wh a small varaon hrogho he geomery. Hgh pressre a pre-senoc regon and sdden declne of he pressre a he hroa are characerzed n S. Agan, fgre able 1 show ha he pressre a mddle pon of cross-secon (), () (l), (m) and (n) for and S are as follows: () mmhg () mmhg (l) mmhg (m) mmhg (n) mmhg Early sysole S Pea sysole S Dasole S Table 1: Pressre a mddle pon of dfferen cross-secon The resls clearly noced ha S can case faal deah by blasng he pre senoc areral regon of he paen of hgh blood pressre. From fgre 8(c), we observe ha, pressre s nform and changng wh a small varaon before he hroa. From able 1, s observed ha, he pressre a before hroa and hroa for are mmhg and mmhg and for S are mmhg and mmhg respecvely. I s well nown ha pressre dfference may provde some nd of era force ha may drve blood hrogh he senoss wh hgh velocy. Snce The pressre dfference n S s hgher han ha of for all case. So a every case S creaes more force han a her pre-senoc regon. Fgre 8:Comparson beween Carod arery and Senosed carod arery pressre conor a (a) early sysole (b) pea sysole and (c) dasole. ( ) S S S rr rr r R DOI: Page

8 pressre(pa) pressre(pa) pressre(pa) pressre(pa) pressre(pa) pressre(pa) pressre(pa) pressre(pa) pressre(pa) pressre(pa) pressre(pa) pressre(pa) ( ) (l ) ( m ) ( n ) rr S S rr S rr S rr (a) (b) (c) Fgre 9: Comparson beween Carod arery and Senosed carod arery pressre dsrbon a cross-secon along dameer (a) early sysole (b) pea sysole and (c) dasole a dsance from velocy-nle () 0.072mm, () 0.08mm (mddle of senoss), (l) 0.088mm (bfrcaon pon), (m) 0.103mm(I) and (n) 0.103(E). 4.4Sreamlne conors Varaons of sreamlnes conors have been nvesgaed for and S o sdy he flow paern of blood hrogh arery. Sreamlne dsrbon for he models a early-sysole, pea-sysole, and dasole s shown n fgre10. A early sysole here s no vore n and sreamlnes rn smoohly hrogho he geomeres. B n S a very small vore s noced afer he hroa and bfrcaon regon, ecep hs mos of he sreamlnes rn smoohly hrogho he geomery. A pea all areres creae vore. Two large vorees a he pos senoc regon near he hroa are seen n S and a lle rblence s seen a I n. B here are dfferences n vore sze and shape. A dasole all areres creae sgnfcan sze of vore. Lengh of vore s ncreased n S. Sgnfcan sreamlne dfferences are observed among he areres. S creae large dfference n sreamlne from. S rr S rr S rr S rr S rr S rr rr S S rr DOI: Page

9 Fgre 10: Comparson beween Carod arery and Senosed carod arery velocy sreamlne a (a) early sysole (b) pea sysole and (c) dasole. V. Conclson Nmercal smlaons are performed on plsale blood flow hrogh carod arery and 75% severy senosed carod arery. Snsodal velocy profles are mposed a nle for he smlaons. Wall shear sress dsrbon shows ha as severy ncreases, shearng of flow also ncreases for all cases. And mamm sress s eered n hroa regon a pea sysole. Pressre dsrbon demonsrae ha a every case senosed carod arery creaes more force han normal carod arery a her pre-senoc regon and 75% areral senoss can case faal deah by blasng he pre senoc areral regon of he paen of hgh blood pressre. De o senoc consrcon, he flow behavor changes abrply n he downsream of he senoss. Drng pea sysole he flow ncreases a he hroa regon formng a e and laer dsrps sddenly formng eddes de o pressre drop n downsream sde. Analyss of sreamlne shows ha a recrclaon regon s vsble a he pos senoc regon a 75% senoc carod arery for all cases and recrclaon regon ncreases wh ncreasng mnmm velocy a he nle flow. Fnally, sgnfcan flow changes are noced n senosed carod arery. So senoss s very harmfl condon o any paen. References [1]. Lopez AD, Mahers CD, Ezza M, Jamson DT and Mrray CJ., Global and regonal brden of dsease and rs facors, 2001:sysemac analyss of poplaon healh daa, Lance 367: , (2006). [2]. Lloyd-Jones D, Adams RJ, Brown TM, Carnehon M, Da S, e al., Hear dsease and sroe sascs-2010 pdae: a repor from he Amercan Hear Assocaon Crclaon 121:e46-e215, (2010). [3]. Haney GJ, Jamroz K, Broadhrs RJ, Forbes S and Anderson CS., Long-erm dsably afer frs-ever sroe and relaed prognosc facors n he Perh Commny Sroe Sdy, Sroe 33: , (2002). [4]. Mohan KM, Crchon SL, Greve AP, Rdd AG, Wolfe CD, e al., Freqency and predcors for he rs of sroe recrrence p o 10 years afer sroe, he Soh London Sroe Regser. J Nerol Nerosrg, Psychary 80: (2009). [5]. Paala J, Haapanem E, Meso AJ, Meso TM, Aro V, e al., Recrren schemc evens n yong adls afer frs-ever schemc sroe., Ann Nerol 68: , (2010). [6]. O'Donnell MJ, Xaver D, L L, Zhang H, Chn SL, e al., Rsfacors for schaemc and nracerebral haemorrhagc sroe n 22 conres(he INTERSTROKE sdy): a case-conrol sdy., Lance 376: , (2010). [7]. S. Khader and B. Shenoy., Effec of ncreased severy n paen specfc senoss,world ornal of Modellng and Smlaon, Vol. 7 No. 2,pp ,(2011). [8]. M.X. L, J.J. Beech-Brand, L.R. John,P.R.Hosns and W.J. Easson., Nmercal analyss of plsale blood flow and vessel wall mechancs n dfferen degrees of senoss,j. Bomech, 40, ,(2007). [9]. C, Cha and G. Shread., Changes n flow and wall sresses hrogh areral consrcon offse from he cenre of he vessel, The Anzam Jornal, 50: C744-C759,(2009). [10]. Ags, A.D., B. Arff, S. A. G. McG. Thom, X.Y. X and A. D. Hghes., Analyss of comple flow and relaonshp beween blood pressre, wall shear sress and nma-meda hcness n he hmen carod arery,am. J. physol. Hear Cr. Physol.293: ,(2007). [11]. Bale-Glcman, J.K. Selby, D. Saloner and O. Savas., Epermenal flow sdes n eac-replca phanoms of aheroscleroc carod bfrcaonsnder seady np condons, J. Bomech. Eng. 125(1):38-48,(2003). DOI: Page

10 [12]. Kohler, U. I. Marshall, M.B. Roberson, Q. Long and X. Y. X., MRI measremen of wall shear sress vecors n bfrcaon models and comparson wh CFD predcons, J. magn. Reson. Imagng 14(5): ,(2001). [13]. Fscher, P. F., F. Loh, S. E. Lee, S.W. Lee, D. Smh and H. Bassony., : Smlaon of hgh Reynolds nmber vasclar flows,comp. Mehods Appl. Mech. Eng. 196: ,(2007). [14]. Senman, D.A., T.L. Poeppng, M. Tambasco, R. N. Rann and D. W. Holdsworh., Flow paerns a he senosed carod bfrcaon: effec of concenrc verss eccenrc senoss, Ann. Bomed. Eng. 28(4): ,(2000). [15]. Srod, J. S. A. Berger and D. Saloner., Nmercal analyss of flow hrogh a severely senoc carod arery bfrcaon, J. Bomech. Eng.124:9-20,(2002). [16]. S. Lee and S. Lee., Drec nmercal smlaon of ransonal flow n a senose carod bfrcaon, Jornal of Bomechancs,41: ,(2008). [17]. V. Yong, A. Paerson M. Graves, Z-Y LI, V. Tavan, T. Tang, and J. H. Gllard, The Brsh Jornal of Radology, S39-S45 (2009). [18]. J. Pno, K.L. Bessa, D.F. Legendre, and R.H. Prado, ABCM Symposm n Boengneerng, vol. 01, [19]. S.A. Ahmed, and D.P. Gddens, Jornal of Bomechancs, 17, (1984). [20]. G. Lorenzn, and E. Casalena, Jornal of Bomechancs, 41, (2008). [21]. K Mamn, M.N.Aher, M.S.H.Mollah,M.A.N.Sheh, and M.Al, The 6h BSME Inernaonal Conference on Thermal Engneerng, Dhaa(2014), eded by A.K.M. Sadrl Islam, M. Rhl Amn, Mohammad Al and Tosha Seogch (Proceda Engneerng, 2015), vol. 150, pp [22]. S.S. Varghese and S.H. Franel, Jornal of Bomechancs, 125, (2003). Mohammad Mar Rahman Comparave Sdy of Neonan Snsodal Blood Flow hrogh Normal and Senosed Carod Arery. IOSR Jornal of Mahemacs (IOSR-JM), vol. 13, no. 6, 2017, pp DOI: Page

OPTIMIZATION OF A NONCONVENTIONAL ENGINE EVAPORATOR

OPTIMIZATION OF A NONCONVENTIONAL ENGINE EVAPORATOR Jornal of KONES Powerran and Transpor, Vol. 17, No. 010 OPTIMIZATION OF A NONCONVENTIONAL ENGINE EVAPORATOR Andre Kovalí, Eml Toporcer Unversy of Žlna, Facly of Mechancal Engneerng Deparmen of Aomove Technology

More information

Calculation of the Resistance of a Ship Mathematical Formulation. Calculation of the Resistance of a Ship Mathematical Formulation

Calculation of the Resistance of a Ship Mathematical Formulation. Calculation of the Resistance of a Ship Mathematical Formulation Ressance s obaned from he sm of he frcon and pressre ressance arables o deermne: - eloc ecor, (3) = (,, ) = (,, ) - Pressre, p () ( - Dens, ρ, s defned b he eqaon of sae Ressance and Proplson Lecre 0 4

More information

Prediction of Wing Downwash Using CFD

Prediction of Wing Downwash Using CFD Predcon of Wng Downwash Usng CFD Mohammed MAHDI* *Correspondng ahor Aeronacal Research Cener-Sdan P.O. Bo 334 momahad7@homal.com DOI:.3/66-8.5.7.. 3 rd Inernaonal Worshop on Nmercal Modellng n Aerospace

More information

by Lauren DeDieu Advisor: George Chen

by Lauren DeDieu Advisor: George Chen b Laren DeDe Advsor: George Chen Are one of he mos powerfl mehods o nmercall solve me dependen paral dfferenal eqaons PDE wh some knd of snglar shock waves & blow-p problems. Fed nmber of mesh pons Moves

More information

Numerical simulation of flow reattachment length in a stilling basin with a step-down floor

Numerical simulation of flow reattachment length in a stilling basin with a step-down floor 5 h Inernaonal Symposm on Hydralc Srcres Brsbane, Asrala, 5-7 Jne 04 Hydralc Srcres and Socey: Engneerng hallenges and Eremes ISBN 97874756 - DOI: 0.464/ql.04.3 Nmercal smlaon of flow reaachmen lengh n

More information

Turbulence Modelling (CFD course)

Turbulence Modelling (CFD course) Trblence Modellng (CFD corse) Sławomr Kbac slawomr.bac@mel.pw.ed.pl 14.11.016 Copyrgh 016, Sławomr Kbac Trblence Modellng Sławomr Kbac Conens 1. Reynolds-averaged Naver-Soes eqaons... 3. Closre of he modelled

More information

Improvement of Two-Equation Turbulence Model with Anisotropic Eddy-Viscosity for Hybrid Rocket Research

Improvement of Two-Equation Turbulence Model with Anisotropic Eddy-Viscosity for Hybrid Rocket Research evenh Inernaonal onference on ompaonal Fld Dynamcs (IFD7), Bg Island, awa, Jly 9-, IFD7-9 Improvemen of Two-Eqaon Trblence Model wh Ansoropc Eddy-Vscosy for ybrd oce esearch M. Mro * and T. hmada ** orrespondng

More information

Outline. Review Solution Approaches. Review Basic Equations. Nature of Turbulence. Review Fluent Exercise. Turbulence Models

Outline. Review Solution Approaches. Review Basic Equations. Nature of Turbulence. Review Fluent Exercise. Turbulence Models Trblence Models Larry areo Mechancal Engneerng 69 ompaonal Fld Dynamcs Febrary, Olne Revew las lecre Nare of rblence Reynolds-average Naver-Soes (RNS) Mng lengh heory Models sng one dfferenal eqaon Two-eqaon

More information

CONSISTENT EARTHQUAKE ACCELERATION AND DISPLACEMENT RECORDS

CONSISTENT EARTHQUAKE ACCELERATION AND DISPLACEMENT RECORDS APPENDX J CONSSTENT EARTHQUAKE ACCEERATON AND DSPACEMENT RECORDS Earhqake Acceleraons can be Measred. However, Srcres are Sbjeced o Earhqake Dsplacemens J. NTRODUCTON { XE "Acceleraon Records" }A he presen

More information

VI. Computational Fluid Dynamics 1. Examples of numerical simulation

VI. Computational Fluid Dynamics 1. Examples of numerical simulation VI. Comaonal Fld Dnamcs 1. Eamles of nmercal smlaon Eermenal Fas Breeder Reacor, JOYO, wh rmar of coolan sodm. Uer nner srcre Uer lenm Flow aern and emerare feld n reacor essel n flow coas down Core Hh

More information

Numerical Simulation on Wind Flow over Step-shaped Cliff Topography with Rough Surface

Numerical Simulation on Wind Flow over Step-shaped Cliff Topography with Rough Surface In. J. Envron. Res., 7(1):173-186, Wner 013 ISSN: 1735-6865 Nmercal Smlaon on Wnd Flow over Sep-shaped Clff Topography wh Rogh Srface Yassn, M.F. 1,* and Al-Harb, M. 1 1 Deparmen of Envronmenal Technology

More information

, t 1. Transitions - this one was easy, but in general the hardest part is choosing the which variables are state and control variables

, t 1. Transitions - this one was easy, but in general the hardest part is choosing the which variables are state and control variables Opmal Conrol Why Use I - verss calcls of varaons, opmal conrol More generaly More convenen wh consrans (e.g., can p consrans on he dervaves More nsghs no problem (a leas more apparen han hrogh calcls of

More information

Solution of a diffusion problem in a non-homogeneous flow and diffusion field by the integral representation method (IRM)

Solution of a diffusion problem in a non-homogeneous flow and diffusion field by the integral representation method (IRM) Appled and ompaonal Mahemacs 4; 3: 5-6 Pblshed onlne Febrary 4 hp://www.scencepblshnggrop.com//acm do:.648/.acm.43.3 olon of a dffson problem n a non-homogeneos flow and dffson feld by he negral represenaon

More information

Dynamic Model of the Axially Moving Viscoelastic Belt System with Tensioner Pulley Yanqi Liu1, a, Hongyu Wang2, b, Dongxing Cao3, c, Xiaoling Gai1, d

Dynamic Model of the Axially Moving Viscoelastic Belt System with Tensioner Pulley Yanqi Liu1, a, Hongyu Wang2, b, Dongxing Cao3, c, Xiaoling Gai1, d Inernaonal Indsral Informacs and Comper Engneerng Conference (IIICEC 5) Dynamc Model of he Aally Movng Vscoelasc Bel Sysem wh Tensoner Plley Yanq L, a, Hongy Wang, b, Dongng Cao, c, Xaolng Ga, d Bejng

More information

Different kind of oscillation

Different kind of oscillation PhO 98 Theorecal Qeson.Elecrcy Problem (8 pons) Deren knd o oscllaon e s consder he elecrc crc n he gre, or whch mh, mh, nf, nf and kω. The swch K beng closed he crc s copled wh a sorce o alernang crren.

More information

Chapter 1 Introduction of boundary layer phenomena

Chapter 1 Introduction of boundary layer phenomena Chaper 1 Inrodcon of bondary layer phenomena T-S Le Jan. 13, 018 Man Topcs Hsory of Fld Mechancs Developmen Idea of Bondary Layer Bondary Layer Eqaons 1 Fld Mechancs Developmen Hsory Ideal fld: Invscd

More information

Solving Parabolic Partial Delay Differential. Equations Using The Explicit Method And Higher. Order Differences

Solving Parabolic Partial Delay Differential. Equations Using The Explicit Method And Higher. Order Differences Jornal of Kfa for Maemacs and Compe Vol. No.7 Dec pp 77-5 Solvng Parabolc Paral Delay Dfferenal Eqaons Usng e Eplc Meod And Hger Order Dfferences Asss. Prof. Amal Kalaf Haydar Kfa Unversy College of Edcaon

More information

The Elastic Wave Equation. The elastic wave equation

The Elastic Wave Equation. The elastic wave equation The Elasc Wave Eqaon Elasc waves n nfne homogeneos soropc meda Nmercal smlaons for smple sorces Plane wave propagaon n nfne meda Freqency, wavenmber, wavelengh Condons a maeral dsconnes nell s Law Reflecon

More information

Cartesian tensors. Order (rank) Scalar. Vector. 3x3 matrix

Cartesian tensors. Order (rank) Scalar. Vector. 3x3 matrix Caresan ensors Order (rank) 0 1 3 a b c d k Scalar ecor 33 mar Caresan ensors Kronecker dela δ = 1 f = 0 f Le- Ca epslon ε k = 1 f,, k are cclc 1 f,, k are ancclc 0 oherse Smmaon conenon (o eqal ncces

More information

Observer Design for Nonlinear Systems using Linear Approximations

Observer Design for Nonlinear Systems using Linear Approximations Observer Desgn for Nonlnear Ssems sng Lnear Appromaons C. Navarro Hernandez, S.P. Banks and M. Aldeen Deparmen of Aomac Conrol and Ssems Engneerng, Unvers of Sheffeld, Mappn Sree, Sheffeld S 3JD. e-mal:

More information

ON THE ACCURACY OF NUMERICAL PREDICTION IN TRANSONIC-SUPERSONIC FLOW ARROUND MISSILES

ON THE ACCURACY OF NUMERICAL PREDICTION IN TRANSONIC-SUPERSONIC FLOW ARROUND MISSILES U.P.B. Sc. Bll., Seres D, Vol. 7, Iss. 3, ISSN 454-358 ON THE ACCURACY OF NUMERICAL PREDICTION IN TRANSONIC-SUPERSONIC FLOW ARROUND MISSILES Crsna MIHAILESCU, Teodor Vorel CHELARU, Seran DANAILA 3, Cornel

More information

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5 TPG460 Reservor Smulaon 08 page of 5 DISCRETIZATIO OF THE FOW EQUATIOS As we already have seen, fne dfference appromaons of he paral dervaves appearng n he flow equaons may be obaned from Taylor seres

More information

Research Article Cubic B-spline for the Numerical Solution of Parabolic Integro-differential Equation with a Weakly Singular Kernel

Research Article Cubic B-spline for the Numerical Solution of Parabolic Integro-differential Equation with a Weakly Singular Kernel Researc Jornal of Appled Scences, Engneerng and Tecnology 7(): 65-7, 4 DOI:.96/afs.7.5 ISS: 4-7459; e-iss: 4-7467 4 Mawell Scenfc Pblcaon Corp. Sbmed: Jne 8, Acceped: Jly 9, Pblsed: Marc 5, 4 Researc Arcle

More information

Solution in semi infinite diffusion couples (error function analysis)

Solution in semi infinite diffusion couples (error function analysis) Soluon n sem nfne dffuson couples (error funcon analyss) Le us consder now he sem nfne dffuson couple of wo blocks wh concenraon of and I means ha, n a A- bnary sysem, s bondng beween wo blocks made of

More information

Stochastic Programming handling CVAR in objective and constraint

Stochastic Programming handling CVAR in objective and constraint Sochasc Programmng handlng CVAR n obecve and consran Leondas Sakalaskas VU Inse of Mahemacs and Informacs Lhana ICSP XIII Jly 8-2 23 Bergamo Ialy Olne Inrodcon Lagrangan & KKT condons Mone-Carlo samplng

More information

Chapter 6 DETECTION AND ESTIMATION: Model of digital communication system. Fundamental issues in digital communications are

Chapter 6 DETECTION AND ESTIMATION: Model of digital communication system. Fundamental issues in digital communications are Chaper 6 DCIO AD IMAIO: Fndaenal sses n dgal concaons are. Deecon and. saon Deecon heory: I deals wh he desgn and evalaon of decson ang processor ha observes he receved sgnal and gesses whch parclar sybol

More information

Mathematical Analysis of non-newtonian Fluid Flow through Multiple Stenotic Artery in the Presence of Catheter-A Pulsatile Flow.

Mathematical Analysis of non-newtonian Fluid Flow through Multiple Stenotic Artery in the Presence of Catheter-A Pulsatile Flow. ISSN 1749-3889 (prn), 1749-3897 (onlne) Inernaonal Journal of Nonlnear Scence Vol.13(1) No.1,pp.15-7 Mahemacal Analyss of non-newonan Flud Flow hrough Mulple Senoc Arery n he Presence of Caheer-A Pulsale

More information

Variational method to the second-order impulsive partial differential equations with inconstant coefficients (I)

Variational method to the second-order impulsive partial differential equations with inconstant coefficients (I) Avalable onlne a www.scencedrec.com Proceda Engneerng 6 ( 5 4 Inernaonal Worksho on Aomoble, Power and Energy Engneerng Varaonal mehod o he second-order mlsve aral dfferenal eqaons wh nconsan coeffcens

More information

Comparison between two solar tower receivers of different geometry

Comparison between two solar tower receivers of different geometry Reve des Energes Renovelables Vol. 20 N 4 (2017) 713-720 Comparson beween wo solar ower recevers of dfferen geomery M. Hazmone 1.2 *, B. Aor 2, M.M. Hada 1 and A. Male 1 1 Cenre de Développemen des Energes

More information

Robustness Experiments with Two Variance Components

Robustness Experiments with Two Variance Components Naonal Insue of Sandards and Technology (NIST) Informaon Technology Laboraory (ITL) Sascal Engneerng Dvson (SED) Robusness Expermens wh Two Varance Componens by Ana Ivelsse Avlés avles@ns.gov Conference

More information

CFD MODELING FOR HELIUM RELEASES IN A PRIVATE GARAGE WITHOUT FORCED VENTILATION

CFD MODELING FOR HELIUM RELEASES IN A PRIVATE GARAGE WITHOUT FORCED VENTILATION CFD MODELING FOR HELIUM RELEASES IN A PRIVATE GARAGE WITHOUT FORCED VENTILATION Papankolao, E.A. 1 and Venesanos, A.G. 1 1 Envronmenal Research Laboraory, NCSR Demokros, Agha Paraskev, Aks, 15310, Greece,

More information

Numerical Simulation of the Dispersion of a Plume of Exhaust Gases from Diesel and Petrol Engine Vehicles

Numerical Simulation of the Dispersion of a Plume of Exhaust Gases from Diesel and Petrol Engine Vehicles World Academy of Scence, Engneerng and Technology 67 01 Numercal Smulaon of he Dsperson of a Plume of Exhaus Gases from Desel and Perol Engne Vehcles H. ZAHLOUL, and M. MERIEM-BENZIANE Absrac The obecve

More information

Block 5 Transport of solutes in rivers

Block 5 Transport of solutes in rivers Nmeral Hydrals Blok 5 Transpor of soles n rvers Marks Holzner Conens of he orse Blok 1 The eqaons Blok Compaon of pressre srges Blok 3 Open hannel flow flow n rvers Blok 4 Nmeral solon of open hannel flow

More information

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria IOSR Journal of Mahemacs (IOSR-JM e-issn: 78-578, p-issn: 9-765X. Volume 0, Issue 4 Ver. IV (Jul-Aug. 04, PP 40-44 Mulple SolonSoluons for a (+-dmensonalhroa-sasuma shallow waer wave equaon UsngPanlevé-Bӓclund

More information

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS R&RATA # Vol.) 8, March FURTHER AALYSIS OF COFIDECE ITERVALS FOR LARGE CLIET/SERVER COMPUTER ETWORKS Vyacheslav Abramov School of Mahemacal Scences, Monash Unversy, Buldng 8, Level 4, Clayon Campus, Wellngon

More information

Variants of Pegasos. December 11, 2009

Variants of Pegasos. December 11, 2009 Inroducon Varans of Pegasos SooWoong Ryu bshboy@sanford.edu December, 009 Youngsoo Cho yc344@sanford.edu Developng a new SVM algorhm s ongong research opc. Among many exng SVM algorhms, we wll focus on

More information

XIII International PhD Workshop OWD 2011, October Three Phase DC/DC Boost Converter With High Energy Efficiency

XIII International PhD Workshop OWD 2011, October Three Phase DC/DC Boost Converter With High Energy Efficiency X nernaonal Ph Workshop OW, Ocober Three Phase C/C Boos Converer Wh Hgh Energy Effcency Ján Perdľak, Techncal nversy of Košce Absrac Ths paper presens a novel opology of mlphase boos converer wh hgh energy

More information

10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current :

10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current : . A. IUITS Synopss : GOWTH OF UNT IN IUIT : d. When swch S s closed a =; = d. A me, curren = e 3. The consan / has dmensons of me and s called he nducve me consan ( τ ) of he crcu. 4. = τ; =.63, n one

More information

ESTIMATION OF HYDRAULIC JUMP LOCATION USING NUMERICAL SIMULATION

ESTIMATION OF HYDRAULIC JUMP LOCATION USING NUMERICAL SIMULATION Twel Inernaonal Waer Tecnology Conerence, IWTC 8, Aleandra, gyp STIMATION OF HYDRAULIC JUMP LOCATION USING NUMRICAL SIMULATION M. T. Samaa,, M. Hasem and H. M. Karr, Cvl ngneerng Deparmen, Facly o ngneerng,

More information

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4 CS434a/54a: Paern Recognon Prof. Olga Veksler Lecure 4 Oulne Normal Random Varable Properes Dscrmnan funcons Why Normal Random Varables? Analycally racable Works well when observaon comes form a corruped

More information

Method of Characteristics for Pure Advection By Gilberto E. Urroz, September 2004

Method of Characteristics for Pure Advection By Gilberto E. Urroz, September 2004 Mehod of Charaerss for Pre Adveon By Glbero E Urroz Sepember 004 Noe: The followng noes are based on lass noes for he lass COMPUTATIONAL HYDAULICS as agh by Dr Forres Holly n he Sprng Semeser 985 a he

More information

TSS = SST + SSE An orthogonal partition of the total SS

TSS = SST + SSE An orthogonal partition of the total SS ANOVA: Topc 4. Orhogonal conrass [ST&D p. 183] H 0 : µ 1 = µ =... = µ H 1 : The mean of a leas one reamen group s dfferen To es hs hypohess, a basc ANOVA allocaes he varaon among reamen means (SST) equally

More information

CFD Analysis of Aerodynamic Drag Effects on Vacuum Tube Trains

CFD Analysis of Aerodynamic Drag Effects on Vacuum Tube Trains Jornal of Appled Fld Mechancs, ol. 1, No. 1, pp. 303-309, 019. Aalable onlne a.afmonlne.ne, ISSN 1735-357, EISSN 1735-3645. DOI: 10.95/afm.75.53.9091 CFD Analss of Aerodnamc Drag Effecs on acm Tbe Trans

More information

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim Korean J. Mah. 19 (2011), No. 3, pp. 263 272 GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS Youngwoo Ahn and Kae Km Absrac. In he paper [1], an explc correspondence beween ceran

More information

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!) i+1,q - [(! ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL The frs hng o es n wo-way ANOVA: Is here neracon? "No neracon" means: The man effecs model would f. Ths n urn means: In he neracon plo (wh A on he horzonal

More information

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany Herarchcal Markov Normal Mxure models wh Applcaons o Fnancal Asse Reurns Appendx: Proofs of Theorems and Condonal Poseror Dsrbuons John Geweke a and Gann Amsano b a Deparmens of Economcs and Sascs, Unversy

More information

Modelling of Diffusion Process in Porous Bricks

Modelling of Diffusion Process in Porous Bricks hp://www.aras.org/aras/ornals/mcm Modellng of Dffson Process n Poros Brcs KRNIATI ORNAM, MASYKR KIMSAN, LA ODE NGKOIMANI 3, EDI CAHYONO 4 Deparmen of Archecre, Hal Oleo nversy, Jl.H.E.A Moodomp Kamps Ha

More information

Wronskian Determinant Solutions for the (3 + 1)-Dimensional Boiti-Leon-Manna-Pempinelli Equation

Wronskian Determinant Solutions for the (3 + 1)-Dimensional Boiti-Leon-Manna-Pempinelli Equation Jornal of Appled Mahemacs and Physcs 0 8-4 Pblshed Onlne ovember 0 (hp://www.scrp.org/jornal/jamp) hp://d.do.org/0.46/jamp.0.5004 Wronskan Deermnan Solons for he ( + )-Dmensonal Bo-Leon-Manna-Pempnell

More information

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study) Inernaonal Mahemacal Forum, Vol. 8, 3, no., 7 - HIKARI Ld, www.m-hkar.com hp://dx.do.org/.988/mf.3.3488 New M-Esmaor Objecve Funcon n Smulaneous Equaons Model (A Comparave Sudy) Ahmed H. Youssef Professor

More information

Study on the Unexpected Wave of the Experiment of Ultrasonic Guided Wave in Square Steel Bar based on 2D Equivalent Simulation

Study on the Unexpected Wave of the Experiment of Ultrasonic Guided Wave in Square Steel Bar based on 2D Equivalent Simulation Sdy on he Unepeced Wave of he Epermen of Ulrasonc Gded Wave n Sqare Seel Bar based on D Eqvalen Smlaon Le Zhang,, Yan Yang,*, Xaoyan We, Wenqng Yao School of aomaon and nformaon engneerng, X'an Unversy

More information

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations. Soluons o Ordnary Derenal Equaons An ordnary derenal equaon has only one ndependen varable. A sysem o ordnary derenal equaons consss o several derenal equaons each wh he same ndependen varable. An eample

More information

ELASTIC MODULUS ESTIMATION OF CHOPPED CARBON FIBER TAPE REINFORCED THERMOPLASTICS USING THE MONTE CARLO SIMULATION

ELASTIC MODULUS ESTIMATION OF CHOPPED CARBON FIBER TAPE REINFORCED THERMOPLASTICS USING THE MONTE CARLO SIMULATION THE 19 TH INTERNATIONAL ONFERENE ON OMPOSITE MATERIALS ELASTI MODULUS ESTIMATION OF HOPPED ARBON FIBER TAPE REINFORED THERMOPLASTIS USING THE MONTE ARLO SIMULATION Y. Sao 1*, J. Takahash 1, T. Masuo 1,

More information

Modelling of test case particle-laden jet with NEPTUNE_CFD

Modelling of test case particle-laden jet with NEPTUNE_CFD Modelln of es case arcle-laden e wh NEPTNE_CFD H. D. Le 1 J-M. Lacome 1 A. Vnes 1 B. Debray 1 B. Trcho 1 P. Fede 3 E. Clmen 3 1 INERIS Parc echnoloe ALATA B.P. FR-60550 Vernel-en-Halae nversé de Tolose

More information

Polymerization Technology Laboratory Course

Polymerization Technology Laboratory Course Prakkum Polymer Scence/Polymersaonsechnk Versuch Resdence Tme Dsrbuon Polymerzaon Technology Laboraory Course Resdence Tme Dsrbuon of Chemcal Reacors If molecules or elemens of a flud are akng dfferen

More information

FTCS Solution to the Heat Equation

FTCS Solution to the Heat Equation FTCS Soluon o he Hea Equaon ME 448/548 Noes Gerald Reckenwald Porland Sae Unversy Deparmen of Mechancal Engneerng gerry@pdxedu ME 448/548: FTCS Soluon o he Hea Equaon Overvew Use he forward fne d erence

More information

Anisotropic Behaviors and Its Application on Sheet Metal Stamping Processes

Anisotropic Behaviors and Its Application on Sheet Metal Stamping Processes Ansoropc Behavors and Is Applcaon on Shee Meal Sampng Processes Welong Hu ETA-Engneerng Technology Assocaes, Inc. 33 E. Maple oad, Sue 00 Troy, MI 48083 USA 48-79-300 whu@ea.com Jeanne He ETA-Engneerng

More information

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC CH.3. COMPATIBILITY EQUATIONS Connuum Mechancs Course (MMC) - ETSECCPB - UPC Overvew Compably Condons Compably Equaons of a Poenal Vecor Feld Compably Condons for Infnesmal Srans Inegraon of he Infnesmal

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Lnear Response Theory: The connecon beween QFT and expermens 3.1. Basc conceps and deas Q: ow do we measure he conducvy of a meal? A: we frs nroduce a weak elecrc feld E, and hen measure

More information

On One Analytic Method of. Constructing Program Controls

On One Analytic Method of. Constructing Program Controls Appled Mahemacal Scences, Vol. 9, 05, no. 8, 409-407 HIKARI Ld, www.m-hkar.com hp://dx.do.org/0.988/ams.05.54349 On One Analyc Mehod of Consrucng Program Conrols A. N. Kvko, S. V. Chsyakov and Yu. E. Balyna

More information

2.1 Constitutive Theory

2.1 Constitutive Theory Secon.. Consuve Theory.. Consuve Equaons Governng Equaons The equaons governng he behavour of maerals are (n he spaal form) dρ v & ρ + ρdv v = + ρ = Conservaon of Mass (..a) d x σ j dv dvσ + b = ρ v& +

More information

Is it necessary to seasonally adjust business and consumer surveys. Emmanuelle Guidetti

Is it necessary to seasonally adjust business and consumer surveys. Emmanuelle Guidetti Is necessar o seasonall adjs bsness and consmer srves Emmanelle Gde Olne 1 BTS feares 2 Smlaon eercse 3 Seasonal ARIMA modellng 4 Conclsons Jan-85 Jan-87 Jan-89 Jan-91 Jan-93 Jan-95 Jan-97 Jan-99 Jan-01

More information

Tools for Analysis of Accelerated Life and Degradation Test Data

Tools for Analysis of Accelerated Life and Degradation Test Data Acceleraed Sress Tesng and Relably Tools for Analyss of Acceleraed Lfe and Degradaon Tes Daa Presened by: Reuel Smh Unversy of Maryland College Park smhrc@umd.edu Sepember-5-6 Sepember 28-30 206, Pensacola

More information

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy Arcle Inernaonal Journal of Modern Mahemacal Scences, 4, (): - Inernaonal Journal of Modern Mahemacal Scences Journal homepage: www.modernscenfcpress.com/journals/jmms.aspx ISSN: 66-86X Florda, USA Approxmae

More information

Notes on the stability of dynamic systems and the use of Eigen Values.

Notes on the stability of dynamic systems and the use of Eigen Values. Noes on he sabl of dnamc ssems and he use of Egen Values. Source: Macro II course noes, Dr. Davd Bessler s Tme Seres course noes, zarads (999) Ineremporal Macroeconomcs chaper 4 & Techncal ppend, and Hamlon

More information

Background and Motivation: Importance of Pressure Measurements

Background and Motivation: Importance of Pressure Measurements Imornce of Pressre Mesremens: Pressre s rmry concern for mny engneerng lcons e.g. lf nd form drg. Cvon : Pressre s of fndmenl mornce n ndersndng nd modelng cvon. Trblence: Velocy-Pressre-Grden ensor whch

More information

EVALUATION OF FORCE COEFFICIENTS FOR A 2-D ANGLE SECTION USING REALIZABLE k-ε TURBULENCE MODEL

EVALUATION OF FORCE COEFFICIENTS FOR A 2-D ANGLE SECTION USING REALIZABLE k-ε TURBULENCE MODEL The Sevenh Asa-Pacfc Conference on Wnd Engneerng, November 8-, 009, Tape, Tawan EVALUATION OF FORCE COEFFICIENTS FOR A -D ANGLE SECTION USING REALIZABLE k-ε TURBULENCE MODEL S. Chra Ganapah, P. Harkrshna,

More information

Available online at Physics Procedia 32 (2012 )

Available online at   Physics Procedia 32 (2012 ) Avalable onlne a www.scencedrec.com Physcs Proceda 32 (2012 ) 614 622 18h Inernaonal Vacuum Congress, Beng of P. R. Chna, Augus 2010 Numercal research abou he nernal flow of seam-e vacuum pump: evaluaon

More information

Hierarchical Sliding Mode Control for Series Double Inverted Pendulums System

Hierarchical Sliding Mode Control for Series Double Inverted Pendulums System Herarchcal Sldng Mode Conrol for Seres Doble Invered Pendlms Sysem Danwe Qan, Janqang Y, Dongbn Zhao, and Ynxng Hao Laboraory of Complex Sysems and Inellgence Scence Inse of Aomaon, Chnese Academy of Scences

More information

Separated Turbulent Flow Simulations Using a Reynolds Stress Model and Unstructured Meshes

Separated Turbulent Flow Simulations Using a Reynolds Stress Model and Unstructured Meshes AIAA-5-194 AIAA 43rd Aerospace Scences Meeng & Exhb Reno, NV, Janary, 5 Separaed rblen Flo Smlaons Usng a Reynolds Sress Model and Unsrcred Meshes Emre Alpman and Lyle N. Long Deparmen of Aerospace Engneerng

More information

Real-Time Trajectory Generation and Tracking for Cooperative Control Systems

Real-Time Trajectory Generation and Tracking for Cooperative Control Systems Real-Tme Trajecor Generaon and Trackng for Cooperave Conrol Ssems Rchard Mrra Jason Hcke Calforna Inse of Technolog MURI Kckoff Meeng 14 Ma 2001 Olne I. Revew of prevos work n rajecor generaon and rackng

More information

II. Light is a Ray (Geometrical Optics)

II. Light is a Ray (Geometrical Optics) II Lgh s a Ray (Geomercal Opcs) IIB Reflecon and Refracon Hero s Prncple of Leas Dsance Law of Reflecon Hero of Aleandra, who lved n he 2 nd cenury BC, posulaed he followng prncple: Prncple of Leas Dsance:

More information

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes. umercal negraon of he dffuson equaon (I) Fne dfference mehod. Spaal screaon. Inernal nodes. R L V For hermal conducon le s dscree he spaal doman no small fne spans, =,,: Balance of parcles for an nernal

More information

Numerical simulation of a solar chimney power plant in the southern region of Iran

Numerical simulation of a solar chimney power plant in the southern region of Iran Energy Equp. Sys./ Vol. 5/No.4/December 2017/ 431-437 Energy Equpmen and Sysems hp://energyequpsys.u.ac.r www.energyequpsys.com Numercal smulaon of a solar chmney power plan n he souhern regon of Iran

More information

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model BGC1: Survval and even hsory analyss Oslo, March-May 212 Monday May 7h and Tuesday May 8h The addve regresson model Ørnulf Borgan Deparmen of Mahemacs Unversy of Oslo Oulne of program: Recapulaon Counng

More information

Comparison of Differences between Power Means 1

Comparison of Differences between Power Means 1 In. Journal of Mah. Analyss, Vol. 7, 203, no., 5-55 Comparson of Dfferences beween Power Means Chang-An Tan, Guanghua Sh and Fe Zuo College of Mahemacs and Informaon Scence Henan Normal Unversy, 453007,

More information

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction ECOOMICS 35* -- OTE 9 ECO 35* -- OTE 9 F-Tess and Analyss of Varance (AOVA n he Smple Lnear Regresson Model Inroducon The smple lnear regresson model s gven by he followng populaon regresson equaon, or

More information

Department of Economics University of Toronto

Department of Economics University of Toronto Deparmen of Economcs Unversy of Torono ECO408F M.A. Economercs Lecure Noes on Heeroskedascy Heeroskedascy o Ths lecure nvolves lookng a modfcaons we need o make o deal wh he regresson model when some of

More information

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS THE PREICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS INTROUCTION The wo dmensonal paral dfferenal equaons of second order can be used for he smulaon of compeve envronmen n busness The arcle presens he

More information

Cubic Bezier Homotopy Function for Solving Exponential Equations

Cubic Bezier Homotopy Function for Solving Exponential Equations Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 omoopy Funcon for Solvng Eponenal Equaons S. S. Raml *,,. Mohamad Nor,a, N. S. Saharzan,b and M.

More information

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Journal of Appled Mahemacs and Compuaonal Mechancs 3, (), 45-5 HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Sansław Kukla, Urszula Sedlecka Insue of Mahemacs,

More information

A Functional-Link-Based Fuzzy Neural Network for Temperature Control

A Functional-Link-Based Fuzzy Neural Network for Temperature Control Proceedngs of he 7 IEEE Symposm on Fondaons of Compaonal Inellgence (FOCI 7) A Fnconal-Ln-Based Fzzy Neral Neor for emperare Conrol Cheng-Hng Chen *, Chn-eng Ln, Fello, IEEE, and Cheng-Jan Ln, ember, IEEE

More information

Numerical Studies on Lip Shock Flow Behaviors over Backward Facing Sharp Edge Step with Hybrid RANS-LES

Numerical Studies on Lip Shock Flow Behaviors over Backward Facing Sharp Edge Step with Hybrid RANS-LES Numercal Sudes on Lp Shock Flow Behavors over Backward Facng Sharp Edge Sep wh Hybrd RANS-LES Dr. Nrmal Kumar Kund 1 1 Deparmen of Producon Engneerng 1 Veer Surendra Sa Unversy of Technology, Burla, Odsha,

More information

Existence of Periodic Solution for a Non-Autonomous Stage-Structured Predator-Prey System with Impulsive Effects

Existence of Periodic Solution for a Non-Autonomous Stage-Structured Predator-Prey System with Impulsive Effects Appled ahemacs 55-6 do:.6/am.. Pblshed Onlne arch (hp://www.scrp.org/jornal/am) Exsence o Perodc Solon or a Non-Aonomos Sage-Srcred Predaor-Prey Sysem wh Implsve Eecs Absrac eng W Zolang Xong Ypng Deng

More information

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Ths documen s downloaded from DR-NTU, Nanyang Technologcal Unversy Lbrary, Sngapore. Tle A smplfed verb machng algorhm for word paron n vsual speech processng( Acceped verson ) Auhor(s) Foo, Say We; Yong,

More information

EFFECT OF HEAT FLUX RATIO FROM BOTH SIDE-WALLS ON THERMAL- FLUID FLOW IN CHANNEL

EFFECT OF HEAT FLUX RATIO FROM BOTH SIDE-WALLS ON THERMAL- FLUID FLOW IN CHANNEL 8h AIAA/ASME Jon Thermophyscs and Hea Transfer Conference 4-6 June 00, S. Lous, Mssour AIAA-00-873 00-873 EFFECT OF HEAT FLUX RATIO FROM BOTH SIDE-WALLS ON THERMAL- FLUID FLOW IN CHANNEL SHUICHI TORII

More information

Appendix H: Rarefaction and extrapolation of Hill numbers for incidence data

Appendix H: Rarefaction and extrapolation of Hill numbers for incidence data Anne Chao Ncholas J Goell C seh lzabeh L ander K Ma Rober K Colwell and Aaron M llson 03 Rarefacon and erapolaon wh ll numbers: a framewor for samplng and esmaon n speces dversy sudes cology Monographs

More information

RELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA

RELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA RELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA Mchaela Chocholaá Unversy of Economcs Braslava, Slovaka Inroducon (1) one of he characersc feaures of sock reurns

More information

ABSTRACT. Keywords: Finite Element, Active Optics, Adaptive Optics, WaveFront Error, System Analysis, Optomechanical 1.

ABSTRACT. Keywords: Finite Element, Active Optics, Adaptive Optics, WaveFront Error, System Analysis, Optomechanical 1. Analyss echne or conrollng Sysem Waveron Error wh Acve/Adapve Opcs Vcor L. Genberg*, Gregory J. Mchels Sgmadyne, 83 Wes Ave, ocheser, NY 14611 *genberg@sgmadyne.com (585)35-746 ABSTACT The lmae goal o

More information

Learning Objectives. Self Organization Map. Hamming Distance(1/5) Introduction. Hamming Distance(3/5) Hamming Distance(2/5) 15/04/2015

Learning Objectives. Self Organization Map. Hamming Distance(1/5) Introduction. Hamming Distance(3/5) Hamming Distance(2/5) 15/04/2015 /4/ Learnng Objecves Self Organzaon Map Learnng whou Exaples. Inroducon. MAXNET 3. Cluserng 4. Feaure Map. Self-organzng Feaure Map 6. Concluson 38 Inroducon. Learnng whou exaples. Daa are npu o he syse

More information

Reconstruction of Variational Iterative Method for Solving Fifth Order Caudrey-Dodd-Gibbon (CDG) Equation

Reconstruction of Variational Iterative Method for Solving Fifth Order Caudrey-Dodd-Gibbon (CDG) Equation Shraz Unvery of Technology From he SelecedWor of Habbolla Lafzadeh Reconrcon of Varaonal Ierave Mehod for Solvng Ffh Order Cadrey-Dodd-Gbbon (CDG Eqaon Habbolla Lafzadeh, Shraz Unvery of Technology Avalable

More information

Chapters 2 Kinematics. Position, Distance, Displacement

Chapters 2 Kinematics. Position, Distance, Displacement Chapers Knemacs Poson, Dsance, Dsplacemen Mechancs: Knemacs and Dynamcs. Knemacs deals wh moon, bu s no concerned wh he cause o moon. Dynamcs deals wh he relaonshp beween orce and moon. The word dsplacemen

More information

Qi Kang*, Lei Wang and Qidi Wu

Qi Kang*, Lei Wang and Qidi Wu In. J. Bo-Inspred Compaon, Vol., Nos. /, 009 6 Swarm-based approxmae dynamc opmzaon process for dscree parcle swarm opmzaon sysem Q Kang*, Le Wang and Qd W Deparmen of Conrol Scence and Engneerng, ongj

More information

Minimum Mean Squared Error Estimation of the Noise in Unobserved Component Models

Minimum Mean Squared Error Estimation of the Noise in Unobserved Component Models Mnmm Mean Sqared Error Esmaon of he Nose n Unobserved Componen Models Agsín Maravall ( Jornal of Bsness and Economc Sascs, 5, pp. 115-10) Absrac In model-based esmaon of nobserved componens, he mnmm mean

More information

Macroscopic Traffic Flow Model in Line with Kerner Theory with Consideration of Driver Aggressiveness

Macroscopic Traffic Flow Model in Line with Kerner Theory with Consideration of Driver Aggressiveness Inernaonal Jornal of Engneerng Innovaon & esearch Volme 5, Isse, ISSN: 77 5668 acroscopc Traffc Flow odel n ne wh Kerner Theory wh Consderaon of Drver Aggressveness warana Amendeo*. Kenya Unversy, Nkb

More information

Bernoulli process with 282 ky periodicity is detected in the R-N reversals of the earth s magnetic field

Bernoulli process with 282 ky periodicity is detected in the R-N reversals of the earth s magnetic field Submed o: Suden Essay Awards n Magnecs Bernoull process wh 8 ky perodcy s deeced n he R-N reversals of he earh s magnec feld Jozsef Gara Deparmen of Earh Scences Florda Inernaonal Unversy Unversy Park,

More information

Numerical Simulation on Supersonic Turbulent Flow past Backward Facing Rounded Step Utilizing Hybrid RANS-LES

Numerical Simulation on Supersonic Turbulent Flow past Backward Facing Rounded Step Utilizing Hybrid RANS-LES Numercal Smulaon on Supersonc Turbulen Flow pas Backward Facng Rounded Sep Ulzng Hybrd RANS-LES Absrac Dr. Nrmal Kumar Kund Assocae Professor, Deparmen of Producon Engneerng Veer Surendra Sa Unversy of

More information

Multi-Fuel and Mixed-Mode IC Engine Combustion Simulation with a Detailed Chemistry Based Progress Variable Library Approach

Multi-Fuel and Mixed-Mode IC Engine Combustion Simulation with a Detailed Chemistry Based Progress Variable Library Approach Mul-Fuel and Med-Mode IC Engne Combuson Smulaon wh a Dealed Chemsry Based Progress Varable Lbrary Approach Conens Inroducon Approach Resuls Conclusons 2 Inroducon New Combuson Model- PVM-MF New Legslaons

More information

PHYS 1443 Section 001 Lecture #4

PHYS 1443 Section 001 Lecture #4 PHYS 1443 Secon 001 Lecure #4 Monda, June 5, 006 Moon n Two Dmensons Moon under consan acceleraon Projecle Moon Mamum ranges and heghs Reerence Frames and relae moon Newon s Laws o Moon Force Newon s Law

More information

Math 128b Project. Jude Yuen

Math 128b Project. Jude Yuen Mah 8b Proec Jude Yuen . Inroducon Le { Z } be a sequence of observed ndependen vecor varables. If he elemens of Z have a on normal dsrbuon hen { Z } has a mean vecor Z and a varancecovarance marx z. Geomercally

More information

Real Time Hybrid Simulation using Shaking Tabels

Real Time Hybrid Simulation using Shaking Tabels Real Tme Hybrd Smlaon sng Shakng Tabels Deparmen o Cvl and Envronmenal Engneerng Unversy o Kassel, Germany Olne Inrodcon A ndamenal sbsrcre algorhm wh sb seppng Applcaons o he algorhm Conclsons Inrodcon

More information