SEARCH · Search NASA

Results for “Statistical 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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1984

Theory and modeling of atmospheric turbulence, part 1

The cascade transfer which is the only function to describe the mode coupling as the result of the nonlinear hydrodynamic state of turbulence is discussed. A kinetic theory combined with a scaling procedure was developed. The transfer function governs the non-linear mode coupling in strong turbulence. The master equation is consistent with the hydrodynamical system that describes the microdynamic state of turbulence and has the advantages to be homogeneous and have fewer nonlinear terms. The modes are scaled into groups to decipher the governing transport processes and statistical characteristics. An equation of vorticity transport describes the microdynamic state of two dimensional, isotropic and homogeneous, geostrophic turbulence. The equation of evolution of the macrovorticity is derived from group scaling in the form of the Fokker-Planck equation with memory. The microdynamic state of turbulence is transformed into the Liouville equation to derive the kinetic equation of the singlet distribution in turbulence. The collision integral contains a memory, which is analyzed with pair collision and the multiple collision. Two other kinetic equations are developed in parallel for the propagator and the transition probability for the interaction among the groups.

Source record↗

NASA NTRS1978

On the bulk viscosity of relativistic matter

It is shown that the Kubo formulation leads to an expression for bulk viscosity involving the differences of the trace equilibrium and nonequilibrium hydrodynamic tensors just as in the Boltzmann formulation. It is argued that, if a dense system is represented by a quantum field, the field-theoretic energy-momentum tensor and, therefore, the hydrodynamic tensor should always have a vanishing trace in the limit of high energies. This bulk viscosity should vanish in the same limit. Finally, it is noted that the explicit use of well-known statistical mechanical formulation can be considered as a justification of Weinberg's intuitive argument from a fundamental point of view.

Canuto, V.↗

NASA NTRS1977

Atmospheric study relating to pad lift-off and entry landing

A relationship between the atmospheric general circulation and geophysical hydrodynamic experiments was sought by attempting to find a relationship between wave number and temperature gradient at mid-latitudes at 500 mb. To this end data were gathered from four winter seasons and analyzed. The statistical analysis failed to provide convincing support for the hypothesis of a direct relationship between wave number and temperature gradient, although an indication that the transient waves may be so related was noted.

King, R. L.↗

NASA NTRS1998

Statistical Properties of Line Centroid Velocity Increments in the rho Ophiuchi Cloud

We present a comparison of histograms of CO (2-1) line centroid velocity increments in the rho Ophiuchi molecular cloud with those computed for spectra synthesized from a three-dimensional, compressible, but non-starforming and non-gravitating hydrodynamic simulation. Histograms of centroid velocity increments in the rho Ophiuchi cloud show clearly non-Gaussian wings, similar to those found in histograms of velocity increments and derivatives in experimental studies of laboratory and atmospheric flows, as well as numerical simulations of turbulence. The magnitude of these wings increases monotonically with decreasing separation, down to the angular resolution of the data. This behavior is consistent with that found in the phase of the simulation which has most of the properties of incompressible turbulence. The time evolution of the magnitude of the non-Gaussian wings in the histograms of centroid velocity increments in the simulation is consistent with the evolution of the vorticity in the flow. However, we cannot exclude the possibility that the wings are associated with the shock interaction regions. Moreover, in an active starforming region like the rho Ophiuchi cloud, the effects of shocks may be more important than in the simulation. However, being able to identify shock interaction regions in the interstellar medium is also important, since numerical simulations show that vorticity is generated in shock interactions.

Lis, D. C.↗

NASA NTRS2021

Atacama Cosmology Telescope: Combined kinematic and thermal Sunyaev-Zel’dovich measurements from BOSS CMASS and LOWZ halos

The scattering of cosmic microwave background (CMB) photons off the free-electron gas in galaxies and clusters leaves detectable imprints on high resolution CMB maps: the thermal and kinematic Sunyaev-Zel’dovich effects (tSZ and kSZ respectively). We use combined microwave maps from the Atacama Cosmology Telescope DR5 and Planck in combination with the CMASS (mean redshifthzi¼0.55and host halo masshMviri¼3×1013M⊙) and LOWZ (hzi¼0.31,hMviri¼5×1013M⊙) galaxy catalogs from the Baryon Oscillation Spectroscopic Survey (BOSS DR10 and DR12), to study the gas associated with these galaxy groups. Using individual reconstructed velocities, we perform a stacking analysis and reject the no-kSZ hypothes is at 6.5σ, the highest significance to date. This directly translates into a measurement of the electron number density profile, and thus of the gas density profile. Despite the limited signal to noise, the measurement shows at high significance that the gas density profile is more extended than the dark matter density profile, for any reasonable baryon abundance (formally>90σfor the cosmic baryon abundance). We simultaneously measure the tSZ signal, i.e., the electron thermal pressure profile of the same CMASS objects, and reject theno-tSZ hypothesis at10σ. We combine tSZ and kSZ measurements to estimate the electron temperature to20% precision in several aperture bins, and find it comparable to the virial temperature. In a companion paper, we analyze these measurements to constrain the gas thermodynamics and the properties of feedback inside galaxy groups. We present the corresponding LOWZ measurements in this paper, ruling out a null kSZ (tSZ)signal at 2.9ð13.9Þσ, and leave their interpretation to future work. This paper and the companion paper demonstrate that current CMB experiments can detect and resolve gas profiles in low mass halos and at high redshifts, which are the most sensitive to feedback in galaxy formation and the most difficult to measure any other way. They will be a crucial input to cosmological hydrodynamical simulations, thus improving our understanding of galaxy formation. These precise gas profiles are already sufficient to reduce the main limiting theoretical systematic in galaxy-galaxy lensing: baryonic uncertainties. Future such measurements will thus unleash the statistical power of weak lensing from the Rubin, Euclid and Roman observatories. Our stacking software Thumb Stack is publicly available and directly applicable to future Simons Observatory andCMB-S4 data.

Emmanuel Schaan↗

DOE OSTI2024

Comparative study of machine learning techniques for post-combustion carbon capture systems

Computational analysis of countercurrent flows in packed absorption columns, often used in solvent-based post-combustion carbon capture systems (CCSs), is challenging. Typically, computational fluid dynamics (CFD) approaches are used to simulate the interactions between a solvent, gas, and column's packing geometry while accounting for the thermodynamics, kinetics, heat, and mass transfer effects of the absorption process. These simulations can then be used explain a column's hydrodynamic characteristics and evaluate its CO 2 -capture efficiency. However, these approaches are computationally expensive, making it difficult to evaluate numerous designs and operating conditions to improve efficiency at industrial scales. In this work, we comprehensively explore the application of statistical ML methods, convolutional neural networks (CNNs), and graph neural networks (GNNs) to aid and accelerate the scale-up and design optimization of solvent-based post-combustion CCSs. We apply these methods to CFD datasets of countercurrent flows in absorption columns with structured packings characterized by several geometric parameters. We train models to use these parameters, inlet velocity conditions, and other model-specific representations of the column to estimate key determinants of CO 2 -capture efficiency without having to simulate additional CFD datasets. We also evaluate the impact of different input types on the accuracy and generalizability of each model. We discuss the strengths and limitations of each approach to further elucidate the role of CNNs, GNNs, and other machine learning approaches for CO 2 -capture property prediction and design optimization.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2021

Accuracy Assessment of Global Internal-Tide Models Using Satellite Altimetry Models Using Satellite Altimetry

Altimeter measurements are corrected for several geophysical parameters in order to access ocean signals of interest, like mesoscale or sub-mesoscale variability. The ocean tide is one of the most critical corrections due to the amplitude of the tidal elevations and to the aliasing phenomena of high-frequency signals into the lower-frequency band, but the internal-tide signatures at the ocean surface are not yet corrected globally. Internal tides can have a signature of several centimeters at the surface with wavelengths of about 50–250 km for the first mode and even smaller scales for higher-order modes. The goals of the upcoming Surface Water Ocean Topography (SWOT) mission and other high-resolution ocean measurements make the correction of these small-scale signals a challenge, as the correction of all tidal variability becomes mandatory to access accurate measurements of other oceanic signals. In this context, several scientific teams are working on the development of new internal-tide models, taking advantage of the very long altimeter time series now available, which represent an unprecedented and valuable global ocean database. The internal-tide models presented here focus on the coherent internal-tide signal and they are of three types: empirical models based upon analysis of existing altimeter missions, an assimilative model and a three-dimensional hydrodynamic model. A detailed comparison and validation of these internal-tide models is proposed using existing satellite altimeter databases. The analysis focuses on the four main tidal constituents: M2, K1, O1 and S2. The validation process is based on a statistical analysis of multi-mission altimetry including Jason-2 and Cryosphere Satellite-2 data. The results show a significant altimeter variance reduction when using internal-tide corrections in all ocean regions where internal tides are generating or propagating. A complementary spectral analysis also gives some estimation of the performance of each model as a function of wavelength and some insight into the residual non-stationary part of internal tides in the different regions of interest. This work led to the implementation of a new internal-tide correction (ZARON'one) in the next geophysical data records version-F (GDR-F) standards.

Loren Carrere↗

NASA NTRS2015

Hazard Assessment from Storm Tides and Rainfall on a Tidal River Estuary

Here, we report on methods and results for a model-based flood hazard assessment we have conducted for the Hudson River from New York City to Troy/Albany at the head of tide. Our recent work showed that neglecting freshwater flows leads to underestimation of peak water levels at up-river sites and neglecting stratification (typical with two-dimensional modeling) leads to underestimation all along the Hudson. As a result, we use a three-dimensional hydrodynamic model and merge streamflows and storm tides from tropical and extratropical cyclones (TCs, ETCs), as well as wet extratropical cyclone (WETC) floods (e.g. freshets, rain-on-snow events). We validate the modeled flood levels and quantify error with comparisons to 76 historical events. A Bayesian statistical method is developed for tropical cyclone streamflows using historical data and consisting in the evaluation of (1) the peak discharge and its pdf as a function of TC characteristics, and (2) the temporal trend of the hydrograph as a function of temporal evolution of the cyclone track, its intensity and the response characteristics of the specific basin. A k-nearest-neighbors method is employed to determine the hydrograph shape. Out of sample validation tests demonstrate the effectiveness of the method. Thus, the combined effects of storm surge and runoff produced by tropical cyclones hitting the New York area can be included in flood hazard assessment. Results for the upper Hudson (Albany) suggest a dominance of WETCs, for the lower Hudson (at New York Harbor) a case where ETCs are dominant for shorter return periods and TCs are more important for longer return periods (over 150 years), and for the middle-Hudson (Poughkeepsie) a mix of all three flood events types is important. However, a possible low-bias for TC flood levels is inferred from a lower importance in the assessment results, versus historical event top-20 lists, and this will be further evaluated as these preliminary methods and results are finalized. Future funded work will quantify the influences of sea level rise and flood adaptation plans (e.g. surge barriers). It would also be valuable to examine how streamflows from tropical cyclones and wet cool-season storms will change, as this factor will dominate at upriver locations.

Hazard assessment↗

DOE OSTI2024

Hydrodynamic fluctuations near a Hopf bifurcation: Stochastic onset of vortex shedding behind a circular cylinder

Here, we investigate hydrodynamic fluctuations in the flow past a circular cylinder near the critical Reynolds number Re c for the onset of vortex shedding. Starting from the fluctuating Navier-Stokes equations, we perform a perturbation expansion around Re c to derive analytical expressions for the statistics of the fluctuating lift force. Molecular-level simulations using the direct simulation Monte Carlo method support the theoretical predictions of the lift power spectrum and amplitude distribution. Notably, we have been able to collect sufficient statistics at distances Re ⁡/ Re c – 1 = O ⁡(10 –3 ) from the instability that confirm the appearance of non-Gaussian fluctuations, and we observe that they are associated with intermittent vortex shedding. These results emphasize how unavoidable thermal-noise-induced fluctuations become dramatically amplified in the vicinity of oscillatory flow instabilities and that their onset is fundamentally stochastic.

42 ENGINEERING↗

DOE OSTI2024

Do not forget the electrons: Extending moderately-sized nuclear networks for multidimensional hydrodynamic codes

Context.Nuclear networks are widely used coupled with hydrodynamical simulations of explosive scenarios to account for the change of nuclear species and energy generation rate due to nuclear reactions. In this way, there is a feedback mechanism between the hydrodynamical state and the nuclear processes. Unfortunately, the timescale of nuclear reactions is orders of magnitude smaller than the dynamical timescale that drives hydrodynamical simulations. Therefore, these nuclear networks are usually very small, reduced in most cases to a dozen elements, especially when simulations are carried out in more than one dimension. Aims.We present here an extended nuclear network, with 90 species, designed for being coupled with hydrodynamic simulations, which includes neutrons, protons, electrons, positrons, and the corresponding neutrino and anti-neutrino emission. This network is also coupled with temperature, making it extremely robust and, together with its size, unique of its kind. The inclusion of electron captures on free protons makes the network very appropriate for multidimensional studies of Type Ia supernova explosions, especially when the exploding object is a massive white dwarf. Methods.We perform several tests that are relevant to simulate explosive scenarios, such as Type Ia supernovae and core-collapse supernovae. We compare the results of the 90 nuclei network with a standardα-chain network with 14 elements to evaluate the differences in the energy generation rate. We also evaluate the relevance of including the electrons in the network in terms of generated yields and how it affects the pressure of a degenerate fluid such as that of white dwarfs. The results obtained with the 90-nuclei network have been verified with a much larger 2000-nuclei network built from REACLIB (WinNet), in terms of nuclear energy generation rate, pressure, and produced yields. Results.The results obtained with the proposed medium-sized network compare fairly well, to a few percent, with those computed withWinNetin scenarios reproducing the gross physical conditions of current Type Ia supernova explosion models. In those cases where the carbon and oxygen fuel ignites at high density, the high-temperature plateau typical of the nuclear statistical equilibrium regime is well defined and stable, allowing large integration time steps. We show that the inclusion of electron captures on free protons substantially improves the estimation of the electron fraction of the mixture. Therefore, the pressure is better determined than in networks where electron captures are excluded, which will ultimately lead to more reliable hydrodynamic models. Explosive combustion of helium at low density, occurring near the surface layer of a white dwarf, is also better described with the proposed network, which gives nuclear energy generation rates much closer toWinNetthan typical reduced alpha networks. Conclusions.A nuclear network withN= 90 species, including electrons, aimed at multidimensional calculations of supernova explosions is described and verified. The proposed network is suitable for the study of Type Ia supernova explosions because it provides better values of pressure and electron abundance than other existing networks with smaller or even a similar size but without including electron capture processes.

Astronomy & Astrophysics↗

DOE OSTI2024

On the transferability of residence time distributions in two 10-km long river sections with similar hydromorphic units

Quantifying hydrologic exchange fluxes (HEFs) at the stream-groundwater interface and their residence time distributions (RTDs) in the subsurface are important for managing the water quality and ecosystem health in dynamic river corridors. However, direct simulating high-spatial resolution HEFs and RTDs can be time-consuming, especially for watershed-scale modeling. Efficient surrogate models linking RTDs to hydromorphic units (HUs) can be alternatives for simulating RTDs in large-scale models. A common concern of these surrogate models, though, is the transferability of the relationship between the RTDs and HUs from one river corridor to another. To address this issue, this work evaluates the HEFs and resulting RTD-HU relationships for two 10-km long river corridors along the Columbia River leveraging a one-way coupled three-dimensional transient surface-subsurface water transport modeling framework we previously developed. Applying such a framework at the two river corridors with similar HUs allows for quantitative comparisons of HEFs and RTDs using both statistical tests and machine learning classification models. Finally, our comparison shows that the similarity and transferability of the RTD-HU relationship is very low for the two investigated river sections, which suggests that devising a general algorithm to estimate RTDs based solely on surface water hydrodynamics and short-distance river channel topography data, as well as HU classification, might be nearly impossible.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1990

Design optimization and probabilistic analysis of a hydrodynamic journal bearing

A nonlinear constrained optimization of a hydrodynamic bearing was performed yielding three main variables: radial clearance, bearing length to diameter ratio, and lubricating oil viscosity. As an objective function a combined model of temperature rise and oil supply has been adopted. The optimized model of the bearing has been simulated for population of 1000 cases using Monte Carlo statistical method. It appeared that the so called 'optimal solution' generated more than 50 percent of failed bearings, because their minimum oil film thickness violated stipulated minimum constraint value. As a remedy change of oil viscosity is suggested after several sensitivities of variables have been investigated.

Liniecki, Alexander G.↗

DOE OSTI2024

The shape of convection in 2D and 3D global simulations of stellar interiors

Context. Theoretical descriptions of convective overshooting in stellar interiors often rely on a basic one-dimensional parameterization of the flow called the filling factor for convection. Several different definitions of the filling factor have been developed for this purpose, based on: (1) the percentage of the volume, (2) the mass flux, and (3) the convective flux that moves through the boundary. Aims. We examine these definitions of the filling factor with the goal of establishing their ability to explain differences between 2D and 3D global simulations of stellar interiors that include fully compressible hydrodynamics and realistic microphysics for stars. Methods. We study convection and overshooting in pairs of identical two-dimensional (2D) and three-dimensional (3D) global simulations of stars produced with MUSIC , a fully compressible, time-implicit hydrodynamics code. We examine pairs of simulations for (1) a 3 M ⊙ red giant star near the first dredge-up point, (2) a 1 M ⊙ pre-main-sequence star with a large convection zone, (3) the current sun, and (4) a 20 M ⊙ main-sequence star with a large convective core. Results. Our calculations of the filling factor based on the volume percentage and the mass flux indicate asymmetrical convection near the surface for each star with an outer convection zone. However, near the convective boundary, convective flows achieve inward-outward symmetry for each star that we study; for 2D and 3D simulations, these filling factors are indistinguishable. A filling factor based on the convective flux is contaminated by boundary-layer-like flows, making a theoretical interpretation difficult. We present two possible new alternatives to these frequently used definitions of a filling factor, which instead compare flows at two different radial points. The first alternative is the penetration parameter of Anders et al. (2022, ApJ, 926, 169). The second alternative is a new statistic that we call the plume interaction parameter. We demonstrate that both of these parameters captures systematic differences between 2D and 3D simulations around the convective boundary.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Modeling the Cosmological Lyman-𝛼 Forest at the Field Level

The distribution of absorption lines in the spectra of distant quasars, called the Lyman-𝛼 (Ly-𝛼) forest, is a unique probe of cosmology and the intergalactic medium at high redshifts and small scales. The statistical power of ongoing redshift surveys demands precise theoretical tools to model the Ly-𝛼 forest. We address this challenge by developing an analytic, perturbative forward model to predict the Ly-𝛼 forest at the field level for a given set of cosmological initial conditions. Our model shows a remarkable performance when compared with the Sherwood hydrodynamic simulations: it reproduces the Ly-𝛼 forest flux power spectrum, its cross-correlation with dark matter halos, and the one-point probability distribution function of both fields at the percent level down to scales of a few Mpc. Our work provides crucial tools that bridge analytic modeling on large scales with simulations on small scales, enabling field-level inference from Ly-𝛼 forest data and simulation-based priors for cosmological analyses. Furthermore, this is especially timely for realizing the full scientific potential of the Ly-𝛼 forest measurements by the dark energy spectroscopic instrument.

Cosmological parameters↗

DOE OSTI2026

ELG×LRG Distribution through Dark Matter Halo Dynamics

We investigate the clustering and halo occupation distribution (HOD) of DESI Y1 emission-line (ELGs) and luminous red (LRGs) galaxies at 0.8 < z < 1.1, including their cross-correlation (ELG×LRG), using the A BACUS S UMMIT suite and a new Halo Occupation Model (H OME ) for galaxy multitracers. This integrates intrahalo dynamics, halo exclusion, and quenching, bridging insights from hydrodynamical, HOD, abundance-matching, and semianalytic studies. Leveraging full phase-space information from the Uchuu N-body simulation, and sampling satellites from dark-matter particle positions via physically motivated prescriptions, Home reproduces the anisotropic clustering down to s = 200 h −1 kpc with unprecedented accuracy. Model parameters are inferred solely from two-point statistics using a two-level Bayesian framework, yielding high-fidelity ELG, LRG, and cross-reference catalogs. We find that satellite ELGs behave as incoherent flows within their parent halos, dominating the clustering below 4 h −1 Mpc. The HOD from the best-fit Home has the following properties: (i) 90.50% (85.91%) of ELGs (LRGs) are central galaxies without satellites, residing in halos of M vir ∼ 6.6 × 10 11 (1.2 × 10 13 ) h −1 M ⊙ ; (ii) the ELG×LRG cross-correlation is governed by central-central pairs and shaped by halo exclusion on 2–5 h −1 Mpc scales; (iii) 9.50% (14.09%) of ELGs (LRGs) are satellites, of which 1.09% (3.52%) inhabit halos with a central galaxy of the same species in a maximally conformal configuration, 7.02% (0.005%) orbit complementary hosts in a minimally conformal state, and 0.58% (10.57%) are orphans. The high sensitivity of Home precisely captures the dynamics of satellites in different host environments, opening a promising avenue for understanding systematics and the dynamical nature of dark matter, potentially distinguishing gravity models.

Favole, Ginevra [Universidad de La Laguna (Spain);↗

DOE OSTI2025

Emergent hydrodynamic mode on SU(2) plaquette chains and quantum simulation

We search for emergent hydrodynamic modes in real-time Hamiltonian dynamics of 2+1-dimensional SU(2) lattice gauge theory on a quasi-one-dimensional plaquette chain, by numerically computing symmetric correlation functions of energy densities on lattice sizes of about 20 with the local Hilbert space truncated at 𝑗 max = $\frac{1}{2}$. Because of the Umklapp processes, we only find a mode for energy diffusion. The symmetric correlator exhibits transport peak near zero frequency with a width approximately proportional to momentum squared at small momentum, when the system is fully quantum ergodic, as indicated by the eigenenergy level statistics. This transport peak leads to a power-law 𝑡 −$\frac{1}{2}$ decay of the symmetric correlator at late time, also known as the long-time tail, as well as diffusionlike spreading in position space. We also introduce a quantum algorithm for computing the symmetric correlator on a quantum computer and find it gives results consistent with exact diagonalization when tested on the IBM emulator. Finally we discuss the future prospect of searching for the sound modes.

Hamiltonian systems↗

NASA NTRS1975

On the lower mass limit for the carbon detonation scenario

The lower mass limit for carbon detonation (about 4 solar masses) is due to an argument by Paczynski (1970) that stars in the range between 1.4 and 3.5 solar masses will lose their envelopes prior to carbon detonation. This argument is based on two criteria for mass loss: (1) the envelope is unstable to adiabatic oscillations, and (2) the total energy of the envelope is positive. Recent hydrodynamic calculations of red-giant envelopes have shown that the envelopes are highly nonadiabatic and that the second criterion is not a sufficient condition for mass loss. Thus, a calculation which takes into account the nonadiabatic nature of the envelope is required in order to determine the lower mass limit for carbon detonation. This lower mass limit is directly related to the formation mechanisms and statistics of supernovae and pulsars.

Endal, A. S.↗