Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ADIABATIC PROCESSES”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Liquid Piston with Spray Cooling Near-Isothermal Compressor

The goal of this project was to prototype and characterize the performance of a liquid-piston spray-cooled gas compressor. The working principle of the compressor enables optimized high-efficiency operation over a very wide range of operating conditions, unlike conventional compressors that are optimized for a narrow range of operating conditions. The compressor technology is suitable for many applications, such as gas pipeline transport, gas storage, and commercial and residential heat pumps. Both physical testing and computational fluid dynamics (CFD) modeling of the processes using the Oak Ridge National Laboratory high-performance computing center were completed. The experimental and CFD studies focused on a near-isothermal liquid-piston compressor (LPC) that uses propylene glycol to compress CO 2 . The first prototype demonstrated isothermal operation during several sequentially executed cycles of CO 2 compression and raised the temperature of the compressed CO 2 by only 2 K, compared with approximately 6 K when the gas was compressed non-isothermally. Isothermal operation was demonstrated at CO 2 flow rates of up to 2 L/min. The second prototype was designed with two compression chambers to allow continuous flow of high-pressure CO 2 . However, the design of the valve train to direct flow between the compression chambers was not sufficient to allow demonstration of CO 2 compression. Numerical simulations of the LPC in which the compression chamber was filled with propylene glycol injected from the bottom inlet were performed using large eddy simulation (LES) with the wall-adapting local eddy-viscosity subgrid-scale model coupled with the multiphase volume of fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfers within the compression chamber. In this effort, the effects of boundary conditions applied to the LES-VOF calculations (i.e., no wall, an adiabatic wall, and a wall with a heat flux subscribed) on the overall pressure and temperature of the CO 2 gas as well as the transient evolution of flow and heat transfer within the compression chamber were investigated and are discussed in this report. The LES calculation with no wall showed no dynamical flow patterns, and the volume-averaged temperature of CO 2 increased from 305 to 392.7 K, whereas LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO 2 temperatures. The results of the LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO 2 gas layers and circulation flow patterns within the chamber that contributed to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor. Though an industrial partnership for commercializing the compressor was not secured, the technology was attractive for an industrial partner to use in two research proposals in response to US Department of Energy funding opportunity announcements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

One-to-one mapping between the electromagnetic modes of cylindrical and coaxial half-wave cavities

Design of radio frequency (RF) couplers and diagnostics require a good understanding of the electromagnetic mode patterns of RF cavities. This study investigates the adiabatic transformation of transverse magnetic (TM) modes in a cylindrical cavity into transverse electromagnetic (TEM) modes of a coaxial cavity by gradually introducing an inner conductor. Using CST Studio Suite, we simulate the eigenmode evolution as the geometry transforms from a pure cylindrical to a coaxial configuration. We track the behavior of TM010 through TM014 modes to observe the continuous evolution into the corresponding TEM0 through TEM4 modes of the coaxial cavity. The process is governed by the evolution of the electric field orientation as the geometry shifts, enabling the axial TM fields to reorient into the radial electric field configuration of TEM modes. Field patterns, eigen-frequencies, and mode indentities are analyzed throughtout the transition. The results provide simulation-based evidence that TM to TEM conversion occurs without generation of newer eigenmodes, offering a valuable insight into the design of transition regions in superconducting RF (SRF) systems and provides a foundation for experimental validation.

Accelerator Physics↗
DOE OSTI2018

Theoretical Analysis of the Entrainment–Mixing Process at Cloud Boundaries. Part I: Droplet Size Distributions and Humidity within the Interface Zone

Abstract The problem of a complex entrainment–mixing process is analyzed by solving a diffusion–evaporation equation for an open region in the vicinity of the cloud–dry air interface. Upon normalization the problem is reduced to a one-parametric one, the governing parameter being the potential evaporation parameter R proportional to the ratio of saturation deficit in the dry air to the available liquid water content in the cloud air. As distinct from previous multiple studies analyzing mixing within closed adiabatic volumes, we consider a principally nonstationary problem that never leads to a homogeneous equilibrium state. It is shown that at R < −1 the cloud edge shifts toward the cloud; that is, the cloud dissipates due to mixing with dry air, and the cloud volume decreases. If R > −1, the cloud edge shifts outside; that is, the mixing leads to an increase in the cloud volume. The time evolution of droplet size distribution and its moments, as well as the relative humidity within the expanding cloud–dry air interface, are calculated and analyzed. It is shown that the values of the mean volume radii rapidly decrease within the interface zone in the direction away from the cloud, indicating significant changes in the cloud edge microstructure. Scattering diagrams plotted for the cloud edge agree well with high-frequency in situ measurements, corroborating the reliability of the proposed approach. It is shown that the humidity front moves toward dry air faster than the front of liquid water content. As a result, the mixing leads to formation of a humid air shell around the cloud. The widths of the interface zone and humid shell are evaluated.

Pinsky, Mark↗
NASA NTRS1986

Possible explanation of the solar-neutrino puzzle

A new derivation of the Mikheyev and Smirnov (1985) mechanism for the conversion of electron neutrinos into mu neutrinos when traversing the sun is presented, and various hypotheses set forth. It is assumed that this process is responsible for the detection of fewer solar neutrinos than expected, with neutrinos below a minimum energy, E(m), being undetectable. E(m) is found to be about 6 MeV, and the difference of the squares of the respective neutrino masses is calculated to be 6 X 10 to the - 5th sq eV. A restriction on the neutrino mixing angle is assumed such that the change of density near the crossing point is adiabatic. It is predicted that no resonance conversion of neutrinos will occur in the dense core of supernovae, but conversion of electron neutrinos to mu neutrinos will occur as they escape outward through a density region around 100.

Bethe, H. A.↗
NASA NTRS2003

Parameterization of Cloud Droplet Formation in Global Climate Models

An aerosol activation parameterization has been developed based on a generalized representation of aerosol size and composition within the framework of an ascending adiabatic parcel; this allows for parameterizing the activation of chemically complex aerosol with an arbitrary size distribution and mixing state. The new parameterization introduces the concept of"population splitting", in which the cloud condensation nuclei (CCN) that form droplets are treated as two separate populations; those that have a size close to their critical diameter and those that do not.Explicit consideration of kinetic limitations of droplet growth is introduced. Our treatment of the activation process unravels much of its complexity. As a result of this, a substantial number of conditions of droplet formation can be treated completely free of empirical information or correlations; there are, however, some conditions of droplet activation for which an empirically derived correlation is utilized. Predictions of the parameterization are compared against extensive cloud parcel model simu;lations for a variety of aerosol activation conditions that cover a wide range of chemical variability and CCN concentrations. The parameterization tracks the parcel model simulations closely and robustly. The parameterization presented here is intended to allow for a comprehensive assessment of the aerosol indirect effect in general circulation models.

Nenes, A.↗
DOE OSTI2016

Theoretical simulations on the antioxidant mechanism of naturally occurring flavonoid: A DFT approach

Synthetic antioxidants such as butylated hydroxyanisole (BHA) and butylated hydroxytoluene (BHT) are found to be toxic, hence non-carcinogenic naturally occurring radical scavengers especially flavonoids have gained considerable importance in the past two decades. In the present investigation, the radical scavenging activity of C-glycosyl flavonoids is evaluated using theoretical approach which could broaden its scope in therapeutic applications. Gas and solvent phase studies of structural and molecular characteristics of C-glycosyl flavonoid, isovitexin is investigated through hydrogen atom transfer mechanism (HAT), Electron transfer-proton transfer (ET–PT) and Sequential proton loss electron transfer (SPLET) by Density functional theory (DFT) using hybrid parameters. The computed values of the adiabatic ionization potential, electron affinity, hardness, softness, electronegativity and electrophilic index indicate that isovitexin possess good radical scavenging activity. The behavior of different –OH groups in polyphenolic compounds is assessed by considering electronic effects of the neighbouring groups and the overall geometry of molecule which in turn helps in analyzing the antioxidant capacity of the polyphenolic molecule. The studies indicate that the H–atom abstraction from 4’–OH site is preferred during the radical scavenging process. From Mulliken spin density analysis and FMOs, B–ring is found to be more delocalized center and capable of electron donation. Comparison of antioxidant activity of vitexin and isovitexin leads to the conclusion that isovitexin acts as a better radical scavenger. This is an evidence for the importance of position of glucose unit in the flavonoid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

The steady growth of the high-energy spectral cut-off in relativistic magnetic reconnection

Magnetic reconnection is invoked as an efficient particle accelerator in a variety of astrophysical sources of non-thermal high-energy radiation. With large-scale two-dimensional particle-in-cell simulations of relativistic reconnection (i.e. with magnetization σ >> 1) in pair plasmas, we study the long-term evolution of the power-law slope and high-energy cut-off of the spectrum of accelerated particles. We find that the high-energy spectral cut-off does not saturate at γcut ~ 4σ, as claimed by earlier studies, but it steadily grows with time as long as the reconnection process stays active. At late times, the cut-off scales approximately as γ cut ∝√t, regardless of the flow magnetization and initial temperature. We show that the particles dominating the high-energy spectral cut-off reside in plasmoids, and in particular in a strongly magnetized ring around the plasmoid core. The growth of their energy is driven by the increase in the local field strength, coupled with the conservation of the first adiabatic invariant. We also find that the power-law slope of the spectrum (p = -d log N/d log γ) evolves with time. For σ ≳ 10, the spectrum is hard at early times (p ≲ 2), but it tends to asymptote to p ~ 2; the steepening of the power-law slope allows the spectral cut-off to extend to higher and higher energies, without violating the fixed energy budget of the system. Our results demonstrate that relativistic reconnection is a viable candidate for accelerating the high-energy particles emitting in relativistic astrophysical sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Theoretical study of Raman chirped adiabatic passage by X-ray absorption spectroscopy: Highly excited electronic states and rotational effects

Raman Chirped Adiabatic Passage (RCAP) is an efficient method to climb the vibrational ladder of molecules. It was shown on the example of fixed-in-space HCl molecule that selective vibrational excitation can thus be achieved by RCAP and that population transfer can be followed by X-ray Photoelectron spectroscopy [S. Engin, N. Sisourat, P. Selles, R. Taïeb, and S. Carniato, Chem. Phys. Lett. 535, 192–195 (2012)]. Here, in a more detailed analysis of the process, we investigate the effects of highly excited electronic states and of molecular rotation on the efficiency of RCAP. Furthermore, we propose an alternative spectroscopic way to monitor the transfer by means of X-ray absorption spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2005

Soot Formation in Purely-Curved Premixed Flames and Laminar Flame Speeds of Soot-Forming Flames

The research addressed here is a collaborative project between University of Delaware and Case Western Reserve University. There are two basic and related scientific objectives. First, we wish to demonstrate the suitability of spherical/cylindrical, laminar, premixed flames in the fundamental study of the chemical and physical processes of soot formation. Our reasoning is that the flame standoff distance in spherical/cylindrical flames under microgravity can be substantially larger than that in a flat burner-stabilized flame. Therefore the spherical/cylindrical flame is expected to give better spatial resolution to probe the soot inception and growth chemistry than flat flames. Second, we wish to examine the feasibility of determining the laminar flame speed of soot forming flames. Our basic assumption is that under the adiabatic condition (in the absence of conductive heat loss), the amount and dynamics of soot formed in the flame is unique for a given fuel/air mixture. The laminar flame speed can be rigorously defined as long as the radiative heat loss can be determined. This laminar flame speed characterizes the flame soot formation and dynamics in addition to the heat release rate. The research involves two integral parts: experiments of spherical and cylindrical sooting flames in microgravity (CWRU), and the computational counterpart (UD) that aims to simulate sooting laminar flames, and the sooting limits of near adiabatic flames. The computations work is described in this report, followed by a summary of the accomplishments achieved to date. Details of the microgra+ experiments will be discussed in a separate, final report prepared by the co-PI, Professor C-J. Sung of CWRU. Here only a brief discussion of these experiments will be given.

Buchanan, Thomas↗
DOE OSTI2020

On the elastocaloric effect in CuAlBe shape memory alloys: A quantitative phase-field modeling approach

We report the reversible stress-induced phase transformation in shape memory alloys (SMAs) is a dissipative process during which heat is absorbed or released. The inherent temperature variations inside the material has an elastocaloric effect (eCE) with appealing applications in solid-state cooling technology such as compact and efficient on-board refrigeration system for eletronic devices. In this manuscript, we conduct the first study of eCE of CuAlBe SMAs utilizing phase-field modeling. For an applied stress of 500 MPa, the results for polycrystalline Cu-Al11-2Be (at. %) show a minimum adiabatic unloading temperature change of -10 K over a pseudoelastic window of 40 K. In the absence of plastic deformation, the material demonstrates good reproducibility of the eCE over a few loading–unloading cycles. The presence of plastic deformation is found to cause functional fatigue that deteriorates the cooling capacity; however, the coefficient of performance only decreases from 9.04 to 8.03, which is still a very good value. These results place CuAlBe as a frontrunner SMA for solid-state cooling compared to the expensive NiTi.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Enhanced strength and toughness in ultra-fine grained 99.9% copper obtained by cryo-hydrostatic extrusion

Highlights: • Cryo-HE can increase the yield strength of Cu by more than threefold. • In cryo-HE of Cu, the grain refinement is stronger by 30% than in room-HE. • Cryo-HE with post-annealing effectively creates a bimodal microstructure in Cu. • Short-term post-annealing applied to cryo-HE-treated Cu can double its ductility. • Grain refinement with recrystallization results in a slight increase of Cu toughness. - Abstract: Ultrafine grained 99.9% Cu with enhanced strength and toughness was produced by hydrostatic extrusion conducted at the temperature of liquid nitrogen. The cryo-extruded Cu has a tensile yield strength of 490 MPa, hardness of 145 HV0.2, and tensile fracture toughness of 39.5 MPa m{sup 1/2}, which means that their values increased by 380%, 185% and ~10%, respectively, compared to those of the starting material. The significant strengthening is due to the strong grain refinement to 320 nm and the high energy of defects stored during the processing at liquid nitrogen temperature. The short term post deformation annealing doubles the ductility of material even by above 10% of the elongation to fracture at a slight expense of the yield strength (485 MPa). The observed behavior of Cu can be explained in terms of its low stacking fault energy, the susceptibility to dynamic recrystallization within the working zone of the extrusion die, and the triaxial strain components with the very large strain gradients, combined with the significant adiabatic heating effect active during the hydrostatic extrusion. The better homogeneity of the grain refinement and the higher efficiency of the distribution of defects and of the energy storage accompanied by the decreased dislocation mobility during the cryo hydrostatic extrusion contributed to the higher strength.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Molecular events in the light of strong fields: A light‐induced potential scenario

A new perspective on how to manipulate molecules by means of very strong laser pulses is emerging with insights from the so‐called light‐induced potentials, which are the adiabatic potential energy surfaces of molecules severely distorted by the effect of the strong field. Different effects appear depending on how the laser frequency is tuned, to a certain electronic transition, creating light‐induced avoided crossings, or very off‐resonant, generating Stark shifts. In the former case it is possible to induce dramatic changes in the geometry and redistribution of charges in the molecule while the lasers are acting and to fully control photodissociation reactions as well as other photochemical processes. Several theoretical proposals taken from the work of the authors are reviewed and analyzed showing the unique features that the strong‐laser chemistry opens to control the transient properties and the dynamics of molecules. © 2016 Wiley Periodicals, Inc.

Chang, Bo Y.↗
DOE OSTI2024

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Pressure balance in a lower collisionality, attached tokamak scrape-off layer

Previously the gyrokinetic neoclassical code XGCa found that the pressure balance (or momentum balance) in the diverted scrape-off layer (SOL) does not follow that from the fluid description based on the Chew–Goldberger–Low (CGL) theory. In this paper a gyrocenter parallel momentum moment equation is derived for pressure balance (or momentum balance) in a tokamak SOL, for use in interpreting this difference. The new gyrocenter-fluid pressure balance equation allows identifying from the XGCa code results which physical processes dominate the setting of pressure variation in the scrape-off layer. This pressure balance equation is applied to the simulation of a DIII-D H-mode discharge, with a lower ion collisionality in the SOL, and the Coulomb and atomic collisions are not strong enough to yield a detached divertor plasma. It is found that the total pressure balance is much better matched using the gyrocenter parallel momentum moment equation. Electrons are shown to be dominantly adiabatic, while ions have multiple contributions to pressure balance, including terms originating from particle drifts. Furthermore, these results show that in strong gradient, low collisionality regions of the SOL, the typical fluid reductions miss important effects captured in the gyrokinetic model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Low Temperature Formation of Carbonaceous Grains from Hydrocarbons and PAHs with the COSmIC Facility

Complex carbon molecules and ions are ubiquitous in space and form the building blocks of the carbonaceous components of cosmic dust grains, ultimately contributing to the formation of planets. We report experimental investigations of the low temperature chemical pathways leading to the formation of cosmic grain analogs from gas phase molecular precursors to better understand the evolution of cosmic carbon. The study of the formation of dust is essential to understand and to quantify the budget of extraterrestrial organic molecules. Although dust plays an important role in the evolution of interstellar chemistry and in the formation of organic molecules, little is known on the formation and destruction processes of carbonaceous dust. The experiments were performed using the Cosmic Simulation Chamber (COSmIC) to generate and characterize solid grains formed from gas phase precursors under controlled conditions representative of astrophysical environments [1]. Using COSmIC, it is possible to investigate the evolution of cosmic carbon from the formation of neutral and ionized hydrocarbons and PAHs in the gas phase [2], to solid carbon grains [3, 4, 5]. This is achieved by using a pulsed slit discharge nozzle to produce an adiabatic jet expansion and cool down Ar-hydrocarbon and PAH gas mixtures to astrophysically relevant temperature before inducing chemistry by generating a plasma discharge in the stream of the expansion. This plasma induced chemistry results in the formation of complex molecules and solid particles, analogs of cosmic grains in situ in the plasma expansion (i.e., without wall effects) that are carried by the accelerated gas in the expansion, and collected on substrates placed a few centimeters downstream of the electrodes. The results of a preliminary solid phase ex-situ analysis of cosmic grain analogs produced from a series of gas mixtures including hydrocarbons such as CH4 and C2H2 [3] and PAHs [4, 5] found in circumstellar ejecta of late C stars will be discussed. SEM imaging was used to provide insight on the impact of the precursors on the morphology and growth structure of the grains produced [3]. Laser desorption mass spectrometry was used to identify the molecules making up the main structures within the condensed grains [4] while X-ray irradiation and X-ray photoemission spectroscopy were used to investigate the impact the impact of stellar X-rays on cosmic dust.

Carbonaceous Grains↗
DOE OSTI2014

Anomalous grain refinement trends during mechanical milling of Bi 2 Te 3

The structural evolution of nanocrystalline bismuth telluride (Bi 2 Te 3 ) during mechanical milling is investigated here under different milling energies and temperatures. After prolonged milling, the compound evolves toward a steady-state nanostructure that is found to be unusually strongly dependent on the processing conditions. In contrast to most literature on mechanical milling, in Bi 2 Te 3 we find that the smallest steady-state grain sizes are attained under the lowest energy milling conditions. An analysis based on the balance between refinement and recovery in the steady state shows that two regimes of behavior are expected based on the thermo-physical properties of the milled powder. Bi 2 Te 3 lies in a relatively unusual regime where greater impact energy promotes adiabatic heating and recovery more than it does defect accumulation; hence more intense milling leads to larger steady-state grain sizes. Implications for other materials are discussed with reference to a “milling intensity map” that delineates the set of material properties for which this behavior will be observed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations

There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.

79 ASTRONOMY AND ASTROPHYSICS↗