Asymptotical Analysis of Internal Gravity Wave. Dynamics in Stratified Medium

Size: px
Start display at page:

Download "Asymptotical Analysis of Internal Gravity Wave. Dynamics in Stratified Medium"

Transcription

1 Alie Mahemaical ciences Vol 8 4 no HIKARI L wwwm-hikaricom h://oiorg/988/ams43637 Asmoical Analsis of Inernal Gravi ave namics in raifie Meium V V Bulaov an Yu V Vlaimirov Insiue for Problems in Mechanics RA Pr Vernaskogo - Moscow 956 Russia fa: Corigh 4 V V Bulaov an Yu V Vlaimirov This is an oen access aricle isribue uner he Creaive Commons Aribuion License which ermis unresrice use isribuion an reroucion in an meium rovie he original work is roerl cie Absrac In aer funamenal roblems of inernal gravi waves namics are consiere ave namics of sraifie meium ocean amoshere is highl eenen on ensi isribuion an boom oograh The eac analical soluion is obaine onl if he waer isribuion ensi an boom shae escribe b some moel funcions The soluion of his roblem is eresse in erms of he Green s funcion an he asmoic reresenaions of he soluions are consiere The uniform asmoic forms of he inernal gravi waves in horionall inhomogeneous an non-saionar sraifie ocean are obaine A moifie saio-emoral ra meho is roose which belongs o he class of geomerical oics mehos KBJ meho The soluion is roose in erms of wave moes roagaing ineenenl a he aiabaic aroimaion an escribe as a non-inegral egree series of a small arameer characeriing he sraifie meium roeries Analical an numerical algorihms of inernal gravi wave calculaions for he real ocean arameers are resene Kewors: KBJ aroimaions wave roagaion sraifie meium inernal gravi waves

2 8 V V Bulaov an Yu V Vlaimirov Inroucion Inusrial aciviies on he coninenal shelf connece wih oil gas an oher minerals eracion became one of he imoran reasons o begin researches of namic inernal gravi waves [-3] his an laforms bus wih rilling an consrucion a he eh use long ubes joining hem wih he sea boom Builers of unerwaer consrucions in euaorial isrics eerience he influence of huge unerwaer inernal waves an srong surface flows which can have he form of see waerfalls ome ime ago when he henomenon of inernal waves an heir srengh were no known i haene ha he builers los heir euimen uch eensive losses mae hem hink ha securi of unerwaer euimen an he influence of inernal gravi waves shoul be conrolle The inernal waves characerisics are use for areciaion of heir influence on he environmen an unerwaer laforms of oil an gas eosis a he shelf Arcic basin China an Yellow eas ec [4-] ecial ineres o research of inernal gravi waves is connece wih also inensive eloiaion of Arcic an is naural wealh Inernal gravi waves in Arcic are oorl suie as he move uner ice an racicall invisible from above bu accessible informaion abou unerwaer objecs movemen show heir eisence omeimes here are eclusions when inernal gravi waves reach ice an ulif an lower i wih efinie erioici which can be fie wih he hel of raiolocaion souning Influence of all kins of waves can be he reason of he ice cover sli in he Arcic Inernal waves make for he movemen of icebergs an ifferen kins of olluion o he research of inernal gravi waves namics is an imoran funamenal scienific an racical roblem aime a ensuring securi while [89] The inernal gravi waves are he oscillaions of a sraifie meium in he gravi force fiel The sraifie meium is such a meium where he ensi increases wih he eh uose ha a volume elemen of he meium is no a he euilibrium for eamle i coul be islace uwar hen i will be heavier han he surrouning meium an herefore Archimeean forces will make i move back o he euilibrium The essenial arameer of an oscillaing ssem is he freuenc I is eermine b he correlaion of wo facors: reurning forces which reurn he erurbe ssem owars is euilibrium an he inerial forces For he inernal gravi waves he reurning forces are roorional o he verical graien of he flui s ensi an he inerial ones are roorional o he ensi iself For he characerisic freuenc of he gravi waves oscillaions we have g he following eression: ρ This freuenc is usuall ρ calle b he Brun-Vaisala freuenc or he buoanc freuenc Here ρ is he ensi consiere as a funcion of he eh g is he acceleraion in he

3 Asmoical analsis of inernal gravi wave namics 9 gravi force fiel he sign - originaes from he increase of he ensi wih he ρ eh an herefore < [-3] The eac soluions of he essenial euaions escribing he inernal gravi waves are onl obaine for secial cases Tha is he reason wh he aroimae asmoical mehos are ssemaicall use for he invesigaion of he inernal gravi wave fiels in sraifie ocean The inernal gravi waves are usuall reresene in he following inegral form: J e λf F λ >> where f an F are analic funcions of γ [ ] he comle variable ; γ is a conour of inegraion on he comle lane The universal wa o consruc he asmoic forms of such eressions is he meho of ealon inegrals [3-7] This aer is evoe o he ssemaical escriion of a generaliaion of he geomerical oics meho KBJ meho ie we iscuss he saioemoral ra meho of ealon funcions This meho allows one o solve he roblem of asmoic moeling of he inharmonic wave acke s namics for he inernal gravi waves in sraifie meia wih slowl varing arameers The main reasons o use he ra mehos are he following: he ra reresenaions are well correlae wih he inuiion an wih he emirical maerial for he roagaion of he inernal gravi waves in naural sraifie meia ocean amoshere These mehos are universal an ver ofen one can use onl hem for he aroimae comuaions of he wave fiels in slowl changing nonhomogeneous sraifie meia [367] As is well known an essenial influence on he roagaion of inernal gravi waves in sraifie naural meia ocean amoshere is cause b he horional inhomogenei an non-saionari of hese meia To he mos ical horional inhomogenei of a real ocean one can refer he moificaion of he relief of he boom an inhomogenei of he ensi fiel an he variabili of he mean flows One can obain an eac analic soluion of his roblem for insance b using he meho of searaion of variables onl if he isribuion of ensi an he shae of he boom are escribe b raher simle moel funcions If he shae of he boom an he sraificaion are arbirar hen one can consruc onl asmoic reresenaions of he soluion in he near an far ones; however o escribe he fiel of inernal waves beween hese ones one nees an accurae numerical soluion of he roblem [-3] Using asmoic mehos one can consier a wie class of ineresing hsical roblems incluing roblems concerning he roagaion of nonharmonic wave ackes of inernal gravi waves in iverse nonhomogeneous sraifie meia uner he assumion ha he moificaion of he arameers of a vericall sraifie meium are slow in he horional irecion From he general oin of view roblems of his kin can be suie in he framework of a combinaion of he aiabaic an semi-classical aroimaions or b using close aroach for eamle ra eansions In aricular he asmoic soluions of

4 V V Bulaov an Yu V Vlaimirov iverse namical roblems can be escribe b using he Maslov canonical oeraor which eermines he asmoic behavior of he soluions incluing he case of neighborhoos of singular ses comose of focal oins causics ec [347] The secific form of he wave acke can be finall eresse b using some secial funcions sa in erms of oscillaing eonenials Air funcion Fresnel inegral Pearce-e inegrals ec The above aroaches are uie general an in rincile enable one o solve a broa secrum of roblems from he mahemaical oin of view; however he roblem of heir racical alicaions an in aricular of he visualiaion of he corresoning asmoic formulas base on he Maslov canonical oeraor is sill far from comleion an in some secific roblems one can use oher schemes o fin he asmoic behavior whose comuer realiaion using sofware of Mahemaica e is raher simle In his aer using he aroaches eveloe in [5] we consruc asmoic soluions of he roblem which is formulae as follows Problem formulaion: horionall inhomogeneous vericall sraifie meium This aer covers he conceual issues of he saio-emoral meho he geomerical oics meho he KBJ aroimaion aking ino consieraion he secifics of inernal gravi waves If we eamine he inernal gravi waves for he case when he unisurbe ensi fiel ρ eens no onl on he eh bu on he horional coorinaes an hen in general erms if he unisurbe ensi is a funcion of horional coorinaes such a isribuion of ensi inuces a fiel of horional flows These flows however are eremel slow an in he firs aroimaion can be neglece o i is commonl suose ha he fiel ρ is efine a riori hus i is assume ha here eis cerain eernal sources or he eamine ssem is non-conservaive I is also evien ha if he inernal gravi waves are roagaing above an non-uniform boom here is no such a roblem because he "inernal wave non-uniform boom" ssem is conservaive an here is no eernal energ flush [-3] e we shall wrie ou he linearie ssem of euaions of hronamics [-35] ~ ~ U U ρ ρ ~ ρ gρ ~ ~ ~ U U ρ ~ ρ ~ ρ ~ ρ U U

5 Asmoical analsis of inernal gravi wave namics ~ ~ ~ Here U U is he veloci vecor of inernal gravi waves an ρ are he ressure an ensi erurbaions g is he acceleraion of gravi - ais is irece ownwars Using he Boussines aroimaion which means he ensi ρ in he firs hree euaions of he ssem is assume a consan value he ssem b aling he cross-iffereniaing will be given as 4 ~ ~ g ~ ρ ~ ρ ~ ρ Δ Δ U U ρ ~ ~ ΔU ~ ~ ΔU Δ / / As he bounar coniions we ake he "rigi-li" coniion a H 3 ~ ~ ~ Consier he harmonic waves U U e i U U Inrouce he non-imensional variable accoring o he formulas: where L is he ical scale of he L L h horional variaions ρ h is he ical scale of he verical variaions ρ for eamle he hermocline wih [35] In non-imension coorinaes he euaion ssem will be wrien as ine is omie hereafer g ρ ρ ρ ε Δ ε εu εu ρ ε ΔU ε ΔU 4 h g where ε << g L h Asmoic soluions The asmoic soluion 4 shall be foun in he form usual for he geomeric oics meho [378-] V iε m Vm e 5 m iε V U U

6 V V Bulaov an Yu V Vlaimirov Funcions an V m m are subjec o efiniion From here on we shall resric ourselves o fining onl he ominan member of he eansion 5 for he verical veloci comonen a ha from he las wo euaions 4 we have i / U 6 U i / ubsiue 5 ino he firs euaion of he ssem 4 an se eual he members of he orer O 7 H g ρ where is he Vaisala-Brun freuenc eening of he ρ horional coorinaes The Vaisala-Brun freuenc buoanc freuenc is a main arameer han eermines inernal gravi wave roeries in real ocean [- 3] The bounar roblem 7 has a calculaion seu of eigenfuncions n an eigenvalues Kn n which are assume o be known [5] From here on he ine n will be omie while assuming ha furher calculaions belong o an iniviuall aken moe For he funcion we have he eikonal euaion [48-] K 8 Iniial coniions for he eikonal for he horional case are efine on he line L : : [8-] For solving he eikonal euaion we L consruc he ras ha is he euaion 8 wih characerisics K 9 K K K where / / is he lengh elemen of he ra The iniial coniions an shall be efine from he ssem

7 Asmoical analsis of inernal gravi wave namics 3 K The euaions 9 an iniial coniions efine he ra Afer he ras are foun he eikonal is efine b inegraion along he ra [569] K ow we are going o fin he eigenfuncion oe ha from 7 we can eermine onl he verical eenenc of he funcion In oher wors he funcion is efine o he accurac of mulilicaion b he arbirar funcion A e shall fin given as ˆ A where ˆ is he soluion of he roblem 7 normalie as follows H ˆ Afer subsiuing 5 ino 4 se eual he members of he orer ε O K g Δ ρ ρ 3 H e we use he orhogonali coniion of he righ-han member of he euaion 3 wih resec o funcion Muliling 3 b an inegraing over from o H we obain H H H g ln ρ 4 Conver he secon erm ino 4 using he inegraion b ars an ero bounar coniions for H H g ln ρ 5 Conver he firs erm ino 4 while accouning for H H A 6 In orer o ransform he hir member ino 4 we al he graien oeraor o he euaion 7 seing Y

8 4 V V Bulaov an Yu V Vlaimirov Y Y K K 7 Muliling 7 b an inegraing over from o H while accouning for we ge ln K A H 8 Finall we wrie over he euaion 4 using 5 6 an 8 ln 3 Δ K A A 9 The ransfer euaion 9 will be solve in characerisics of he eikonal euaion 9 Using he formula for Δ along he ras: JK J Δ where J is he geomeric ra ivergence we reuce he ransfer euaion 9 o he following conservaion law along he ras [38-] ln K J A The conservaion law can be wrien as well in he form suiable for fining he funcion A [58-] a K A a K A where J a is he uni ra ube wih oe ha he wave energ flash is roorional o a K A hus from i follows ha in his case here survives he value eual o he wave energ flash ivie b he wave vecor moulus [ 57] on-harmonic inernal gravi wave ackes To rocee o suing he roblem of non-harmonic wave ackes evoluion in a smoohl non-uniform horionall sraifie meium we resuose he choice of Ana Ana is he German a soluion e [34] which efine he roagaion of Air an Fresnel inernal waves wih cerain heurisic argumens

9 Asmoical analsis of inernal gravi wave namics 5 Air wave Le s inrouce he slow variables ε ε ε no slowness is suose over he ine is omie hereafer where ε λ / L << is he small arameer ha characeries he sofness of ambien horional changes λ is he ical wave lengh L is he scale of a horional non-uniformi Then he ssem in slow variables will be wrien as follows: 4 g ρ ρ ρ ε Δ εu εu ρ ε ΔU ε ΔU e we eamine he suerimosiion of harmonic waves in slow variables m i iε m e [ m ] 3 m ε ih resec o funcions m i is assume ha he are o-numbere on an min / is reache a for all an ubsiuing 3 ino we can easil have i rove ha he funcion m has a a ole of he m-h orer Therefore as he moel inegral R m for iniviual erms in 3 will serve he following formulas: m 3 i R i m e π where he inegraion conour is going 3 aroun he oin from overhea which enables he funcions R m o eoneniall eca a >> [-] The funcions R m have he following feaure: Rm R m a ha R Ai R Ai R Ai u u ec I is evien consiering cerain roeries of Air inegrals ha he funcions R m relae wih each oher as: R R R 3 R R [] Fresnel wave As he moel inegrals R m ha escribe he roagaion of Fresnel waves aking ino accoun he soluion srucure for he islacemen in he horionall uniform case we use he following formulas: * R Re e i i ReΦ Φ Is eas o see ha funcion * Φ has he following feaure * Φ i e i i

10 6 V V Bulaov an Yu V Vlaimirov e e * Φ i i i i i i i i From which we can obain for insance: i i * * 3 * 3 Φ Φ Φ or in erms of funcions m R : R i R 3 R i R i R [] Base on he above аn as well on he firs member srucure of he Air an Fresnel uniform wave asmoics for a horionall uniform meium he soluion of he ssem in can be foun for insance in he form for an iniviuall aken moe n n U furher omiing he ine n [5] K ε ε ε R R R a a 4 K 3 ε ε ε R R R a a U U U U where he argumen a a a ε is assume o be of he orer of uni Eansion Eression 4 agrees wih a common aroach of he geomeric oics meho saial ime ra meho KBJ meho [-3367] oe also ha from such a soluion srucure i follows ha he soluion for a horionall non-uniform an non-saionar meium shall een on boh he "fas" verical coorinae an "slow" ime an horional coorinaes variables e we generall are going o fin a soluion in "slow" variables a ha he soluions srucural elemens which een on he "fas" variables aear in he form of inegrals of some slowl varing funcions along he sace-ime ras This soluion choice allows us o efine he uniform asmoics for inernal gravi wave fiels roagaing wihin sraifie laer wih slowl varing arameers which hols rue eiher near or far awa from he wave frons of a single wave moe If we nee onl o efine he behavior of a fiel near he wave fron hen we can use one of he geomeric oics mehos he "rogressing wave" meho an a weakl isersive aroimaion in he form of aroriae local asmoics an fin he reresenaion for he hase funcions argumen in he form: ε ε a where he funcion efines he wave fron osiion an is eermine from he eikonal euaion soluion: c [35] A funcion c is he maimum grou veloci of a resecive wave moe ie he firs member of he isersion curve eansion in ero [5] The funcion he secon member of he eane isersion curve escribes he sace-ime imulse wih evoluion of Air or Fresnel non-harmonic waves an hen i will be efine from some arbirar laws of conservaion along he eikonal euaion characerisics wih heir acual form o be eermine b he roblem hsical coniions [-3 4]

11 Asmoical analsis of inernal gravi wave namics 7 Eikonal euaion: ras an iniial coniions formulaion e we focus on he eikonal euaion an formulaion of iniial seings for i In his secion we shall eamine generall wo basic eikonal euaions an relaive characerisic ssems for solving hese euaions The firs e eikonal euaion is obaine b consrucing he asmoics wihin a limie sace omain ie near he wave fron of a single moe ha is o sa we use a "weak isersion" aroimaion [5] In his case he isersion curve having is roeries efine b he wave fiel near he fron is aroimae b aroriae Talor eansion limie b he hir member of he asmoic series Then as i will be illusrae in [5] o efine he momen of he wave fron arrival we have he following eikonal euaion c 5 where c is he maimum grou veloci of roagaing iniviual inernal gravi waves moe The soluion of he verical secrum roblem wih aroriae bounar coniions A A c is assume o be known [-3] Then for solving he euaion of his eikonal 5 he ssem of characerisic euaions will be wrien as [8-] c c c 6 c c c Hence we ge for eamle he relaionshi as follows: c c c 7 which inicaes ha for his case he eikonal is he linear funcion of ime As he iniial seings for solving he ssem 6 we ma use he aa on some line The roblem uner invesigaion is wo-imensional because i is assume ha he meium arameers are slowl varing in ime an horionall an no slowness in he verical coorinae is reuire Therefore le he lane have a aramericall efine line: l λ λ on which we can lo he eikonal iniial isribuion: λ The hsical meaning of seing he iniial coniions in such a form is raher evien: he horional inernal wave falls on he bounar beween horionall uniform an non-uniform meiums

12 8 V V Bulaov an Yu V Vlaimirov he inerface beween he meiums is efine b his line all wave fiel arameers a his borer line are known Each value of he arameer λ ienicall eermines is value λ an he Caresian orinaes λ λ from which he characerisics is rouce ha is he ssem 6 soluion In he eamine ssem 6 b reason of is efiniion he arameer ha is varing along he characerisics is he ime elece as he iniial value is he ime of ei of raiaion Le here be a oin on he line l from which here originaes a characerisics ra wih arameer λ for which i akes ime o reach an arbirar oin Evienl we ma assume ha hese values are he funcions of he arameers λ an ha is λ λ If he ransformaion Jacobian from ra coorinaes λ o he Caresian iffers from ero hen hese euaions are solvable wih resec o variables λ : λ λ The funcion a λ λ is given as λ an as he iniial value we ake he momen ime of ei for he characerisics ra Hence we have he ssem λ a λ λ Then using he relaion 7 we obain ha The secon e eikonal euaion resuls from consrucing he roorional asmoics of he far fiels of inernal gravi waves The eikonal euaion for eermining he hase of non-harmonic wave ackes ha escribe hese roorional asmoics shall be wrien as K 8 k / olving he verical secral roblem wih resecive bounar coniions A k A is also assume o be known The isersion eenenc enoe as K is efine from he soluion of his verical secral roblem is he secral arameer The euaion 8 is he Hamilon-Jacobi euaion wih Hamilonian oeraor k K [8-] an he characerisic ssem of his euaion is given b K K K K K K K K Consier he characerisic ssem as follows

13 Asmoical analsis of inernal gravi wave namics 9 9 where he freuenc an he ime of ei of raiaion are he ra coorinaes The soluions of his ssem are: The Jacobian of ransiion from he Caresian orinaes o ra coorinaes is wrien as: e we calculae: Using he characerisic ssem 9 we can obain Then we have e we show ha here is he relaionshi Hence he following euali will be realie ivc c 3

14 3 V V Bulaov an Yu V Vlaimirov c ivc In he case uner invesigaion he oeraor ivc accouns for he elici eenence of c c on he variables an he imlici eenence on hese variables since i is evien ha he freuenc is also a funcion of hese variables ow we calculae on he characerisics 9 he values of funcions iffereniaing he euaion of characerisics on variables we have conseuenl From which we obain he following relaionshis / / / / Then we have: / / B his means he relaionshi 3 has been rove I is evien ha in he sace such roving is no necessar whereas he euali ivc ln in he sace is evien b virue of he oeraor efiniion iv wihin his sace The reuiremen of roving he relaionshi 3 comes abou from he fac ha he characerisics ras are no viewe in he four-imension sace bu in he wo-imension Caresian sace [ 3 8-] Take a ifferen seing of he iniial aa for solving he characerisic ssem in 9 such as he iniial aa seing a he momen of ime on he aramericall efine line are he ra coorinaes

15 Asmoical analsis of inernal gravi wave namics 3 Then we have: Ω a On his line he eikonal is he well-known funcion Then he characerisics euaion can be wrien as Ω a 3 a For he eamine case we wrie ou a formula for he Jacobian ransiion from he Caresian orinaes o he ra coorinaes : e we calculae 3 Using he characerisic ssem 9 an iniial aa in he form 3 we obain he following relaionshis: Ω Ω Ω Ω Using hese relaionshis we can calculae Ω Ω Ω Ω Laer iffereniaing he iniial coniions 3 on resecivel an b reeaing he above algorihm we can ge / / Ω Ω Ω Ω / / Ω Ω Hence resuling from 3 we also come o he formula 3 e have o make a few oins wih regar o he Jacobian geomeric meaning of he iscusse roblem in relaion o non-harmonic ackes of inernal gravi waves generae b isurbing sources in moion [7] Consier a characerisic ssem for he horionall uniform saionar scenario an a oin isurbing source in moion wih veloci V The euaions of characerisics ras are efining he sraigh lines

16 3 V V Bulaov an Yu V Vlaimirov V VK K ν K ν K K K V The Jacobian ransiion from he ra coorinaes o Caresian orinaes is given b ν ν V VK K K K VK K VK K / Fig The geomerical meaning of he Jacobian L is he source moion rajecor wihin he Caresian sace; А is he characerisic ra wih a consan value of ra coorinaes ; C is he characerisic ra wih a consan value of ra coorinaes Δ ; B is he characerisic ra wih a consan value of ra coorinaes Δ ; is he characerisic ra wih a consan value of ra coorinaes Δ Δ Then on accoun of he geomeric consrucion shown in Fig i is evien ha in he firs aroimaion he area Σ of he surface elemen ABC euals: Σ Δ Δ Here lines AC an B are he consan value lines ha is he conour lines of ra variables an in he Caresian sace I is evien ha in non-isersive meiums because of no ivergence of characerisics ras he area sie Σ of he surface elemen ABC b virue of is efiniion is eual o ero If aar from he isersion he sraifie meium is

17 Asmoical analsis of inernal gravi wave namics 33 also horionall non-uniform hen he characerisics ras are no sraigh lines he surface elemen ABC onl in he firs aroimaion an a uie small Δ an Δ is aroimae b a resecive arallelogram For he case of a horionall uniform meium he characerisics ras are alwas sraigh lines an he elemen ABC is alwas a arallelogram Thus he Jacobian of ra coorinaes ransiing o he Caresian orinaes oulines he geomeric ivergence in he Caresian sace of he characerisics ras for he resecive eikonal euaion [38-] ow we go back o he wo rincial mehos of efining he iniial seings for solving he characerisic ssem of an arbirar eikonal euaion: F / which eermines he hase funcion [358-] The firs meho consiss in ha we efine an arbirar iniial sace-ime omain : λ λ λ where λ are some ra variables in which all iniial aa is efine incluing he hase funcion: λ The characerisics ras are «le ou» of his sace-ime omain which is he beginning of furher namics for he characerisics ras I is evien ha he following relaionshis are holing for his omain λ λ λ λ λ λ λ Using he eikonal euaion for he iniial values F 33 we obain hree euaions o eermine he remaining unknown iniial values wihin his iniial sace-ime omain Us a rule he hsical saemen of a roblem relae o non-harmonic ackes of inernal gravi waves generae in real-worl naural environmens involves he following There are wo meiums: a horionall uniform an a horionall non-uniform which have a cerain inerface borer line in beween he form of which is well-known [5] On his borer line base he roblem soluion for a uniform meium we assume as known all roeries of he eamine wave fiel ha is he hase funcion is consiere as known for all momens of he ime Accouning for ha we have o efine he evoluion of non-harmonic wave ackes in a horionall non-uniform meium To his en i seems more convenien o se he iniial aa for solving a ssem of characerisics in oher forma ha is in he form of arbirar "bounar" coniions as efine on his borer line beween wo meiums ha his means is ha we efine he line a over which he hase funcion is well-known a an momen of ime : Is

18 34 V V Bulaov an Yu V Vlaimirov worh o menion ha he bounar line beween wo meiums is no ime eenen On his line evienl we can b virue of efiniion realie he following ar: Then o eermine he iniial values on his line we have which ogeher wih he relaionshi 33 allows us o efine he iniial aa values for hese wo variables In wha follows we shall al his meho o formulae he iniial aa for solving a characerisic ssem [578-] Asmoical analsis of sraifie non-saionar meium wave namics Uner he real oceanic coniions he Vaisala-Brun freuenc g ln ρ / efines he basic characerisics of inernal gravi wave namic an i shall no een solel on saial variables bu also on he ime The mos characerisic es of ime-o-ime variabili are he hermocline going u or own an changing is wih ec [479] There is a number of ime scales for variaions of hro-hsical fiels in he oceans an seas: a small-scale wih erios of abou minues a meso-scale wih erios of abou a a wen-four hours аs well as snoical an global variaions wih erios of a few monhs o a few ears [5] In wha follows we shall anale he inernal gravi fiel roagaion in non-saionar meiums wih arameer variaion erios of a a an over which allows us o use he geomeric oics aroimaion because he erio of inernal gravi waves is ens of minues an less The ssem of linearie euaions of hronamics when he non-erurbae ensi ρ eens on variables an reuces o a single euaion for eamle for he verical veloci comonen: ln ρ ln ρ ln ρ g Δ Δ Δ If ln ρ / is neglece we obain an euaion in he Boussines aroimaion: ln ρ Δ Δ I aears naural o neglec as well he member wih ln ρ / which woul correson o a conseuen alicaion of he Boussines hohesis [-3] I means he ensi characeriing he fuis s iner mass can be assume consan Then we have: Δ Δ The resuling euaion iffers from a sanar

19 Asmoical analsis of inernal gravi wave namics 35 euaion of inernal gravi waves in a saionar sraifie meium jus b he ime aramerical inclusion ino he Vaisala-Brun freuenc The asmoic soluion is foun in he form of a sum of moes wih ever one of hem roagaing ineenenl of each oher he aiabaic aroimaion e are going o eamine a single iniviuall aken moe while omiing is ine e we focus solel on he sace region near he wave fron which means ha we consier he ime as being close o he arrival ime of he wave fron henceforh enoe b τ ie we use a weakl isersive aroimaion Consier he wave roagaion in a laer of sraifie meium H < < wih he Vaisala-Brun freuenc e shall seek he soluion wih bounar coniions H in he form O ε A ε ε τ A ε ε τ ε τ L F ϕ τ F m ϕ Fm ϕ τ τ ε ε B ε ε τ B ε ε τ ε τ L F ϕ τ where / 3 F ϕ Ai is he Air erivaive having is argumen ϕ ε ε ε τ ε ε of he orer of uni The funcion τ efines he ϕ ε wave fron osiion funcion escribes he evoluion of he Air wave wih he small arameer ε secifies "slow variables" ince our focus is onl on slow imes ε being close o he ime of he wave fron arrival τ hen all funcions receing funcions F m are given in he form of Talor series bε τ ε owers Le ε be wrien as: τ ε τ ε τ O ε e consier he far inernal τ gravi waves generae b moving source Afer comlicae analical calculaions we can obain an eression for he erm in he form: 5 / c τ f τ f τ τ Ai 3 / / / 4 / 3 c τ R V c τ ε where c f - eigenvalue an eigenfuncion of inernal gravi verical secral roblem V- source see - eh of source moion R - some funcion ha are eermine b he arameers of he roblem [57] The figures emonsrae he numerical resuls of inernal gravi wave calculaions for ical oceanic arameers [569] The Fig shows a ssem of ras hin line causics bol line generae b source moving in a nonsaionar sraifie ocean I is a general rule ha causic of a famil of ras single ou an area in sace so ha ras of ha famil canno aear in he marke area There is also anoher area an each oin of ha area has wo ras ha ass hrough his oin One of hose ras has alrea asse his oin an anoher is

20 36 V V Bulaov an Yu V Vlaimirov going o ass he oin Formal aroimaion of geomerical oics or KBJ aroimaion canno be alie near he causic ha is because ras merge ogeher in ha area afer he were reflece b causic If we wan o fin wave fiel near he causic hen i is necessar o use secial aroimaion of he soluion an in he aer a moifie ra meho is roose in orer o buil uniform asmoic eansion of inegral forms of he inernal gravi wave fiel Afer he ras are reflece b he causic here aears a hase shif I is clear ha he hase shif can onl haen in he area where mehos of geomerical oics which were use in revious secions can be alie If he ras ouch he causic several imes hen aiional hase shifs will be ae Phase shif which was creae b he causic is raher small in comarison wih he change in hase along he ra bu his shif can consierabl affec inerference aern of he wave fiel The Fig3 emonsraes an evoluion of inernal gravi wave acke in a nonsaionar sraifie ocean wihin he ssem of coorinaes ha is in moion ogeher wih he isurbing source Time inerval for calculaions is eual hours A ha i s evien ha if here were no Vaisala-Brun ime-o-ime freuenc variaions such a wave coorinae woul be saionar umerical resuls show ha inernal gravi waves namic in he real ocean is subsaniall influence b non-saionari of hro-hsical fiels The obaine asmoic soluions are uniform an allow far inernal gravi wave fiels o be escribe boh near an far from urning oins an wave frons Fig Ras an causics in sraifie non-saionar meium

21 Asmoical analsis of inernal gravi wave namics 37 Fig 3 Evoluion of inernal gravi wave acke in non-saionar sraifie meium iscussion Thus we have he following scheme for far inernal gravi waves asmoic calculaing in heerogeneous an non-saionar sraifie meium: for an arbirar ensi isribuion inernal gravi waves verical secral roblems is solve numericall b a shooing meho an he corresoning normalie eigenfuncions an eigenvalues are obaine [5] eikonal euaion is solve numericall along hese characerisics ras wih aroriae iniial coniions [3-5] 3 afer eerminaions of characerisics ras eikonal hase value hase funcions is calculae b numerical inegraion along hese ras [-3] 4 geomeric ivergence ra ubes is eermine for eamle numerical iffereniaion of closel sace ras [3 5-7] 5 inernal gravi waves amliue is calculae from he corresoning conservaion laws energ conservaion along he ras characerisics [ ] In his aer a moifie sace-ime ra meho is roose which belongs o he class of geomerical oics mehos KBJ aroimaion The ke oin of he roose echniue is he ossibili o erive he asmoic reresenaion of he soluion in erms of a non-ineger ower series of he small arameer

22 38 V V Bulaov an Yu V Vlaimirov ε Λ / L where Λ is he characerisic wave lengh an L is he characerisic scale of he horional heerogenei The elici form of he asmoic soluion was eermine base on he rinciles of locali an asmoic behavior of he soluion in he case of a saionar an horionall homogeneous meium The wave acke amliues are eermine from he energ conservaion laws along he characerisic curves A ical assumion mae in suies on he inernal wave evoluion in sraifie meia is ha he wave ackes are locall harmonic A moificaion of he geomerical oics meho base on an eansion of he soluion in moel funcions allows one o escribe he wave fiel srucure boh far from an a he vicini of he wave fron The universal characer of he asmoic meho roose for moeling inernal gravi fiels makes i ossible o effecivel calculae wave fiels an in aiion ualiaivel anale he obaine soluions This meho offers broa ooruniies for he analsis of wave fiels on a large scale which is imoran for eveloing correc mahemaical moels of wave namics an for assessing in siu measuremens of wave fiels in he ocean The aricular role of he roose asmoic mehos is eermine b he fac ha he arameers of naural sraifie meia are usuall known aroimael an aems a heir aeuae numerical soluion using he iniial euaions of hronamics an such arameers ma resul in a noable loss of accurac for he resuls obaine In aiion o heir funamenal significance he obaine asmoic moels are also imoran for alie invesigaions since he roose meho of geomerical oics allows soluion of a wie secrum of roblems relae o moeling wave fiels In such a siuaion he escriion an analsis of wave namics ma be realie hrough eveloing asmoic moels an using analical mehos for heir soluion base on he roose KBJ moifie meho The resuls of his work reresen a significan ineres for hsics mahemaics an engineers Besies ha ineres analical asmoic an numerical soluions which were obaine in his aer can resen significan imorance for engineering alicaions since resene meho which were o calculae he inernal gravi waves fiel make i ossible o calculae ifferen wave fiels in he raher big class of anoher roblems Acknowlegmens The resuls resene in he aer have been obaine b research erforme uner rojecs suore b he Russian Founaion for Basic Research o Program of he Russian Acaem of ciences Funamenal Problems of Oceanolog: Phsics Geolog Biolog Ecolog

23 Asmoical analsis of inernal gravi wave namics 39 References [] VV Bulaov YuV Vlaimirov ave namics of sraifie meiums Moscow auka Publishers [] J Pelosk aves in he ocean an amoshere: inroucion o wave namics Berlin- Heielberg ringer [3] BR uherlan Inernal gravi waves Cambrige Cambrige Universi Press [4] C Garre E Kune Inernal ie generaion in he ee ocean Rev Flui Mech [5] J Grue JKveen A scaling law of inernal run-u uraion Ocean namics [6] ZJong BYGou L Lua ZM hi YXiao Y Qu Comarisons of inernal soliar wave an surface wave acions on marine srucures an heir resonses Alie Ocean Research 33-9 [7] MK Hsu AK Liu C Liu A su of inernal waves in he China eas an Yellow ea using AR Coninenal helf Research [8] J Grue A Jensen Orbial veloci an breaking in see ranom gravi waves J Geohs Res 7 C7-C3 [9] J Grue RH Gabrielsen Es Marine ransor in he High orh Oslo ovus Forlag e orske Vienskas-Akaemi orges ekniske Vienskasakaemi [] J oong BC Lee CC Kao ave measuremen using GP veloci signal ensors [] Yang-Yih Chen GYChen Chia-Hao Lin Chiu-Long Chou Progressive waves in real fluis over a rigi ermeable boom Coas Eng J [3] VM Babich V Bulrev Asmoic mehos in shor-wavelengh iffracion heor Ofor Alha cience 7 [4] VIArnol Caasrohe heor Berlin Heielberg ringer 99

24 4 V V Bulaov an Yu V Vlaimirov [5] VVBulaov YuV Vlaimirov The uniform asmoic form of he inernal gravi wave fiel generae b a source moving above a smoohl varing boom J Eng Mah [6] J Asen J Vannese Ineria-gravi-wave generaion: a geomeric-oics aroach II/IUTAM roceeings 9 [7] Yu obrokhoov A Lohnikov CA Vargas Asmoics of waves on he shallow waer generae b saiall-localie sources an rae b unerwaer riges Rus J Mah Phs 3-4 [8] RB hie Asmoic analsis of ifferenial euaions Lonon: Imerial College Press 5 [9] FC Hoensea Quasi-sea sae analsis of ifferenial ifference inegral an graien ssems Couran Lecure oes Amer Mah oc [] JAaners F Verhuls JMurock Averaging mehos in nonlinear namical ssems n e ew York ringer-verlag 7 [] M Abramowi IA egun Hanbook of mahemaical funcions Rerin of he 97 e ew York over Publicaions Inc 99 [] G ason A reaise on he heor of Bessel funcions Rerin of he n e Cambrige Cambrige Universi Press 995 Receive: ovember 3

Lecture 6 - Testing Restrictions on the Disturbance Process (References Sections 2.7 and 2.10, Hayashi)

Lecture 6 - Testing Restrictions on the Disturbance Process (References Sections 2.7 and 2.10, Hayashi) Lecure 6 - esing Resricions on he Disurbance Process (References Secions 2.7 an 2.0, Hayashi) We have eveloe sufficien coniions for he consisency an asymoic normaliy of he OLS esimaor ha allow for coniionally

More information

Check in: 1 If m = 2(x + 1) and n = find y when. b y = 2m n 2

Check in: 1 If m = 2(x + 1) and n = find y when. b y = 2m n 2 7 Parameric equaions This chaer will show ou how o skech curves using heir arameric equaions conver arameric equaions o Caresian equaions find oins of inersecion of curves and lines using arameric equaions

More information

Design flood calculation for ungauged small basins in China

Design flood calculation for ungauged small basins in China 236 Preicions in Ungauge Basins: PUB Kick-off (Proceeings of he PUB Kick-off meeing hel in Brasilia, 20 22 November 2002). IAHS Publ. 309, 2007. Design floo calculaion for ungauge small basins in China

More information

Chapter Three Systems of Linear Differential Equations

Chapter Three Systems of Linear Differential Equations Chaper Three Sysems of Linear Differenial Equaions In his chaper we are going o consier sysems of firs orer orinary ifferenial equaions. These are sysems of he form x a x a x a n x n x a x a x a n x n

More information

Chapter 3 Common Families of Distributions

Chapter 3 Common Families of Distributions Chaer 3 Common Families of Disribuions Secion 31 - Inroducion Purose of his Chaer: Caalog many of common saisical disribuions (families of disribuions ha are indeed by one or more arameers) Wha we should

More information

Suyash Narayan Mishra, Piyush Kumar Tripathi & Alok Agrawal

Suyash Narayan Mishra, Piyush Kumar Tripathi & Alok Agrawal IOSR Journal o Mahemaics IOSR-JM e-issn: 78-578 -ISSN: 39-765X. Volume Issue Ver. VI Mar - Ar. 5 PP 43-5 www.iosrjournals.org A auberian heorem or C α β- Convergence o Cesaro Means o Orer o Funcions Suash

More information

Chapter 2 The Derivative Applied Calculus 107. We ll need a rule for finding the derivative of a product so we don t have to multiply everything out.

Chapter 2 The Derivative Applied Calculus 107. We ll need a rule for finding the derivative of a product so we don t have to multiply everything out. Chaper The Derivaive Applie Calculus 107 Secion 4: Prouc an Quoien Rules The basic rules will le us ackle simple funcions. Bu wha happens if we nee he erivaive of a combinaion of hese funcions? Eample

More information

Homework 2 Solutions

Homework 2 Solutions Mah 308 Differenial Equaions Fall 2002 & 2. See he las page. Hoework 2 Soluions 3a). Newon s secon law of oion says ha a = F, an we know a =, so we have = F. One par of he force is graviy, g. However,

More information

Module II, Part C. More Insight into Fiber Dispersion

Module II, Part C. More Insight into Fiber Dispersion Moule II Par C More Insigh ino Fiber Dispersion . Polariaion Moe Dispersion Fiber Birefringence: Imperfec cylinrical symmery leas o wha is known as birefringence. Recall he HE moe an is E x componen which

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

The Fundamental Theorems of Calculus

The Fundamental Theorems of Calculus FunamenalTheorems.nb 1 The Funamenal Theorems of Calculus You have now been inrouce o he wo main branches of calculus: ifferenial calculus (which we inrouce wih he angen line problem) an inegral calculus

More information

( ) ( ) ( ) ( u) ( u) = are shown in Figure =, it is reasonable to speculate that. = cos u ) and the inside function ( ( t) du

( ) ( ) ( ) ( u) ( u) = are shown in Figure =, it is reasonable to speculate that. = cos u ) and the inside function ( ( t) du Porlan Communiy College MTH 51 Lab Manual The Chain Rule Aciviy 38 The funcions f ( = sin ( an k( sin( 3 38.1. Since f ( cos( k ( = cos( 3. Bu his woul imply ha k ( f ( = are shown in Figure =, i is reasonable

More information

Lecture 2: Telegrapher Equations For Transmission Lines. Power Flow.

Lecture 2: Telegrapher Equations For Transmission Lines. Power Flow. Whies, EE 481/581 Lecure 2 Page 1 of 13 Lecure 2: Telegraher Equaions For Transmission Lines. Power Flow. Microsri is one mehod for making elecrical connecions in a microwae circui. I is consruced wih

More information

Math 2214 Sol Test 2B Spring 2015

Math 2214 Sol Test 2B Spring 2015 Mah 14 Sol Tes B Sring 015 roblem 1: An objec weighing ounds sreches a verical sring 8 fee beond i naural lengh before coming o res a equilibrium The objec is ushed u 6 fee from i s equilibrium osiion

More information

Physics 218 Exam 1. with Solutions Fall 2010, Sections Part 1 (15) Part 2 (20) Part 3 (20) Part 4 (20) Bonus (5)

Physics 218 Exam 1. with Solutions Fall 2010, Sections Part 1 (15) Part 2 (20) Part 3 (20) Part 4 (20) Bonus (5) Physics 18 Exam 1 wih Soluions Fall 1, Secions 51-54 Fill ou he informaion below bu o no open he exam unil insruce o o so! Name Signaure Suen ID E-mail Secion # ules of he exam: 1. You have he full class

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

Seminar 5 Sustainability

Seminar 5 Sustainability Seminar 5 Susainabiliy Soluions Quesion : Hyperbolic Discouning -. Suppose a faher inheris a family forune of 0 million NOK an he wans o use some of i for himself (o be precise, he share ) bu also o beques

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

BOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at

BOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at BOX-JENKINS MODEL NOAION he Box-Jenkins ARMA(,q) model is denoed b he equaion + + L+ + a θ a L θ a 0 q q. () he auoregressive (AR) ar of he model is + L+ while he moving average (MA) ar of he model is

More information

D.I. Survival models and copulas

D.I. Survival models and copulas D- D. SURVIVAL COPULA D.I. Survival moels an copulas Definiions, relaionships wih mulivariae survival isribuion funcions an relaionships beween copulas an survival copulas. D.II. Fraily moels Use of a

More information

Dirac s hole theory and the Pauli principle: clearing up the confusion.

Dirac s hole theory and the Pauli principle: clearing up the confusion. Dirac s hole heory and he Pauli rincile: clearing u he conusion. Dan Solomon Rauland-Borg Cororaion 8 W. Cenral Road Moun Prosec IL 656 USA Email: dan.solomon@rauland.com Absrac. In Dirac s hole heory

More information

VS203B Lecture Notes Spring, 2011 Topic: Diffraction

VS203B Lecture Notes Spring, 2011 Topic: Diffraction Diffracion Diffracion escribes he enency for ligh o ben aroun corners. Huygens principle All poins on a wavefron can be consiere as poin sources for he proucion of seconary waveles, an a a laer ime he

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

Integral representations and new generating functions of Chebyshev polynomials

Integral representations and new generating functions of Chebyshev polynomials Inegral represenaions an new generaing funcions of Chebyshev polynomials Clemene Cesarano Faculy of Engineering, Inernaional Telemaic Universiy UNINETTUNO Corso Viorio Emanuele II, 39 186 Roma, Ialy email:

More information

Computer Vision. Motion Extraction

Computer Vision. Motion Extraction Comuer Moion Eracion Comuer Alicaions of moion eracion Change / sho cu deecion Surveillance / raffic monioring Moion caure / gesure analsis HC image sabilisaion Moion comensaion e.g. medical roboics Feaure

More information

Vector autoregression VAR. Case 1

Vector autoregression VAR. Case 1 Vecor auoregression VAR So far we have focused mosl on models where deends onl on as. More generall we migh wan o consider oin models ha involve more han one variable. There are wo reasons: Firs, we migh

More information

Rise-Time Distortion of Signal without Carrying Signal

Rise-Time Distortion of Signal without Carrying Signal Journal of Physics: Conference Series PAPER OPEN ACCESS Rise-Time Disorion of Signal wihou Carrying Signal To cie his aricle: N S Bukhman 6 J. Phys.: Conf. Ser. 738 8 View he aricle online for udaes and

More information

total distance cov ered time int erval v = average speed (m/s)

total distance cov ered time int erval v = average speed (m/s) Physics Suy Noes Lesson 4 Linear Moion 1 Change an Moion a. A propery common o eeryhing in he unierse is change. b. Change is so imporan ha he funamenal concep of ime woul be meaningless wihou i. c. Since

More information

ES2A7 - Fluid Mechanics Example Classes Example Questions (Set IV)

ES2A7 - Fluid Mechanics Example Classes Example Questions (Set IV) ESA7 - Flui Mechanics Eample Classes Eample Quesions (Se IV) Quesion : Dimensional analysis a) I is observe ha he velociy V of a liqui leaving a nozzle epens upon he pressure rop P an is ensiy ρ. Show

More information

The Natural Logarithm

The Natural Logarithm The Naural Logarihm 5-4-007 The Power Rule says n = n + n+ + C provie ha n. The formula oes no apply o. An anierivaive F( of woul have o saisfy F( =. Bu he Funamenal Theorem implies ha if > 0, hen Thus,

More information

ON DETERMINATION OF SOME CHARACTERISTICS OF SEMI-MARKOV PROCESS FOR DIFFERENT DISTRIBUTIONS OF TRANSIENT PROBABILITIES ABSTRACT

ON DETERMINATION OF SOME CHARACTERISTICS OF SEMI-MARKOV PROCESS FOR DIFFERENT DISTRIBUTIONS OF TRANSIENT PROBABILITIES ABSTRACT Zajac, Budny ON DETERMINATION O SOME CHARACTERISTICS O SEMI MARKOV PROCESS OR DIERENT DISTRIBUTIONS O R&RATA # 2(3 ar 2 (Vol. 2 29, June ON DETERMINATION O SOME CHARACTERISTICS O SEMI-MARKOV PROCESS OR

More information

Practice Problems - Week #4 Higher-Order DEs, Applications Solutions

Practice Problems - Week #4 Higher-Order DEs, Applications Solutions Pracice Probles - Wee #4 Higher-Orer DEs, Applicaions Soluions 1. Solve he iniial value proble where y y = 0, y0 = 0, y 0 = 1, an y 0 =. r r = rr 1 = rr 1r + 1, so he general soluion is C 1 + C e x + C

More information

IB Physics Kinematics Worksheet

IB Physics Kinematics Worksheet IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?

More information

Section 5: Chain Rule

Section 5: Chain Rule Chaper The Derivaive Applie Calculus 11 Secion 5: Chain Rule There is one more ype of complicae funcion ha we will wan o know how o iffereniae: composiion. The Chain Rule will le us fin he erivaive of

More information

Advanced Integration Techniques: Integration by Parts We may differentiate the product of two functions by using the product rule:

Advanced Integration Techniques: Integration by Parts We may differentiate the product of two functions by using the product rule: Avance Inegraion Techniques: Inegraion by Pars We may iffereniae he prouc of wo funcions by using he prouc rule: x f(x)g(x) = f (x)g(x) + f(x)g (x). Unforunaely, fining an anierivaive of a prouc is no

More information

THEORY OF CUMULATIVE FUEL CONSUMPTION BY LPG POWERED CARS

THEORY OF CUMULATIVE FUEL CONSUMPTION BY LPG POWERED CARS Journal of KONES Powerrain an Transor, Vol. 22, No. 4 205 THEORY OF CUMULATIVE FUEL CONSUMPTION BY LPG POWERED CARS Lech Jerzy Sinik Wroclaw Universiy of Technology Faculy of Mechanical Engineering Wysianskiego

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

Sub Module 2.6. Measurement of transient temperature

Sub Module 2.6. Measurement of transient temperature Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,

More information

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11. 1 Mah 334 Tes 1 KEY Spring 21 Secion: 1 Insrucor: Sco Glasgow Daes: Ma 1 and 11. Do NOT wrie on his problem saemen bookle, excep for our indicaion of following he honor code jus below. No credi will be

More information

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion

More information

Classical Thermodynamics as a Consequence of Spacetime Geometry

Classical Thermodynamics as a Consequence of Spacetime Geometry Classical Thermoynamics as a Consequence of Spaceime Geomery Jay R. Yablon 9 Norhumberlan Drive Schenecay, New York 239-284 jyablon@nycap.rr.com January 6, 25 Absrac: Jus as Maxwell s magneic equaions

More information

Chapter 12: Velocity, acceleration, and forces

Chapter 12: Velocity, acceleration, and forces To Feel a Force Chaper Spring, Chaper : A. Saes of moion For moion on or near he surface of he earh, i is naural o measure moion wih respec o objecs fixed o he earh. The 4 hr. roaion of he earh has a measurable

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Lecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance

Lecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance Lecure 5 Time series: ECM Bernardina Algieri Deparmen Economics, Saisics and Finance Conens Time Series Modelling Coinegraion Error Correcion Model Two Seps, Engle-Granger procedure Error Correcion Model

More information

COMPETITIVE GROWTH MODEL

COMPETITIVE GROWTH MODEL COMPETITIVE GROWTH MODEL I Assumpions We are going o now solve he compeiive version of he opimal growh moel. Alhough he allocaions are he same as in he social planning problem, i will be useful o compare

More information

Topics in Combinatorial Optimization May 11, Lecture 22

Topics in Combinatorial Optimization May 11, Lecture 22 8.997 Topics in Combinaorial Opimizaion May, 004 Lecure Lecurer: Michel X. Goemans Scribe: Alanha Newman Muliflows an Disjoin Pahs Le G = (V,E) be a graph an le s,,s,,...s, V be erminals. Our goal is o

More information

Math Spring 2015 PRACTICE FINAL EXAM (modified from Math 2280 final exam, April 29, 2011)

Math Spring 2015 PRACTICE FINAL EXAM (modified from Math 2280 final exam, April 29, 2011) ame ID number Mah 8- Sring 5 PRACTICE FIAL EXAM (modified from Mah 8 final exam, Aril 9, ) This exam is closed-book and closed-noe You may use a scienific calculaor, bu no one which is caable of grahing

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

Introduction to Physical Oceanography Homework 5 - Solutions

Introduction to Physical Oceanography Homework 5 - Solutions Laure Zanna //5 Inroducion o Phsical Oceanograph Homework 5 - Soluions. Inerial oscillaions wih boom fricion non-selecive scale: The governing equaions for his problem are This ssem can be wrien as where

More information

15. Bicycle Wheel. Graph of height y (cm) above the axle against time t (s) over a 6-second interval. 15 bike wheel

15. Bicycle Wheel. Graph of height y (cm) above the axle against time t (s) over a 6-second interval. 15 bike wheel 15. Biccle Wheel The graph We moun a biccle wheel so ha i is free o roae in a verical plane. In fac, wha works easil is o pu an exension on one of he axles, and ge a suden o sand on one side and hold he

More information

CH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

CH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC CH.7. PLANE LINEAR ELASTICITY Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Overview Plane Linear Elasici Theor Plane Sress Simplifing Hpohesis Srain Field Consiuive Equaion Displacemen Field The Linear

More information

1 birth rate γ (number of births per time interval) 2 death rate δ proportional to size of population

1 birth rate γ (number of births per time interval) 2 death rate δ proportional to size of population Scienific Comuing I Module : Poulaion Modelling Coninuous Models Michael Bader Par I ODE Models Lehrsuhl Informaik V Winer 7/ Discree vs. Coniuous Models d d = F,,...) ) =? discree model: coninuous model:

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

Classical Thermodynamics Entropy Laws as a Consequence of Spacetime Geometry

Classical Thermodynamics Entropy Laws as a Consequence of Spacetime Geometry Classical Thermoynamics Enropy Laws as a Consequence of Spaceime Geomery Jay R. Yablon 9 Norhumberlan Drive Schenecay, New York 239-284 jyablon@nycap.rr.com January 7, 25 Absrac: Jus as Maxwell s magneic

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Second-Order Differential Equations

Second-Order Differential Equations WWW Problems and Soluions 3.1 Chaper 3 Second-Order Differenial Equaions Secion 3.1 Springs: Linear and Nonlinear Models www m Problem 3. (NonlinearSprings). A bod of mass m is aached o a wall b means

More information

Computing with diode model

Computing with diode model ECE 570 Session 5 C 752E Comuer Aided Engineering for negraed Circuis Comuing wih diode model Objecie: nroduce conces in numerical circui analsis Ouline: 1. Model of an examle circui wih a diode 2. Ouline

More information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

When to Sell an Asset Where Its Drift Drops from a High Value to a Smaller One

When to Sell an Asset Where Its Drift Drops from a High Value to a Smaller One American Journal of Operaions Research, 5, 5, 54-55 Publishe Online November 5 in SciRes. hp://www.scirp.org/journal/ajor hp://x.oi.org/.436/ajor.5.564 When o Sell an Asse Where Is Drif Drops from a High

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

1. (16 points) Answer the following derivative-related questions. dx tan sec x. dx tan u = du d. dx du tan u. du tan u d v.

1. (16 points) Answer the following derivative-related questions. dx tan sec x. dx tan u = du d. dx du tan u. du tan u d v. Exam #2 Soluions. (6 poins) Answer he following eriaie-relae quesions. (a) (8 poins) If y an sec x, fin. This is an applicaion of he chain rule in wo sages, in which we shall le u sec x, an sec x: an sec

More information

( ) = b n ( t) n " (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2.

( ) = b n ( t) n  (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2. Andrei Tokmakoff, MIT Deparmen of Chemisry, 3/14/007-6.4 PERTURBATION THEORY Given a Hamilonian H = H 0 + V where we know he eigenkes for H 0 : H 0 n = E n n, we can calculae he evoluion of he wavefuncion

More information

Exam I. Name. Answer: a. W B > W A if the volume of the ice cubes is greater than the volume of the water.

Exam I. Name. Answer: a. W B > W A if the volume of the ice cubes is greater than the volume of the water. Name Exam I 1) A hole is punched in a full milk caron, 10 cm below he op. Wha is he iniial veloci of ouflow? a. 1.4 m/s b. 2.0 m/s c. 2.8 m/s d. 3.9 m/s e. 2.8 m/s Answer: a 2) In a wind unnel he pressure

More information

On a problem of Graham By E. ERDŐS and E. SZEMERÉDI (Budapest) GRAHAM stated the following conjecture : Let p be a prime and a 1,..., ap p non-zero re

On a problem of Graham By E. ERDŐS and E. SZEMERÉDI (Budapest) GRAHAM stated the following conjecture : Let p be a prime and a 1,..., ap p non-zero re On a roblem of Graham By E. ERDŐS and E. SZEMERÉDI (Budaes) GRAHAM saed he following conjecure : Le be a rime and a 1,..., a non-zero residues (mod ). Assume ha if ' a i a i, ei=0 or 1 (no all e i=0) is

More information

) were both constant and we brought them from under the integral.

) were both constant and we brought them from under the integral. YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha

More information

Vector autoregression VAR

Vector autoregression VAR Vecor auoregression VAR So far we have focused mosly on models where y deends on as y. More generally we migh wan o consider models for more han on variable. If we only care abou forecasing one series

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

C4: QUESTIONS FROM PAST PAPERS - PARAMETRICS

C4: QUESTIONS FROM PAST PAPERS - PARAMETRICS C: QUESTIONS FROM PAST PAPERS - PARAMETRICS.. The iagram above shows a skech of he curve C wih parameric equaions = 5, = (9 ) The curve C cus he -ais a he poins A an B. (a) Fin he -coorinae a he poin A

More information

Linear Dynamic Models

Linear Dynamic Models Linear Dnamic Models and Forecasing Reference aricle: Ineracions beween he muliplier analsis and he principle of acceleraion Ouline. The sae space ssem as an approach o working wih ssems of difference

More information

THE IMPACT OF SEASONAL ADJUSTMENT ON TIME SERIES PREDICTION

THE IMPACT OF SEASONAL ADJUSTMENT ON TIME SERIES PREDICTION Review of he Air Force Acaemy No 3 (27) 2014 THE IMPACT OF SEASONAL ADJUSTMENT ON TIME SERIES PREDICTION Aela SASU Transilvania Universiy of Braşov Absrac: In his aricle we inen o comare he erformance

More information

Section 2.6 Derivatives of products and quotients

Section 2.6 Derivatives of products and quotients Secion 2.6 Derivaives of proucs an quoiens (3/19/08) Overview: In his secion, we erive formulas for erivaives of funcions ha are consruce by aking proucs an quoiens of oher funcions, an we use hese formulas

More information

And the solution to the PDE problem must be of the form Π 1

And the solution to the PDE problem must be of the form Π 1 5. Self-Similar Soluions b Dimensional Analsis Consider he diffusion problem from las secion, wih poinwise release (Ref: Bluman & Cole, 2.3): c = D 2 c x + Q 0δ(x)δ() 2 c(x,0) = 0, c(±,) = 0 Iniial release

More information

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations Concourse Mah 80 Spring 0 Worked Examples: Marix Mehods for Solving Sysems of s Order Linear Differenial Equaions The Main Idea: Given a sysem of s order linear differenial equaions d x d Ax wih iniial

More information

Vorticity equation 2. Why did Charney call it PV?

Vorticity equation 2. Why did Charney call it PV? Vorici eqaion Wh i Charne call i PV? The Vorici Eqaion Wan o nersan he rocesses ha roce changes in orici. So erie an eression ha incles he ime eriaie o orici: Sm o orces in irecion Recall ha he momenm

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

This is an example to show you how SMath can calculate the movement of kinematic mechanisms.

This is an example to show you how SMath can calculate the movement of kinematic mechanisms. Dec :5:6 - Kinemaics model of Simple Arm.sm This file is provided for educaional purposes as guidance for he use of he sofware ool. I is no guaraeed o be free from errors or ommissions. The mehods and

More information

t 2 B F x,t n dsdt t u x,t dxdt

t 2 B F x,t n dsdt t u x,t dxdt Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.

More information

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8. Kinemaics Vocabulary Kinemaics and One Dimensional Moion 8.1 WD1 Kinema means movemen Mahemaical descripion of moion Posiion Time Inerval Displacemen Velociy; absolue value: speed Acceleraion Averages

More information

The law of conservation of mass: Mass can be neither created nor destroyed. It can only be transported or stored.

The law of conservation of mass: Mass can be neither created nor destroyed. It can only be transported or stored. UDMETL COCEPTS OR LOW LYSIS We covere mehos of analysis of nonflowing fluis in he previous chaper. In his chaper, we evelop he funamenal conceps of flow analysis, incluing he way o escribe flui flow, naural

More information

ψ(t) = V x (0)V x (t)

ψ(t) = V x (0)V x (t) .93 Home Work Se No. (Professor Sow-Hsin Chen Spring Term 5. Due March 7, 5. This problem concerns calculaions of analyical expressions for he self-inermediae scaering funcion (ISF of he es paricle in

More information

Math 10B: Mock Mid II. April 13, 2016

Math 10B: Mock Mid II. April 13, 2016 Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.

More information

MOMENTUM CONSERVATION LAW

MOMENTUM CONSERVATION LAW 1 AAST/AEDT AP PHYSICS B: Impulse and Momenum Le us run an experimen: The ball is moving wih a velociy of V o and a force of F is applied on i for he ime inerval of. As he resul he ball s velociy changes

More information

Intermediate Macroeconomics: Mid-term exam May 30 th, 2016 Makoto Saito

Intermediate Macroeconomics: Mid-term exam May 30 th, 2016 Makoto Saito 1 Inermediae Macroeconomics: Mid-erm exam May 30 h, 2016 Makoo Saio Try he following hree roblems, and submi your answer in handwrien A4 aers. You are execed o dro your aers ino he mailbox assigned for

More information

3rd Int. Conference on Inverse Problems in Engineering

3rd Int. Conference on Inverse Problems in Engineering Inverse Problems in Engineering: Theor and Pracice 3rd In. Conference on Inverse Problems in Engineering June 13-18, 1999, Por Ludlow, WA, USA EXP06 OPTIMIZATION METHOD FOR ONE- AND TWO-DIMENSIONAL INVERSE

More information

NUMERICAL CALIBRATION FOR 3-AXIS ACCELEROMETERS AND MAGNETOMETERS

NUMERICAL CALIBRATION FOR 3-AXIS ACCELEROMETERS AND MAGNETOMETERS NUMERICAL CALIBRATION FOR 3-AXIS ACCELEROMETERS AND MAGNETOMETERS Frédéric Cams Sofware Engineer Sébasien Harasse Posdocoral fellow CNRS/LAAS - Universi of Toulouse France 7 ave colonel Roche, 377 Toulouse

More information

4.5. Embedded Closed Curves in the Plane

4.5. Embedded Closed Curves in the Plane 88 4. TYPE II SINGUARITIES 4.5. Embee Close Curves in he Plane In he secial case of he evoluion of an embee close curve in he lane, i is ossible o exclue a all he ye II singulaies. This, ogeher wih he

More information

Online Learning with Partial Feedback. 1 Online Mirror Descent with Estimated Gradient

Online Learning with Partial Feedback. 1 Online Mirror Descent with Estimated Gradient Avance Course in Machine Learning Spring 2010 Online Learning wih Parial Feeback Hanous are joinly prepare by Shie Mannor an Shai Shalev-Shwarz In previous lecures we alke abou he general framework of

More information

2. The following diagram shows a circular loop of wire in a uniform magnetic field that points out of the page.

2. The following diagram shows a circular loop of wire in a uniform magnetic field that points out of the page. 1. Two elecromagneic waves ravel hrough emp space. Wave A as a wavelengh of 700 nm (red ligh), while Wave B has a wavelengh of 400 nm (blue ligh). Which saemen is rue? A) Wave A ravels faser because i

More information

Linear Surface Gravity Waves 3., Dispersion, Group Velocity, and Energy Propagation

Linear Surface Gravity Waves 3., Dispersion, Group Velocity, and Energy Propagation Chaper 4 Linear Surface Graviy Waves 3., Dispersion, Group Velociy, and Energy Propagaion 4. Descripion In many aspecs of wave evoluion, he concep of group velociy plays a cenral role. Mos people now i

More information

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time. Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike

More information

Solutions of Ecological models by Homotopy-Perturbation Method (HPM)

Solutions of Ecological models by Homotopy-Perturbation Method (HPM) Ques Journals Journal of Research in Alied Mahemaics Volume ~ Issue 0 (0) : 7- ISSN(Online) : 9-07 ISSN (Prin):9-075 www.quesjournals.org Research Paer Soluions of Ecological models by Homooy-Perurbaion

More information

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters?

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters? of Topics wih MATLAB Shan He School for Compuaional Science Universi of Birmingham Module 6-3836: Compuaional Modelling wih MATLAB Wha will we learn from his lecure Aim of analsis: Aim of analsis. Some

More information

Second Order Linear Differential Equations

Second Order Linear Differential Equations Second Order Linear Differenial Equaions Second order linear equaions wih consan coefficiens; Fundamenal soluions; Wronskian; Exisence and Uniqueness of soluions; he characerisic equaion; soluions of homogeneous

More information

Chapter 3 Kinematics in Two Dimensions

Chapter 3 Kinematics in Two Dimensions Chaper 3 KINEMATICS IN TWO DIMENSIONS PREVIEW Two-dimensional moion includes objecs which are moing in wo direcions a he same ime, such as a projecile, which has boh horizonal and erical moion. These wo

More information

Settling Time Design and Parameter Tuning Methods for Finite-Time P-PI Control

Settling Time Design and Parameter Tuning Methods for Finite-Time P-PI Control Journal of Conrol Science an Engineering (6) - oi:.765/8-/6.. D DAVID PUBLISHING Seling ime Design an Parameer uning Mehos for Finie-ime P-PI Conrol Keigo Hiruma, Hisaazu Naamura an Yasuyui Saoh. Deparmen

More information

Material Resistance and Friction in Cold Rolling

Material Resistance and Friction in Cold Rolling h World Congresses of Srucural and Mulidiscilinary Oimizaion Rio de Janeiro, 30 May - 03 June 200, Brazil Maerial Resisance and Fricion in Cold Rolling A.K. Tieu, C. You, H.T. Zhu, C. Lu, Z.Y. Jiang and

More information

Introduction to continuum mechanics-ii

Introduction to continuum mechanics-ii Ioan R. IONESCU LSPM, Universiy Paris 3 (Sorbonne-Paris-Cié) July 3, 27 Coninuum kinemaics B uneforme/reference configuraion, B eforme/acual configuraion Moion: x(, X ) = X + u(, X ) wih x(, B ) = B. Maerial/Lagrangian

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

Position predictive measurement method for time grating CNC rotary table

Position predictive measurement method for time grating CNC rotary table Posiion redicive measuremen mehod for ime graing CC roary able Liu Xiaokang a, Peng Donglin a, Yang Wei a and Fei Yeai b a Engineering Research Cener of Mechanical Tesing Technology and Equimen, Minisry

More information