Vraag:
Middelpuntvliedende kracht in de Navier Stokes-vergelijkingen
BarocliniCplusplus
2017-09-22 22:35:35 UTC
view on stackexchange narkive permalink

De Navier-Stokes-vergelijkingen zijn een reeks niet-lineaire differentiaalvergelijkingen die de windsnelheid en -richting diagnosticeren. Ze worden (ongeveer) uitgedrukt als $$ \ frac {d \ vec {u}} {dt} = - \ frac {1} {\ rho} \ nabla P + f \ hat {k} \ times \ vec {u} -g \ hat {k} + \ nu \ nabla ^ 2 \ vec {u} $$ waarbij $ \ frac {d} {dt} = \ frac {\ partieel} {\ partieel t} + u \ frac {\ partieel } {\ partiële x} + v \ frac {\ partiële} {\ partiële y} + w \ frac {\ partiële} {\ partiële z} $, $ f $ is de coriolisparameter, $ \ rho $ is dichtheid, $ P $ is druk, $ g $ is zwaartekracht, $ \ nu $ is de kinematische viscositeit en $ \ vec {u} $ is de windvector.

In synoptische meteorologie wordt geleerd dat in jet strepen beïnvloedt de kromming van de jetstreak de versnelling door de middelpuntvliedende kracht. In de bovenstaande vergelijking zie ik echter geen centrifugale component.

Ontbreekt de op kromming gebaseerde centrifugale beweging gewoon of is deze verborgen? Als het ontbreekt, hoe zou het dan worden uitgedrukt (in cartesiaanse coördinaten)? Als het verborgen is, in welke term is het dan verborgen? Kunnen we het ontleden in centrifugale en niet-centrifugale versnelling? Een afleiding van cyclostrofisch zou voldoende zijn, aangezien het expliciet zal zeggen uit welk deel van de vergelijking kromming is.

@BaroclinicCplusplus eigenlijk heb je een geopotentieel dat is samengesteld uit zwaartekracht + centrifugaal. Je moet de krul nemen van het vectorproduct van de omloopsnelheid van de aarde en de rotatiesnelheid van de aarde. Dit blijkt dan volgens de Grassmann-regel nul te zijn. Dus als de krul van het vectorveld nul is, kan dit worden uitgedrukt als een scalair potentieel. Mijn leerboek bevat het middelpuntvliedende potentieel op de RHS of NS-vergelijking en combineert het met het zwaartekrachtpotentieel om het geopotentieel te krijgen.
Dus hier heb je $$ - g \ hat k $$. Ik neem aan dat dit betekent dat de zwaartekracht naar beneden is gericht. Andere afleidingen hebben een andere benadering door de middelpuntvliedende kracht en de zwaartekracht te combineren en het geopotentiaal te noemen.
http://glossary.ametsoc.org/wiki/Apparent_gravity
@gansub Ik had het niet over een zwaartekrachtaanpassing vanwege de veranderende coördinaten. Ik had het over kromming in stralen die tot een middelpuntvliedende kracht leiden. Kan op dezelfde manier het cyclostrofische evenwicht worden afgeleid uit de N-S-vergelijkingen?
@gansub bijvoorbeeld: http://journals.ametsoc.org/doi/pdf/10.1175/1520-0493%281992%29120%3C2429%3ATEOJSC%3E2.0.CO%3B2
@BaroclinicCplusplus - Ik heb opnieuw gelezen wat ik hierboven heb geschreven. In mijn opmerking hierboven suggereer ik nergens een verandering van coördinaat en bijbehorende zwaartekrachtaanpassing. De eerste opmerking van mij is ook mogelijk in Cartesiaanse coördinaten.
@gansub Schijnbare zwaartekracht is de combinatie van de zwaartekracht van de aarde + de middelpuntvliedende kracht van de aarde. Wat ik bedoelde met "veranderende coördinaten" was de verandering van inertiële naar niet-intertiële referentiekaders. Dit omvat schijnbare zwaartekracht, Coriolis en de andere metrische termen verminderd door schaalanalyse.
@BaroclinicCplusplus - Dan heb ik de vraag die je hebt gesteld beantwoord. In cartesische coodinaten in een niet-inertiaal referentiekader is de middelpuntvliedende kracht aanwezig als onderdeel van de schijnbare zwaartekracht, op voorwaarde dat je $$ - g \ hat k $$ vervangt door $$ - \ nabla \ phi_a $$ waarbij $$ \ phi_a $$ is het geopotentieel.
@gansub Mijn vraag had geen betrekking op inertiële en niet-inertiële referentiekaders. Mijn vraag had betrekking op de kromming van de stroming, die aanzienlijk kleiner van schaal is. Zie http://www.theweatherprediction.com/advanced/habyextra7/.
@BaroclinicCplusplus - dus je hebt een coördinatensysteem nodig om de kromming uit te drukken. Vraagt ​​u zich af of de kromming van de stroom ook in cartesiaanse coördinaten kan worden uitgedrukt?
Herschik de momentumvergelijking van Cartesiaans (x, y) naar polair (r, theta) coördinatensysteem. Centrifugale versnelling zal een van de termen zijn in du / dt.
@BarocliniCplusplus Ik kan een cylostrofisch evenwicht afleiden uit de NS-vergelijking. Zou je dat als antwoord willen zien?
@BarocliniCplusplus kunt u de vraag bewerken om dat feit weer te geven? Ik zou niet willen antwoorden zonder de expliciete vermelding in de vraag
@BarocliniCplusplus - de afleiding is in dit boek - https://books.google.co.in/books?id=7IsmskcU3msC&pg=PA93&lpg=PA93&dq=material+derivative+natural+coordinate&source=bl&ots=5u9HHDBiQ_&sig=_OYZYvbZRHJa8tvyeGZWPn2cHNE&hl=en&sa=X&ved=0ahUKEwi4nJbv_JLXAhVLuo8KHYK_BaQQ6AEITzAH # v = onepage & q = material% 20derivative% 20natural% 20coordinate & f = false. Weet je toevallig wat het ESSE-beleid is met betrekking tot afleidingen uit een tekstboek, kunnen we verwijzen naar auteursrechtelijk beschermd materiaal? Zo ja, dan kan ik een afleiding laten zien, maar de bron is het boek.
@gansub Hoewel ik niet zeker ben van het auteursrechtbeleid (kan een goede metavraag vormen), ben ik er vrij zeker van dat u kunt verwijzen naar auteursrechtelijk beschermd materiaal, zolang u maar een citaat maakt, aangezien verwijzingen naar gepubliceerde boeken regelmatig voorkomen in de academische wereld.
Twee antwoorden:
AtmosphericPrisonEscape
2019-04-17 05:14:45 UTC
view on stackexchange narkive permalink

Op basis van de stand van de discussie in de commentaren wist ik niet zeker of dit nu een antwoord heeft, maar ik struikelde erover en ik denk dat ik een kort antwoord kan geven waarvan men veel kan leren.

Ik zal mijn discussie baseren op polaire $ (\ rho, \ phi) $ coördinaten, voor de duidelijkheid. In de klassieke mechanica gebruiken we voor de positievector $$ \ vec r = r \ vec {e} _r = r \; (\ cos \ phi, \ sin \ phi), $$ omdat $ \ punt {\ vec {e}} _ r = \ punt {\ phi} \; (- \ sin \ phi, \ cos \ phi) = \ punt {\ phi} \; \ vec {e} _ {\ phi} $ . In 3D sferische coördinaten hebben we nog een paar van die identiteiten nodig, maar de algebra blijft hetzelfde. Hieruit volgt $$ \ punt {\ vec r} = \ punt {r} \ vec {e} _r + r \ punt {\ phi} \ vec {e} _ {\ phi} $$ en tot slot $$ \ ddot {\ vec r} = (\ ddot {r} -r \ dot {\ phi} ^ 2) \ vec { e} _r + (r \ ddot {\ phi} + 2 \ punt {r} \ punt {\ phi}) \ vec {e} _ {\ phi} $$ waarbij we al rekening hebben gehouden met de verandering van coördinaten, evenals eenheidsvectoren. De volgende stap zou $ \ dot \ phi $ ontleden als $ \ dot \ phi = \ Omega + u / r $ , terwijl $$ r \ dot \ phi ^ 2 = r \ Omega ^ 2 + 2 \ Omega u + u ^ 2 / r $$ waar we identificeren (in volgorde van de termen) 1.) de middelpuntvliedende term, 2.) de coriolis-term en 3.) de "kromming / metrische" term. Deze procedure werkt ook wonderwel in sferische coördinaten en men herstelt perfect de termen die bijvoorbeeld Holton of Vallis in hun boeken krijgen via langdurige en ondoorzichtige 'proof-by-picture'-methoden.

Ik zeg het laatste omdat nu trek ik conclusies voor de volledige vergelijkingen, uit mijn afleidingen die werden gedaan uit de onvolledige: Conclusies:

  1. Centrifugale termen zijn niet wat meteorologen 'kromming / metrische' termen noemen. Ze kunnen hetzelfde lijken, als men slordig is en $ \ Omega = vr $ schrijft. Of wat vaak wordt gedaan, is om de baan van een satelliet te berekenen om de middelpuntvliedende kracht en de zwaartekracht gelijk te stellen, terwijl men $ u ^ 2 / r $ gebruikt, wat in feite is de 'kromming' niet de middelpuntvliedende term. De fysiek metrische term is ook al $ r \ dot {\ phi} ^ 2 $ , aangezien deze voortkomt uit een gebogen coördinatensysteem.
  2. Zoals aangegeven in de commentaren, is de middelpuntvliedende term $ r \ Omega ^ 2 $ gewoonlijk verborgen in de meteorologie met de zwaartekracht. We vermoeden dat dit ook in uw verzameling gebeurt, omdat de Coriolis-term verschijnt, wat erop duidt dat de bovenstaande berekening is uitgevoerd volgens de eerste orde en dat de middelpuntvliedende term dus ergens moet zweven.
  3. De krommingstermen zijn nergens te vinden in uw set vergelijkingen. Dit komt waarschijnlijk doordat de termen simpelweg worden verwaarloosd. Ook als ze in hun volledige vorm in sferische coördinaten worden bekeken, lijken ze nogal rommelig en ik vond geen gemakkelijke manier om ze uit te drukken door een mooie vector- of matrixbewerking. Daarom heb ik het sterke vermoeden dat mensen bij het geven van de vectoriële vorm van de NS-vergelijkingen deze meestal verwaarlozen en ze alleen toevoegen als ze de NS-vergelijkingen in componentvorm presenteren.
  4. De bovenstaande minimalistische afleiding maakt duidelijk dat de kromming , zijn centrifugale en Coriolis-termen afkomstig van een gekromd coördinatensysteem. Elk gekromd coördinatensysteem is noodzakelijkerwijs niet-traag. Dus in Cartesiaanse coördinaten, die traag zijn, moeten die termen noodzakelijkerwijs verdwijnen.
Ik leidde eigenlijk de middelpuntvliedende kracht of de middelpuntzoekende kracht af in een http://tornado.sfsu.edu/Geosciences/classes/e260/CoordinateSystems/NaturalCoordinateSystem.html
@gansub: Welnu, in natuurlijke coördinaten laat men vaak de bespreking van de mate van verandering van de coördinatenbasis achterwege, of hoe daardoor nieuwe termen worden geïntroduceerd. En ik ben niet duidelijk op welke, middelpuntvliedende of kromming, in deze afleiding overeenkomt met de middelpuntvliedende kracht in natuurlijke coördinaten.
Nee. Er is een versie van het afleiden van de NS-vergelijkingen met behulp van de https://en.wikipedia.org/wiki/Frenet%E2%80%93Serret_formulas waarin de mate van verandering van de coördinatenbasis wordt beschouwd in differentiële meetkunde. Dus als je transformeert van cartesiaanse naar natuurlijke coördinaten, krijg je de middelpuntzoekende kracht. Bekijk Dynamics of the Atmosphere door Bott Chapter 8.
@gansub: "in welke mate van verandering van de coördinatenbasis", welnu, wat Frenet de middelpuntvliedende krachten $ V ^ 2 / R $ noemt, is gewoon onze $ r \ punt {\ phi} ^ 2 $. Ik zie geen probleem? Of weet u niet zeker waar uw opmerking over gaat?
Ik voeg een parallel antwoord toe. We kunnen daarna discussiëren. Daar heb je gelijk in. Overigens ben ik het niet oneens met uw antwoord. Ik heb erover gestemd.
gansub
2019-04-17 21:43:54 UTC
view on stackexchange narkive permalink

Ik heb beide referenties - The effect of jet streak curvature on kinematic fields and the background reference in that paper - Isolatie van de traagheidszwaartekrachtcomponent in een niet-lineair atmosferisch model en geen van deze verwijzingen vermeldt enige verbinding met een op Cartesiaans gebaseerde NS en het definiëren van kromming in die coördinatenruimte. Voor mij kan het krommingsprobleem in de meteorologie het best worden bestudeerd (voor gelokaliseerde luchtstroom) in een natuurlijk coördinatenstelsel, wat zijn natuurlijke coördinaten?

Gezien het feit dat uitgangspunt Ik wil het hebben over de formule van Frenet Serret. In deze context heeft het natuurlijke coördinatensysteem drie orthogonale basisvectoren ( $ e_s, e_n, e_z) $ waarbij $ e_s $ is de richting van de horizontale windvector, dwz $$ v_h = | v_h | * e_s $$ en $ e_s $ is gewoon de bekende Streamline en $ e_n $ is normaal voor de stroomlijn. We moeten ook de kromtestraal introduceren die wordt gebruikt in de Frenet Serret-formules. Opgemerkt moet worden dat de kromming waarover in dit antwoord wordt gesproken de extrinsieke kromming

is. Vanwege de orthogonaliteit $ e_s \ times e_n = e_z $

Ik zal eerst de vraag van AtmosphericPrisonEscape behandelen en vervolgens de vraag van OP beantwoorden.

De formule van Frenet Serret omvat de veranderingssnelheid van basisvectoren. Kun je vanuit het perspectief van een leek uitleggen wat de veranderingssnelheid van basisvectoren is? Nou, meteorologische problemen kunnen worden gedefinieerd in een cartesiaanse basis (x, y, z) die "vast" is ten opzichte van de vaste sterren. Dit systeem kan dus worden beschouwd als een traagheidssysteem voor een aardgebonden waarnemer. In dit systeem zijn er geen "fictieve" krachten die ontstaan ​​door de specifieke keuze van het coördinatensysteem. Details over dit "vaste" traagheidsframe worden gegeven in de referentie die aan het einde van het antwoord in hoofdstuk 1 wordt aangehaald.

Aangezien meteorologische waarnemingen worden uitgevoerd in een roterend coördinatensysteem, kunnen ze het beste worden beschreven in een coördinatensysteem die meedraait met de aarde. Dus het Cartesiaans systeem beweegt van punt naar punt, de basisvector blijft vast en varieert niet. OTOH in een aardachtig systeem (zoals een bol) beschrijven we over het algemeen in termen van kromlijnige coördinaten waarvan de basisvectoren van punt tot punt variëren.

Wiskundig in een stationair Cartesiaans systeem $ \ frac {de_x} {dx} = 0 $ omdat de basisvectoren onafhankelijk zijn van de positie.

In een algemeen kromlijnig systeem laten we deze aannames los en introduceren we het idee van kromming. Dus vanuit een Frenet Serret-perspectief

$ \ frac {\ partiële e_s} {\ partiële n} = \ frac {e_n} {R_n} $ en soortgelijke andere termen.

Nu om de vraag van OP te beantwoorden.

Als we deze vergelijking nemen en het differentiëren naar tijd $$ v_h = | v_h | * e_s $$

we krijgen

$$ \ frac {dv_h} {dt} = e_s \ frac {d | v_h |} {dt} + | v_h | \ frac {de_s} {dt} $$

Zoals je hier kunt zien, maak ik expliciet onderscheid tussen de basisvector $ e_s $ tov tijd. Dus de eerste term op de RHS staat bekend als de tangentiële versnelling (omdat deze in de richting van de stroomlijn is) en de tweede term na een gebruik van de formule van Frenet Serret (en het aanpassen van enkele termen) kan worden gelijkgesteld aan

$$ | v_h | \ frac {de_s} {dt} = e_n \ frac {v ^ 2_h} {R_t} $$ en $ R_t $ is de straal van Traject

De tweede term staat bekend als de centripetale versnelling en staat loodrecht op de stroomlijn.

Nu gewapend met deze informatie kunnen we proberen alle vragen van OP te beantwoorden.

1) De kromming van jetstreaks kan worden gebruikt om gebieden met samenvloeiing te bestuderen. en diffluence. Het enige wat je hoeft te doen is de divergentie van de horizontale windvector nemen. Houd er rekening mee dat de helling en de windvector moeten worden uitgedrukt in natuurlijke coördinaten. In het diffluentiegebied is de kromtestraal positief en in het samenvloeiingsgebied negatief.

2) Eindelijk kan de kwestie van Cyclostrofisch evenwicht gemakkelijk worden afgeleid door de NS-vergelijking om te zetten in natuurlijke coördinaten, zoals weergegeven in dit ESSE-antwoord Oog van een tornado

Ik ga het niet afleiden, maar geef OP voldoende aanwijzingen zodat OP het kan afleiden. De normale component van de NS-vergelijking wordt gelijkgesteld aan de normale component van de drukgradiënt. Ik heb de vergelijking voor de normale component van de bovenstaande NS-vergelijking al gegeven. Als we die twee gelijk stellen, krijgen we de cyclostrofische balans zoals hieronder weergegeven

$$ | v_h | \ frac {de_s} {dt} = e_n \ frac {v ^ 2_h} {R_t} $$

$$ \ frac {v ^ 2_h} {R_t} = - \ frac {1} {\ rho } \ frac {\ partiële p} {\ partiële n} $$

Zoals men kan zien, beïnvloedt de kromming van de jetstreak de versnelling door de middelpuntzoekende kracht.

Referenties - Dynamiek van de atmosfeer: een cursus in theoretische meteorologie Hoofdstuk 8.

Goed, ik begrijp wat je nu bedoelt. Je denkt dat de 'kromming' simpelweg verwijst naar de straal van de baankromming. Ja, dat hangt af van waar OP's naar vroeg. 2 problemen die ik heb met dit antwoord: 1.) ijk-coördinaten zijn een al lokaal meedraaiend systeem, daarom zijn ze niet traag en niet identiek aan het traagheidssysteem $ \ vec {e} _x, \ vec {e} _y , \ vec {e} _z $. Volgens Newton zijn alleen die coördinaten traag, wat kan worden uitgedrukt als $ \ vec e_ {whatevr} = x_0 + a_0 t $ (de beroemde lineaire beweging), want dan $ \ ddot {\ vec e} _ {whatevr} = 0 $ en Newtons F = ma geldt in al die traagheidsframes.
2.) Het andere probleem dat ik heb, is dat volgens ten minste de auteursrechtregels in wetenschappelijke artikelen, en daarom ben ik er zeker van dat SE iets soortgelijks redelijk heeft, dat als je iets alleen terugkrijgt, je het altijd kunt presenteren. Er is geen copyright op afleidingen, in het bijzonder zou dit ene boek niet het copyright op de afleiding hebben, maar degene die het het eerst heeft gedaan, die allang dood is, en daarom zou het auteursrecht hoe dan ook vervallen.
@AtmosphericPrisonEscape Ik zal uw bezorgdheid over inertiële en niet-inertiële referentiekaders in een update behandelen. Tweede punt Ik kan de afleiding voor de NS-vergelijking presenteren, maar ik denk dat OP al genoeg materiaal heeft om zelf op te zoeken en af ​​te leiden. Het is heel eenvoudig.


Deze Q&A is automatisch vertaald vanuit de Engelse taal.De originele inhoud is beschikbaar op stackexchange, waarvoor we bedanken voor de cc by-sa 3.0-licentie waaronder het wordt gedistribueerd.
Loading...