Jump to content

User:GuntarsKeks/sandbox

From Wikipedia, the free encyclopedia
File:Https://www.researchgate.net/profile/Heng-Xiao-7/publication/325929714/figure/fig19/AS:675547747192833@1538074542320/A-schematic-representation-of-the-hierarchy-of-turbulence-modeling-approaches-based-on W640.jpg
Boeing supercomputing

Turbulences modelēšana ir šķidrumu un gāzes mehānikas skaitlisko aprēķinu (Computational Fluid Dynamics, CFD) apakšnozare kas ļauj aprakstīt turbulences parādības plūsmā. CFD nozares pamatā ir Navje-Stoksa (NS) vienādojums, kas raksturo plašu plūsmu klāstu ka arī, principā, arī satur informāciju par turbulenci. Lai izšķirtu smalkas turbulences paradības, skaitliska modeļa režģim jābūt pietiekami smalkam un tadu metodi sauc par Direct Numerical Simulation (DNS), līdz ar ko var pateikt ka DNS tipa aprēķiniem turbulences modelēšana nav nepieciešama. Līdz ar modeļa augstu izšķirtspēju palielinās dator resursi un uz simulāciju patērētais laiks, padarot šada tipa apreķinus nepraktiskus. Tas fakts mudina meklēt jaunas metodes turbulences parādību aprakstīšanai, netieši, izmantojot mazāk resursu.

Reinoldsa Videjotais Navje Stoksa vienādojums

[edit]

Turbulences parādības var noverot mums visapkārt, piemēram ka virpuļi upes straumē un ceļošos dūmos. Atšķirīgas pazīmes- plūsma ir laika mainīga, neregulāra un haotiska. Analizējot tādu plūsmu laboratorijā var pamanīt ka kādā punktā ātrums ir vidēji regulārs, bet papildus pastāv šķietami patvaļīgas fluktuācijas.

Reinoldsa dekompozīcija

[edit]

Šīs veids ātrumu (ātruma i-to komponenti) apskatīt ka kādu vidēju vērtību un statistiska rakstura fluktuāciju sauc par Reinoldsa dekompozīciju (Reynolds decomposition)

un balstoties uz tā, var izvest Reinoldsa Videjoto Navje Stoksa (Reynolds Avereged Navier Stokes, RANS) vienādojumu. Primāri apskata tādas fluktuācijas lai to laika vidējas vertības būtu nulles, kas novēd pie sekojošas RANS vienādojuma formas, pierakstīto izmantojot Einšteina indeksus, Dekartā koordināšu sistēmā:

Plūsmas apraksta pilnvertigums (noslegšanas problēma)

[edit]

RANS ir vienādojums vidējām ātruma komponentēm un spiedienam. Papildus plūsmas aprakstam izmanto vai nu nepartrauktības vienādojumu, vai Pusasona vienādojumu spiedienam. Tas veido 4 (3+1) vienādojumus, bet RANS satūr ari papildus nelineāru termu (pēdējais sakaitamais) , kas ir atkarīgs no nezināmam fluktuācijām. Šo saskaitāmo sauc par Reinodlsa stressu un izvest no apakšprincipien to nav iespējams. Tā kā vienādojumu skaits ir mazaks par nezināmo skaitu, vienozīmīgo atrisinājumu atrast nav iespējams. Tāpēc pastāv vairāki modeļi (turbulences modeļi) kas mēģina Reinodlsa stressu saistīt to ar citiem nezināmiem (noslēgt sistēmu).

Virpuļu viskozitāte

[edit]

Virpuļu viskozitātes modeļu klasse balstās uz Buseneska (Joseph Valentin Boussinesq) hipotēzi[1] par Reinoldsa stresa terma saistību ar ātruma gradientiem.

Kur

ir "turbulenta viskozitāte"
ir turbulenta kinetiska enerģija
and ir Kronekera simbols.

Šo izteiksmi var pierakstīt īsāk sadalot izmantojot deviatorisko daļu no bīdes ātruma tenzora

ka

Saskaņa ar šo pieņemumu, sakotnēji nelineāru stressa termu var aizvietot ar lineāru vienādojumu ar proporcionalitātes koeficientu un turbulēto kinētisko enerģiju (TKE) . Pēc savas uzbūves jaunais saskaitamais atgādina NS vienādojuma difūzijas termu ar viskozitātes koeficientu , līdz ar ko koeficientu tipiski sauc par “turbulentu viskozitāti” vai “virpuļu viskozitāti”. vērtība var mainīties telpā atkarībā no plūsmas parametriem, robežnosacījumiem un ģeometrijas. Lai to noteiktu tiek izstrādāti modeļi, kuros tipiski izmantojot empīriskas sakarības ar koeficientiem noteiktiem no eksperimentiem. Viskozitātes modeļu grupa sadalas 3 tipos, pēc nepieciešamu papildus vienadijumu skaita:

  1. Algebriski (nulles vienādojumu modeļi)- turbulentā viskozitāte tiek modificēta ņemot vērā tikai plūsmas parametrus. Plūsmas vesture nespele lomu. Piemeram The Baldwin-Lomax un Cebeci-Smith modeļi.
  2. Viena vienādojuma modeļi- tiek iekļauts kāda lauka transporta vienādojums, piemēram TKE. Popular piemērs ir Spalart-Allmaras turbulences modelis.
  3. Divu vienādojumu modeļi- papildus TKE transportam tiek risināts turbulentas dissipacijas (transporta) vienādojums. Piemeram turbulences modelis un turbulences modelis.

Common models

[edit]

The following is a brief overview of commonly employed models in modern engineering applications.

The Spalart–Allmaras model[2] is a one-equation model that solves a modelled transport equation for the kinematic eddy turbulent viscosity. The Spalart–Allmaras model was designed specifically for aerospace applications involving wall-bounded flows and has been shown to give good results for boundary layers subjected to adverse pressure gradients. It is also gaining popularity in turbomachinery applications.[citation needed]

K-epsilon (k-ε) turbulence model[3] is the most common model used in computational fluid dynamics (CFD) to simulate mean flow characteristics for turbulent flow conditions. It is a two-equation model which gives a general description of turbulence by means of two transport equations (PDEs). The original impetus for the K-epsilon model was to improve the mixing-length model, as well as to find an alternative to algebraically prescribing turbulent length scales in moderate to high complexity flows.

In computational fluid dynamics, the k–omega (k–ω) turbulence model[4] is a common two-equation turbulence model that is used as a closure for the Reynolds-averaged Navier–Stokes equations (RANS equations). The model attempts to predict turbulence by two partial differential equations for two variables, k and ω, with the first variable being the turbulence kinetic energy (k) while the second (ω) is the specific rate of dissipation (of the turbulence kinetic energy k into internal thermal energy).

SST (Menter's shear stress transport) turbulence model[5] is a widely used and robust two-equation eddy-viscosity turbulence model used in computational fluid dynamics. The model combines the k-omega turbulence model and K-epsilon turbulence model such that the k-omega is used in the inner region of the boundary layer and switches to the k-epsilon in the free shear flow.

The Reynolds stress equation model (RSM), also referred to as second moment closure model,[6] is the most complete classical turbulence modelling approach. Popular eddy-viscosity based models like the k–ε (k–epsilon) model and the k–ω (k–omega) models have significant shortcomings in complex engineering flows. This arises due to the use of the eddy-viscosity hypothesis in their formulation. For instance, in flows with high degrees of anisotropy, significant streamline curvature, flow separation, zones of recirculating flow or flows influenced by rotational effects, the performance of such models is unsatisfactory.[7] In such flows, Reynolds stress equation models offer much better accuracy.[8]

Eddy viscosity based closures cannot account for the return to isotropy of turbulence,[9] observed in decaying turbulent flows. Eddy-viscosity based models cannot replicate the behaviour of turbulent flows in the Rapid Distortion limit,[10] where the turbulent flow essentially behaves like an elastic medium.[11]

Sienas funkcijas

[edit]
“Sienas likums”: ātruma sadalījums pie sienas un aproksimācijas funkcijas

Pielietojot plaši izmantotu “pielipšanas robežnosacījumu” (vides nulles ātrumu pie sienas, ), gadījumos ja vides ātrums ārpus sienas nav nulle, var sagaidīt lielus ātruma gradientus sienas normāles virzienā. Kas savukārt izraisa turbulences parādības. Lai izšķirtu parādības pie sienas skaitliskam modeļa režģim šeit jābūt pietiekami smalkam, kas palielina slodzi uz resursiem. Viens risinajums būtu, zinot ātruma sadalījumu attiecīgas plūsmas veidam, modificēt robežslāņa izteiksmes lai procesi tiktu nevis izšķirti telpisi, bet modelēti sienas tuvāka elementa ietvaros. Aktuāli pieminēt ka Launder un Spalding [12] 1970jos gados piedāvāja izmantot daļēji logaritmisko ātruma sadalījumu (“Sienas likumu”) turbulences modeli. Implementācija ietver vairākas nezināmas konstantes, kuras piemeklētas labākai liela klāsta problēmu atbilstībai un joprojam tiek izmantotas programmas ka OpenFOAM[13] un Ansys FLUENT. CFD ietvaros “Sienas likuma” sadalījums, izteikts bezdimensionālās vienībās , var būt aproksimēts ar funkciju/-am. vērtības nav zināmas apriori un var mainīties sienas garumā un mainoties plūsmai. Tapēc ideāli ja aproksimācijas funkcija labi atbilstu visā garumā. Piemēram Spaldinga sienas funckija. Tomēr izplatītākais paņēmiens ir izmantot divas funkcijas kas labi apraksta sadalījumu viskoza slānī () un logaritmiska slānī (). Elementārākais paņēmiens ir parslēgties starp abam funkcijam (vai citu vērtību pieņemtu par piemērotu parēju) apgabalā. Šeit neviena no funkcijām neapraksta sadalījumu pietiekami precīzi. Lai novērstu aproksimācijas neatbilstību bufera slāņi, divu funkciju gadījumā, šajā intervālā tiek pielietota abu funkciju “sajaukšana” (blending).

Prandtl's mixing-length concept

[edit]

Later, Ludwig Prandtl introduced the additional concept of the mixing length,[14] along with the idea of a boundary layer. For wall-bounded turbulent flows, the eddy viscosity must vary with distance from the wall, hence the addition of the concept of a 'mixing length'. In the simplest wall-bounded flow model, the eddy viscosity is given by the equation:

where:
is the partial derivative of the streamwise velocity (u) with respect to the wall normal direction (y);
is the mixing length.

This simple model is the basis for the "law of the wall", which is a surprisingly accurate model for wall-bounded, attached (not separated) flow fields with small pressure gradients.

More general turbulence models have evolved over time, with most modern turbulence models given by field equations similar to the Navier–Stokes equations.

Smagorinsky model for the sub-grid scale eddy viscosity

[edit]

Joseph Smagorinsky was the first who proposed a formula for the eddy viscosity in Large Eddy Simulation models,[15] based on the local derivatives of the velocity field and the local grid size:

In the context of Large Eddy Simulation, turbulence modeling refers to the need to parameterize the subgrid scale stress in terms of features of the filtered velocity field. This field is called subgrid-scale modeling.

Spalart–Allmaras, k–ε and k–ω models

[edit]

The Boussinesq hypothesis is employed in the Spalart–Allmaras (S–A), k–ε (k–epsilon), and k–ω (k–omega) models and offers a relatively low cost computation for the turbulence viscosity . The S–A model uses only one additional equation to model turbulence viscosity transport, while the k–ε and k–ω models use two.


References

[edit]

Notes

[edit]
  1. ^ "CFD Direct". Boussinesq eddy viscosity assumption.{{cite web}}: CS1 maint: url-status (link)
  2. ^ Spalart, P.; Allmaras, S. (1992). "A one-equation turbulence model for aerodynamic flows". 30th Aerospace Sciences Meeting and Exhibit, AIAA. doi:10.2514/6.1992-439.
  3. ^ Hanjalic, K.; Launder, B. (1972). "A Reynolds stress model of turbulence and its application to thin shear flows". Journal of Fluid Mechanics. 52 (4): 609–638. doi:10.1017/S002211207200268X.
  4. ^ Wilcox, D. C. (2008). "Formulation of the k-omega Turbulence Model Revisited". AIAA Journal. 46: 2823–2838. doi:10.2514/1.36541.
  5. ^ Menter, F. R. (1994). "Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications" (PDF). AIAA Journal. 32 (8): 1598–1605. doi:10.2514/3.12149.
  6. ^ Hanjalić, Hanjalić; Launder, Brian (2011). Modelling Turbulence in Engineering and the Environment: Second-Moment Routes to Closure.
  7. ^ Mishra, Aashwin; Girimaji, Sharath (2013). "Intercomponent energy transfer in incompressible homogeneous turbulence: multi-point physics and amenability to one-point closures". Journal of Fluid Mechanics. 731: 639–681. Bibcode:2013JFM...731..639M. doi:10.1017/jfm.2013.343.
  8. ^ Pope, Stephen. "Turbulent Flows". Cambridge University Press, 2000.
  9. ^ Lumley, John; Newman, Gary (1977). "The return to isotropy of homogeneous turbulence". Journal of Fluid Mechanics. 82: 161–178. Bibcode:1977JFM....82..161L. doi:10.1017/s0022112077000585.
  10. ^ Mishra, Aashwin; Girimaji, Sharath (2013). "Intercomponent energy transfer in incompressible homogeneous turbulence: multi-point physics and amenability to one-point closures". Journal of Fluid Mechanics. 731: 639–681. Bibcode:2013JFM...731..639M. doi:10.1017/jfm.2013.343.
  11. ^ Sagaut, Pierre; Cambon, Claude (2008). Homogeneous Turbulence Dynamics.
  12. ^ Launder, B.E.; Spalding, D.B. (1974). "The numerical computation of turbulent flows". Computer Methods in Applied Mechanics and Engineering. 3 (2). Elsevier BV: 269–289. doi:10.1016/0045-7825(74)90029-2. ISSN 0045-7825.
  13. ^ "OpenFOAM: User Guide k-epsilon". OpenFOAM.{{cite web}}: CS1 maint: url-status (link)
  14. ^ Prandtl, Ludwig (1925). "Bericht uber Untersuchungen zur ausgebildeten Turbulenz". Zs. Angew. Math. Mech. 2.
  15. ^ Smagorinsky, Joseph (1963). "Smagorinsky, Joseph. "General circulation experiments with the primitive equations: I. The basic experiment". Monthly Weather Review. 91 (3): 99–164. doi:10.1175/1520-0493(1963)091<0099:GCEWTP>2.3.CO;2.

Other

[edit]
  • Absi, R. (2019) "Eddy Viscosity and Velocity Profiles in Fully-Developed Turbulent Channel Flows" Fluid Dyn (2019) 54: 137. https://doi.org/10.1134/S0015462819010014
  • Townsend, A.A. (1980) "The Structure of Turbulent Shear Flow" 2nd Edition (Cambridge Monographs on Mechanics), ISBN 0521298199
  • Bradshaw, P. (1971) "An introduction to turbulence and its measurement" (Pergamon Press), ISBN 0080166210
  • Wilcox C. D., (1998), "Turbulence Modeling for CFD" 2nd Ed., (DCW Industries, La Cañada), ISBN 0963605100

Category:Turbulence Category:Turbulence models