Effect of Structure of Internal Waves on Energy Transfer and Flow Control Between South and Middle Basins of Caspian Sea

Size: px
Start display at page:

Download "Effect of Structure of Internal Waves on Energy Transfer and Flow Control Between South and Middle Basins of Caspian Sea"

Transcription

1 Modern Appled Scence; Vol. 10, No. 3; 016 ISSN E-ISSN Publshed by Canadan Center of Scence and Educaton Effect of Structure of Internal Waves on Energy Transfer and Flow Control Between South and Mddle Basns of Caspan Sea Ashrafoalsadat Shekarbaghan 1 1 Organzaton for Educatonal Research and plannng(oerp), Iran Correspondence: Ashrafoalsadat Shekarbaghan, Assstant Professor n Physcs Department of Scence & Mathematcs Educaton, of Educaton Studes Research Center, Organzaton for Educatonal Research and plannng(oerp), Iran, Tehran. E-mal: a0.baghan@gmal.com Receved: December 10, 015 Accepted: January 10, 016 Onlne Publshed: January 13, 016 do: /mas.v10n3p33 URL: Abstract In ths study, role of nternal waves n controllng stratfed exchange flows between south and mddle basns of Caspan Sea and energy transfer have been studed. The studes have been conducted based on feld and theoretcal studes based on numercal calculatons. In ths nvestgaton, Caspan Sea has been consdered as stratfed and Boussnesq approxmaton has been consdered. In equatons governng the movement, coheson and compressblty have been neglected. In the feld studes, nvestgaton of densty profles ndcate sgns of stratfed structure whch so far has been attrbuted to dual dstrbuton convecton phenomenon, but by consderng densty rato profle and values of densty rato whch are mostly negatve numbers, the factor causng stratfed structure of Caspan Sea can be attrbuted to exstence of nternal waves. Densty feld n these two basns ndcate that the mddle basn has hgher average densty compared to the south basn, ths fact can cause creaton of a gravty exchange flow. Ths exchange flow causes creaton of nternal waves n ths sea that can have a modal structure towards vertcal drecton. In theoretcal studes, the two mddle and south basns of Caspan Sea are networked n two dmenson x-z systems and equatons governng movement are solved through a numercal method by usng fnte dfference method and the vertcal velocty values obtaned n 90 ponts of the network was almost 10-6 to 10-4 m/s, meanng that vertcal movements are much smaller. Values obtaned for horzontal velocty are almost 10-3 to 10-1 m/s. Vertcal and horzontal velocty profles obtaned from soluton of these equatons and also densty profle obtaned n four ntegraton perods ndcate that as a result of passage of gravty flow between the two basns, densty n varous ponts of the two basns change. The structure created from nternal waves on the crest causes creaton of shear strata observed. Ths structure can be effectve n exchange of the two basns. Therefore, n transfer and energy converson, dentfcaton of the created strata s of a hgh mportance. Keywords: nternal waves, normal modes, dstrbuton convecton, gravty flow, Caspan Sea 1. Introducton Observatons ndcate exstence of stratfed structure n seas; so that the way of creaton of stratfed structure has always been a queston. Dfferent mechansms such as dual dstrbuton, nternal waves, and turbulent mxng have been proposed for creaton of them []. Ruddck [199] beleves that exstence of stratfed structure depends on dual dstrbuton [3], but n the Caspan Sea the condtons of creaton of dual dstrbuton are less observed. Hence, exstence of nternal waves should cause creaton of stratfed structure. Therefore, mechansms of nternal waves are more acceptable for creaton of a stratfed structure. Exstence of nternal waves n stratfed envronments such as sea s very common. These waves have a more effectve role n transferrng the momentum n sea [4, 5]. Normal modes of these waves cause creaton of shear layers whch have also been observed n expermental medums [6]. The stratfed structure related to normal modes of these waves can be mportant n vertcal and horzontal transfer n marne envronments or dffuson of sound waves whch are usually used for searchng n seas. Wong et al. (001) ndcated that n closed envronments, flows due to plumes that move n the bottom of the closed basn have caused nternal waves that ther phase velocty s towards down and ther group velocty s towards up drectons. Equalty of ther phase velocty wth advecton velocty due to plume wth fllng box 33

2 Modern Appled Scence Vol. 10, No. 3; 016 mechansm, cause creaton of sem-statonary modal structure. Also ths structure has caused creaton of shear strata (5 to 7) whch are effectve n creaton of stratfed structure of physcal quanttes such as salnty and temperature [6]. Vertcal structure of physcal propertes such as salnty and temperature n lakes also ndcate these strata. For nstance, Boehrer (000) by nvestgatng vertcal structure of salnty and temperature n Constance Lake ndcated that regular strata are observed n them so that nternal waves can be analyzed nto modes that each can transfer ndependently, consderng ths fact that ndependent waves n lakes can be observed wth hgher probablty. Therefore, the hypothess that n lakes, all water motons end n nternal waves s a reasonable hypothess. Hence, he attrbuted the profle of flow velocty n Constance Lake to the velocty obtaned from solvng equatons of nternal waves and by proposng an analytc model of normal modes of nternal waves, justfed them [7]. Hogg et al. (000) nvestgated role of nternal waves n controllng stratfed exchange flows between two closed marne basns by solvng equatons governng movement wth numercal methods by determnng Egen value matrx. At frst, created nternal waves n an expermental smulaton n order to create exchange flow n a reservor ncludng the two basns, then by examnng waves n terms of locaton and tme, calculated propagaton velocty. After ths operaton, based on theory of hydraulc control n stratfed current conservaton, flud velocty n defned and separate strata that ther nteracton occurs only n lmted condtons was calculated and then compared wth velocty of nternal waves n the water column [8]. Bans and Turner (1969) ndcated that gravty flows cause modal movements (together wth nternal waves) [9]. Hydrographcal structure of Caspan Sea whch s a closed lake s n such a way that can be consdered n two mddle and south basns. Exstence of a mound between mddle and south basns of Caspan Sea whch s almost located besde Abshouran pennsula, causes that flow exchange between the two basns from top of the mound, cause creaton of nternal waves. Meanng that the gravty flow between the two basns especally on the bottom towards the south basn from top of the mound n the stratfed area, causes creaton of such waves. On the other hand, modal structure of these waves should be effectve n the way of flow exchange between these two basns. In ths research an analytc (numercal) model has been proposed and compared wth feld observatons and the obtaned results are dscussed. These results can be compared wth those obtaned from studes conducted by Boehrer and Hogg et al., and t s found that there s an approprate consstency between these results. Several studes have been conducted on coastal areas and waves caused by wnds n the Caspan Sea, however, few studes have been conducted on marne physcs of varous parts of Caspan Sea. Therefore, t seems that the present study s a new approach n terms of ssues of marne physcs related to nternal parts of the Caspan Sea especally n Abshouran regon.. Feld Study In ths research, the obtaned data by CTD measurements have been used for provdng σ t profle. Plan of the Caspan Sea and CTD measurement statons have been ndcated n the frst research tour of the Internatonal Atomc Energy Agency n fgure 1 [1]. Feld observatons ndcate σ t profle of stratfed structure of the Caspan Sea n fgure 3. σ t values n these waters are often negatve and less than one. 3. Analytc Model Internal waves n an area wth densty stratfcaton can be obtaned by analyss of the movement equatons. In small movements, ampltude of pressure and densty changes are as follows, respectvely: PP(z)+P (x,y,z,t) ρρ(z)+ρ (x,y,z,t) Where, P(Z) s hydrostatc dstrbuton of pressure n water column and P' represents fluctuatons of pressure, ρ(z) s vertcal dstrbuton of densty feld and ρ s also changes related to densty feld. Horzontal components of flud's movement are as follows: ρ(ut-fv)-p'x (1) ρ(vt-fu)-p'y And ts vertcal component s: ρ(wt) -P z - gρ () 34

3 Modern Appled Scence Vol. 10, No. 3; 016 And the ndexes are ndcatng the partal dervatve (for example Ut u t ). Also f Ω snφ s the Corols parameter so that φs the lattude and Ωs the angular velocty of the ground. Mass contnuty equaton for ncompressble condtons s also as follows: u x +v y +w z 0 (3) Also buoyancy conservaton equaton s: 0 d ρ t +w ρ (4) dz The equaton ndcates local changes of densty due to vertcal movements of the flud n densty feld wth vertcal gradent. By combnng these equatons an equaton based on vertcal component of the velocty of flud movement can be obtaned [10]: w N H w f w t z Where N s frequency of buoyancy of the envronment and s: N ρ ρ z, g H + x y Equaton 5 has a wave soluton. For nstance, for small ampltude movements n a flud wth stable densty stratfcaton, ths soluton can be n the followng form: kx ( + ly ωt) w wˆ( z) e (6) Where wz ˆ( ) s vertcal structure of ampltude of vertcal movements that can nclude multple modes of Wn (n0, 1,, 3, 4, ). These modes are orthogonal functons and are known as normal modes. K and L are horzontal components of vector of turbulence wave number andω s turbulence frequency and exponental part of ths term ndcates tme horzontal structure of the turbulence. Ampltude of the vertcal movement of partcles of the flud proportonal to w n s h n and n case where Boussnesq approxmaton s applcable, by replacng equaton 6 nto equaton 5, t can be ndcated that vertcal structure of normal modes s as the followng relaton [10]: d w ˆ ( N ω )( k + l + ) w ˆ 0 (7) dz ω So that for fxed and postve values of N has trgonometrc soluton as follows: Also we have: ( N ω )( k + l ) wˆ sn( ( z+ H) ω (5) n d h dz N + hn 0 Cn (8) Where C n s phase velocty of mode "n" and we have C n ghe n whch e H s the equvalent heght [10]. For n0, movement mode s the same barotropc mode whch ts phase velocty s gh and H s heght of the stratum. For n>0, modes are related to baroclnc movements and 1 C n s the especal value accordng to 35

4 Modern Appled Scence Vol. 10, No. 3; 016 especal functons of h ( z ). For a medum wth fxed N, soluton of equaton 8 wth boundary condtons h 0 n z 0, H s as follows: n n n d h N + hn 0 (9) dz Cn Boundary condtons of h n 0 n surface of the flud stratum ndcates rgd ld approxmaton. Structure of these functons for varous modes s n a way that for n1, the functon has one peak and for n, the functon has two peaks, and etc. By havng w n n accordance wth h n, horzontal veloctes of u n can also be obtaned by usng conservaton equaton [6]. Of course, n the model proposed by Wong et al., effect of molecular coheson has also been consdered whch ndcates that vertcal structure of normal modes have attenuaton. Ths effect (attenuaton) can also be observed n feld and expermental observatons (fgure 6) [11]. Due to exstence of densty dfference between two ponts of the same heght from two basns of the Caspan Sea and also densty dfference between two ponts n the same dstance from coast and of dfferent heghts, value of N s varable. Therefore, there s no analytc soluton for equaton 7 and numercal method should be appled. For ths purpose, frst the two basns n the Caspan Sea are networked and then the equaton governng nternal waves on ponts of the network are solved. Equaton of nternal waves wll be wrtten n drectons of x (along north-south drecton of the Caspan Sea) and z (along vertcal drecton from surface to bottom of the sea) n D coordnate system accordng to equaton 5 by consderng assumptons of ncompressblty and non-cohesveness of the flud and usng Boussnesq approxmaton as follows [10]: w N w t x z x After nettng the basn, orgn of the coordnate s assumed on the southern coast of the Caspan Sea (Nowshahar Cty), then general answer of the equaton s defned as follows [7]: wxzt (,, ) { aw j j( z).cos( kx ω jt) } (11) Consderng boundary condtons for all "J", the followng relaton exsts: j W (0) W ( H) 0 j By replacng equaton 11 n equaton 10, the followng equaton s obtaned: { cos( kx ωjt) a j( W j ωj + W jk ( N ω j) } 0 (1) j j j ωj j ωj (10) W + W k ( N ) 0 (13) u w Consderng conservaton equaton of + 0 x z, we have: 1 uxzt (,, ) Σj{ aw j j ( z) sn( kx ω jt)} k 1 U j Wj (15) k Frst, Wj, vertcal velocty of nternal waves, s obtaned n the J defned n equaton 13, then Uj whch s horzontal velocty of nternal waves n the J defned n equaton 15 s obtaned. (14) For calculatng Wj, equaton 13 s solved by usng fnte dfference method whch means the term z w to the defned and lmted value. For ths purpose, backward dfference method s used whch means the dervatve s 36

5 Modern Appled Scence Vol. 10, No. 3; 016 stuated as functon n varous ponts, n fact dervatve s made separated by Taylor Seres: W W + W j, j,( 1) j,( ) Δz j, 1 dw dz (16) By replacng equaton 16 n equaton 13, the followng relaton s obtaned: Wj, Wj,( 1) + W j,( ) N j,( 1) + k 1 W j,( 1) 0 Δz ω j, K s wavenumber and L s wavelength. If wavelength s assumed to be one eghth of length of the regon, by usng relaton of the wavenumber; (17) π k (18) L wavenumber value for calculatons obtaned as equal to 5 1 π 10 m. N s the buoyancy frequency, ω s wave frequency, C s wave propagaton velocty, g' s the reduced gravty ntensty, and h s the heght of water from sea level, then the followng relatons are obtaned: N Δρ. g Δ z ρ N j,( 1) ρ ρ ρ g + ρ Δz j,( 1) j, j,( 1) j, (19) ω kc 1/ ω j ( 1) kgh ( ) (0) ω j,( 1) k ρ ρ ρ j,( 1) j, + ρ j,( 1) j, gh N ( j, 1) 1 ω ( j, 1) k Δzh (1) It s observed that for solvng these equatons, frst j and should be defned n the regon beng studed. Hence, networkng method should be used. W j, ndcates vertcal velocty n specfc and j whch are stuated n a specfc pont, Wj, 1 ndcates vertcal velocty n a pont wth dstance Δx from the coast and wth heght dfference equal to Δ z before the pont (,j), and Wj,( ) ndcates vertcal velocty n a pont wth Δx dstance from the coast and heght dfference equal to Δ z before the pont (,j). 3.1 Networkng Method In order to calculate W, the regon of the study should be networked n terms of and j.for ths purpose, accordng to fgure, we dvde X axs nto 1 equal dvsons wth 70 km dstance from each other, ndcate j on the X axs so that we would have ponts from 1 to 13 and we dvde Z axs nto 40 equal dvsons wth 0 km dstance from each other, ndcate on Z axs so that we would have from 1 to 40. In ths case, j1 and 1 are related to orgn ofd coordnate system, so for and j we wll have: j1,, 3, 4,,13 37

6 Modern Appled Scence Vol. 10, No. 3; 016 1,, 3, 4,,40 Therefore, dstance between two successve j parameter s 70km and the dstance between two successve parameters s 0km. Dstance between two successve parameters s equal to Δz between the two pontsand h value whch s heght from sea level s equal to ( 1) Δ z: h ( 1) Δ z 0( 1) () 3. Calculaton of Vertcal Velocty Accordng to the networkng method, frst we compare varous ponts of the network wth frst tour research statons of the Internatonal Atomc Energy Agency: j1 s located n margn of south coast of the Caspan Sea and border of the south basn s northern coast of Iran, n ths pont velocty s equal to zero, whch means: W(1, ) 0 j s located between the south coast of Caspan Sea and staton 1 of frst marne tour of Internatonal Atomc Energy Agency. For calculatng data related to ths pont, nterpolaton method s used and the data related to j whch extends n 70km dstance from coast of Caspan Sea s obtaned from surface to bottom of the sea. j3 s almost located on staton 1 of frst marne tour of Internatonal Atomc Energy Agency. The data of CTD of the staton can be used for j3 whch extends n 140km dstance from coast of Caspan Sea from surface to bottom of the sea. Accordngly, the pont j4 s almost located on staton, the pont j5 s almost located on staton 3, the pont j6 s almost located on staton 4, the pont j7 s almost located on staton 5, the pont j8 s almost located on staton 6, the pont j9 s almost located on staton 7, the pont j10 s almost located on staton 8, the pont j11 s almost located on staton 9, and the pont j13 s almost located on staton 11, so that we can apply the measurements related to marne frst tour statons as the mentoned order for ponts wth specfc J. But j1 s almost located on staton 11 of frst marne tour of Internatonal Atomc Energy Agency so that was deleted because of low depth of measurements of frst tour of Internatonal Atomc Energy Agency, so the requred data for j1 whch have been obtaned by nterpolaton are used. Therefore, ponts j1 and j13 are located n northern part of the Caspan Sea and these regons have low depth, then because of low depth we can assume the velocty n pont 13 as zero: W(13, ) 0 So for solvng equaton 13 we can apply the followng hypothess as one of the boundary condtons: W(1, ) W(13, ) 0 Also vertcal velocty n bottom of the sea s equal to zero. Therefore, for each specfc j, consderng topography of bottom of the Caspan Sea, when reachng bottom of the sea, vertcal velocty of the sad pont becomes zero. Equaton 13 after replacement and smplfcaton s wrtten as below: 1 8 Wj, π 10 Wj,( 1) + Wj,( ) 0 1 (3) We solve equaton 3 by usng MATLAB.6 software, we obtan vertcal velocty n each node, meanng that we obtan W j n certan j wth dfferent parameters and hence, frst we obtan W j n 1,,. We follow the same process for all dfferent j parameters so that W j profles are respectvely depcted. Profles of vertcal velocty n statons 4, 5, and 6 are observed n fgure 4. In fact these profles ndcate vertcal modes of nternal waves. 3.3 Calculaton of Horzontal Velocty After calculatng W j, we fnd values of U j, and snce U j, ndcates horzontal velocty n pont (,j), for calculatng U j, we use the followng relatons: U U 1 W k W, k z j, j, j, 1 j 38

7 Modern Appled Scence Vol. 10, No. 3; 016 W j, ku j, z Hence, by applyng fnte dfference method through backward Euler method, equaton 4 s solved: W U j, j, 1 j, (4) W ku j, Δz ( Wj, Wj, 1) (5) 0k Therefore, by havng values of Wj, and Wj, 1, we can obtan U j, n each pont. Lke calculaton of W j,, we also fnd U j, n the certan j for parameters exsted, then U j profle lke W j profle s drawn n varous j parameters. Profles of horzontal velocty n statons 4, 5, and 6 are ndcated n fgure 5. In fact these profles ndcate vertcal and horzontal modes of nternal waves. 3.4 Densty Calculaton Method After calculatng U and W, densty value wll be consdered. Accordng to mass conservaton equaton we have: ρ u ρ + + w ρ 0 (6) t x z We use fnte dfference method by usng backward Euler method for smplfyng equaton 6, so that we have: ρ x j, ρ ρ j, ( j 1), Δx (7) ρ z j, ρ ρ j, j,( 1) Δz (8) We replace equatons 7 and 8 n equaton 6 and obtan equaton 9: ρ ρj, ρ( j 1), ρj, ρj,( 1) U j, Wj, t Δx Δz (9) 4 Δ x 7 10 m Δ z 0m Hence, equaton 9 s solved and value of ρ t for certan and j can be obtaned. Also t can be noted that: ρ Δρ t Δt ρ Δρ Δt t (30) when we have Δ t assumed as regulaton tme or collapse tme of mxed water column after mxng. Consderng the dstance between the frst two ponts symmetry to each other n two sdes of the crest between south and mddle basns whch have notceable densty dfference, we obtan value of L from equaton 31, then by consderng propertes related to these ponts whch one belongs to j6 and the other one belongs to j8 (j7 s above the stuated obstacle), we can obtan Δ t whch s collapse tme [1]. L Δ t 1/ ( gh ) (31) Value of Δt s obtaned almost as equal to 9300 sec, value of Δρ s obtaned n each staton for less values 39

8 Modern Appled Scence Vol. 10, No. 3; 016 than ths tme as 800sec, and more values of tme as 45600sec and sec and sec. Also we have: ρ ρ0 + ρ That also can be wrtten as follows: ρ ρ0 +Δ ρ (3) Therefore, by havng Δρ and ρ0 we obtan value of ρ for the sad tmes so that we can compare densty profles obtaned wth the prmary densty profle (drawn by data resulted from fst research tour measurements of UNESCO). Accordng to the comparson t can be observed that the calculated densty profle has an approprate coordnaton wth the prmary densty profle (fgures 6). Meanng that nternal waves are effectve n structure of strata of waters of south and mddle bass of Caspan Sea. 4. Effect of Internal Waves on Energy Transfer Accordng to the observatons from feld studes (drawng densty rato profle), t was found that dual dstrbuton convecton s not often created n Caspan Sea and t doesn t seems that ths convecton has a man role n stratfed structure n the Caspan Sea, but also waves have a man role n these waters. Feld studes confrm exstence of surface flows from north to south and also creaton of flow from surface water of the mddle basn and upwellng of ths flow from depth. Transfer of ths flow occurs together wth water stratfcaton and creaton of nternal waves. Also analytc studes confrm exstence of such waves and water stratfcaton and fgure 4 ndcates exstence of vertcal modes of nternal waves n Caspan Sea. Fgure 6 ndcates approprate consstency exstng between calculaton results and results obtaned from frst tour research measurements. Ths fact mples exstence of nternal waves whch cause creaton of stratfed structure. Generally, n the Caspan Sea, normal modes of nternal waves are produced wth long wavelength by entrance of rvers n an envronment wth stratfed layers, and fnally result n complete stratfcaton of Caspan Sea, so that wthn long perod of tme, the flow transfers from one stratum to another, movement of ths flow causes exctaton of the nternal waves whch propagate energy towards vertcal drecton. Internal waves, n envronments wth stable densty stratfcaton lke Caspan Sea waters, can transfer energy. Meanng that ths flow, together wth changes n velocty due to passage of varous strata of water wth dfferent denstes, reaches water surface and then the phenomenon of water rotaton on surface of the sea can be observed. Sometmes when nternal waves break, turbulence caused by ths event, creates stratfed structure [10]. The man reason of creaton of exchange flow between the two basns of the Caspan Sea s horzontal gradent of densty, that n ths phenomenon, nternal waves are created and have effect on ths exchange, and sgn of exstence of such nternal waves s stratfed structure observed. The Stratfed structure n Caspan Sea manly exsts n the part under the mxed stratum. Thckness of the strata n Caspan Sea waters s almost 10m to 0m and the number of them vares from 3 to 7. Exstence of stratfed structure also has effect on thermal dstrbuton coeffcent and cause transfer of heat n Caspan Sea waters and s also effectve n extent of exchange of waters. Also exstence of varous water strata wth dfferent denstes s effectve n propagaton of sound waves whch s consdered as mechancal waves, so that the more an stratum s denser, the more velocty of sound wave ncreases n that stratum. Hence, n varous water strata, sound wave s also dfferent and at the tme of passage of sound from one stratum to another, falure phenomenon for sound waves occurs. For convertng knematc energy of nternal waves created along vertcal drecton nsde Caspan Sea waters to the electrcal energy, buoyancy systems whch functon wth fluctuatons of nternal waves for drvng hydraulc pumps, are used. Ths system s nstalled on a buoyant rod n the sea and s fxed on a devce n bottom of the sea and a seres of buoyant devces attached to t or restraned by t fluctuate wth that wave. Ths movement strkes electrc generator and creates electrcty and electrc energy produced s transferred by electrcty cable under the sea to the coast, and for changng knematc energy of nternal waves created along vertcal drecton of Caspan Sea waters to electrc energy, some devces whch work wth water fluctuaton are used, so that by wave moton of water, and fluctuaton of water, water column s flled and empted. In ths process, water s entered the compressed column and creates energy by workng, then ths energy s taken under control and transferred to the coast by electrcty cable. 5. Concluson Consderng numercal calculatons, value of vertcal velocty related to movement of nternal waves s obtaned as 10-6 to 10-3 m/s, meanng that vertcal movements of nternal waves occurs slowly. The calculated horzontal 40

9 Modern Appled Scence Vol. 10, No. 3; 016 velocty s also obtaned as 10-3 to 10-1 m/s. Profle of vertcal and horzontal veloctes obtaned from solvng equatons and also profle of densty s drawn after pass of four dfferent tme ntegraton steps. These profles ndcate that due to passage of gravtatonal flow between the two basns, densty n varous ponts of the two basns change, and snce south and mddle basns of Caspan Sea are related to each other through a mound (Abshouran pont), n the conducted studes t s observed that flow of the mddle basn moves towards the south basn and overflows from top of the mound between the two basns, then from bottom of the mound, the flow returns towards the mddle basn. And ths fact s consstent wth measurements conducted by UNESCO (1995). Snce for creaton of dual dstrbuton convecton densty rato should be almost equal to 1, ths phenomenon does not often occur n ths sea (fgure 3), hence we cannot attrbute the observed strata to ths phenomenon. But the modal structure created from nternal waves on the mound causes creaton of shear strata observed. The created strata, n addton to dstrbutng organc substance and also oxygen exstng n Caspan Sea waters, causes energy transfer that by usng energy convertor of waves n ths place, electrc energy s produced. References Bans, W. D., & Turner, I. S. (1969). Turbulent buoyant Convecton from a source n a confned regon. J. Flud Mech., 34, Bdokht, A. A., & Grffths, R. W. (001). Internal waves as an addtonal mechansm for layerng n the outflows from sem enclosed seas.14th Australasan Flud Mechancs Conference. Bdokht, A. A., & Norooz, M. (004). A Physcal Model for the layered structure of a densty drven flow over a slope. Tenth Asan Congress of flud Mechancs, Ser Lanka. Bdokht, A. A., & Shekarbaghan, A. (011). The layered structure n exchange flows between two basns (Mddle and Southern basns of the Caspan Sea). Int. J. Mar. Sc. Eng., 1(1), 13-. Boehrer, B. (000). Modal response of a deep stratfed Lake: Western lake Constance. J. Geophys. Res., 105, Fedrov, K. N. (013). The Physcal Nature and Structure of Oceanc Fronts. Gll, A. E. (198). Atmoshperc and ocean dynamcs. Academc press N. Grffths, R. E. (008). Interleavng ntrusons produced by nternal waves: a laboratory experment. J. Flud Mech., 60, Hogg, A., Mc, C., Wnters, K. B., & Ivey, G. N. (001). Lnear nternal Waves and the control of stratfed exchange flows. J. Flud Mech., 447, Ruddck, B. R. (199). Intrusve mxng n Medterranean salt lens ntruson slopes and dynamcal mechansms. J. Phys. Oceanogr,, UNESCO, IHP-IOC-IAEA. (1995).Workshop on sea level rse and multdscplnary studes of envronmental processes n the Caspan regon 9-1 May. Pars, Farance IOC workshop, 108. Wong, B. D., Grffths, R. W., & Hughes, G. O. (001). Shear Layers drven by turbulent plumes. J. Flud Mech, 434, Appendx A 41

10 Modern Appled Scence Vol. 10, No. 3; 016 Fgure A1. Map of the Caspan Sea and the statons of measurement Fgure A. Bottom topography of the Caspan Sea s ndcated by dstance from southern shore n southern northern drecton 4

11 Modern Appled Scence Vol. 10, No. 3; 016 Fgure A3. Vertcal profle of densty raton n statons 4, 5, and 6 Fgure A4. Profles of vertcal velocty calculated by numercal method n statons 4, 5, and 6 (x dstances are from Caspan Sea south coast) Fgure A5. Profles of horzontal velocty calculated by numercal method n statons 4, 5, and6 (x dstances are from Caspan Sea south coast) 43

12 Modern Appled Scence Vol. 10, No. 3; 016 depth (m ) potental densty (kg/m*3) [1] [] depth (m ) potental densty 0 (kg/m^3) [1] [] s4 s5 depth (m ) potental densty (kg/m*3) [1] [] depth (m ) potental densty (kg/m^3) [1] [] s6 s7 Fgure A6. Comparson between the prmary potental densty profle obtaned from frst research tour measurements of UNESCO [1] and potental densty profle calculated by numercal method [] n statons 4, 5, 6, and 7 Fgure A7. Horzontal velocty made dmensonless by nternal velocty of mxng obtaned n two successve tmes created n the same stratfed structure by a plume n a closed envronment wth densty stratfcaton Copyrghts Copyrght for ths artcle s retaned by the author(s), wth frst publcaton rghts granted to the journal. Ths s an open-access artcle dstrbuted under the terms and condtons of the Creatve Commons Attrbuton lcense ( 44

Physics 181. Particle Systems

Physics 181. Particle Systems Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,

More information

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy Comparatve Studes of Law of Conservaton of Energy and Law Clusters of Conservaton of Generalzed Energy No.3 of Comparatve Physcs Seres Papers Fu Yuhua (CNOOC Research Insttute, E-mal:fuyh1945@sna.com)

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

The Governing Equations

The Governing Equations The Governng Equatons L. Goodman General Physcal Oceanography MAR 555 School for Marne Scences and Technology Umass-Dartmouth Dynamcs of Oceanography The Governng Equatons- (IPO-7) Mass Conservaton and

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Mathematical Preparations

Mathematical Preparations 1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1

More information

The Tangential Force Distribution on Inner Cylinder of Power Law Fluid Flowing in Eccentric Annuli with the Inner Cylinder Reciprocating Axially

The Tangential Force Distribution on Inner Cylinder of Power Law Fluid Flowing in Eccentric Annuli with the Inner Cylinder Reciprocating Axially Open Journal of Flud Dynamcs, 2015, 5, 183-187 Publshed Onlne June 2015 n ScRes. http://www.scrp.org/journal/ojfd http://dx.do.org/10.4236/ojfd.2015.52020 The Tangental Force Dstrbuton on Inner Cylnder

More information

Spin-rotation coupling of the angularly accelerated rigid body

Spin-rotation coupling of the angularly accelerated rigid body Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s

More information

Grid Generation around a Cylinder by Complex Potential Functions

Grid Generation around a Cylinder by Complex Potential Functions Research Journal of Appled Scences, Engneerng and Technolog 4(): 53-535, 0 ISSN: 040-7467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around

More information

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850) hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity 1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum

More information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LIX, 013, f.1 DOI: 10.478/v10157-01-00-y ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION BY ION

More information

Tensor Smooth Length for SPH Modelling of High Speed Impact

Tensor Smooth Length for SPH Modelling of High Speed Impact Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices Amplfcaton and Relaxaton of Electron Spn Polarzaton n Semconductor Devces Yury V. Pershn and Vladmr Prvman Center for Quantum Devce Technology, Clarkson Unversty, Potsdam, New York 13699-570, USA Spn Relaxaton

More information

Research Article Green s Theorem for Sign Data

Research Article Green s Theorem for Sign Data Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of

More information

The Symmetries of Kibble s Gauge Theory of Gravitational Field, Conservation Laws of Energy-Momentum Tensor Density and the

The Symmetries of Kibble s Gauge Theory of Gravitational Field, Conservation Laws of Energy-Momentum Tensor Density and the The Symmetres of Kbble s Gauge Theory of Gravtatonal Feld, Conservaton aws of Energy-Momentum Tensor Densty and the Problems about Orgn of Matter Feld Fangpe Chen School of Physcs and Opto-electronc Technology,Dalan

More information

FREQUENCY DISTRIBUTIONS Page 1 of The idea of a frequency distribution for sets of observations will be introduced,

FREQUENCY DISTRIBUTIONS Page 1 of The idea of a frequency distribution for sets of observations will be introduced, FREQUENCY DISTRIBUTIONS Page 1 of 6 I. Introducton 1. The dea of a frequency dstrbuton for sets of observatons wll be ntroduced, together wth some of the mechancs for constructng dstrbutons of data. Then

More information

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

More information

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

Perfect Fluid Cosmological Model in the Frame Work Lyra s Manifold

Perfect Fluid Cosmological Model in the Frame Work Lyra s Manifold Prespacetme Journal December 06 Volume 7 Issue 6 pp. 095-099 Pund, A. M. & Avachar, G.., Perfect Flud Cosmologcal Model n the Frame Work Lyra s Manfold Perfect Flud Cosmologcal Model n the Frame Work Lyra

More information

Implicit Integration Henyey Method

Implicit Integration Henyey Method Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure

More information

Normally, in one phase reservoir simulation we would deal with one of the following fluid systems:

Normally, in one phase reservoir simulation we would deal with one of the following fluid systems: TPG4160 Reservor Smulaton 2017 page 1 of 9 ONE-DIMENSIONAL, ONE-PHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases

More information

2 Finite difference basics

2 Finite difference basics Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T

More information

11. Dynamics in Rotating Frames of Reference

11. Dynamics in Rotating Frames of Reference Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons

More information

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

Electrical double layer: revisit based on boundary conditions

Electrical double layer: revisit based on boundary conditions Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer

More information

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent

More information

1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave

1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave 014 4th Internatonal Conference on Future nvronment and nergy IPCB vol.61 (014) (014) IACSIT Press, Sngapore I: 10.776/IPCB. 014. V61. 6 1-mensonal Advecton-ffuson Fnte fference Model ue to a Flow under

More information

Publication 2006/01. Transport Equations in Incompressible. Lars Davidson

Publication 2006/01. Transport Equations in Incompressible. Lars Davidson Publcaton 2006/01 Transport Equatons n Incompressble URANS and LES Lars Davdson Dvson of Flud Dynamcs Department of Appled Mechancs Chalmers Unversty of Technology Göteborg, Sweden, May 2006 Transport

More information

Temperature. Chapter Heat Engine

Temperature. Chapter Heat Engine Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the

More information

FORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO

FORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO ISTP-,, PRAGUE TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA FORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO Mohammad Rahnama*, Seyed-Mad Hasheman*, Mousa Farhad**

More information

Linear Momentum. Center of Mass.

Linear Momentum. Center of Mass. Lecture 6 Chapter 9 Physcs I 03.3.04 Lnear omentum. Center of ass. Course webste: http://faculty.uml.edu/ndry_danylov/teachng/physcsi Lecture Capture: http://echo360.uml.edu/danylov03/physcssprng.html

More information

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11) Gravtatonal Acceleraton: A case of constant acceleraton (approx. hr.) (6/7/11) Introducton The gravtatonal force s one of the fundamental forces of nature. Under the nfluence of ths force all objects havng

More information

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH Turbulence classfcaton of load data by the frequency and severty of wnd gusts Introducton Oscar Moñux, DEWI GmbH Kevn Blebler, DEWI GmbH Durng the wnd turbne developng process, one of the most mportant

More information

Title: Radiative transitions and spectral broadening

Title: Radiative transitions and spectral broadening Lecture 6 Ttle: Radatve transtons and spectral broadenng Objectves The spectral lnes emtted by atomc vapors at moderate temperature and pressure show the wavelength spread around the central frequency.

More information

Solutions to Exercises in Astrophysical Gas Dynamics

Solutions to Exercises in Astrophysical Gas Dynamics 1 Solutons to Exercses n Astrophyscal Gas Dynamcs 1. (a). Snce u 1, v are vectors then, under an orthogonal transformaton, u = a j u j v = a k u k Therefore, u v = a j a k u j v k = δ jk u j v k = u j

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

The Analysis of Convection Experiment

The Analysis of Convection Experiment Internatonal Conference on Appled Scence and Engneerng Innovaton (ASEI 5) The Analyss of Convecton Experment Zlong Zhang School of North Chna Electrc Power Unversty, Baodng 7, Chna 469567@qq.com Keywords:

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

Air Age Equation Parameterized by Ventilation Grouped Time WU Wen-zhong

Air Age Equation Parameterized by Ventilation Grouped Time WU Wen-zhong Appled Mechancs and Materals Submtted: 2014-05-07 ISSN: 1662-7482, Vols. 587-589, pp 449-452 Accepted: 2014-05-10 do:10.4028/www.scentfc.net/amm.587-589.449 Onlne: 2014-07-04 2014 Trans Tech Publcatons,

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

A large scale tsunami run-up simulation and numerical evaluation of fluid force during tsunami by using a particle method

A large scale tsunami run-up simulation and numerical evaluation of fluid force during tsunami by using a particle method A large scale tsunam run-up smulaton and numercal evaluaton of flud force durng tsunam by usng a partcle method *Mtsuteru Asa 1), Shoch Tanabe 2) and Masaharu Isshk 3) 1), 2) Department of Cvl Engneerng,

More information

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth

More information

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI

More information

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

COMPLEX NUMBERS AND QUADRATIC EQUATIONS COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

An identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites

An identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites IOP Conference Seres: Materals Scence and Engneerng PAPER OPE ACCESS An dentfcaton algorthm of model knetc parameters of the nterfacal layer growth n fber compostes o cte ths artcle: V Zubov et al 216

More information

The effect of a background shear current on large amplitude internal solitary waves

The effect of a background shear current on large amplitude internal solitary waves PHYSICS OF FLUIDS 8, 03660 006 The effect of a background shear current on large ampltude nternal soltary waves Wooyoung Cho a Department of Mathematcal Scences, Center for Appled Mathematcs and Statstcs,

More information

HYBRID LBM-FVM AND LBM-MCM METHODS FOR FLUID FLOW AND HEAT TRANSFER SIMULATION

HYBRID LBM-FVM AND LBM-MCM METHODS FOR FLUID FLOW AND HEAT TRANSFER SIMULATION HYBRID LBM-FVM AND LBM-MCM METHODS FOR FLUID FLOW AND HEAT TRANSFER SIMULATION Zheng L a,b, Mo Yang b and Yuwen Zhang a* a Department of Mechancal and Aerospace Engneerng, Unversty of Mssour, Columba,

More information

Lecture 5.8 Flux Vector Splitting

Lecture 5.8 Flux Vector Splitting Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

Uncertainty and auto-correlation in. Measurement

Uncertainty and auto-correlation in. Measurement Uncertanty and auto-correlaton n arxv:1707.03276v2 [physcs.data-an] 30 Dec 2017 Measurement Markus Schebl Federal Offce of Metrology and Surveyng (BEV), 1160 Venna, Austra E-mal: markus.schebl@bev.gv.at

More information

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information

The Finite Element Method

The Finite Element Method The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of

More information

Supplementary Information for Observation of Parity-Time Symmetry in. Optically Induced Atomic Lattices

Supplementary Information for Observation of Parity-Time Symmetry in. Optically Induced Atomic Lattices Supplementary Informaton for Observaton of Party-Tme Symmetry n Optcally Induced Atomc attces Zhaoyang Zhang 1,, Yq Zhang, Jteng Sheng 3, u Yang 1, 4, Mohammad-Al Mr 5, Demetros N. Chrstodouldes 5, Bng

More information

1. Governing Equations

1. Governing Equations 1. Governng Equatons 1a. Governng Equatons for Mean Varables The governng equatons descrbe the varaton n space and tme of the zonal, merdonal and vertcal wnd components, densty, temperature, specfc humdty

More information

829. An adaptive method for inertia force identification in cantilever under moving mass

829. An adaptive method for inertia force identification in cantilever under moving mass 89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,

More information

Numerical Solution of Boussinesq Equations as a Model of Interfacial-wave Propagation

Numerical Solution of Boussinesq Equations as a Model of Interfacial-wave Propagation BULLETIN of the Malaysan Mathematcal Scences Socety http://math.usm.my/bulletn Bull. Malays. Math. Sc. Soc. (2) 28(2) (2005), 163 172 Numercal Soluton of Boussnesq Equatons as a Model of Interfacal-wave

More information

NUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD

NUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD THERMAL SCIENCE: Year 2018, Vol. 22, No. 5, pp. 1955-1962 1955 NUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD Introducton by Tomok IZUMI a* and

More information

Advanced Quantum Mechanics

Advanced Quantum Mechanics Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

MEASUREMENT OF MOMENT OF INERTIA

MEASUREMENT OF MOMENT OF INERTIA 1. measurement MESUREMENT OF MOMENT OF INERTI The am of ths measurement s to determne the moment of nerta of the rotor of an electrc motor. 1. General relatons Rotatng moton and moment of nerta Let us

More information

Frequency dependence of the permittivity

Frequency dependence of the permittivity Frequency dependence of the permttvty February 7, 016 In materals, the delectrc constant and permeablty are actually frequency dependent. Ths does not affect our results for sngle frequency modes, but

More information

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15 NGN 40 ynamcs and Vbratons Homework # 7 ue: Frday, Aprl 15 1. Consder a concal hostng drum used n the mnng ndustry to host a mass up/down. A cable of dameter d has the mass connected at one end and s wound/unwound

More information

/ n ) are compared. The logic is: if the two

/ n ) are compared. The logic is: if the two STAT C141, Sprng 2005 Lecture 13 Two sample tests One sample tests: examples of goodness of ft tests, where we are testng whether our data supports predctons. Two sample tests: called as tests of ndependence

More information

Research & Reviews: Journal of Engineering and Technology

Research & Reviews: Journal of Engineering and Technology Research & Revews: Journal of Engneerng and Technology Case Study to Smulate Convectve Flows and Heat Transfer n Arcondtoned Spaces Hussen JA 1 *, Mazlan AW 1 and Hasanen MH 2 1 Department of Mechancal

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton

More information

Thermodynamics General

Thermodynamics General Thermodynamcs General Lecture 1 Lecture 1 s devoted to establshng buldng blocks for dscussng thermodynamcs. In addton, the equaton of state wll be establshed. I. Buldng blocks for thermodynamcs A. Dmensons,

More information

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5 - HIKARI Ltd, www.m-har.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for One-Characterstc Dfferental Equaton

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices Internatonal Mathematcal Forum, Vol 11, 2016, no 11, 513-520 HIKARI Ltd, wwwm-hkarcom http://dxdoorg/1012988/mf20166442 The Jacobsthal and Jacobsthal-Lucas Numbers va Square Roots of Matrces Saadet Arslan

More information

Turbulent Flow. Turbulent Flow

Turbulent Flow. Turbulent Flow http://www.youtube.com/watch?v=xoll2kedog&feature=related http://br.youtube.com/watch?v=7kkftgx2any http://br.youtube.com/watch?v=vqhxihpvcvu 1. Caothc fluctuatons wth a wde range of frequences and

More information

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force. The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,

More information

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on

More information

Chapter 3. Estimation of Earthquake Load Effects

Chapter 3. Estimation of Earthquake Load Effects Chapter 3. Estmaton of Earthquake Load Effects 3.1 Introducton Sesmc acton on chmneys forms an addtonal source of natural loads on the chmney. Sesmc acton or the earthquake s a short and strong upheaval

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

and Statistical Mechanics Material Properties

and Statistical Mechanics Material Properties Statstcal Mechancs and Materal Propertes By Kuno TAKAHASHI Tokyo Insttute of Technology, Tokyo 15-855, JAPA Phone/Fax +81-3-5734-3915 takahak@de.ttech.ac.jp http://www.de.ttech.ac.jp/~kt-lab/ Only for

More information

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018 MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

Higher Order Wall Boundary Conditions for Incompressible Flow Simulations

Higher Order Wall Boundary Conditions for Incompressible Flow Simulations THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN KOREA OCTOBER 7-30 003 Hgher Order Wall Boundary Condtons for Incompressble Flow Smulatons Hdetosh Nshda. Department of Mechancal and System Engneerng

More information

The Solution of the Two-Dimensional Gross-Pitaevskii Equation Using Lattice-Boltzmann and He s Semi-Inverse Method

The Solution of the Two-Dimensional Gross-Pitaevskii Equation Using Lattice-Boltzmann and He s Semi-Inverse Method Internatonal Journal of Mathematcal Analyss Vol., 7, no., 69-77 HIKARI Ltd, www.m-hkar.com https://do.org/.988/jma.7.634 The Soluton of the Two-Dmensonal Gross-Ptaevsk Equaton Usng Lattce-Boltzmann and

More information

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA 14 th Internatonal Users Conference Sesson: ALE-FSI Statstcal Energy Analyss for Hgh Frequency Acoustc Analyss wth Zhe Cu 1, Yun Huang 1, Mhamed Soul 2, Tayeb Zeguar 3 1 Lvermore Software Technology Corporaton

More information

THEOREMS OF QUANTUM MECHANICS

THEOREMS OF QUANTUM MECHANICS THEOREMS OF QUANTUM MECHANICS In order to develop methods to treat many-electron systems (atoms & molecules), many of the theorems of quantum mechancs are useful. Useful Notaton The matrx element A mn

More information