Английская Википедия:Averaged Lagrangian

High-altitude wave cloud formed over the Hampton area at Burra, South Australia on 16 January 2007.

In continuum mechanics, Whitham's averaged Lagrangian method – or in short Whitham's method – is used to study the Lagrangian dynamics of slowly-varying wave trains in an inhomogeneous (moving) medium. The method is applicable to both linear and non-linear systems. As a direct consequence of the averaging used in the method, wave action is a conserved property of the wave motion. In contrast, the wave energy is not necessarily conserved, due to the exchange of energy with the mean motion. However the total energy, the sum of the energies in the wave motion and the mean motion, will be conserved for a time-invariant Lagrangian. Further, the averaged Lagrangian has a strong relation to the dispersion relation of the system.

The method is due to Gerald Whitham, who developed it in the 1960s. It is for instance used in the modelling of surface gravity waves on fluid interfaces,^[1]^[2] and in plasma physics.^[3]^[4]

Resulting equations for pure wave motion

In case a Lagrangian formulation of a continuum mechanics system is available, the averaged Lagrangian methodology can be used to find approximations for the average dynamics of wave motion – and (eventually) for the interaction between the wave motion and the mean motion – assuming the envelope dynamics of the carrier waves is slowly varying. Phase averaging of the Lagrangian results in an averaged Lagrangian, which is always independent of the wave phase itself (but depends on slowly varying wave quantities like wave amplitude, frequency and wavenumber). By Noether's theorem, variation of the averaged Lagrangian <math>\mathcal{L}</math> with respect to the invariant wave phase <math>\theta(\boldsymbol{x},t)</math> then gives rise to a conservation law:^[5]

Шаблон:Equation box 1 = 0.</math> | ref = 1 | cellpadding = 5 }}

This equation states the conservation of wave action – a generalization of the concept of an adiabatic invariant to continuum mechanics – with^[6]

<math display="block">\mathcal{A} \equiv -\frac{\partial\mathcal{L}}{\partial(\partial_t \theta)} = +\frac{\partial\mathcal{L}}{\partial\omega}</math> <math display="block">\boldsymbol{\mathcal{B}} \equiv -\frac{\partial\mathcal{L}}{\partial(\boldsymbol{\nabla}\theta)} = -\frac{\partial\mathcal{L}}{\partial\boldsymbol{k}}</math>

being the wave action <math>\mathcal{A}</math> and wave action flux <math>\boldsymbol{\mathcal{B}}</math> respectively. Further <math>\boldsymbol{x}</math> and <math>t</math> denote space and time respectively, while <math>\boldsymbol{\nabla}</math> is the gradient operator. The angular frequency <math>\omega(\boldsymbol{x},t)</math> and wavenumber <math>\boldsymbol{k}(\boldsymbol{x},t)</math> are defined as^[7]

Шаблон:Equation box 1

and both are assumed to be slowly varying. Due to this definition, <math>\omega(\boldsymbol{x},t)</math> and <math>\boldsymbol{k}(\boldsymbol{x},t)</math> have to satisfy the consistency relations:

Шаблон:Equation box 1

The first consistency equation is known as the conservation of wave crests, and the second states that the wavenumber field <math>\boldsymbol{k}(\boldsymbol{x},t)</math> is irrotational (i.e. has zero curl).

Method

The averaged Lagrangian approach applies to wave motion – possibly superposed on a mean motion – that can be described in a Lagrangian formulation. Using an ansatz on the form of the wave part of the motion, the Lagrangian is phase averaged. Since the Lagrangian is associated with the kinetic energy and potential energy of the motion, the oscillations contribute to the Lagrangian, although the mean value of the wave's oscillatory excursion is zero (or very small).

The resulting averaged Lagrangian contains wave characteristics like the wavenumber, angular frequency and amplitude (or equivalently the wave's energy density or wave action). But the wave phase itself is absent due to the phase averaging. Consequently, through Noether's theorem, there is a conservation law called the conservation of wave action.

Originally the averaged Lagrangian method was developed by Whitham for slowly-varying dispersive wave trains.^[8] Several extensions have been made, e.g. to interacting wave components,^[9]^[10] Hamiltonian mechanics,^[8]^[11] higher-order modulational effects,^[12] dissipation effects.^[13]

Variational formulation

The averaged Lagrangian method requires the existence of a Lagrangian describing the wave motion. For instance for a field <math>\varphi(\boldsymbol{x},t)</math>, described by a Lagrangian density <math>L\left(\partial_t\varphi,\boldsymbol{\nabla}\varphi,\varphi\right),</math> the principle of stationary action is:^[14]

 \delta \left( \iiiint\, 
   L\left( \partial_t\varphi(\boldsymbol{x},t),
           \boldsymbol{\nabla}\varphi(\boldsymbol{x},t),
           \varphi(\boldsymbol{x},t)
    \right)\, 
   \text{d}\boldsymbol{x}\, \text{d}t \right) = 0,

</math>

with <math>\boldsymbol{\nabla}</math> the gradient operator and <math>\partial_t</math> the time derivative operator. This action principle results in the Euler–Lagrange equation:^[14]

 \partial_t \left( \frac{\partial L}{\partial \left( \partial_t \varphi \right)} \right)
 +\boldsymbol{\nabla} \cdot
   \left( \frac{\partial L}{\partial \left( \boldsymbol{\nabla} \varphi \right)} \right)
 - \frac{\partial L}{\partial \varphi}
 = 0,

</math>

which is the second-order partial differential equation describing the dynamics of <math>\varphi.</math> Higher-order partial differential equations require the inclusion of higher than first-order derivatives in the Lagrangian.^[14]

Example

For example, consider a non-dimensional and non-linear Klein–Gordon equation in one space dimension <math>x</math>:^[15] Шаблон:NumBlk

This Euler–Lagrange equation emerges from the Lagrangian density:^[15] Шаблон:NumBlk

The small-amplitude approximation for the Sine–Gordon equation corresponds with the value <math display="inline">\sigma = -\tfrac{1}{24}.</math>^[16] For <math>\sigma = 0</math> the system is linear and the classical one-dimensional Klein–Gordon equation is obtained.

Slowly-varying waves

Slowly-varying linear waves

Whitham developed several approaches to obtain an averaged Lagrangian method.^[14]^[17] The simplest one is for slowly-varying linear wavetrains, which method will be applied here.^[14]

The slowly-varying wavetrain –without mean motion– in a linear dispersive system is described as:^[18]

 \varphi \sim \Re\left\{ A(\boldsymbol{x},t)\, e^{i\theta(\boldsymbol{x},t)} \right\} 
 = a(\boldsymbol{x},t)\, \cos \left( \theta(\boldsymbol{x},t) + \alpha \right),

</math> with <math display="block">a = \left| A \right|</math> and <math display="block">\alpha = \arg\left\{ A \right\},</math>

where <math>\theta(\boldsymbol{x},t)</math> is the real-valued wave phase, <math>|A|</math> denotes the absolute value of the complex-valued amplitude <math>A(\boldsymbol{x},t),</math> while <math>\arg\{A\}</math> is its argument and <math>\Re\{A\}</math> denotes its real part. The real-valued amplitude and phase shift are denoted by <math>a</math> and <math>\alpha</math> respectively.

Now, by definition, the angular frequency <math>\omega</math> and wavenumber vector <math>\boldsymbol{k}</math> are expressed as the time derivative and gradient of the wave phase <math>\theta(\boldsymbol{x},t)</math> as:^[7]

<math display="block">\omega \equiv -\partial_t \theta</math> and <math display="block">\boldsymbol{k} \equiv +\boldsymbol{\nabla} \theta.</math>

As a consequence, <math>\omega(\boldsymbol{x},t)</math> and <math>\boldsymbol{k}(\boldsymbol{x},t)</math> have to satisfy the consistency relations:

<math display="block">\partial_t \boldsymbol{k} + \boldsymbol{\nabla} \omega = \boldsymbol{0}</math> and <math display="block">\boldsymbol{\nabla} \times \boldsymbol{k}=\boldsymbol{0}.</math>

These two consistency relations denote the "conservation of wave crests", and the irrotationality of the wavenumber field.

Because of the assumption of slow variations in the wave train – as well as in a possible inhomogeneous medium and mean motion – the quantities <math>A,</math> <math>a,</math> <math>\omega,</math> <math>\boldsymbol{k}</math> and <math>\alpha</math> all vary slowly in space <math>\boldsymbol{x}</math> and time <math>t</math> – but the wave phase <math>\theta</math> itself does not vary slowly. Consequently, derivatives of <math>a,</math> <math>\omega,</math> <math>\boldsymbol{k}</math> and <math>\alpha</math> are neglected in the determination of the derivatives of <math>\varphi(\boldsymbol{x},t)</math> for use in the averaged Lagrangian:^[14]

<math display="block">\partial_t\varphi \approx +\omega\, a\, \sin (\theta+\alpha)</math> and <math display="block">\boldsymbol{\nabla} \varphi \approx -\boldsymbol{k}\, a\, \sin (\theta + \alpha).</math>

Next these assumptions on <math>\varphi(\boldsymbol{x},t)</math> and its derivatives are applied to the Lagrangian density <math>L\left(\partial_t\varphi,\boldsymbol{\nabla}\varphi,\varphi\right).</math>

Slowly-varying non-linear waves

Several approaches to slowly-varying non-linear wavetrains are possible. One is by the use of Stokes expansions,^[19] used by Whitham to analyse slowly-varying Stokes waves.^[20] A Stokes expansion of the field <math>\varphi(\boldsymbol{x},t)</math> can be written as:^[19]

<math display="block">\varphi = a\, \cos \left( \theta + \alpha \right) + a_2\, \cos \left( 2 \theta + \alpha_2 \right) + a_3\, \cos \left( 3 \theta + \alpha_3 \right) + \cdots,</math>

where the amplitudes <math>a,</math> <math>a_2,</math> etc. are slowly varying, as are the phases <math>\alpha,</math> <math>\alpha_2,</math> etc. As for the linear wave case, in lowest order (as far as modulational effects are concerned) derivatives of amplitudes and phases are neglected, except for derivatives <math>\omega</math> and <math>\boldsymbol{k}</math> of the fast phase <math>\theta</math>:

<math display="block">\partial_t \varphi \approx +\omega a\, \sin \left( \theta + \alpha \right) + 2\omega a_2\, \sin \left( 2 \theta + \alpha_2 \right) + 3\omega a_3\, \sin \left( 3 \theta + \alpha_3 \right) + \cdots,</math> and <math display="block">\boldsymbol{\nabla} \varphi \approx -\boldsymbol{k} a\, \sin \left( \theta + \alpha \right) - 2\boldsymbol{k} a_2\, \sin \left( 2 \theta + \alpha_2 \right) - 3\boldsymbol{k} a_3\, \sin \left( 3 \theta + \alpha_3 \right) + \cdots.</math>

These approximations are to be applied in the Lagrangian density <math>L</math>, and its phase average <math>\bar{L}.</math>

Averaged Lagrangian for slowly-varying waves

For pure wave motion the Lagrangian <math>L\left(\partial_t\varphi,\boldsymbol{\nabla}\varphi,\varphi\right)</math> is expressed in terms of the field <math>\varphi(\boldsymbol{x},t)</math> and its derivatives.^[14]^[17] In the averaged Lagrangian method, the above-given assumptions on the field <math>\varphi(\boldsymbol{x},t)</math> –and its derivatives– are applied to calculate the Lagrangian. The Lagrangian is thereafter averaged over the wave phase <math>\theta</math>:^[14]

<math display="block"> \bar{L} = \frac{1}{2\pi} \int_0^{2\pi} L\left(\partial_t\varphi,\boldsymbol{\nabla}\varphi,\varphi\right) \text{d}\theta.</math>

As a last step, this averaging result <math>\bar{L}</math> can be expressed as the averaged Lagrangian density <math>\mathcal{L}(\omega,\boldsymbol{k},a)</math> – which is a function of the slowly varying parameters <math>\omega,</math> <math>\boldsymbol{k}</math> and <math>a</math> and independent of the wave phase <math>\theta</math> itself.^[14]

The averaged Lagrangian density <math>\mathcal{L}</math> is now proposed by Whitham to follow the average variational principle:^[14]

Шаблон:Equation box 1

From the variations of <math>\mathcal{L}</math> follow the dynamical equations for the slowly-varying wave properties.

Example

Continuing on the example of the nonlinear Klein–Gordon equation, see equations Шаблон:EquationNote and Шаблон:EquationNote, and applying the above approximations for <math>\varphi,</math> <math>\partial_t \varphi</math> and <math>\partial_x \varphi</math> (for this 1D example) in the Lagrangian density, the result after averaging over <math>\theta</math> is:

<math display="block"> \bar{L} = \tfrac 1 4 \left(\omega^2 - k^2 - 1 \right) a^2

       - \tfrac{3}{32} \sigma a^4
       + \left(\omega^2 - k^2 - \tfrac 1 4\right) a_2^2
       + \mathcal{O}{\left(a^6\right)},

</math> where it has been assumed that, in big-O notation, <math>a_2 = \mathcal{O}(a^2)</math> and <math>a_3 = \mathcal{O}(a^3)</math>. Variation of <math>\bar{L}</math> with respect to <math>a_2</math> leads to <math>a_2=0.</math> So the averaged Lagrangian is: Шаблон:NumBlk

For linear wave motion the averaged Lagrangian is obtained by setting <math>\sigma</math> equal to zero.

Set of equations emerging from the averaged Lagrangian

Applying the averaged Lagrangian principle, variation with respect to the wave phase <math>\theta</math> leads to the conservation of wave action:

 \partial_t \left( + \frac{\partial \mathcal{L}}{\partial \omega} \right)
 + \boldsymbol{\nabla} \cdot \left( - \frac{\partial \mathcal{L}}{\partial \boldsymbol{k}} \right) = 0,

</math>

since <math>\omega = -\partial_t\theta</math> and <math>\boldsymbol{k} = \boldsymbol{\nabla}\theta</math> while the wave phase <math>\theta</math> does not appear in the averaged Lagrangian density <math>\mathcal{L}</math> due to the phase averaging. Defining the wave action as <math>\mathcal{A}\equiv+\partial\mathcal{L}/\partial\omega</math> and the wave action flux as <math>\boldsymbol{\mathcal{B}}\equiv-\partial\mathcal{L}/\partial\boldsymbol{k}</math> the result is:

<math display="block">\partial_t \mathcal{A} + \boldsymbol{\nabla} \cdot \boldsymbol{\mathcal{B}} = 0.</math>

The wave action equation is accompanied with the consistency equations for <math>\omega</math> and <math>\boldsymbol{k}</math> which are:

<math display="block">\partial_t \boldsymbol{k} + \boldsymbol{\nabla} \omega = \boldsymbol{0}</math> and <math display="block">\boldsymbol{\nabla} \times \boldsymbol{k} = \boldsymbol{0}.</math>

Variation with respect to the amplitude <math>a</math> leads to the dispersion relation <math>\partial\mathcal{L}/\partial a = 0.</math>

Example

Continuing with the nonlinear Klein–Gordon equation, using the average variational principle on equation Шаблон:EquationNote, the wave action equation becomes by variation with respect to the wave phase <math>\theta:</math> <math display="block">

 \partial_t \left( \tfrac 1 2 \omega a^2 \right) + \partial_x \left( \tfrac 1 2 k a^2 \right) = 0,

</math> and the nonlinear dispersion relation follows from variation with respect to the amplitude <math>a:</math> <math display="block">\omega^2 = k^2 + 1 + \tfrac 3 4 \sigma a^2.</math>

So the wave action is <math display="inline">\mathcal{A} = \tfrac 1 2 \omega a^2</math> and the wave action flux <math display="inline">\mathcal{B} = \tfrac 1 2 k a^2.</math> The group velocity <math>v_g</math> is <math>v_g \equiv \mathcal{B}/\mathcal{A} = k/\omega.</math>

Mean motion and pseudo-phase

Шаблон:Empty section

Conservation of wave action

The averaged Lagrangian is obtained by integration of the Lagrangian over the wave phase. As a result, the averaged Lagrangian only contains the derivatives of the wave phase <math>\theta</math> (these derivatives being, by definition, the angular frequency and wavenumber) and does not depend on the wave phase itself. So the solutions will be independent of the choice of the zero level for the wave phase. Consequently – by Noether's theorem – variation of the averaged Lagrangian <math>\bar{\mathcal{L}}</math> with respect to the wave phase results in a conservation law:

Шаблон:Equation box 1 = 0,</math> }} where

<math display="block"> \mathcal{A}\equiv\frac{\delta \bar{\mathcal{L}}}{\delta \omega} = -\frac{\delta \bar{\mathcal{L}}}{\delta \left( \partial_t \theta \right)},</math>
<math display="block"> \boldsymbol{\mathcal{B}}\equiv-\frac{\delta \bar{\mathcal{L}}}{\delta \boldsymbol{k}} = -\frac{\delta \bar{\mathcal{L}}}{\delta \left( \boldsymbol{\nabla} \theta \right)},</math>

with <math>\mathcal{A}</math> the wave action and <math>\boldsymbol{\mathcal{B}}</math> the wave action flux. Further <math>\partial_t</math> denotes the partial derivative with respect to time, and <math>\boldsymbol{\nabla}</math> is the gradient operator. By definition, the group velocity <math>\boldsymbol{v}_g</math> is given by:

<math display="block">\boldsymbol{\mathcal{B}} \equiv \boldsymbol{v}_g \mathcal{A}.</math>

Note that in general the energy of the wave motion does not need to be conserved, since there can be an energy exchange with a mean flow. The total energy – the sum of the energies of the wave motion and the mean flow – is conserved (when there is no work by external forces and no energy dissipation).

Conservation of wave action is also found by applying the generalized Lagrangian mean (GLM) method to the equations of the combined flow of waves and mean motion, using Newtonian mechanics instead of a variational approach.^[21]

Conservation of energy and momentum

Шаблон:Empty section

Connection to the dispersion relation

Pure wave motion by linear models always leads to an averaged Lagrangian density of the form:^[14] <math display="block">\mathcal{L} = G(\omega,\boldsymbol{k}) a^2.</math>

Consequently, the variation with respect to amplitude: <math>\partial \mathcal{L}/\partial a = 0</math> gives <math display="block">G(\omega,\boldsymbol{k}) = 0.</math>

So this turns out to be the dispersion relation for the linear waves, and the averaged Lagrangian for linear waves is always the dispersion function <math>G(\omega,\boldsymbol{k})</math> times the amplitude squared.

More generally, for weakly nonlinear and slowly modulated waves propagating in one space dimension and including higher-order dispersion effects – not neglecting the time and space derivatives <math>\partial_t a</math> and <math>\partial_x a</math> of the amplitude <math>a(\mu x,\mu t)</math> when taking derivatives, where <math>\mu\ll 1</math> is a small modulation parameter – the averaged Lagrangian density is of the form:^[22] <math display="block">\mathcal{L} =

 G(\omega,k) a^2 + G_2(\omega,k) a^4 
 + \tfrac 1 2 \mu^2 \left( 
         G_{\omega\omega} (\partial_T a)^2 
     + 2 G_{\omega k}     (\partial_T a) (\partial_X a) 
     +   G_{kk}           (\partial_X a)^2 
   \right),

</math> with the slow variables <math>X=\mu x</math> and <math>T=\mu t.</math>

References

Notes

Шаблон:Reflist

Publications by Whitham on the method

Шаблон:Ref begin An overview can be found in the book:

Шаблон:Citation

Some publications by Whitham on the method are:

Шаблон:Ref end

Партнерские ресурсы
Криптовалюты	Обмен криптовалют - www.bestchange.ru Криптовалютная биржа CoinEx Криптовалютная биржа Binance HIVE OS - операционная система для майнинга e4pool - Мультивалютный пул для майнинга.
Магазины	AliExpress — глобальная виртуальная (в Интернете) торговая площадка, предоставляющая возможность покупать товары производителей из КНР; computeruniverse.net - Интернет-магазин компьютеров(Промо код 5 Евро на первую покупку:FWWC3ZKQ);
Хостинг	DigitalOcean - американский провайдер облачных инфраструктур, с главным офисом в Нью-Йорке и с центрами обработки данных по всему миру;
Разное	Викиум - Онлайн-тренажер для мозга Like Центр - Центр поддержки и развития предпринимательства. Gamersbay - лучший магазин по бустингу для World of Warcraft. Ноотропы OmniMind N°1 - Усиливает мозговую активность. Повышает мотивацию. Улучшает память. Санкт-Петербургская школа телевидения - это федеральная сеть образовательных центров, которая имеет филиалы в 37 городах России. Lingualeo.com — интерактивный онлайн-сервис для изучения и практики английского языка в увлекательной игровой форме. Junyschool (Джунискул) – международная школа программирования и дизайна для детей и подростков от 5 до 17 лет, где ученики осваивают компьютерную грамотность, развивают алгоритмическое и креативное мышление, изучают основы программирования и компьютерной графики, создают собственные проекты: игры, сайты, программы, приложения, анимации, 3D-модели, монтируют видео. Умназия - Интерактивные онлайн-курсы и тренажеры для развития мышления детей 6-13 лет SkillBox - это один из лидеров российского рынка онлайн-образования. Среди партнеров Skillbox ведущий разработчик сервисного дизайна AIC, медиа-компания Yoola, первое и самое крупное русскоязычное аналитическое агентство Tagline, онлайн-школа дизайна и иллюстрации Bang! Bang! Education, оператор PR-рынка PACO, студия рисования Draw&Go, агентство performance-маркетинга Ingate, scrum-студия Sibirix, имидж-лаборатория Персона. «Нетология» — это университет по подготовке и дополнительному обучению специалистов в области интернет-маркетинга, управления проектами и продуктами, дизайна, Data Science и разработки. В рамках Нетологии студенты получают ценные теоретические знания от лучших экспертов Рунета, выполняют практические задания на отработку полученных навыков, общаются с экспертами и единомышленниками. Познакомиться со всеми продуктами подробнее можно на сайте https://netology.ru, линейка курсов и профессий постоянно обновляется. StudyBay Brazil – это онлайн биржа для португалоговорящих студентов и авторов! Студент получает уникальную работу любого уровня сложности и больше свободного времени, в то время как у автора появляется дополнительный заработок и бесценный опыт. Автор24 — самая большая в России площадка по написанию учебных работ: контрольные и курсовые работы, дипломы, рефераты, решение задач, отчеты по практике, а так же любой другой вид работы. Сервис сотрудничает с более 70 000 авторов. Более 1 000 000 работ уже выполнено. StudyBay – это онлайн биржа для англоязычных студентов и авторов! Студент получает уникальную работу любого уровня сложности и больше свободного времени, в то время как у автора появляется дополнительный заработок и бесценный опыт.

Английская Википедия:Averaged Lagrangian

Содержание

Resulting equations for pure wave motion

Method

Variational formulation

Example

Slowly-varying waves

Slowly-varying linear waves

Slowly-varying non-linear waves

Averaged Lagrangian for slowly-varying waves

Example

Set of equations emerging from the averaged Lagrangian

Example

Mean motion and pseudo-phase

Conservation of wave action

Conservation of energy and momentum

Connection to the dispersion relation

References

Notes

Publications by Whitham on the method

Further reading

Навигация

Действия на странице

Действия на странице

Персональные инструменты

Навигация

Поиск

Инструменты