SEARCH · Search NASA

Results for “Mixed Compression”

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 1,099 records · Page 61

All sources NASA NTRS DOE OSTI

NASA NTRS1993

The compositional structure and the effects of exothermicity in a nonpremixed planar jet flame

Results are presented of direct numerical simulation (DNS) of a randomly perturbed compressible, spatially developing two-dimensional planar jet under the influence of a finite rate chemical reaction of the type F + O yields product. The objectives of the simulations are to assess the compositional structure of the flame and to determine the influence of reaction exothermicity by means of statistical sampling of the DNS generated data. It is shown that even with this idealized kinetics model the simulated results exhibit features in accord with experimental data. These results indicate that the Damkohler number is an important parameter in determining the statistical composition of the reacting field and that the results are not very sensitive to the mechanism by which this parameter is varied. It is demonstrated that as the intensity of mixing is increased and the effect of finite rate chemistry is more pronounced, the magnitudes of the ensemble mean and variance of the product mass fraction decrease and those of the reactants' mass fraction increase. Also, at higher mixing rates the joint probability density functions of the reactants' mass fractions shift towards higher values within the composition domain, indicating a lower reactedness.

Steinberger, C. J.↗

NASA NTRS1988

Data compression using adaptive transform coding. Appendix 1: Item 1

Adaptive low-rate source coders are described in this dissertation. These coders adapt by adjusting the complexity of the coder to match the local coding difficulty of the image. This is accomplished by using a threshold driven maximum distortion criterion to select the specific coder used. The different coders are built using variable blocksized transform techniques, and the threshold criterion selects small transform blocks to code the more difficult regions and larger blocks to code the less complex regions. A theoretical framework is constructed from which the study of these coders can be explored. An algorithm for selecting the optimal bit allocation for the quantization of transform coefficients is developed. The bit allocation algorithm is more fully developed, and can be used to achieve more accurate bit assignments than the algorithms currently used in the literature. Some upper and lower bounds for the bit-allocation distortion-rate function are developed. An obtainable distortion-rate function is developed for a particular scalar quantizer mixing method that can be used to code transform coefficients at any rate.

Rost, Martin Christopher↗

NASA NTRS1973

The behavior of a compressible turbulent boundary layer in a shock-wave-induced adverse pressure gradient

The results of an experimental investigation of the mean- and fluctuating-flow properties of a compressible turbulent boundary layer in a shock-wave-induced adverse pressure gradient are presented. The turbulent boundary layer developed on the wall of an axially symmetric nozzle and test section whose nominal free-stream Mach number and boundary-layer thickness Reynolds number were 4 and 100,000, respectively. The adverse pressure gradient was induced by an externally generated conical shock wave. Mean and time-averaged fluctuating-flow data, including the complete experimental Reynolds stress tensor and experimental turbulent mass- and heat-transfer rates are presented for the boundary layer and external flow, upstream, within and downstream of the pressure gradient. The mean-flow data include distributions of total temperature throughout the region of interest. The turbulent mixing properties of the flow were determined experimentally with a hot-wire anemometer. The calibration of the wires and the interpretation of the data are discussed. From the results of the investigation, it is concluded that the shock-wave - boundary-layer interaction significantly alters the turbulent mixing characteristics of the boundary layer.

Rose, W. C.↗

DOE OSTI2018

Final Scientific / Technical Report: Enabling Technologies for Heavy-Duty Vehicles – Cummins 55BTE

The Cummins 55% BTE (55BTE) program has completed the planned technical work on the project. This work includes the planned engine system demonstration in pursuit of the goal of demonstrating a peak system brake thermal efficiency (BTE) of 55%. The engine system included a high efficiency diesel engine integrated with a state-of-the-art waste heat recovery (WHR) system and an advanced aftertreatment system capable of meeting the current emissions standards. While the ultimate program goal of 55% BTE was not fully achieved due to hardware issues during the final testing phase, the program demonstrated a significant increase in reported engine system BTE for a heavy duty sized engine. The previous demonstrations in the Department of Energy funded SuperTruck 1 programs ranged between 50% - 51% BTE. The Cummins 55% BTE program demonstrated 54% BTE. Additionally, the program established a revised path-to-target showing how the system could be improved to reach the ultimate program goal of 55% BTE with some minor modification to the engine system. The changes in the revised path-to-target were unable to be completed during the course of this program due to time and money constraints placed on the program. The program’s goals were challenging in both scope and timing. Although the program’s goal can be simply stated in terms of demonstrating a system efficiency of 55%, the achievement of a heavy-duty diesel engine capable of full torque curve operation, adequate transient performance, low emissions and high brake thermal efficiency (>51% BTE) had not been previously demonstrated in a heavy-duty engine. The closest efforts to date included the SuperTruck 1 programs referenced earlier. The program made several advancements in all system areas. The combustion system was re-designed for a shorter combustion duration and lower in-cylinder heat loss. This was achieved through optimization of the fuel injection rate shape, number of spray holes, piston bowl shape, compression ratio, piston oil cooling, heat flow through the piston and in-cylinder charge motion. The air handling system was re-designed to provide cooled EGR at virtually no pumping penalty. This was achieved through implementation of a dual loop EGR system, reduction of EGR system pressure drop, and implementation of advanced turbocharger efficiency technologies. The engine friction and parasitic signature was dramatically reduced by the program through adoption of variable flow pumps, advanced rings and coatings, rollerized valvetrain, adoption of low viscosity lubricants and reduction in cylinder line bore distortion. The aftertreatment system was optimized through use of low dP substrates and an ammonia gas injection system. The WHR system development included the use of a dual-entry turbine and a mixed charger cooler. Finding and implementing these solutions in a short two-year program at the budgeted funding provided the greatest challenge for the program.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1993

A survey of low frequency waves at Jupiter: The Ulysses encounter

We report the results of a survey of low-frequency (LF) plasma waves detected during the Ulysses Jupiter flyby. In the Jovian foreshock, two predominant wave periods are detected: 10(exp 2)-s and 5-s, as measured in the spacecraft frame. The 10(exp 2)-s waves are highly nonlinear propagate at large angles to vector-B(sub 0) (typically 50 deg), are steepened, and sometimes have attached whistler packets. For the interval analyzed the 10(exp 2)-s waves had mixed right-and left-hand polarizations. We argue that these are all consistent with being right-hand magnetosonic waves in the solar wind frame. The 10(exp 2)-s waves with attached whistler are similar to cometary waves. The trailing portions are linearly polaraized and the whistler portions circularly polarized with amplitudes decreasing linearly with time. The emissions are generated by approximately 2-keV protons flowing from the Jovian bow shock/magnetosheath into the upstream region. The instability is the ion beam instability. Higher Z ions were considered as a source of the waves but have been ruled out because of the low sunward velocities needed for their resonance. The 5-s waves have delta vector-B/B(sub 0 approximately = 0.5, are compressive and are left-hand polarized in the spacecraft frame. Local generation by three different resonant interactions were considered and have been ruled out. One possibility is that these waves are whistler mode by-products of the steepened lower-frequency magnetosonic waves. Mirror mode structures were detected throughout the outbound magnetosheath passes. For these structures, the theta(sub kB) values were consistently in the range of 80 deg to 90 deg, exceptionally high values.

Tsurutani, Bruce T.↗

DOE OSTI2025

Uncertainty Visualization of Critical Points of 2D Scalar Fields for Parametric and Nonparametric Probabilistic Models

This paper presents a novel end-to-end framework for closed-form computation and visualization of critical point uncertainty in 2D uncertain scalar fields. Critical points are fundamental topological descriptors used in the visualization and analysis of scalar fields. The uncertainty inherent in data (e.g., observational and experimental data, approximations in simulations, and compression), however, creates uncertainty regarding critical point positions. Uncertainty in critical point positions, therefore, cannot be ignored, given their impact on downstream data analysis tasks. Here, in this work, we study uncertainty in critical points as a function of uncertainty in data modeled with probability distributions. Although Monte Carlo (MC) sampling techniques have been used in prior studies to quantify critical point uncertainty, they are often expensive and are infrequently used in production-quality visualization software. We, therefore, propose a new end-to-end framework to address these challenges that comprises a threefold contribution. First, we derive the critical point uncertainty in closed form, which is more accurate and efficient than the conventional MC sampling methods. Specifically, we provide the closed-form and semianalytical (a mix of closed-form and MC methods) solutions for parametric (e.g., uniform, Epanechnikov) and nonparametric models (e.g., histograms) with finite support. Second, we accelerate critical point probability computations using a parallel implementation with the VTK-m library, which is platform portable. Finally, we demonstrate the integration of our implementation with the ParaView software system to demonstrate near-real-time results for real datasets.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1995

The microscopic state of the solar wind: Links between composition, velocity distributions, waves and turbulence

An overview is given of the microscopic state of the solar wind with emphasis on recent Ulysses high-latitude observations and previous Helios in-ecliptic observations. The possible links between composition, ionization state. velocity distribution functions of electrons, protons and heavy ions. kinetic plasma waves and MHD-scale turbulence are elaborated. Emphasis is placed on a connection of interplanetary kinetic-scale phenomena with their generating microscopic processes in the corona. The fast streams seem to consist of mesoscale pressure-balanced plasma filaments and magnetic flux tubes, reminiscent of the supergranular-size structures building the open corona, from which copious Alfven waves emanate. The wind from the magnetically structured and active corona shows developed compressive turbulence and considerable abundance and ionization state variations. Some modelling attempts to explain the observed element fractionation are briefly reviewed. The causes of the nonthermal particle features. such as proton-ion differential streaming, ion beams. temperature anisotropies, and skewed distributions associated with collisionless heat conduction, are ultimately to be searched in the fact, that the corona is never quiescent but fundamentally variable in space and time. Consequently, the radial evolution of the internal state of the wind resembles at all latitudes a complicated relaxation process, in the course of which the free (in comparison with LTE conditions) particle kinetic energy is converted into plasma waves and MHD turbulence on a wide range of scales. This leads to intermittent wave-particle interactions and unsteady anomalous transport, mixed with the weak effects of the rare Coulomb collisions. Spherical expansion and large-scale inhomogeneity forces the wind to attain microscopically a complex internal state of dynamic statistical equilibrium.

Marsch, E.↗

DOE OSTI2022

Mechanical Response and Processability of Wet-Laid Recycled Carbon Fiber PE, PA66 and PET Thermoplastic Composites

The interest in recycled carbon fiber (rCF) is growing rapidly and the supply chain for these materials is gradually being established. However, the processing routes, material intermediates and properties of rCF composites are less understood for designers to adopt them into practice. This paper provides a practical pathway for rCFs in conjunction with low cost and, for the most part, commodity thermoplastic resins, namely polyethylene (PE), polyamide 66 (PA66) and polyethylene terephthalate (PET). Industrially relevant wet-laid (WL) process routes have been adopted to produce mats using two variants of WL mats, namely (a) high speed wet-laid inclined wire to produce broad good ‘roll’ forms and (b) 3DEPTM process patented by Materials Innovation Technologies (MIT)-recycled carbon fiber (RCF), now Carbon Conversions, which involves mixing fibers and water and depositing the fibers on a water-immersed mold. These are referred to as ‘sheet’ forms. The produced mats were evaluated for their processing into composites as ‘fully consolidated mats’ and ‘non-consolidated’ as-produced mats. Comprehensive mechanical data in terms of tensile strength, tensile modulus and impact toughness for rCF C/PE, C/PA66 and C/PET are presented. The work is of high value to sustainable composite designers and modelers.

compression molding↗

NASA NTRS1996

Turbulent Compressible Convection with Rotation: Flow Structure and Evolution - Part 1

The effects of Coriolis forces on compressible convection are studied using three-dimensional numerical simulations carried out within a local modified f-plane model. The physics is simplified by considering a perfect gas occupying a rectilinear domain placed tangentially to a rotating sphere at various latitudes, through which a destabilizing heat flux is driven. The resulting convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers, evaluating conditions where the influence of rotation is both weak and strong. Given the computational demands of these high-resolution simulations, the parameter space is explored sparsely to ascertain the differences between laminar and turbulent rotating convection. The first paper in this series examines the effects of rotation on the flow structure within the convection, its evolution, and some consequences for mixing. Subsequent papers consider the large-scale mean shear flows that are generated by the convection, and the effects of rotation on the convective energetics and transport properties. It is found here that the structure of rotating turbulent convection is similar to earlier nonrotating studies, with a laminar, cellular surface network disguising a fully turbulent interior punctuated by vertically coherent structures. However, the temporal signature of the surface flows is modified by inertial motions to yield new cellular evolution patterns and an overall increase in the mobility of the network. The turbulent convection contains vortex tubes of many scales, including large-scale coherent structures spanning the full vertical extent of the domain involving multiple density scale heights. Remarkably, such structures align with the rotation vector via the influence of Coriolis forces on turbulent motions, in contrast with the zonal tilting of streamlines found in laminar flows. Such novel turbulent mechanisms alter the correlations which drive mean shearing flows and affect the convective transport properties. In contrast to this large-scale anisotropy, small-scale vortex tubes at greater depths are randomly orientated by the rotational mixing of momentum, leading to an increased degree of isotropy on the medium to small scales of motion there. Rotation also influences the thermodynamic mixing properties of the convection. In particular, interaction of the larger coherent vortices causes a loss of correlation between the vertical velocity and the temperature leaving a mean stratification which is not isentropic.

Brummell, Nicholas H.↗

DOE OSTI2018

The effects of boost pressure on stratification and burn duration of gasoline homogeneous charge compression ignition combustion

This article investigates the effects of intake pressure (boost) on the pre-ignition stratification and burn duration of homogeneous charge compression ignition combustion. Full cycle computational fluid dynamics simulations are performed with gasoline kinetics. An intake pressure sweep is performed while maintaining the same combustion timing and mean composition. The burn duration reduces with increasing boost, even though intake temperature is reduced to hold combustion timing constant. It is shown that the compositional stratification increases with boost whereas thermal stratification decreases. A quasi-dimensional model is employed to assess the effect of compositional stratification, pressure, mean temperature and isolate the effect of thermal stratification on burn duration. The analysis reveals that reducing charge temperature neutralizes the effect of increased boost on reactivity and the shorter burn durations at higher boost are primarily due to the lower thermal stratification. It is shown that higher pressures do not significantly increase the mixing and the lower thermal stratification is due to lower wall heat losses per unit charge mass. A follow-up set of non-reacting simulations with adiabatic walls corroborate this claim by revealing a constant magnitude of thermal stratification across the boost sweep.

Shingne, Prasad S.↗

DOE OSTI2023

Ultrafine-grained Fe-TiB 2 high-modulus nanocomposite steel with high strength and isotropic mechanical properties by laser powder bed fusion

Fe-TiB 2 metal matrix composite, also called high-modulus steels (HMSs), are of great interest for applications in fuel-efficient transportation infrastructure, aerospace, and wear industries due to their high specific stiffness and yield strength. However, conventional cast Fe-TiB 2 HMSs often contain coarse and sharp-edged TiB 2 particles which easily trigger premature cracking during loading. Here, we synthesized a Fe-TiB 2 nanocomposite HMS via laser powder bed fusion (LPBF) additive manufacturing of mixed micro-sized powders of Fe, Ti, and Fe 2 B. We investigated the microstructure formation and mechanical behavior of the Fe-TiB 2 HMS. We found that in situ chemical reaction of Ti and Fe 2 B enables the formation of TiB 2 particles at nanoscale during rapid solidification of LPBF. These nanoscale TiB 2 particles can serve as heterogeneous nucleation sites and promote the formation of ultrafine and equiaxed α-Fe grains with random crystallographic textures, which differ from many other additively manufactured (AM) metal alloys characteristic of strong crystallographic textures. As such, isotropic mechanical properties were achieved in the AM Fe-TiB 2 nanocomposite HMS with a high elastic modulus of ~ 240 GPa, an exceptional yield strength of ~ 1450 MPa, and a large plasticity of ~ 20% under compression. Quantitative analysis reveals that the high yield strength primarily originates from strengthening contributions of the ultrafine grains with an average grain size of ~450 nm, the nanoscale TiB 2 reinforcing particles of 20–180 nm, and a high density of printing-induced dislocations of the order of 10 15 m –2 . In situ synchrotron high-energy X-ray diffraction unveils the load partitioning from the softer α-Fe matrix to the stiffer and stronger TiB 2 nanoparticles, contributing to the sustained strain hardening during compression. Our work not only provides a general pathway for achieving high-performance metal matrix nanocomposites by in situ chemical reaction and precipitation of ceramic nanoparticles during additive manufacturing, but also offers mechanistic insights into the deformation mechanism of nanoparticle-reinforced HMS composites.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Solar Thermochemical Fuel Production via a Novel Low Pressure, Magnetically Stabilized, Non-volatile Iron Oxide Looping Process (University of Florida)

The UF Solar Fuels team has developed a modular reactor at the 10kW scale for the thermochemical conversion of H 2 O and CO 2 to H 2 and CO, respectively, using concentrated solar energy to drive high-temperature redox reactions. Various materials, including ferrites, doped and mixed metal ferrites, and ceria, were investigated and reported on. Ceria was chosen as the reactive material, with thermogravimetric and bench-scale analysis concurring with literature review and pointing towards reasonable kinetics and stability at temperatures on the order of 1500°C. A combined radiation, conduction, and species reaction/mass transport model was developed, utilizing Monte Carlo ray tracing, Lattice Boltzmann, and random walk particle tracking for the respective components to direct development and optimization of the reactor. With experimental data used for reaction rates, the model showed a path to efficiencies >20%, which could be market competitive with current PV-hydrolysis hydrogen generation systems. Economic analysis of a larger-scale plant (100kW modular system with cost performed at 500 modules) has been performed, assuming 0.5% and 10% internal rate of return. Without compression costs, H 2 production cost using this is technology at 12 USD/kg H 2 and 17 USD/kg H 2 for realized efficiencies of 20% and 15%, respectively.

14 SOLAR ENERGY↗

DOE OSTI2018

Magnetosheath Reconnection Before Magnetopause Reconnection Driven by Interplanetary Tangential Discontinuity: A Three‐Dimensional Global Hybrid Simulation With Oblique Interplanetary Magnetic Field

Abstract Terrestrial dayside dynamics associated with a southward turning, oblique interplanetary magnetic field (IMF) carried by an interplanetary tangential discontinuity (TD) is investigated by performing a three‐dimensional global‐scale hybrid simulation systematically for cases in which the incoming solar wind TD possesses various magnetic field rotation angles ΔΦ = 90° to 180° and half widths w = 2 d i 0 to w = 30 d i 0 , where d i 0 is the ion inertial length in the solar wind. Overall, the TD is compressed while being transmitted into the magnetosheath, with different compression processes downstream of the quasi‐parallel (Q‐||) and quasi‐perpendicular (Q‐ ⊥ ) shocks. It is found that magnetosheath reconnection may take place downstream of both the Q‐|| and Q‐ ⊥ shocks due to interaction of the directional TD with the bow shock and magnetopause, but the existence of magnetosheath reconnection depends on w and ΔΦ . Magnetosheath flux ropes are formed through three‐dimensional patchy reconnection in the thinned current sheet, with a longer rope lengths under a larger ΔΦ . There exists a dawn‐dusk asymmetry in the spatial extent of the flux ropes, which becomes more significant as ΔΦ decreases. When ΔΦ decreases to 90° , no reconnection flux ropes are found. Magnetopause reconnection is initiated when the magnetic fluxes with a southward turning IMF (on the sunward side) reach the magnetopause, and the magnetopause flux ropes can be mixed with the magnetosheath ones. Our simulation demonstrates that the effects of a southward turning of the IMF may not be a simple field direction change that leads to reconnection only at the magnetopause.

Guo, Zhifang↗

DOE OSTI2024

Dark energy survey year 3 results: likelihood-free, simulation-based w CDM inference with neural compression of weak-lensing map statistics

We present simulation-based cosmological wcold dark matter (wCDM) inference using dark energy survey year 3 weak-lensing maps, via neural data compression of weak-lensing map summary statistics: power spectra, peak counts, and direct map-level compression/inference with convolutional neural networks (CNN). Using simulation-based inference, also known as likelihood-free or implicit inference, we use forward-modelled mock data to estimate posterior probability distributions of unknown parameters. This approach allows all statistical assumptions and uncertainties to be propagated through the forward-modelled mock data; these include sky masks, non-Gaussian shape noise, shape measurement bias, source galaxy clustering, photometric redshift uncertainty, intrinsic galaxy alignments, non-Gaussian density fields, neutrinos, and non-linear summary statistics. We include a series of tests to validate our inference results. This paper also describes the Gower Street simulation suite: 791 full-sky pkdgrav3 dark matter simulations, with cosmological model parameters sampled with a mixed active-learning strategy, from which we construct over 3000 mock dark energy survey lensing data sets. For wCDM inference, for which we allow –1 < w < –$\frac{1}{3}$⁠, our most constraining result uses power spectra combined with map-level (CNN) inference. Using gravitational lensing data only, this map-level combination gives Ω m = 0.283$^{+0.020}_{–0.027}$⁠, S 8 = 0.804$^{+0.025}_{–0.017⁠}$, and w < –0.80 (with a 68 per cent credible interval); compared to the power spectrum inference, this is more than a factor of two improvement in dark energy parameter (Ω⁠ DE , w⁠) precision.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Pressure-enhanced interplay between lattice, spin, and charge in the mixed perovskite La 2 FeMnO 6

Spin crossover plays a central role in the structural instability, net magnetic moment modification, metallization, and even in superconductivity in corresponding materials. Most reports on the pressure-induced spin crossover with a large volume collapse have so far focused on compounds with a single transition metal. Here we report a comprehensive high-pressure investigation of a mixed Fe-Mn perovskite La 2 FeMnO 6 . Under pressure, the strong coupling between Fe and Mn leads to a combined valence/spin transition: Fe 3+ (S=5/2) → Fe 2 + (S = 0) and Mn 3+ (S = 2) → Mn 4+ (S = 3/2), with an isostructural phase transition. The spin transitions of both Fe and Mn are offset by ~20 GPa of the onset pressure, and the lattice collapse occurs in between. Interestingly, Fe 3+ ion shows an abnormal behavior when it reaches a lower valence state (Fe 2+ ) accompanied by a +0.5 eV energy shift in the Fe K -absorption edge at 15 GPa. This process is associated with the charge-spin-orbital state transition from high spin Fe 3+ to low spin Fe 2+ , caused by the significantly enhanced t 2 g -e g crystal field splitting in the compressed lattice under high pressure. Density functional theory calculations confirm the energy preference of the high-pressure state with charge redistribution accompanied by spin state transition of Fe ions. Moreover, La 2 FeMnO 6 maintains semiconductor behaviors even when the pressure reached 144.5 GPa as evidenced by the electrical transport measurements, despite the huge resistivity decreasing seven orders of magnitude compared with that at ambient pressure. The investigation carried out here demonstrates high flexibility of double perovskites and their good potentials for optimizing the functionality of these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1971

A numerical solution of the Navier-Stokes equations for supercritical fluid thermodynamic analysis

An explicit numerical solution of the compressible Navier-Stokes equations is applied to the thermodynamic analysis of supercritical oxygen in the Apollo cryogenic storage system. The wave character is retained in the conservation equations which are written in the basic fluid variables for a two-dimensional Cartesian coordinate system. Control-volume cells are employed to simplify imposition of boundary conditions and to ensure strict observance of local and global conservation principles. Non-linear real-gas thermodynamic properties responsible for the pressure collapse phenomonon in supercritical fluids are represented by tabular and empirical functions relating pressure and temperature to density and internal energy. Wall boundary conditions are adjusted at one cell face to emit a prescribed mass flowrate. Scaling principles are invoked to achieve acceptable computer execution times for very low Mach number convection problems. Detailed simulations of thermal stratification and fluid mixing occurring under low acceleration in the Apollo 12 supercritical oxygen tank are presented which model the pressure decay associated with de-stratification induced by an ordinary vehicle maneuver and heater cycle operation.

Heinmiller, P. J.↗

NASA NTRS1986

Nonlinear compressible convection penetrating into stable layers and producing internal gravity waves

Penetrative convection spanning multiple scale heights is studied within a simple stellar envelope consisting of three layers: a convectively unstable middle layer bounded above and below by stably stratified polytropes. Two-dimensional numerical simulations are used to investigate the fully compressible nonlinear motions that ensue. The cellular flows display prominent downward-directd plumes surrounded by broader regions of upflow. Such asymmetry arises because pressure fluctuations accentuate buoyancy driving in the concentrated plumes and can even lead to weak buoyancy braking in the surrounding ascending flows. As the plumes plunge downward into a region of stable stratification, they serve to excite a broad spectrum of internal gravity waves there. The induced waves are not passive, for they feed back upon the plumes by deflecting them sideways, thereby modulating the amplitude of the convection in time even in the unstable layer. The penetrative motions that billow upward into the upper stable zone are distinctly weaker, and they cascade back downward toward the unstable zone over a broad horizontal scale. The strong excitation of gravity waves by the convection has implications for gradual mixing deep within a star.

Hurlburt, Neal E.↗

DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗