Utafiti wa kitaaluma, lugha inayoeleweka

Verianla | Akademik Araştırmalardan Türkçe Ekonomi ve Bilim İçerikleri

27 Septemba 2026, Jumapili
VERİANLAUchapishaji huru wa sayansi
Fungua au funga menyu
...
Home / Sayansi Tumizi / Uhandisi / Uhifadhi wa Nishati katika Simuleringi za Mtiririko wa Hypersonic: Visuluhishi vya Riemann Vilivyorekebishwa Vilidumisha Enthalpy Jumla kwa Usahihi wa Mashine
Uhandisi

Uhifadhi wa Nishati katika Simuleringi za Mtiririko wa Hypersonic: Visuluhishi vya Riemann Vilivyorekebishwa Vilidumisha Enthalpy Jumla kwa Usahihi wa Mashine

Mbinu za nambari zinazokokotoa mtiririko unaozunguka vyombo vya hypersonic zinapaswa kukamata mawimbi makali ya mshtuko, kudumisha msongamano na shinikizo kuwa chanya, kufanya kazi kwa uthabiti na kuwakilisha sheria za uhifadhi wa fizikia kwa usahihi kadiri inavyowezekana.

30/07/2026  Veri Anla Imetazamwa mara 25
Uhifadhi wa Nishati katika Simuleringi za Mtiririko wa Hypersonic: Visuluhishi vya Riemann Vilivyorekebishwa Vilidumisha Enthalpy Jumla kwa Usahihi wa Mashine

Mbinu za nambari zinazokokotoa mtiririko unaozunguka vyombo vya hypersonic zinapaswa kukamata mawimbi makali ya mshtuko, kudumisha msongamano na shinikizo kuwa chanya, kufanya kazi kwa uthabiti na kuwakilisha sheria za uhifadhi wa fizikia kwa usahihi kadiri inavyowezekana. Hata hivyo, visuluhishi vya aina ya Godunov vinavyotumika sana kama Rusanov, HLL, HLLC na Roe vinaweza kushindwa kuhifadhi kwa njia ya nambari enthalpy jumla ambayo inapaswa kubaki thabiti kando ya mistari ya mtiririko katika mtiririko thabiti wa Euler. Mabadiliko bandia ya enthalpy yanayotokea mtiririko unapopita kwenye shock yanaweza baadaye kusafirishwa hadi maeneo mengine ya uwanja wa mtiririko na kuvuruga mgawanyo wa nishati na msongamano.

Katika utafiti, sehemu ya nishati ya Simple Riemann solvers zilizotengenezwa na Gallice ilibadilishwa na kuunda mipango ya one- na multidimensional inayohifadhi enthalpy iitwayo MGallice na MGallice-2D. Mabadiliko yanategemea Häenel condition, ambayo huhakikisha numerical diffusion katika energy flux ni sawa na diffusion katika mass flux iliyozidishwa kwa total enthalpy.

Katika jaribio la channel yenye triangular bumps kwa Mach 2, schemes zote nne zilifikia steady solution baada ya takribani 4.000 iterations. Classical Gallice solvers zilitengeneza small but clear jumps katika total enthalpy wakati wa kuvuka stationary shocks. Modified MGallice solvers, kwa upande mwingine, zilidumisha inlet total enthalpy value ya 6,3 katika entire computational domain.

Katika jaribio gumu zaidi la mtiririko kuzunguka half-cylinder kwa Mach 20, one-dimensional solvers zilikumbwa na shock instability inayojulikana kama carbuncle. Multidimensional solvers zilizuia instability hii; hata hivyo classical Gallice-2D ilikokotoa vibaya density distribution katika stagnation point kwa sababu ya enthalpy loss nyuma ya shock. MGallice-2D ilidumisha total enthalpy constant hadi machine precision na kurejesha expected density field.

Study inaonyesha kwamba preservation ya total enthalpy katika steady hypersonic Euler calculations si theoretical property tu; inaweza kubadilisha moja kwa moja energy na density fields nyuma ya strong shocks. Hata hivyo entropy stability ya proposed schemes haijathibitishwa mathematically, na method bado haijajaribiwa kwenye viscous Navier–Stokes problems.

Kwa nini numerical solution ya hypersonic flows ni ngumu?

Mbele ya vehicle inayosafiri kwa hypersonic speed, strong detached bow shock huunda. Across shock, pressure, density na temperature hubadilika sana ndani ya distance fupi. Expansion waves hutokea nyuma ya vehicle, huku velocity na temperature gradients zikijitokeza karibu na surface.

Real hypersonic flows zinaweza kujumuisha viscosity, heat conduction, chemical reactions, thermochemical non-equilibrium effects, turbulence, surface heating na ablation. Study haishughulikii system hii yote, bali inviscid part ya compressible Navier–Stokes equations, yaani Euler equations, inayohusika na shock formation.

Hata numerical scheme ikikamata strong shock stably, inaweza kuzalisha artificial change katika quantity inayopaswa kuhifadhiwa physically. Main problem inayolengwa na research ni numerical non-preservation ya total enthalpy katika steady Euler flows.

Total enthalpy ni nini?

Katika two-dimensional Euler equations, conserved variables ni density, momentum na total energy:

\[ \mathbf{U}= \begin{pmatrix} \rho \\ \rho\mathbf{u} \\ \rho e \end{pmatrix} \]

Hapa ρ ni density, u ni velocity vector na e ni total energy per unit mass. Total energy ni sum ya internal na kinetic energy:

\[ e=\varepsilon+\frac{1}{2}\mathbf{u}\cdot\mathbf{u} \]

Total enthalpy inafafanuliwa kama:

\[ h=e+\frac{p}{\rho} \]

p ni pressure. Study imetumia perfect-gas equation:

\[ p=\rho\varepsilon(\gamma-1) \]

γ ni polytropic au specific-heat ratio. Katika numerical tests, γ = 7/5 imetumika.

Kwa nini enthalpy hubaki constant katika steady Euler flow?

Two-dimensional Euler equations zinaeleza conservation ya mass, momentum na total energy katika form ya jumla:

\[ \frac{\partial\mathbf{U}}{\partial t} + \nabla\cdot\mathbf{F}(\mathbf{U}) =0 \]

Katika steady state, time-dependent term hutoweka. Mass na energy equations zikichunguzwa pamoja, hupatikana:

\[ \mathbf{u}\cdot\nabla h=0 \]

Expression hii inaonyesha kwamba katika steady Euler flow bila viscosity na heat conduction, total enthalpy haibadiliki along streamline.

Ikiwa streamlines zimeunganishwa na inlet boundary, total enthalpy inapaswa kubaki sawa na inlet value h∞ katika entire computational domain. Shock hubadilisha pressure, temperature na density; lakini katika ideal steady Euler solution total enthalpy huhifadhiwa across shock.

Kwa nini numerical schemes zinaweza kupoteza property hii?

Katika finite-volume methods, computational domain hugawanywa katika cells. Mabadiliko ya mass, momentum na energy katika kila cell huamuliwa na numerical fluxes zinazopita kwenye cell faces. Riemann problem huundwa kati ya states pande mbili za face, na numerical flux hupatikana kutoka approximate solution yake.

Classical Godunov-type schemes kama Rusanov, HLL, HLLC na Roe ni conservative, lakini haziongezi numerical diffusion katika mass na energy fluxes kwa namna inayohifadhi exactly same physical relation. Kwa hiyo, across stationary shock, artificial diffusion katika energy flux inaweza kutokuwa sawa na artificial diffusion katika mass flux times total enthalpy.

Matokeo yake total enthalpy iliyokokotolewa upande mmoja wa shock inakuwa tofauti na upande mwingine. Error hii kisha husafirishwa na flow na inaweza kuathiri energy distribution hasa katika low-speed regions nyuma ya shock.

Umuhimu wa enthalpy error katika engineering

Mtafiti anasema kwamba katika baadhi ya Navier–Stokes calculations, criteria zinazotumiwa kufafanua boundary layer zinaweza kutegemea total enthalpy. Artificial enthalpy trace inayotokana na shock ikisafirishwa downstream along wall, inaweza kuwa vigumu kutenganisha physical boundary layer na numerical energy error.

Katika hypersonic vehicles, density, pressure, temperature na heat flux karibu na stagnation point ni muhimu sana. Mach 20 test ya study inaonyesha kwamba total-enthalpy error haivurugi auxiliary variable pekee; inaweza kubadilisha shape ya density field pia.

Simple Riemann solver ni nini?

Simple Riemann solvers hufafanua constant intermediate states kati ya initial left na right states. States hizi hutenganishwa na discontinuities zinazosafiri kwa different speeds:

\[ W(\xi)= \begin{cases} U_l, & \xi\leq\lambda_1 \\ U_k, & \lambda_{k-1}<\xi\leq\lambda_k \\ U_r, & \xi>\lambda_m \end{cases} \]

Hapa ξ = x/t ni self-similar variable, Ul na Ur ni left na right states, na λk ni wave speeds.

Katika Gallice approach, solver kwanza hujengwa katika Lagrange coordinates. Kisha hubadilishwa kuwa Euler form inayofanya kazi kwenye fixed spatial cells kupitia Lagrange–Euler transformation. Method hii hurahisisha analysis ya positive density na internal-energy conditions katika intermediate states.

One-dimensional Gallice scheme

Katika scheme inayoitwa one-dimensional, numerical flux kwenye cell face hukokotolewa kutoka cells mbili tu pande za face. Hata mesh ikiwa two-dimensional, Riemann problem kwenye kila face hushughulikiwa kama one-dimensional katika direction normal to face.

Normal velocity ya intermediate state inafafanuliwa:

\[ u_f= \frac{ \lambda_lu_{n,l} + \lambda_ru_{n,r} - \Delta p }{ \lambda_l+\lambda_r } \]

λl na λr ni Lagrange wave speeds, un ni normal velocity component, na Δp ni pressure difference kati ya cells mbili.

Paper inaeleza kwamba one-dimensional Gallice scheme katika Euler form ni largely equivalent na HLLC solver inayotumia specific wave-speed estimates.

Multidimensional Gallice-2D scheme

Katika multidimensional scheme, flux haiathiriwi tu na face neighbors bali na all cells zinazoshiriki same node. Figure 1 kwenye page 4 inalinganisha classical face fluxes na multidimensional fluxes zilizokokotolewa around a node.

Normal velocity ya intermediate state hufafanuliwa kama projection ya common nodal velocity vector kwenye face normal:

\[ u_n^*=\mathbf{u}_p\cdot\mathbf{n} \]

Nodal velocity up hupatikana kutoka linear system inayotimiza conservation condition ya all faces around node. Multidimensional stencil hii inaweza kuelekeza numerical diffusion kulingana na flow na mesh geometry.

Kulingana na study, multidimensional structure ni more resistant kuliko one-dimensional solvers dhidi ya carbuncle instability katika strong bow shocks.

Positivity na stability conditions

Katika gas-dynamics solution, density na internal energy hazipaswi kuwa negative. Study inaonyesha kwamba kwa one-dimensional solver, intermediate states hubaki positive ikiwa Lagrange wave speeds zinazidi lower bounds fulani.

General condition ni:

\[ \lambda_l> \max \left( \rho_lc_l,\, \sqrt{\rho_l|\Delta p|},\, -\rho_l\Delta u_n \right) \]

\[ \lambda_r> \max \left( \rho_rc_r,\, \sqrt{\rho_r|\Delta p|},\, -\rho_r\Delta u_n \right) \]

c ni speed of sound. Katika calculations, largest slope kutoka left na right conditions ilitumika kwa kila face.

Courant–Friedrichs–Lewy number ya time step ilichaguliwa kuwa 0,5. Local time stepping ilitumika katika steady problems.

Häenel condition

Main condition ya total-enthalpy-preserving schemes ni kwamba numerical diffusion katika energy flux iwe proportional na ile katika mass flux:

\[ D_{\rho e}=hD_\rho \]

Hapa Dρe ni numerical diffusion katika energy flux na Dρ ni numerical diffusion katika mass flux. Katika steady flow, wakati h = h∞, physical link kati ya energy na mass flux huhifadhiwa.

Intermediate energy states za classical Gallice solver hazitimizi equality hii. Energy difference iliyokokotolewa kupitia intermediate states si sawa na density difference times total enthalpy.

MGallice modification

Mtafiti anahifadhi mass na momentum structure ya Gallice solver lakini anafafanua upya intermediate energy states kupitia total enthalpy. Modified Lagrange state vector ni:

\[ \widetilde{\mathbf{V}}_s= \begin{pmatrix} v_s \\ \mathbf{u}_s \\ h_s \end{pmatrix} \]

v ni specific volume, u velocity na h total enthalpy.

Energy component ya modified intermediate state katika Euler frame huundwa kama:

\[ \widetilde{\mathbf{U}}_s^* = \rho_s^* \begin{pmatrix} 1 \\ \mathbf{u}_s+(u_n^*-u_{n,s})\mathbf{n} \\ h_s \end{pmatrix} \]

Kwa njia hii, change katika energy component across each wave inaunganishwa na density change kupitia same total enthalpy.

Katika steady solution, ikiwa left na right total enthalpies ni sawa na h∞:

\[ \rho_l^*e_l^*-\rho_lh_l = (\rho_l^*-\rho_l)h_\infty \]

na equalities zinazofanana hutimia kwa other intermediate waves. Structure hii inatimiza Häenel condition kwa both one- na multidimensional modified Gallice solvers.

Gharama muhimu ya modification

Corresponding energy flux katika Lagrange form inakuwa diffusion-free central flux baada ya modification:

\[ \overline{pu} = \frac{1}{2} \left( p_lu_{n,l} + p_ru_{n,r} \right) \]

Sufficient numerical diffusion katika mass na momentum equations ilitosha kukamata correct solution katika tests. Hata hivyo author hawezi kutoa condition inayohakikisha modified MGallice solvers ni entropy-increasing au entropy-stable.

Kwa hiyo method inapohifadhi total enthalpy, haiinherit automatically full proven entropy property ya classical scheme. Hii ni mojawapo ya most important theoretical limitations za study.

Second-order MUSCL extension

First-order finite-volume schemes ni stable katika strong shocks lakini huzalisha high numerical diffusion katika smooth regions. Study imetumia MUSCL-type reconstruction kwa second-order solution.

  • Density, velocity na pressure gradients zilikokotolewa kwenye nodes kwa weighted linear least-squares method.
  • Cell-center gradients ziliunganishwa kwa weights zinazotegemea inverse norms za nodal gradients.
  • Katika MGallice schemes total enthalpy pia ilireconstructed.
  • Gradients zililimited kwa R3 Nishikawa limiter.
  • Katika limiter p = 3 na cell-scale-dependent εc zilitumika.

Schemes nne zilizolinganishwa

SchemeFlux stencilTotal-enthalpy modification
GalliceOne-dimensional based on face neighborsHakuna
Gallice-2DMultidimensional based on node neighborsHakuna
MGalliceOne-dimensional based on face neighborsIpo
MGallice-2DMultidimensional based on node neighborsIpo

Jaribio la kwanza: Channel yenye triangular bumps kwa Mach 2

Problem ya kwanza ni supersonic flow katika two-dimensional channel yenye small triangular bumps kwenye opposite walls. Domain ina urefu wa 2 meter na height ya 1 meter. Bumps mbili ziliwekwa kwenye lower na upper walls.

Free-stream Mach number2
Computational domain[0, 2] m × [−0,5, 0,5] m
Mesh200 × 100 quadrilateral cells
Initial na inlet stateρ = 1, ux = M∞√γ, uy = 0, p = 1
Specific-heat ratioγ = 7/5
Inlet total enthalpyh∞ = 6,3
OrderSecond-order MUSCL
Convergence criterionDensity residual falling to 10−8

Flow inapogonga bumps, stationary shocks mbili huunda mbele yake. Expansion waves zinazokua nyuma huingiliana na kureflect kutoka upper na lower walls, na kutengeneza diamond-shaped wave pattern.

Figure 3 kwenye page 9 inaonyesha kwamba schemes zote ziliconverge hadi steady solution baada ya takribani 4.000 iterations. Kwa kuwa residual curves ziko close, enthalpy modification haikuathiri noticeably convergence capability.

Density fields zilikuwa karibu sawa

Figure 4 kwenye page 10 inaonyesha density contours za solvers zote nne kati ya 0,63 na 1,77. Solvers zote zinarudisha shocks, expansion waves na diamond interaction region kwa namna inayofanana.

One-dimensional schemes zimeripotiwa kuwa slightly less diffusive. Hata hivyo ukiangalia density field pekee, main difference kati ya classical na modified solvers haionekani wazi.

Enthalpy fields zilikuwa tofauti kwa uwazi

Figure 5 kwenye page 10 inaonyesha total-enthalpy field kati ya 6,2 na 6,44. Classical Gallice scheme huunda small enthalpy jumps across stationary shocks na errors hizi husafirishwa downstream along wall.

Enthalpy distortion katika Gallice-2D huwa more pronounced. Kwa upande mwingine MGallice na MGallice-2D fields zinaonekana uniform; inlet total enthalpy ya 6,3 inaripotiwa kuhifadhiwa katika entire domain.

Test hii inaonyesha kwamba hata density fields zikionekana similar, energy-preservation properties za solvers zinaweza kuwa significantly different.

Jaribio la pili: Half-cylinder kwa Mach 20

Problem ya pili ni steady hypersonic flow around half-cylinder yenye radius 1 meter. Problem hii ni challenging kwa sababu ya both strong detached bow shock na carbuncle instability.

Free-stream Mach number20
Cylinder radius1 m
Computational domainElliptical domain
Mesh32 × 128 quadrilateral cells
Initial na inlet stateρ = 1, ux = M∞√γ, uy = 0, p = 1
Specific-heat ratioγ = 7/5
Inlet total enthalpyh∞ = 283,5
Initial comparisonFirst-order schemes

Reporting inconsistency katika iteration count

Text inasema first-order solutions ziliendeshwa hadi 105 iterations na residual ikapungua kwa three-four orders kufikia point hiyo. Hata hivyo horizontal axis ya Figure 6 kwenye page 10 inaishia karibu 11.000 iterations.

Kwa hiyo kuna internal preprint reporting inconsistency kati ya reported 105 value na plotted graph. Graph inaonyesha kwamba multidimensional schemes huzalisha lower na more regularly decreasing residuals katika displayed interval.

Carbuncle instability

Carbuncle ni nonphysical distortion ya strong shocks aligned with grid lines katika baadhi ya shock-capturing schemes. Stagnation line ya bow shock mbele ya half-cylinder ni classic situation ambako instability hii inaweza kutokea.

Study inasema one-dimensional Gallice na MGallice schemes hazikuweza correctly resolve strong bow shock. Enthalpy-preserving modification pekee haikuondoa carbuncle.

Gallice-2D na MGallice-2D schemes, kwa multidimensional flux structures zao, zilitengeneza carbuncle-free bow shocks. Result hii inaonyesha kwamba shock stability na total-enthalpy preservation ni two separate numerical properties.

Kuhifadhi enthalpy kulirekebisha density field

Figure 7 kwenye page 10 inaonyesha density contours za four schemes kutoka 1 hadi 6,3. Katika classical Gallice-2D, maximum density haitokei katika expected stagnation point. Hii inaonyesha kwamba energy distribution nyuma ya bow shock imebadilika nonphysically.

MGallice-2D hutumia same multidimensional flux stencil, lakini enthalpy-preserving energy modification yake inarudisha density maximum karibu na stagnation point.

Figure 8 kwenye page 10 inaonyesha classical solvers zinatengeneza total-enthalpy decrease across shock. Artificial loss hii hubadilisha energy distribution katika subsonic region nyuma ya shock na hatimaye kuvuruga density field.

Katika MGallice na MGallice-2D solutions total enthalpy inaripotiwa kubaki constant katika entire domain hadi machine precision.

Athari ya second-order solution

Half-cylinder problem ilikokotolewa upya kwa second-order MUSCL extensions za Gallice-2D na MGallice-2D. Figure 9 kwenye page 11 inaonyesha density na total-enthalpy fields side by side.

Second-order reconstruction ilipunguza total-energy diffusion katika classical Gallice-2D. Kwa hiyo Gallice-2D na MGallice-2D density fields zilikuwa closer kuliko first-order results.

Hata hivyo classical Gallice-2D iliendelea kutohifadhi total enthalpy exactly, wakati MGallice-2D iliendelea kuiweka field constant at inlet value. Higher order ilipunguza error magnitude lakini haikuchukua nafasi ya conservation property.

Hitimisho kuu la engineering

Proposed modification haikutoa tu smoother total-enthalpy field kwenye plots. Katika Mach 20 test, ilirekebisha energy distribution nyuma ya shock na hivyo kubadilisha density field ya stagnation region.

Result hii inaonyesha three separate numerical requirements:

  • Scheme inapaswa kukamata strong shocks stably.
  • Inapaswa kuwa na resistance ya kutosha dhidi ya multidimensional shock instabilities.
  • Inapaswa kuhifadhi physical invariants za steady solution, kama total enthalpy, katika discrete level.

MGallice-2D ndiyo ilikuwa structure pekee katika examined Mach 20 test iliyotoa both multidimensional shock stability na total-enthalpy preservation.

Maeneo yanayoweza kutumika kwa Türkiye

Method inaweza kuhamishwa kwenye tathmini ya numerical-flux selection katika domestic au academic CFD software inayotengenezwa kwa high-speed air vehicles na atmospheric-entry problems. Inahusiana hasa na research areas zifuatazo:

  • Hypersonic na supersonic vehicle aerodynamics,
  • External-flow analyses za rockets na high-speed flight vehicles,
  • Bow-shock na stagnation-point calculations,
  • Finite-volume solvers on unstructured meshes,
  • Shock-instability-resistant multidimensional Riemann solvers,
  • Boundary-layer na wall-heat-flux calculations katika Navier–Stokes solutions.

Hata hivyo results za research haziwezi kutumika directly kama vehicle surface temperature, heat-shield thickness au aerodynamic forces. Kwa quantities hizi viscosity, heat conduction, real-gas effects, chemical reactions, turbulence na solid-material behavior lazima zimodeliwe separately.

Matokeo ambayo study inaonyesha

  • Gallice-type Simple Riemann solvers zinaweza kufanywa total-enthalpy-preserving kwa modifying intermediate energy states.
  • Modification inaweza kutumika kwa both one- na multidimensional solvers.
  • Katika Mach 2 triangular-bump channel test, modified schemes zilihifadhi total enthalpy katika entire domain.
  • Katika Mach 20 half-cylinder test, multidimensional schemes zilizuia carbuncle instability.
  • MGallice-2D iliondoa total-enthalpy loss na kurekebisha density distribution katika stagnation region.
  • Second-order reconstruction ilipunguza enthalpy error ya classical scheme lakini haikuiondoa completely.

Matokeo ambayo study haionyeshi

  • Haijaonyeshwa kwamba proposed solvers ndiyo best method kwa all hypersonic geometries.
  • Entropy stability haijathibitishwa mathematically.
  • Validation kwenye viscous Navier–Stokes equations haijafanywa.
  • Wall heat flux au boundary-layer accuracy haijapimwa directly.
  • Real-gas, chemical-reaction au thermochemical non-equilibrium effects hazijachunguzwa.
  • Three-dimensional geometry au three-dimensional Riemann solver haijatestwa.
  • Hakuna comparison na experimental wind-tunnel data.
  • Hakuna detailed performance comparison ya computation time na processor cost.

Mbinu na Matokeo ya Utafiti

Muundo wa utafiti

Study ni computational-engineering research inayotegemea mathematical modification ya one- na multidimensional Godunov-type Simple Riemann solvers kwa steady compressible Euler equations na comparison yao katika two numerical test problems.

Mathematical method

Governing equationsTwo-dimensional compressible Euler equations
Gas modelPerfect gas, γ = 7/5
Spatial discretizationCell-centered finite-volume method
Flux approachGodunov-type Simple Riemann solvers
Coordinate approachSolver construction in Lagrange frame na Lagrange–Euler transformation
One-dimensional solverGallice/HLLC-like scheme based on two cells at a face
Multidimensional solverGallice-2D based on common nodal velocity from cells around a node
Enthalpy modificationRedefinition of intermediate energy states using total enthalpy
Time steppingLocal time stepping, CFL = 0,5
Second-order extensionMUSCL reconstruction na R3 Nishikawa limiter

Comparison ya test results

TestClassical GalliceClassical Gallice-2DMGalliceMGallice-2D
Mach 2 channel convergenceStableStableStableStable
Mach 2 density fieldSimilar overall solutionSimilar overall solutionSimilar overall solutionSimilar overall solution
Mach 2 total enthalpySmall jumps at shocksMore pronounced distortionPreserved across domainPreserved across domain
Mach 20 carbuncle behaviorCarbuncle occursNo carbuncleCarbuncle occursNo carbuncle
Mach 20 total enthalpyDecrease at shockDecrease at shockPreserved to machine precisionPreserved to machine precision
Mach 20 density fieldDistorted due to shock instabilityIncorrect distribution at stagnation pointDistorted due to shock instabilityExpected stagnation-point maximum recovered

Numerical values

  • Mach 2 channel mesh: 200 × 100 quadrilateral cells.
  • Mach 2 inlet total enthalpy: 6,3.
  • Katika Mach 2 solutions, density residual ilishuka hadi 10−8 katika takribani 4.000 iterations.
  • Mach 2 density visualization range: 0,63-1,77.
  • Mach 2 enthalpy visualization range: 6,2-6,44.
  • Mach 20 half-cylinder mesh: 32 × 128 quadrilateral cells.
  • Mach 20 inlet total enthalpy: 283,5.
  • Mach 20 density visualization range: 1-6,3.
  • Mach 20 enthalpy visualization range: 240-330.
  • Mach 20 text inaripoti residual reduction ya three-four orders.

Engineering interpretation ya figures

  • Page 4, Figure 1: Inalinganisha one-dimensional face fluxes na node-based multidimensional flux stencil.
  • Page 6, Figure 2: Inaonyesha four-state Simple Riemann solver katika Lagrange frame na wave speeds.
  • Page 9, Figure 3: Inaonyesha residual curves za schemes nne zikiconverge kwa takribani 4.000 iterations katika Mach 2 channel problem.
  • Page 10, Figure 4: Inaonyesha density-wave patterns ni generally similar kwa schemes nne katika Mach 2 problem.
  • Page 10, Figure 5: Inaonyesha enthalpy error ikisafirishwa downstream kutoka shocks katika classical schemes na kutoweka katika modified schemes.
  • Page 10, Figure 6: Inaonyesha residual history ya Mach 20 test na more regular convergence tendency ya multidimensional schemes.
  • Page 10, Figure 7: Inaonyesha bow shock na density field; carbuncle katika one-dimensional solvers na stagnation-point difference kati ya Gallice-2D na MGallice-2D.
  • Page 10, Figure 8: Inaonyesha shock-induced enthalpy loss katika classical schemes na uniform enthalpy field katika modified schemes.
  • Page 11, Figure 9: Inaonyesha second-order reconstruction ikipunguza error katika classical scheme huku MGallice-2D ikiendelea kuhifadhi total enthalpy exactly.

Nguvu

  • Proposed modification inategemea directly physical steady-flow invariant.
  • One- na multidimensional solvers zimechunguzwa ndani ya same mathematical framework.
  • Scheme haijawasilishwa analytically tu; imetestwa katika two strong-shock problems.
  • Density na enthalpy fields za classical na modified solvers zimelinganishwa separately.
  • Mach 20 test inaonyesha shock stability na enthalpy preservation ni separate requirements.
  • First- na second-order solutions zimelinganishwa.
  • Improvement ya MGallice-2D katika density field inaonyesha conservation property inaweza kubadilisha flow solution.

Mapungufu muhimu

  • Study ni preprint ambayo haijapitia peer review.
  • Only inviscid Euler equations zimetatuliwa.
  • Viscosity, heat conduction na direct wall heat flux hazipo.
  • Real-gas na high-temperature chemistry hazijamodeliwa.
  • Hakuna mathematical guarantee ya entropy stability.
  • Kwa multidimensional solver, explicit slope condition inayohakikisha positivity haijatolewa; slopes za one-dimensional solver zimetumika.
  • Only two numerical test problems zimetolewa.
  • Hakuna three-dimensional solver au three-dimensional vehicle geometry iliyochunguzwa.
  • Hakuna quantitative error analysis dhidi ya physical experiment au independent numerical reference.
  • Mach 2 comparison inategemea kwa kiasi kikubwa field visualizations.
  • Katika Mach 20 problem kuna inconsistency kati ya text statement ya 105 iterations na residual graph axis ya takribani 11.000 iterations.
  • CPU time, memory use na additional computational cost dhidi ya classical solvers hazijaripotiwa.
  • MGallice-2D performance kwa wall heat flux au boundary-layer identification bado haijatestwa.

Dokezo la Chanzo na Mbinu

Utafiti wa asili: “Enthalpy preserving Simple Riemann solvers for steady hypersonic flows.”

Mwandishi: Lucas Tallois.

Corresponding author: Lucas Tallois. Study ina mwandishi mmoja na hakuna equal-contribution statement.

Taasisi: CEA-CESTA, 15 avenue des Sablières, CS 60001, 33116 Le Barp Cedex, France.

Aina ya chanzo: Original computational-fluid-dynamics research na preprint.

Tarehe ya preprint: 12 Haziran 2026.

Hali ya peer review: Study haijapitia peer review. Results hazipaswi kutathminiwa kama findings za final peer-reviewed journal article.

Hali ya jarida: Text ina kauli “Preprint submitted to Elsevier”. Hakuna specific journal name, acceptance date au published final version information.

Platform: SSRN.

DOI:10.2139/ssrn.6948779

Kiungo rasmi:SSRN study record

Ufadhili: Hakuna funding statement katika study text.

Mgongano wa maslahi: Hakuna conflict-of-interest statement katika study text.

Ufikiaji wa data na code: Hakuna open data repository, numerical mesh files au source-code link iliyotolewa.

Scientific method, equations, numerical schemes, test conditions, figure interpretations, findings na limitations katika content hii zinatokana na study iliyochunguzwa. External sources zilitumika tu kwa bibliographic verification ya title, author, institution, platform na DOI; hakuna scientific result ambayo haikuwepo katika study iliyoongezwa.


Shiriki:

Maoni huchapishwa baada ya kukaguliwa.Maoni yako yatapitia mchakato wa idhini na yataonekana yakikubaliwa.

Acha maoni

Anwani yako ya barua pepe haitachapishwa. Sehemu za lazima zimewekewa alama ya *

Your experience on this site will be improved by allowing cookies Cookie Policy