A review on Lamb's atmospheric oscillations using initial value problem approach

Similar documents
MET 4302 Midterm Study Guide 19FEB18

ε(ω,k) =1 ω = ω'+kv (5) ω'= e2 n 2 < 0, where f is the particle distribution function and v p f v p = 0 then f v = 0. For a real f (v) v ω (kv T

Generalized logistic map and its application in chaos based cryptography

Wave Drift Force in a Two-Layer Fluid of Finite Depth

On Gravity Waves on the Surface of Tangential Discontinuity

ME scope Application Note 16

International Journal of Mathematics Trends and Technology- Volume3 Issue4-2012

Free-surface potential flow of an ideal fluid due to a singular sink

Central Force Motion Challenge Problems

Simplifications to Conservation Equations

Notes on pressure coordinates Robert Lindsay Korty October 1, 2002

Dimensional perturbation theory for Regge poles

Highly improved convergence of the coupled-wave method for TM polarization

Statics and dynamics: some elementary concepts

Contribution of the cosmological constant to the relativistic bending of light revisited

Quantum radiation force on a moving mirror for a thermal and a coherent field

AN ALTERNATIVE DESCRIPTION OF SWING-BY TRAJECTORIES IN TWO AND THREE DIMENSIONS

ESCI 342 Atmospheric Dynamics I Lesson 10 Vertical Motion, Pressure Coordinates

AP Calculus Testbank (Chapter 10) (Mr. Surowski)

STRONG TYPE INEQUALITIES AND AN ALMOST-ORTHOGONALITY PRINCIPLE FOR FAMILIES OF MAXIMAL OPERATORS ALONG DIRECTIONS IN R 2

Topological theory of non-hermitian photonic systems

Analytic approximation for the modified Bessel function I 2/3

Session 5: Review of Classical Astrodynamics

Unsteady Flow of a Dusty Conducting Fluid through porous medium between Parallel Porous Plates with Temperature Dependent Viscosity and Heat Source

Solutions 4: Free Quantum Field Theory

Section 4: Electromagnetic Waves 2

Bulletin of the Seismological Society of America

Slow viscous flow in a microchannel with similar and different superhydrophobic walls

Chapter 2 Introductory Concepts of Wave Propagation Analysis in Structures

Introduction to Landau s Fermi Liquid Theory

Numerical modeling of a cutting torch

A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE

Methods for detecting fatigue cracks in gears

References: 1. Cohen Tannoudji Chapter 5 2. Quantum Chemistry Chapter 3

Single and double coincidence nucleon spectra in the weak decay of Λ hypernuclei

Effect of Hall current and rotation on heat transfer in MHD flow of oscillating dusty fluid in a porous channel

You may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.

δq T = nr ln(v B/V A )

Lesson 16 Electric Potential

High speed wind tunnels 2.0 Definition of high speed. 2.1 Types of high speed wind tunnels

First order sensitivity analysis of electron acceleration in dual grating type dielectric laser accelerator structures

NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS

A Note on Massless Quantum Free Scalar Fields. with Negative Energy Density

International Association of Scientific Innovation and Research (IASIR) (An Association Unifying the Sciences, Engineering, and Applied Research)

State Estimation with ARMarkov Models

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 13: PF Design II The Coil Solver

Application of the Shannon-Kotelnik theorem on the vortex structures identification

Probing spin-charge separation in a Tomonaga-Luttinger liquid

Paper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation

Higher order theory for analytic saddle point approximations to the Ρ Ρ and Ρ Ś reflected arrivals at a solid/solid interface

Natural frequency analysis of fluid-conveying pipes in the ADINA system

On the Fluid Dependence of Rock Compressibility: Biot-Gassmann Refined

DIVISIBILITY CRITERIA FOR CLASS NUMBERS OF IMAGINARY QUADRATIC FIELDS

Green s Functions and Topological Configurations

Gravitational radiation from a particle in bound orbit around a black hole; relativistic correction

Wolfgang POESSNECKER and Ulrich GROSS*

Week 8 lectures. ρ t +u ρ+ρ u = 0. where µ and λ are viscosity and second viscosity coefficients, respectively and S is the strain tensor:

Hans-Joachim Herzog, Gerd Vogel DWD, FE14, Potsdam. COSMO meeting -, Milano September 2004

How beaming of gravitational waves compares to the beaming of electromagnetic waves: impacts to gravitational wave detection

A Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition

Dynamic System Eigenvalue Extraction using a Linear Echo State Network for Small-Signal Stability Analysis a Novel Application

The Propagation of Strong Shock Waves in Magnetohydrodynamics

Landau Theory of the Fermi Liquid

SEG Houston 2009 International Exposition and Annual Meeting

1 University of Edinburgh, 2 British Geological Survey, 3 China University of Petroleum

FE FORMULATIONS FOR PLASTICITY

Understanding DPMFoam/MPPICFoam

Classical gas (molecules) Phonon gas Number fixed Population depends on frequency of mode and temperature: 1. For each particle. For an N-particle gas

PASSAGE OF PARTICLES THROUGH A PLASMA

ME 375 System Modeling and Analysis. Homework 11 Solution. Out: 18 November 2011 Due: 30 November 2011 = + +

Anisotropy included in a nanoscale superconductor: Theoretical development

Pretest (Optional) Use as an additional pacing tool to guide instruction. August 21

NONRELATIVISTIC STRONG-FIELD APPROXIMATION (SFA)

Recent flow results at RHIC

Principles of Computed Tomography (CT)

Chapter 6: Sound Wave Equation

Casimir Force Between the Two Moving Conductive Plates.

Low frequency modes in strongly coupled dusty plasmas

Solved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points.

Study of the circulation theory of the cooling system in vertical evaporative cooling generator

Weather and Climate Laboratory Spring 2009

The influence of a semi-infinite atmosphere on solar oscillations

SELF-SIMILAR FLOW OF A MIXTURE OF A NON-IDEAL GAS AND SMALL SOLID PARTICLES WITH INCREASING ENERGY BEHIND A SHOCK WAVE UNDER A GRAVITATIONAL FIELD

4. Score normalization technical details We now discuss the technical details of the score normalization method.

Cybernetic Interpretation of the Riemann Zeta Function

Dufour-Driven Thermosolutal Convection of the Veronis Type

Damage Identification from Power Spectrum Density Transmissibility

arxiv: v1 [nucl-th] 19 Feb 2018

Towards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK

Determining the Critical Point of a Sigmoidal Curve via its Fourier Transform

Physics 582, Problem Set 6 Solutions

Chapter 1 Fundamentals

Beam-Beam Stability in Electron-Positron Storage Rings

arxiv: v2 [cond-mat.stat-mech] 10 Feb 2012

RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES

arxiv: v2 [cond-mat.str-el] 3 Jul 2014

A finite volume method for solving the gravity wave-model equations

Hypersonic flow: introduction

Supplementary Material: Crumpling Damaged Graphene

Salt intrusion response to changes in tidal amplitude during low river flow in the Modaomen Estuary, China

Transcription:

Journal of Physics: Conference Series OPEN ACCESS A review on Lamb's atmosheric oscillations using initial value roblem aroach To cite this article: Ángel De Andrea González 04 J. Phys.: Conf. Ser. 56 004 View the article online for udates and enhancements. This content was downloaded from IP address 48.5.3.83 on 09/0/08 at 9:53

4th International Worksho & Summer School on Plasma Physics 00 Journal of Physics: Conference Series 56 (04 004 doi:0.088/74-6596/56//004 A review on Lamb s atmosheric oscillations using initial value roblem aroach Ángel De Andrea González Deartamento de Física. Escuela Politécnica Suerior. Universidad Carlos III de Madrid. Av. Universidad, 30. 89 Leganés. Sain. E-mail: aandrea@fis.uc3m.es Abstract. Waves at a surface of discontinuity in the atmoshere were analysed in 90 by Lamb, who derived, using normal mode aroach, an analytical disersion relation for a disete mode (surface mode. Lamb examined the case of waves roagated along a horizontal lane where the equilibrium temerature is discontinuous. For simlicity, the uer and the lower regions are considered incomressible. The oscillations are treated in the ideal (dissiationless limit and the uniform gravitational acceleration is taken to be co-aligned with the revailing temerature gradient. In this work, in order to show how the modes aear in the resonse of a surface discontinuity to an initial erturbation, we consider the initial value roblem (IPV. The main difference from the standard analysis is that solutions to the linearized equations of motion which satisfy general conditions are obtained in terms of Fourier-Lalace transform of the hydrodynamics variables. These transforms can be inverted exlicitly to exress the fluid variables as integrals of Green s functions multilied by initial data. In addition to disete mode (surface mode, sets of continuum modes due to branch cuts in the comlex lane, not treated exlicitly in the literature, aears.. Introduction Several authors have urosed solutions of the initial value roblem through the use of a Lalace Transform in time [-4]. The rincial finding of the initial value roblem aroach is that, in addition to the disete eigenvalues linked to the normal modes, there exists a continuous sectrum of eigenvalues. Thus, the modal aroach cannot rovide a comlete solution.. Linear stability analysis We consider two infinitely incomressible lasmas in a steady state which are searated by a lane interface at y0 and subjected to a gravitational field g r directed erendicular to the interface, where q and q are both the inverse density-gradient scale length at each region. By following the aroach develoed by Lamb (90 [], it is ossible to reduce the linearized equations of continuity, momentum and energy for a olytroe to a set of artial differential equations. These equations can be easily combined into a stream-function equation: Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s and the title of the work, journal citation and DOI. Published under licence by Ltd

4th International Worksho & Summer School on Plasma Physics 00 Journal of Physics: Conference Series 56 (04 004 doi:0.088/74-6596/56//004 3 ψ ψ ψ ψ 0 y k y ( Where the dimensionless normalized arameters are: t kg, k k q, t t kg. We define the Lalace transform ψ ( y, ω ω in the comlex ω -lane: i σ i ωt ψ ( y, ω ω ψ ( y, t e dω π ( i σ Then, we consider the initial value roblem for equation (, and introduce the Lalace transform. So, we obtain the following equation: ψ ψ ψ ω Γ y y k y k ω ω ω ( (3 being: i ω F F i ω F F Γ ( y iω F F ω y k y (4 is: where the Lalace transform function ψ ( y, ( y, ψ L ω ψ s i ωf F With the following initial conditions: ψ F ψ ( y, ; F t 0 y ( y, 0 (5 (6 With Lalace contour σ chosen that the conditions at ψ s y at y 0. has no singularities in Re( ω ψ ω σ. Equation is subjected to y that ψ ω 0, the continuity of ψ ω at y 0 and the jum conditions on

4th International Worksho & Summer School on Plasma Physics 00 Journal of Physics: Conference Series 56 (04 004 doi:0.088/74-6596/56//004 We consider the Green s function associated with equation (3. So, the solution of equation ( is: (, (,, Γ( ψ y ω G y ω y y dy (7 s o o 0.. Results and discussion The Green s function for y > 0 : (, ω, y G y o y β y o α k kr ω e r ω α r ω β kr k (8 G ±,, y o 0, the square roots are defined so that Re(α > 0 y Re( β > 0 for Re (ω > σ. It is evident that Green s function has a ole (disete sectra. Exanding Green s function about α 0, and β 0, this function does not turn out to be an In accord with the boundary conditions ( ω even function of α and β. So, the Green s function has branch oints at continuous sectra neglected by Lamb: (associated with ω γ b 4k γ 4k b y en 4kr γ 4k r b (i.e., α 0 y β 0 (9 It is interesting to follow the migration of the ole as k is varied. For k > k, γ > γ b. If k is diminished, then the γ ole moves toward γ b. For k k, then γ γ b. For k < k the ole dros on to the lower sheet of β and the eigenmode disaears. The disaearance of the γ ole was neglected by Lamb too. Thus, there will be some itical wave number k such for k > k k < k eigenmodes exist (do not exist. Lamb obtained the following disersion relation []: ( ( ω ω - - 0 k ( r r r- (0 Being, r q.we show that Lamb s disersion relation, the itical wave number q k is: 3 r k r for r ( 3

4th International Worksho & Summer School on Plasma Physics 00 Journal of Physics: Conference Series 56 (04 004 doi:0.088/74-6596/56//004 w kg.0 0.8 0.6 0.4 0. k 0.0 0 4 6 8 0 q Figure. Dimensionless normalized growth rate versus dimensionless normalized wave number for r4. Solid line corresonds to the Lamb s results; dots corresond to our numerical results. There is an excellent agreement with Lamb s results for k > k. Dashed line and dot dashed line corresond to branch oints. 3. Conclusion Waves at a surface of discontinuity in the atmoshere are analysed using initial value roblem (IPV. So that, sets of continuum modes due to branch cuts in the comlex lane, not treated exlicitly in the literature, exlain the disaearance of the disersion relation. Thus, there will be some itical wave number k such that for k > k ( k < k eigenmodes exist (do not exist. 4. References [] Lamb H 90 Proc. Lond. Math. Soc. A 84 55-7 [] Ott E and Russell D A 978 Physical Review Letters 4 048-5 [3] Russell D A and Ott E 979 Journal of Geohysical Research 84 6573-79. [4] Bogdan T J and Cally P S 997 Proc. R. Soc. Lond. A, Math. Phys. Eng. Sci., 453 943-6 4