Equacions de Navier-Stokes
Mecànica dels medis continus |
---|
![]() |
Lleis
|
|
|
|
Les equacions de Navier-Stokes reben el seu nom de Claude-Louis Navier i George Gabriel Stokes. Es tracta d'un conjunt d'equacions en derivades parcials no lineals que descriuen el moviment d'un fluid. Aquestes equacions governen l'atmosfera terrestre, els corrents oceànics i el flux al voltant de vehicles o projectils i, en general, qualsevol fenomen en el que s'involucrin fluids newtonians.
Aquestes equacions s'obtenen aplicant els principis de conservació de la mecànica i la termodinàmica a un volum fluid. Fent això s'obté l'anomenada formulació integral de les equacions. Per arribar a la seva formulació diferencial es manipulen aplicant certes consideracions, principalment aquella en què els esforços tangencials guarden una relació lineal amb el gradient de velocitat (llei de viscositat de Newton), obtenint d'aquesta manera la formulació diferencial que generalment és més útil per la resolució dels problemes que es plantegen en la mecànica de fluids.
Com ja s'ha dit, les equacions de Navier-Stokes són un conjunt d'equacions en derivades parcials no lineals. No es disposa d'una solució general per a aquest conjunt d'equacions, i excepte certs tipus de flux i situacions molt concretes no és possible trobar una solució analítica; el qual en moltes ocasions hem de recórrer a l'anàlisi numèrica per determinar una solució aproximada. A la branca de la mecànica de fluids que s'ocupa de l'obtenció d'aquestes solucions mitjançant l'ordinador s'anomena dinàmica de fluids computacional (CFD, del seu acrònim anglosaxó Computational Fluid Dynamics ).
Conceptes previs
Derivada substancial o material
Com que generalment adoptem la descripció euleriana, la derivada ordinària ja no representa tota la variació per unitat de temps d'una determinada propietat del fluid seguint la partícula fluida. Això és degut al moviment del fluid. Per reflectir aquesta variació usarem la derivada substancial (o derivada seguint a la partícula fluida). La derivada substancial o derivada material es defineix com l'operador:
On és la velocitat del fluid. El primer terme representa la variació de la propietat en un punt fix de l'espai i per això se l'anomena derivada local, mentre que el segon representa la variació de la propietat associat al canvi de posició de la partícula fluida, i se l'anomena derivada convectiva. Aquest és el procediment que segueix Josep de Echegarai per demostrar la derivada material. Vegeu una demostració de com arribar a una derivada material. Prenent les coordenades d'Euler com:
- .
Calcularem l'acceleració per a aquestes coordenades:
Desenvolupem cada derivada total de cada component, així podrem seguir un desenvolupament fàcil de recordar:
Si sumem terme a termes i traiem factor comú ens adonem que podem factoritzar bastant:
Veiem que la part de les derivades parcials espacials es poden escriure com:
Si ara substituïm velocitat per obtenim formalment l'expressió de la derivada material:
Teorema del transport de Reynolds
Si la derivada substancial permet calcular la variació d'una propietat del fluid seguint una partícula fluida, el teorema del transport de Reynolds permetrà calcular la variació d'una magnitud fluida extensiva lligada a un volum fluid. En la seva forma general, el teorema del transport de Reynolds s'expressa com:
on és una propietat extensiva definida per unitat de volum, és un volum fluid, és un volum de control que coincideix amb en l'instant t, la superfície de control lligada a aquest volum, la velocitat del fluid i la velocitat de la superfície de control.
Expressat en termes col·loquials es pot dir que el teorema del transport de Reynolds ve a dir que la variació d'una propietat extensiva en un volum fluid, és igual a la variació d'aquesta propietat en l'interior d'aquest volum més la quantitat d'aquesta propietat que travessa la superfície del volum.
Teorema de la divergència
El teorema de la divergència (o teorema de Gauss), ens permet transformar integrals de superfície en integrals de volum (i viceversa). En el cas particular de les tres dimensions podem expressar-ho com:
Les equacions de Navier-Stokes
- .
Aquesta expressió representa el principi de conservació del moment lineal aplicada a un fluid general. La llei de la conservació de la massa s'escriu:
En aquestes equacions ρ representa la densitat, ui (i = 1,2,3) les components cartesianes de la velocitat, Fi les forces aplicades sobre el cos, com la gravetat, P la pressió del fluid, i µ la viscositat dinàmica.
on Δ = eij és la divergència del fluid i δij la delta de Kronecker. D/Dt és la derivada total o derivada material temporal seguint el fluid:
La no-linealitat de les equacions es deu precisament al terme relacionat amb la derivada total. Quan µ és uniforme sobretot el fluid les equacions de fluid se simplifiquen de la manera següent:
Casos particulars
Quan la densitat d'un fluid és constant l'anomenem incompressible i l'equació de continuïtat adquireix la forma següent:
Per fluids ideals (viscositat µ nul·la) que siguin alhora incompressibles, les equacions resultants es denominen equacions d'Euler:
En el cas que el fluid sigui incompressible, però a diferència de l'anterior cas, també es tracti d'un fluid viscós, les equacions de Navier-Stokes en coordenades cartesianes són:
Altres consideracions
Les equacion de Navier-Stokes no tenen solució analítica. Per aplicacions pràctiques, aquestes equacions de derivades parcials es ressolen de manera numèrica amb ordinadors i es fan les simulacions corresponents.
Avui dia no existeix cap explicació matemàtica de com un fluid passa de ser un fluid regular a un turbulent. Doncs, hi ha alguna garantia matemàtica que les solucions de les equacions de Navier-Stokes, donades unes condicions inicials del moviment del fluid, existeixin en tot instant futur del temps? La resposta és sí, en un fluid en dues dimensions, però ningú sap respondre a aquesta qüestió en tres dimensions. Aquesta pregunta constitueix un dels Problemes del Mil·lenni que l'Institut de Matemàtiques Clay premia amb 1 milió de dòlars nord-americans a qui pugui resoldre'l.
Ús d'equacions de Navier-Stokes en videojocs
Les equacions de Navier - Stokes s'utilitzen àmpliament en videojocs per modelar una àmplia varietat de fenòmens naturals. Les simulacions de fluids gasosos a petita escala, com foc i fum, sovint es basen en l'article seminal "Real-Time Fluid Dynamics for Games" ("Dinàmica de fluids en temps real per a jocs")[1] de Jos Stam, que elabora un dels mètodes proposats a l'article anterior i més famós de Stam "Stable Fluids" [2] de 1999. Stam proposa una simulació de fluids estable utilitzant un mètode de solució de Navier – Stokes de 1968, juntament amb un esquema de advecció semi-lagrangià incondicionalment estable, com es va proposar per primera vegada el 1992.
Les implementacions més recents basades en aquest treball s'executen a les unitats de processament de gràfics (GPU) dels sistemes de joc en lloc de la unitat central de processament (CPU) i aconsegueixen un grau de rendiment molt més gran.[3][4] S'han proposat moltes millores al treball original de Stam, que pateix inherentment una alta dissipació numèrica tant en velocitat com en massa.
Es pot trobar una introducció a la simulació interactiva de fluids en el curs de l'Association for Computing Machinery SIGGRAPH de 2007, Fluid Simulation for Computer Animation (Simulació de fluids per a animació per ordinador).[5]
Cultura popular
- Aquesta equació és esmentada a la sèrie Young Sheldon. Sheldon diu a la seva àvia que necessita diners per poder comprar un ordinador i així poder resoldre l'equació de Navier-Stokes; ell ho aconsegueix i guanya el premi Nobel a una primerenca edat.
- També és esmentada a la pel·lícula Gifted (2017).
Vegeu també
Referències
- ↑ «Real-Time Fluid Dynamics for Games (Dinámica de fluidos en tiempo real para juegos)» (en inglés). semanticscholar.org, 2003.
- ↑ Stam, Jos. «Stable Fluids (Fluidos estables)» (en inglés). stanford.edu, 1999. Arxivat de l'original el 2019-07-15. [Consulta: 13 març 2022].
- ↑ GPUGems, ed.
- ↑ ShaderX5, ed.
- ↑ «Fluid Simulation for Computer Animation (Simulación de fluidos para animación por computadora)» (en inglés). www.cs.ubc.ca. Arxivat de l'original el 2024-03-30. [Consulta: 27 juny 2024].
Bibliografia
- Acheson, D. J. (1990), Elementary Fluid Dynamics, Oxford Applied Mathematics and Computing Science Series, Oxford University Press, ISBN 978-0-19-859679-0, <https://books.google.com/books?id=IGfDBAAAQBAJ>
- Batchelor, G. K. (1967), An Introduction to Fluid Dynamics, Cambridge University Press, ISBN 978-0-521-66396-0
- Currie, I. G. (1974), Fundamental Mechanics of Fluids, McGraw-Hill, ISBN 978-0-07-015000-3
- V. Girault and P. A. Raviart. Finite Element Methods for Navier–Stokes Equations: Theory and Algorithms. Springer Series in Computational Mathematics. Springer-Verlag, 1986.
- Landau, L. D. & Lifshitz, E. M. (1987), Fluid mechanics, vol. Course of Theoretical Physics Volume 6 (2nd revised ed.), Pergamon Press, ISBN 978-0-08-033932-0, OCLC 15017127
- Polyanin, A. D.; Kutepov, A. M. & Vyazmin, A. V. et al. (2002), Hydrodynamics, Mass and Heat Transfer in Chemical Engineering, Taylor & Francis, London, ISBN 978-0-415-27237-7
- Rhyming, Inge L. (1991), Dynamique des fluides, Presses polytechniques et universitaires romandes
- Smits, Alexander J. (2014), A Physical Introduction to Fluid Mechanics, Wiley, ISBN 0-47-1253499
- Temam, Roger (1984): Navier–Stokes Equations: Theory and Numerical Analysis, ACM Chelsea Publishing, ISBN 978-0-8218-2737-6
- White, Frank M. (2006), Viscous Fluid Flow, McGraw-Hill, ISBN 978-0-07-124493-0
Enllaços externs
- Derivació simplificada de les equacions de Navier-Stokes Arxivat 2017-11-29 a Wayback Machine.
- Pàgina de l'Institut Clay sobre el problema de Navier-Stokes Arxivat 2010-06-13 a Wayback Machine.
- Potter, Merle; Wiggert, David. Mecánica de Fluidos (en espanyol). 3a ed.. Thomson. ISBN 970-686-205-6.