singular definition reform - université grenoble alpes

15
The purpose of this lecture is to explain how boundaries can be taken into account within singular perturbation problems By definition ce singular problem is an equation of Reform ct ne f with EZO where of E is an operator of order d i with the following property if ne à foe f Consone strong sense then à is a solution of JE f where À is an operator of order Lol Asa consequence we cannot enforce the same number of boundary conditions force and for at presence of boundary layers within ne n es à BL part Questions 1 How can we construct boundary layers Systematic construction 2 Justification of the convergence and validity of the approximation

Upload: others

Post on 10-May-2022

1 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: singular definition Reform - Université Grenoble Alpes

The purpose of this lecture is to explain howboundaries can be taken into account within

singular perturbation problems Bydefinitionce singular problem is an equation ofReform

ct ne fwith EZO where of E is an operator oforderd i with the following propertyif ne à foe f Consone strong

sense then à is a solution of

JE fwhere À is an operator of order Lol

Asa consequence we cannot enforce the samenumber of boundary conditions force and

for at presence of boundary layers within nen es à BL part

Questions1 How can we construct boundary layersSystematic construction

2 Justification ofthe convergence and

validity of the approximation

Page 2: singular definition Reform - Université Grenoble Alpes

Remarks D In some cases the BL hason effect on the mean flow2 Unfortunately there is no completelyunified framework to address the two pointsabove Study of two examples one simpleone more complicated

I IntroductionDAtoymodel.im

Enfin usé L x ECO D

ukot.ua O

Existence and uniqueness of a solution n'byLax Milgram theoremA priori estimate x ne exp Ex

SÈME esseulée ÎE Ç nen é fée e fuséeEN

Elway Muette ICul ci init à sel

Ansatzi vieux tube tué f

Page 3: singular definition Reform - Université Grenoble Alpes

1Boundary conditions for à

Need to look at BL close to second set

Plug Ansatz into equation

t LITUBL satisfies dzzubtr dzuf O

u.BE temps f 3

vite satisfies 2g une ciel 0

NE 0

Consequence 7 0 0 no BL ondeleftucxt ffsxueu enpf.EE un

watt satisfies up with small remaindersatisfies the Be at xd

x D with exponentiallysmall remainder

consequence EW uahxlpftlli uaq.seRenard EH Bible E Èeenpf 2 de

OLD

Page 4: singular definition Reform - Université Grenoble Alpes

Thea energyestimate mes the boundarylayer

2 A bit of methodologymm

Ref Gérard Vaut Paul

Simplifications Amine that of l is a linearoperator with constant coefficients and that1 RI IRM xRIfflat boundaryItis legitimate to look for modal solutionsofcffi.ie we lookfor special solutions ofthe form

UE x UE enpfik.ae dxnk'ERw tangential Fourier moded C Q Re d 0UE c pend

Boundary layer Red I

Plug into equation linear systemoftheformA k d UE o ci

where Alisa did matrix with polynomialcoefficients in k d land E in generalCe symbolof et eThere exists a mon trivial solution iff

det A Kkl D 0 2

Page 5: singular definition Reform - Université Grenoble Alpes

Polynomial int with coefficients that are

polynomial in lsland E sayof degree m

This polynomial has ne complex rootsdz du

Complex analysis IggC Q Jj lé such that

Dj E k EH Jj as E 0

NE kD d solution of 2 with Red 0

d dim Vet U polofen de MCE le

d number of boundary conditions thatcourbe lifted by a boundary layerin the case of Dirichlet boundary conditions

Theorem Gérard Vaut Paul If pj is independent of knon degeneracy then WE mm of modal solfurthertechnicalassumptionsGeneralization in an arbitrary domain 1 work withlocal coordinates vs Symbole AEC k DNondegeneracy assumption The gig's are independentof x karcade connected componentofthe boundary

Remark Point of view equivalent to butslightly different from matched asymptoticexpansions

Page 6: singular definition Reform - Université Grenoble Alpes

I Example 1 Boundary layers forthe critical

reflection problem ima stratified fluidJoint work with Roberta Bianchini and LaureSaint Raymond following a previous analysisof Thierry Dauxoir and Bill Young

Starting point Bousainesq model

q v SG F v Vp ni Du bez in rB fb is v F b N'vz j.EEf

in r

b buoyancyGN k s N L inthe following tosimplifyr xs xD EN _x singt nz cosy 0

Conditions on 21 7par 0 dub par 0

Purpose ofAndy by Danois Young studyreflection of an incident wave on a sloping

boundary and retrieve some experimentalobservations trapping of reflected wave dosetothe boundary creation of harmonics etc

Page 7: singular definition Reform - Université Grenoble Alpes

Idea if f v k Les then linear plane waves

ofthe in viscid system are approximate solutions tantthey do not satisfy the boundary conditions

Creation of a reflected wave BL in

orderto match the boundary condition Thenonlinear interaction between the incident wave andthe BL creates hamrmonics

Today focus on boundary layers fr43only Forget about nonlinear kun 8 0

Look for solutions of B with 5 0 intheforme

j espfilkx ut d e Ê2 X singtzzasz e normal variable

mary seeing tangential variable

Plug into linear system

As Ç 0

As 4x4 complex matrix

A io auIf.sihI r igxa

Page 8: singular definition Reform - Université Grenoble Alpes

Nontrivial solution det A 0

Polynomial of order 6 ind

del Agarn kw ico vu ki 5 22k singury

à n'n'y d k'Casey d

Criticality parameter 5 co singstudyof the behavior ofthe roots of det Ajdepending on N K and 5Proposition Amine u s Nek Ikki WILLc det Ay has exactly 3roots with so dalpaitii The rooto.ve inthe following way321 n s 151 Eu aireradjipotold d ouiJ'fuitÊtre pas Lynn 3Mo ER M2 Mfrreflectedwave bz k dz v

dadzvi l'a Superposition of two BtN BD Superpositionof one of size n'k

2dB BL of one of size v s

different sizes

iii the BL can always lift 3 Bc

dim Vert Ypg dee det Asa Remo3

Page 9: singular definition Reform - Université Grenoble Alpes

Construction of an approximate solution ira

weakly nonlinear regime1 Take on invisaid linear wave Corwavepacket2 Couture the associated BL correctsthe defect in B c

Approximate solution of linear systemsatisfying exactly the boundary conditions

3 Deal with NL remainder Out ofthe scope of today's lecture

Justification of the approximationEnergy method case verra 13K n'b

HE Vandchlly templeEst NE Vapp trotteCap C d t s vbulb

Remark Both the incident wave packet andthe BL carry an energyof order OCD

Example 2 Ekman layers

starting point rotatingfluid equationwithsmallvertical viscosity

Page 10: singular definition Reform - Université Grenoble Alpes

zu Et VT n È ez ru Tpe Dune EdéhEoinerdiva 0

es Yara1 T O D T R 12oz

Reminder L P Lez r Coriolis operatorE projection auto 1K

1K net r divu quziz.co aziz D

Proposition There exists ce bilbertion basis of 1Kdenokd eDeeezsyoz.cousistimgofeigenvedorsofL

LNG ideaNkdes lez T

Kal

Ref Chemin Desjardins Gallagher Grenier 2006Lions Temaru Wang 1994Babin MahaloriNicolaenko 1999

Nea eikh.dk men b as CTkzznz k simCTkzz

Following Lecture 2 it is natural to lookfor we in the four

Page 11: singular definition Reform - Université Grenoble Alpes

n'Cttenpf_tel v Ect

In other wordsideaE

yegebellt Nn corrector

Eg for beenProblem Nea pz 0

2 I

from a spectral pointof view consider the operatorAsX H defined by AE f we where

Lu Due fdiva 0

near _0 hïmmpatüThen Al is not a normal operator L and

D associated with Dirichlet BD do notAself adjoint

commute

partSource of difficulty inthe analysisresponsible forthe presence of Ekmanlayers

Page 12: singular definition Reform - Université Grenoble Alpes

Back to the general methodologyLookfor a solution offre linear part fin the forme

E

enplika.xa iaotq.dz BÉE to

anticipatethehigh oscillationsthe system

with Re d 0

linear system

A E É 0

where

AE fixe Kat EM

O I o i kaE

det AE polynomial of degree 6in leCI fact polynomial of degree 3in d

Page 13: singular definition Reform - Université Grenoble Alpes

dette Ê tl 1kf fin élite Kelly

f ICI

Investigate the roots of I focussing ontheroots that do not belong to A pet X 42

Remark Inthe présent case since the BL is

created by the trace ofNeon at the boundary itis sufficient to assume that de fork Cha lez C7f Yo fou some lez EZ

Regime I ka 0 IL 1

2 roots X of order E X Ëwhere 2

M t 2in I 2 0

n i u ti

clamsical Ekman layers

1 root oforder 1 Xo No

where I we no Ibert ce 2 0

Desmet correspond to a BL

Page 14: singular definition Reform - Université Grenoble Alpes

spore

Regime 2 ka 0 1

Abit artificial But corresponds totheasymptotie

Kd1 root of order E

2 X 2in2

2 roots of order E X

whereLiu YI Liu Ikari Ikeda 0

quasi resonant BL much larger than classicalEkman layer

Regime 3 kenzo 101412classical Ekman layersf root in exactly zero

Regime 4 ka 0 lol L1 classical Ekman layer2 roots are exactly zero

In Regime si hourdagou deal withthe not which is OLD

classical construction use dire de no thatuB k BÉ nulle 0 at 0

forget about do

Page 15: singular definition Reform - Université Grenoble Alpes

Since UBL is divergence free and When BCis of sine E u etc Ekman pumpingSome fluid exits the BL and enteurs theinterior

Lifted by a couche vlilt OC.es

oscillating at frequency E D

vlilttfezrvh.lt CD

Source term in the eq.fr vefriction terms