Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hydrodynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1980

Hydrodynamic role of fish squamosal integument as an analog of the surfaces directly formed by the turbulent flow. Report 2: Hydrodynamic function of squamosal integument

The stream flowing round the slowly swimming squama free fish can be laminized with the aid of the external slime coat alone. The slime of the fish with well developed squamae can laminize the stream together with the squamatic integument. Adjustments preventing a loss of the slime during laminization are better developed in the fastest squama free fishes.

Kudryashov, A. F.↗
NASA NTRS1993

Comparison of hydrodynamic and semi-kinetic treatments for plasma flow along closed field lines

Hydrodynamic and semi-kinetic treatments of plasma flow along closed geomagnetic field lines are compared. The hydrodynamic treatment is based on a simplified 16-moment set of transport equations as the equations for the heat flows are not solved; the heat flows are treated heuristically. The semi-kinetic treatment is based on a particle code. The comparison deals with the distributions of the plasma density, flow velocity, and parallel and perpendicular temperatures as obtained from the two treatments during the various stages of the flow. In the kinetic treatment, the appropriate boundary condition is the prescription of the velocity distribution functions for the particles entering the flux tubes at the ionospheric boundaries; those particles leaving the system are determined by the processes occurring in the flux tube. The prescribed distributions are half-Maxwellian with temperature T(sub 0) and density n(sub 0). In the hydrodynamic model, the prescribed boundary conditions are on density (n(sub 0)), flow velocity (V(sub 0)) and temperature (T(sub 0). It was found that results from the hydrodynamic treatment critically depend on V(sub 0); for early stages of the flow this treatment yields results in good agreement with those from the kinetic treatment, when V(sub 0) = square root of (kT(sub 0)/2 (pi)m), which is the average velocity of particles moving in a given direction for a Maxwellian distribution. During this early stage, the flows developing form the conjugate ionospheres show some distinct transitions. For the first hour or so, the flows are highly supersonic and penetrate deep into the opposite hemispheres, and both hydrodynamics and kinetic treatments yield almost similar features. It is found that during this period heatflow effects are negligibly small. When a flow penetrates deep into the opposite hemisphere, the kinetic treatment predicts reflection and setting up of counterstreaming. In contrast, the hydrodynamic treatment yields a shock in the flow. The reasons for this difference in the two treatments is discussed, showing that in view of the relatively warm ions, the coupling of ion beams and the consequent shock formation in the offequatorial region are not likely due to the enhancements in the beam temperatures. The counterstreaming in the kinetic treatment and the shock in the hydrodynamic treatment first advance upward to the equator and then downward to the ionospheric boundary from where the flow originated. The transit time for this advancement is found to be about 1 hour for the respective models. After 2 hours or so, both models predict that the flows from the ionospheric boundaries are generally subsonic with respect to the local ion-sound speed. At late stages of the flow, when a substantial fraction of ions entering the flux tube begin to return back in the kinetic treatment, the hydrodynamic treatment with the boundary condition V(sub 0) = square root of (kT(sub 0)/2(pi)m) yields an over-refilling, and the choice of V(sub 0) becomes uncertain.

Singh, Nagendra↗
NASA NTRS1999

Hydrodynamic Instability and Thermal Coupling in a Dynamic Model of Liquid-Propellant Combustion

For liquid-propellant combustion, the Landau/Levich hydrodynamic models have been combined and extended to account for a dynamic dependence of the burning rate on the local pressure and temperature fields. Analysis of these extended models is greatly facilitated by exploiting the realistic smallness of the gas-to-liquid density ratio rho. Neglecting thermal coupling effects, an asymptotic expression was then derived for the cellular stability boundary A(sub p)(k) where A(sub p) is the pressure sensitivity of the burning rate and k is the disturbance wavenumber. The results explicitly indicate the stabilizing effects of gravity on long-wave disturbances, and those of viscosity and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limit of weak gravity, hydrodynamic instability in liquid-propellant combustion becomes a long-wave, instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumbers. In addition, surface tension and viscosity (both liquid and gas) each produce comparable effects in the large-wavenumber regime, thereby providing important modifications to the previous analyses in which one or more of these effects was neglected. For A(sub p)= O, the Landau/Levich results are recovered in appropriate limiting cases, although this typically corresponds to a hydrodynamically unstable parameter regime for p << 1. In addition to the classical cellular form of hydrodynamic stability, there exists a pulsating form corresponding to the loss of stability of steady, planar burning to time-dependent perturbations. This occurs for negative values of the parameter A(sub p), and is thus absent from the original Landau/Levich models. In the extended model, however, there exists a stable band of negative pressure sensitivities bounded above by the Landau type of instability, and below by this pulsating form of hydrodynamic instability. Indeed, nonsteady modes of combustion have been observed at low pressures in hydroxylammonium nitrate (HAN)-based liquid propellants, which often exhibit negative pressure sensitivities. While nonsteady combustion may correspond to secondary and higher-order bifurcations above the cellular boundary, it may also be a manifestation of this pulsating type of hydrodynamic instability. In the present work, a nonzero temperature sensitivity is incorporated into our previous asymptotic analyses. This entails a coupling of the energy equation to the previous purely hydrodynamic problem, and leads to a significant modification of the pulsating boundary such that, for sufficiently large values of the temperature-sensitivity parameter, liquid-propellant combustion can become intrinsically unstable to this alternative form of hydrodynamic instability. For simplicity, further attention is confined here to the inviscid version of the problem since, despite the fact that viscous and surface-tension effects are comparable, the qualitative nature of the cellular boundary remains preserved in the zero-viscosity limit, as does the existence of the pulsating boundary. The mathematical model adopts the classical assumption that there is no distributed reaction in either the liquid or gas phases, but now the reaction sheet, representing either a pyrolysis reaction or an exothermic decomposition at the liquid/gas interface, is assumed to depend on local conditions there.

Margolis, S. B.↗
DOE OSTI2025

A three-dimensional laser ray-tracing methodology for radiation-hydrodynamics simulations

We report on a methodology for performing laser ray-tracing in three spatial dimensions for radiation-hydrodynamics simulation codes. Our method, which is an extension of that developed in Haines et al., Comput. Fluids 201, 104478 (2020), utilizes an automatically generated separate mesh for the laser ray-tracing from the radiation-hydrodynamics mesh. This enables the laser mesh to be tailored to minimize ray noise with significantly fewer rays than would be required when the ray-tracing is performed on the radiation-hydrodynamics mesh, primarily by allowing the use of high-aspect-ratio cells that are not suitable for hydrodynamics solvers. For a planar target, we show that our method provides a ≈ 100× reduction in computational expense to achieve a fixed level of ray noise relative to ray-tracing directly on the radiation-hydrodynamics mesh. The relatively low ray requirement also enables efficient computation of cross-beam energy transfer. Each cell in the logically cubic laser mesh is a non-convex dodecahedron with triangular sides, and numerical integration of the ray trajectories and inverse bremsstrahlung is performed by mapping each cell to the unit cube. We will describe our methodology in detail as well as its implementation in the xRAGE radiation-hydrodynamics code, discuss performance, and present the results from applying the methodology to test problems with analytic solutions for laser ray-tracing through a quadratic density gradient with an analytic solution as well as for a laser-driven heat front. In 3D radiation-hydrodynamics simulations of laser-driven experiments performed on the National Ignition Facility, laser ray-tracing with our methodology uses less than 1% of total computational time while introducing acceptably low levels of ray noise.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

Design of Fiber Optic Sensors for Measuring Hydrodynamic Parameters

The science of optical hydrodynamics involves relating the optical properties to the fluid dynamic properties of a hydrodynamic system. Fiber-optic sensors are being designed for measuring the hydrodynamic parameters of various systems. As a flowing fluid makes an encounter with a flat surface, it forms a boundary layer near this surface. The region between the boundary layer and the flat plate contains information about parameters such as viscosity, compressibility, pressure, density, and velocity. An analytical model has been developed for examining the hydrodynamic parameters near the surface of a fiber-optic sensor. An analysis of the conservation of momentum, the continuity equation and the Navier-Stokes equation for compressible flow were used to develop expressions for the velocity and the density as a function of the distance along the flow and above the surface. When examining the flow near the surface, these expressions are used to estimate the sensitivity required to perform direct optical measurements and to derive the shear force for indirect optical measurements. The derivation of this result permits the incorporation of better design parameters for other fiber-based sensors. Future work includes analyzing the optical parametric designs of fiber-optic sensors, modeling sensors to utilize the parameters for hydrodynamics and applying different mixtures of hydrodynamic flow. Finally, the fabrication of fiber-optic sensors for hydrodynamic flow applications of the type described in this presentation could enhance aerospace, submarine, and medical technology.

Lyons, Donald R.↗
NASA NTRS1998

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion

The classical Landau/Levich models of liquid propellant combustion, which serve as seminal examples of hydrodynamic instability in reactive systems, have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and/or temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity, surface tension and viscosity on the hydrodynamic instability of the propagating liquid/gas interface. In particular, a composite asymptotic expression, spanning three distinguished wavenumber regimes, is derived for both cellular and pulsating hydrodynamic neutral stability boundaries A(sub p)(k), where A(sub p) is the pressure sensitivity of the burning rate and k is the disturbance wavenumber. For the case of cellular (Landau) instability, the results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limiting case of weak gravity, it is shown that cellular hydrodynamic instability in this context is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(l) wavenumber disturbances. It is also demonstrated that, in the large wavenumber regime, surface tension and both liquid and gas viscosity all produce comparable stabilizing effects in the large-wavenumber regime, thereby providing significant modifications to previous analyses of Landau instability in which one or more of these effects were neglected. In contrast, the pulsating hydrodynamic stability boundary is found to be insensitive to gravitational and surface-tension effects, but is more sensitive to the effects of liquid viscosity, which is a significant stabilizing effect for O(l) and higher wavenumbers. Liquid-propellant combustion is predicted to be stable (i.e., steady and planar) only for a range of negative pressure sensitivities that lie between the two types of hydrodynamic stability boundaries.

Margolis, Stephen B.↗
NASA NTRS2001

Hydrodynamic Stability of Liquid-Propellant Combustion: Landau's Problem Revisited

Hydrodynamic, or Landau, instability in combustion is typically associated with the onset of wrinkling of a flame surface, corresponding to the formation of steady cellular structures as the stability threshold is crossed. As its name suggests, it stems from hydrodynamic effects connected with thermal expansion across the reaction region. In the context of liquid-propellant combustion, the classical models that originally predicted this phenomenon have been extended to include the important effects that arise from a dynamic dependence of the burning rate on the local pressure and temperature fields. Thus, the onset of Landau instability has now been shown to occur for sufficiently small negative values of the pressure sensitivity of the burning rate, significantly generalizing previous classical results for this problem that assumed a constant normal burning rate. It has also been shown that the onset of instability occurs for decreasing values of the disturbance wave number as the gravitational-acceleration parameter decreases. Consequently, in an appropriate weak-gravity limit, Landau instability becomes a long-wave phenomena associated with the formation of large cells on the liquid-propellant surface. Additionally, a pulsating form of hydrodynamic instability has been shown to occur as well, corresponding to the onset of temporal oscillations in the location of the liquid/gas interface. This instability occurs for sufficiently large negative values of the pressure sensitivity, and is enhanced by increasing values of the burning-rate temperature sensitivity. It is further shown that for sufficiently small values of this parameter, there exists a stable range of pressure sensitivities for steady, planar burning such that the classical cellular form of hydrodynamic instability and the more recent pulsating form of hydrodynamic instability can each occur as the corresponding stability threshold is crossed. For larger thermal sensitivities, however, the pulsating stability boundary evolves into a C-shaped curve in the (disturbance-wave number, pressure-sensitivity) plane, indicating loss of stability to pulsating perturbations for all sufficiently large disturbance wavelengths. It is thus concluded, based on characteristic parameter values that an equally likely form of hydrodynamic instability in liquid-propellant combustion is of a nonsteady, long-wave nature, distinct from the steady, cellular form originally predicted by Landau.

Margolis, S. B.↗
NASA NTRS2017

Early Hydrodynamic Escape Limits Rocky Planets to Less Than or Equal to 1.6 Earth Radii

In the past decade thousands of exoplanet candidates and hundreds of confirmed exoplanets have been found. For sub-Neptune-sized planets, those less than approx. 10 Earth masses, we can separate planets into two broad categories: predominantly rocky planets, and gaseous planets with thick volatile sheaths. Observations and subsequent analysis of these planets show that rocky planets are only found with radii less than approx. 1.6 Earth radii. No rocky planet has yet been found that violates this limit. We propose that hydrodynamic escape of hydrogen rich protoatmospheres, accreted by forming planets, explains the limit in rocky planet size. Following the hydrodynamic escape model employed by Luger et al. (2015), we modelled the XUV driven escape from young planets (less than approx.100 Myr in age) around a Sun-like star. With a simple, first-order model we found that the rocky planet radii limit occurs consistently at approx. 1.6 Earth radii across a wide range of plausible parameter spaces. Our model shows that hydrodynamic escape can explain the observed cutoff between rocky and gaseous planets. Fig. 1 shows the results of our model for rocky planets between 0.5 and 10 Earth masses that accrete 3 wt. % H2/He during formation. The simulation was run for 100 Myr, after that time the XUV flux drops off exponentially and hydrodynamic escape drops with it. A cutoff between rocky planets and gaseous ones is clearly seen at approx. 1.5-1.6 Earth radii. We are only interested in the upper size limit for rocky planets. As such, we assumed pure hydrogen atmospheres and the highest possible isothermal atmospheric temperatures, which will produce an upper limit on the hydrodynamic loss rate. Previous work shows that a reasonable approximation for an upper temperature limit in a hydrogen rich protoatmosphere is 2000-3000 K, consistent with our assumptions. From these results, we propose that the observed dichotomy between mini-Neptunes and rocky worlds is simply explained by an early episode of thermally-driven hydrodynamic escape when host stars have saturated XUV fluxes.

Lehmer, O. R.↗
DOE OSTI2025

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗
DOE OSTI2024

High-precision Galaxy Clustering Predictions from Small-volume Hydrodynamical Simulations via Control Variates

Abstract Cosmological simulations of galaxy formation are an invaluable tool for understanding galaxy formation and its impact on cosmological parameter inference from large-scale structures. However, their high computational cost is a significant obstacle for running simulations that probe cosmological volumes comparable to those analyzed by contemporary large-scale structure experiments. In this work, we explore the possibility of obtaining high-precision galaxy clustering predictions from small-volume hydrodynamical simulations such as MillenniumTNG and FLAMINGO via control variates. In this approach, the hydrodynamical full-physics simulation is paired with a matched low-resolution gravity-only simulation. By learning the galaxy–halo connection from the hydrodynamical simulation and applying it to the gravity-only counterpart, one obtains a galaxy population that closely mimics the one in the more expensive simulation. One can then construct an estimator of galaxy clustering that combines the clustering amplitudes in the small-volume hydrodynamical and gravity-only simulations with clustering amplitudes in a large-volume gravity-only simulation. Depending on the galaxy sample, clustering statistic, and scale, this galaxy clustering estimator can have an effective volume of up to around 100 times the volume of the original hydrodynamical simulation in the nonlinear regime. With this approach, we can construct galaxy clustering predictions from existing simulations that are precise enough for mock analyses of next-generation large-scale structure surveys such as the Dark Energy Spectroscopic Instrument and the Legacy Survey of Space and Time.

Doytcheva, Alexandra (ORCID:0009000111254888)↗
NASA NTRS1981

Effect of Surface Roughness on Hydrodynamic Bearings

A theoretical analysis on the performance of hydrodynamic oil bearings is made considering surface roughness effect. The hydrodynamic as well as asperity contact load is found. The contact pressure was calculated with the assumption that the surface height distribution was Gaussian. The average Reynolds equation of partially lubricated surface was used to calculate hydrodynamic load. An analytical expression for average gap was found and was introduced to modify the average Reynolds equation. The resulting boundary value problem was then solved numerically by finite difference methods using the method of successive over relaxation. The pressure distribution and hydrodynamic load capacity of plane slider and journal bearings were calculated for various design data. The effects of attitude and roughness of surface on the bearing performance were shown. The results are compared with similar available solution of rough surface bearings. It is shown that: (1) the contribution of contact load is not significant; and (2) the hydrodynamic and contact load increase with surface roughness.

Majumdar, B. C.↗
NASA NTRS1995

Non-parabolic hydrodynamic formulations for the simulation of inhomogeneous semiconductor devices

Hydrodynamic models are becoming prevalent design tools for small scale devices and other devices in which high energy effects can dominate transport. Most current hydrodynamic models use a parabolic band approximation to obtain fairly simple conservation equations. Interest in accounting for band structure effects in hydrodynamic device simulation has begun to grow since parabolic models can not fully describe the transport in state of the art devices due to the distribution populating non-parabolic states within the band. This paper presents two different non-parabolic formulations of the hydrodynamic model suitable for the simulation of inhomogeneous semiconductor devices. The first formulation uses the Kane dispersion relationship (hk)(exp 2)/2m = W(1 + alpha(W)). The second formulation makes use of a power law ((hk)(exp 2)/2m = xW(sup y)) for the dispersion relation. Hydrodynamic models which use the first formulation rely on the binomial expansion to obtain moment equations with closed form coefficients. This limits the energy range over which the model is valid. The power law formulation readily produces closed form coefficients similar to those obtained using the parabolic band approximation. However, the fitting parameters (x,y) are only valid over a limited energy range. The physical significance of the band non-parabolicity is discussed as well as the advantages/disadvantages and approximations of the two non-parabolic models. A companion paper describes device simulations based on the three dispersion relationships: parabolic, Kane dispersion, and power low dispersion.

Smith, Arlynn W.↗
NASA NTRS1996

Non-Parabolic Hydrodynamic Formulations for the Simulation of Inhomogeneous Semiconductor Devices

Hydrodynamic models are becoming prevalent design tools for small scale devices and other devices in which high energy effects can dominate transport. Most current hydrodynamic models use a parabolic band approximation to obtain fairly simple conservation equations. Interest in accounting for band structure effects in hydrodynamic device simulation has begun to grow since parabolic models cannot fully describe the transport in state of the art devices due to the distribution populating non-parabolic states within the band. This paper presents two different non-parabolic formulations or the hydrodynamic model suitable for the simulation of inhomogeneous semiconductor devices. The first formulation uses the Kane dispersion relationship ((hk)(exp 2)/2m = W(1 + alphaW). The second formulation makes use of a power law ((hk)(exp 2)/2m = xW(exp y)) for the dispersion relation. Hydrodynamic models which use the first formulation rely on the binomial expansion to obtain moment equations with closed form coefficients. This limits the energy range over which the model is valid. The power law formulation readily produces closed form coefficients similar to those obtained using the parabolic band approximation. However, the fitting parameters (x,y) are only valid over a limited energy range. The physical significance of the band non-parabolicity is discussed as well as the advantages/disadvantages and approximations of the two non-parabolic models. A companion paper describes device simulations based on the three dispersion relationships; parabolic, Kane dispersion and power law dispersion.

Smith, A. W.↗
NASA NTRS1999

Pulsating Hydrodynamic Instability and Thermal Coupling in an Extended Landau/Levich Model of Liquid-Propellant Combustion

Hydrodynamic (Landau) instability in combustion is typically associated with the onset of wrinkling of a flame surface, corresponding to the formation of steady cellular structures as the stability threshold is crossed. In the context of liquid-propellant combustion, such instability has recently been shown to occur for critical values of the pressure sensitivity of the burning rate and the disturbance wavenumber, significantly generalizing previous classical results for this problem that assumed a constant normal burning rate. Additionally, however, a pulsating form of hydrodynamic instability has been shown to occur as well, corresponding to the onset of temporal oscillations in the location of the liquid/gas interface. In the present work, we consider the realistic influence of a nonzero temperature sensitivity in the local burning rate on both types of stability thresholds. It is found that for sufficiently small values of this parameter, there exists a stable range of pressure sensitivities for steady, planar burning such that the classical cellular form of hydrodynamic instability and the more recent pulsating form of hydrodynamic instability can each occur as the corresponding stability threshold is crossed. For larger thermal sensitivities, however, the pulsating stability boundary evolves into a C-shaped curve in the disturbance-wavenumber/ pressure-sensitivity plane, indicating loss of stability to pulsating perturbations for all sufficiently large disturbance wavelengths. It is thus concluded, based on characteristic parameter values, that an equally likely form of hydrodynamic instability in liquid-propellant combustion is of a nonsteady, long-wave nature, distinct from the steady, cellular form originally predicted by Landau.

Margolis, S. B.↗
NASA NTRS1999

Pulsating Hydrodynamic Instability in a Dynamic Model of Liquid-Propellant Combustion

Hydrodynamic (Landau) instability in combustion is typically associated with the onset of wrinkling of a flame surface, corresponding to the formation of steady cellular structures as the stability threshold is crossed. In the context of liquid-propellant combustion, such instability has recently been shown to occur for critical values of the pressure sensitivity of the burning rate and the disturbance wavenumber, significantly generalizing previous classical results for this problem that assumed a constant normal burning rate. Additionally, however, a pulsating form of hydrodynamic instability has been shown to occur as well, corresponding to the onset of temporal oscillations in the location of the liquid/gas interface. In the present work, we consider the realistic influence of a nonzero temperature sensitivity in the local burning rate on both types of stability thresholds. It is found that for sufficiently small values of this parameter, there exists a stable range of pressure sensitivities for steady, planar burning such that the classical cellular form of hydrodynamic instability and the more recent pulsating form of hydrodynamic instability can each occur as the corresponding stability threshold is crossed. For larger thermal sensitivities, however, the pulsating stability boundary evolves into a C-shaped curve in the disturbance-wavenumber/ pressure-sensitivity plane, indicating loss of stability to pulsating perturbations for all sufficiently large disturbance wavelengths. It is thus concluded, based on characteristic parameter values, that an equally likely form of hydrodynamic instability in liquid-propellant combustion is of a nonsteady, long-wave nature, distinct from the steady, cellular form originally predicted by Landau.

Margolis, Stephen B.↗
NASA NTRS1999

Pulsating Hydrodynamic Instability and Thermal Coupling in an Extended Landau/Levich Model of Liquid-Propellant Combustion: Inviscid Analysis - 1

Hydrodynamic (Landau) instability in combustion is typically associated with the onset of wrinkling of a flame surface, corresponding to the formation of steady cellular structures as the stability threshold is crossed. In the context of liquid-propellant combustion, such instability has recently been shown to occur for critical values of the pressure sensitivity of the burning rate and the disturbance wavenumber, significantly generalizing previous classical results for this problem that assumed a constant normal burning rate. Additionally, however, a pulsating form of hydrodynamic instability has been shown to occur as well, corresponding to the onset of temporal oscillations in the location of the liquid/gas interface. In the present work, we consider the realistic influence of a non-zero temperature sensitivity in the local burning rate on both types of stability thresholds. It is found that for sufficiently small values of this parameter, there exists a stable range of pressure sensitivities for steady, planar burning such that the classical cellular form of hydrodynamic instability and the more recent pulsating form of hydrodynamic instability can each occur as the corresponding stability threshold is crossed. For larger thermal sensitivities, however, the pulsating stability boundary evolves into a C-shaped curve in the (disturbance-wavenumber, pressure-sensitivity) plane, indicating loss of stability to pulsating perturbations for all sufficiently large disturbance wavelengths. It is thus concluded, based on characteristic parameter values, that an equally likely form of hydrodynamic instability in liquid-propellant combustion is of a non-steady, long-wave nature, distinct from the steady, cellular form originally predicted by Landau.

Margolis, Stephen B.↗