SEARCH · Engineering Papers

Results for “G CODES”

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 685 records · Page 38

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Two-temperature effects in Hall-MHD simulations of the HIT-SI experiment

A two-temperature Hall-magnetohydrodynamic (Hall-MHD) model, which evolves the electron and ion temperatures separately, is implemented in the PSI-Tet 3D MHD code and used to model plasma dynamics in the Helicity Injected Torus–Steady Inductive (HIT-SI) experiment. The two-temperature model is utilized for HIT-SI simulations in both the PSI-Tet and NIMROD codes at a number of different injector frequencies in the 14.5–68.5 kHz range. At all frequencies, the NIMROD two-temperature model results in increased toroidal current, lower chord-averaged density, higher average temperatures, outward radial shift of the current centroid, and axial symmetrization of the current centroid, relative to the single-temperature NIMROD simulations. The two-temperature PSI-Tet model illustrates similar trends, but at high frequency operation, it exhibits lower electron temperature, smaller toroidal current, and decreased axial symmetrization with respect to the single-temperature PSI-Tet model. With all models, average temperatures and toroidal currents increase with the injector frequency. Power balance and heat fluxes to the wall are calculated for the two-temperature PSI-Tet model and illustrate considerable viscous and compressive heating, particularly at high injector frequency. Parameter scans are also presented for artificial diffusivity, wall temperature, and density. Both artificial diffusivity and the density boundary condition significantly modify the plasma density profiles, leading to larger average temperatures, toroidal current, and relative density fluctuations at low densities. As a result, a low density simulation achieves sufficiently high current gain (G > 5) to generate significant volumes of closed flux lasting 1–2 injector periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Design and modeling of indirectly driven magnetized implosions on the NIF

The use of magnetic fields to improve the performance of hohlraum-driven implosions on the National Ignition Facility (NIF) is discussed. The focus is on magnetically insulated inertial confinement fusion, where the primary field effect is to reduce electron-thermal and alpha-particle loss from the compressed hotspot (magnetic pressure is of secondary importance). We summarize the requirements to achieve this state. The design of recent NIF magnetized hohlraum experiments is presented. These are close to earlier shots in the three-shock, high-adiabat (BigFoot) campaign, subject to the constraints that magnetized NIF targets must be fielded at room-temperature, and use ≲1 MJ of laser energy to avoid the risk of optics damage from stimulated Brillouin scattering. We present results from the original magnetized hohlraum platform, as well as a later variant that gives a higher hotspot temperature. In both platforms, imposed fields (at the capsule center) of up to 28 T increase the fusion yield and hotspot temperature. Integrated radiation-magneto-hydrodynamic modeling with the Lasnex code of these shots is shown, where laser power multipliers and a saturation clamp on cross-beam energy transfer are developed to match the time of peak capsule emission and the P2 Legendre moment of the hotspot x-ray image. The resulting fusion yield and ion temperature agree decently with the measured relative effects of the field, although the absolute simulated yields are higher than the data by 2.0−2.7×. The tuned parameters and yield discrepancy are comparable for experiments with and without an imposed field, indicating the model adequately captures the field effects. Self-generated and imposed fields are added sequentially to simulations of one BigFoot NIF shot to understand how they alter target dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

WHONDRS Surface Water Dissolved Organic Carbon and FTICR-MS across Glacial Features in Svalbard 2022

This dataset supports a broader study examining climatic forcing of terrestrial greenhouse gas emission through permafrost and glaciers in Svalbard. The dataset provides geochemistry and organic matter characterization data generated from surface water collected from 11 sites across central Svalbard. Data were also collected in 2021 and can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/1888571. Additional related data were collected and will be published separately in collaboration with Yde and Kleber.This dataset is comprised of one main data folder containing (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) dissolved organic carbon (DOC; measured as non-purgeable organic carbon; NPOC); (5) surface water sampling protocol; (6) readme; (7) methods codes; (8) international geo-sample number (IGSN) mapping file; (9) a folder of field photos; and (10) folder of high resolution characterization of organic matter via 7 Tesla Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) generated through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory). The FTICR folder contains two subfolders, one containing the .xml data files and the other containing instructions for using Formularity (https://omics.pnl.gov/software/formularity) and an R script to process the data based on the user's specific needs. All files are .csv, .pdf, .R, .ref, or .xml.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

A modified cosmic brane proposal for holographic Renyi entropy

We propose a new formula for computing holographic Renyi entropies in the presence of multiple extremal surfaces. Our proposal is based on computing the wave function in the basis of fixed-area states and assuming a diagonal approximation for the Renyi entropy. For Renyi index n ≥ 1, our proposal agrees with the existing cosmic brane proposal for holographic Renyi entropy. For n < 1, however, our proposal predicts a new phase with leading order (in Newton’s constant G) corrections to the cosmic brane proposal, even far from entanglement phase transitions and when bulk quantum corrections are unimportant. Recast in terms of optimization over fixed-area states, the difference between the two proposals can be understood to come from the order of optimization: for n < 1, the cosmic brane proposal is a minimax prescription whereas our proposal is a maximin prescription. We demonstrate the presence of such leading order corrections using illustrative examples. In particular, our proposal reproduces existing results in the literature for the PSSY model and high-energy eigenstates, providing a universal explanation for previously found leading order corrections to the n < 1 Renyi entropies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Conceptual design of ELM control coils for the TCABR tokamak

An upgrade of the Tokamak à Chauffage Alfvén Brésilien (TCABR) is being designed to make it capable of creating a well controlled environment where the impact of resonant magnetic perturbation (RMP) fields on edge localised modes can be addressed over a wide range of (i) plasma shapes, (ii) divertor configurations, (iii) RMP coil geometries and (iv) perturbed magnetic field spectra. To this end, a unique set of in-vessel RMP coils is being designed and, in this work, their conceptual design is presented. This unique set of coils is composed of three toroidal arrays of coils on the low field side and three toroidal arrays of coils on the high field side. Each of these six toroidal arrays is composed of 18 coils thus allowing for the creation of RMP fields with toroidal mode numbers n ≤ 9 and with increased control of the poloidal mode number spectrum. To study dynamical effects of RMP fields of different toroidal mode numbers, all rotating simultaneously with different velocities, each of the 108 RMP coils will be powered independently by power supplies that can provide voltages of up to 4 kV and electric currents of up to 2 kA, with frequencies varying continuously from 0Hz up to 10kHz. A set of physical criteria were used to determine the optimal coil geometry and their respective number of turns to reduce the coil currents and voltages during operation with alternate current. Further, the conceptual design was carried out using both the vacuum approach (no plasma response) and the single-fluid response approach, which accounts for the response of a linear, single-fluid, visco-resistive plasma calculated using the M3D-C 1 code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

DECAF Cross-device characterization of tokamak disruptions indicated by abnormalities in plasma vertical position and current

Abnormal (deviating from target) variations in the plasma vertical position Z and current I p (such as vertical displacements, transient I p 'spikes' and quenches) constitute common elements of a disruption, a phenomenon that is to be mitigated, or ultimately avoided in future reactor-relevant tokamaks. While those abnormalities are generally recognized cross-shot and cross-device, details in terms of appearance (or not) and order of those abnormalities in disruption event chains is bound with the plasma state at the time of the chain initiation. Detection of those abnormalities is thus indicative not only of the onset of the plasma collapse itself, but also of the disruption driving cause that is promoted at the particular plasma state. Here, occurrence of disruptions, explored via detection of a Ip quench, and analysis of disruption event chains constituted by Ip and Z abnormalities, is reported for in total 7 full years of operation of 3 devices (KSTAR, MAST-U and NSTX-U) using the DECAFTM code expanded tools and capabilities. It is shown that the disruption occurrence depends not only on details of the plasma state, but also on (device-dependent) technical elements of the shot exit scenario. A year-to-year change in main disruption causes and a reduction of the disruptivity rate, bound with device and operation upgrades, are reported. Particular trigger instances of disruption event chains (and the full chains, when applicable) are shown to occupy different parts of the operation space diagrams, in accordance with prior expectation. Plasma elongation is identified as an important factor influencing details of the chains and its role will be further explored.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Photosynthetic reaction center variants made via genetic code expansion show Tyr at M210 tunes the initial electron transfer mechanism

Significance Reaction centers (RCs) are critical to photosynthetic energy conversion. RCs in all characterized photosynthetic organisms contain two symmetrically arranged branches of chromophores and enable light-induced electron transfer with high yield. We fine-tune the properties of a key bacterial RC symmetry-breaking tyrosine via its replacement with noncanonical tyrosine analogs and determine kinetic outcomes. Results are interpreted through energetic characterization made possible by resonance Stark spectroscopy. Analysis indicates this tyrosine modulates the mechanism of the initial light-induced electron transfer, affording an alternative functional pathway that maintains the RC’s robust electron transfer. Modern molecular biology, ultrafast spectroscopy, crystallography, and energetic characterization enable the mechanistic model we describe. Our results deepen understanding of RC function and may have implications for other photocatalysts and enzymes.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Radially resolved active charge exchange measurements of the hydrogenic isotope fraction on DIII-D

Radially resolved hydrogenic isotope fraction measurement capabilities have been developed for DIII-D using the main-ion charge exchange recombination (MICER) spectroscopy system in preparation for mixed hydrogen and deuterium experiments. Constraints on the hydrogenic ion temperatures and velocities based on measurements of the impurity ion properties are required to accurately fit the spectrum. Corrections for cross sectional distortions, spatial smearing due to the halo, and a neoclassical offset between the impurity and hydrogenic toroidal rotation are applied to the constraints prior to fitting the MICER spectrum. Extensive atomic physics calculations have been performed using the FIDASIM code, which has recently been improved to allow simulations using mixtures of hydrogenic species. Overall, these results demonstrate that for the same plasma parameters, the Dα emission is 20%–30% brighter than Hα due to differences in rate coefficients associated with the different ion thermal velocities for the same temperature and therefore must be taken into consideration when calculating absolute densities. However, despite these differences, the absolute error when estimating the hydrogen isotope fraction [n H /(n H + n D )] by using the Hα radiance fraction [L Hα /(L Hα + L Dα )] is typically less than 5% due to the way the fraction is formed, making the radiance fraction a reasonably accurate estimate of the isotope fraction for most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Experimental studies of plasma-antenna coupling with the JET Alfvén Eigenmode Active Diagnostic

Here, this paper presents a dedicated study of plasma-antenna (PA) coupling with the Alfvén Eigenmode Active Diagnostic (AEAD) in JET. Stable AEs and their resonant frequencies f, damping rates γ < 0, and toroidal mode numbers n are measured for various PA separations and limiter versus X-point magnetic configurations. Two stable AEs are observed to be resonantly excited at distinct low and high frequencies in limiter plasmas. The values of f and n do not vary with PA separation. However, |γ| increases with PA separation for the low-f, but not high-f, mode, yet this may be due to slightly different edge conditions. The high-f AE is detected throughout the transition from limiter to X-point configuration, though its damping rate increases; the low-f mode, on the other hand, becomes unidentifiable. The linear, resistive MHD code CASTOR is used to simulate the frequency scan of an AEAD-like external antenna. For the limiter pulses, the high-f mode is determined to be an n = 0 GAE, while the low-f mode is likely an n = 2 TAE. During the transition from limiter to X-point configuration, CASTOR indicates that n = 1 and 2 EAEs are excited in the edge gap. These results extend previous experimental studies in JET and Alcator C-Mod; validate the computational work performed by Dvornova et al; and provide guidance for the optimization of PA coupling in upcoming JET energetic particle experiments, for which the AEAD will aim to identify the contribution of alpha particles to AE drive during the DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Mitigating baryonic effects with a theoretical error covariance

ABSTRACT One of the primary sources of uncertainties in modelling the cosmic-shear power spectrum on small scales is the effect of baryonic physics. Accurate cosmology for stage-IV surveys requires knowledge of the matter power spectrum deep in the non-linear regime at the per cent level. Therefore, it is important to develop reliable mitigation techniques to take into account baryonic uncertainties if information from small scales is to be considered in the cosmological analysis. In this work, we develop a new mitigation method for dealing with baryonic physics for the case of the shear angular power spectrum. The method is based on an augmented covariance matrix that incorporates baryonic uncertainties informed by hydrodynamical simulations. We use the results from 13 hydrodynamical simulations and the residual errors arising from a fit to a ΛCDM model using the extended halo model code HMCode to account for baryonic physics. These residual errors are used to model a so-called theoretical error covariance matrix that is added to the original covariance matrix. In order to assess the performance of the method, we use the 2D tomographic shear from four hydrodynamical simulations that have different extremes of baryonic parameters as mock data and run a likelihood analysis comparing the residual bias on Ωm and σ8 of our method and the HMCode for an LSST-like survey. We use different modelling of the theoretical error covariance matrix to test the robustness of the method. We show that it is possible to reduce the bias in the determination of the tested cosmological parameters at the price of a modest decrease in the precision.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Modeling active control of resistive wall mode with power saturation and sensor noise on HL-2M

The resistive wall mode (RWM) control on the HL-2M tokamak is simulated with the MARS-F code (Liu et al 2000 Phys. Plasmas 7 3681), aiming at quantifying control current and voltage requirements when more realistic issues are taken into account, i.e. the control power saturation and the sensor signal noise. The fluid model predicts a narrow stability region for the n = 1 RWM without magnetic feedback, in the 2D parameter space of the plasma pressure versus the toroidal flow speed. Magnetic feedback can fully stabilize the RWM on HL-2M. Without considering the voltage limitation and the sensor signal noise, it is found that plasma flow helps active control of the mode, by reducing the required critical feedback gain for both flux-to-current and flux-to-voltage control schemes. In the absence of the sensor signal noise, the lowest control voltage saturation level, below which the RWM control is lost, is found to roughly satisfy a linear relation to the plasma flow frequency, indicating that subsonic plasma flow is effective in relaxing the control power requirement for the RWM feedback stabilization. The presence of the sensor signal noise substantially modifies the feedback results. A statistical study finds that the sensor signal noise, with the standard deviation of 0.1 G on HL-2M, roughly doubles the required control voltage for successful mode control. The synergistic stabilization effect due to plasma flow is somewhat weakened by the presence of the sensor signal noise. At a given rotation, the tolerable voltage limit generally increases with increasing feedback gain due to the sensor signal noise.

Physics↗

DOE OSTI2025

Charge-exchange losses of beam ions in MAST Upgrade: experiments and modelling

Dedicated experiments were performed on MAST Upgrade to study beam-ion losses caused by charge exchange (CX) with edge neutrals. The fuelling was switched from the high-field side to the low-field side mid-discharge. Direct measurements suggest a strong increase in the neutral density around the plasma and a decrease in the beam-ion density, which is qualitatively explained by CX losses. Measurements by a resistive bolometer have suggested particle bombardment during neutral beam injection, providing a unique opportunity to separate CX from other loss mechanisms. To verify and quantify CX losses, the orbit-following code ASCOT, which accounts both for CX neutralization and reionization, was used to simulate beam-particle power loads on the bolometer. Simulations reproduce measured bolometer power loads during high-field-side fuelling, verifying CX losses of approximately 10% of the off-axis beam power. Toroidally symmetric simulations overestimate power loads on the bolometer during low-field-side fuelling, which is explained by toroidal asymmetry in the neutral density distribution, as is demonstrated by toroidally asymmetric simulations. Results suggest significantly higher CX losses during low-field-side fuelling, up to about 50% of off-axis beam power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

SICAS, a new integrated divertor-SOL-pedestal-core framework for self-consistent modeling of ion and impurity transport in fusion devices

We present a new framework for core-edge integration studies named SICAS (SOLPS-ITER coupled to ASTRA-STRAHL) which enables high fidelity simulations of the core, edge, and divertor regions encompassing the transport of ions as well as the impurities through the entire plasma domain. SICAS handles the exchanging of the particle and power fluxes as well as transport coefficients to ensure consistency through the codes. An overlapping region is defined from the inside of the separatrix of the plasma to the SOLPS-ITER core boundary allowing for matching profiles and fluxes between the two codes for a self-consistent approach. The results presented here demonstrated the flexibility of SICAS to simulate different configurations, scenarios, divertor geometries, and plasma species with good agreement with DIII-D experimental data. This tool opens new possibilities in integrated modeling of fusion devices integrating all relevant phenomena in the core and the divertor plasmas. These capabilities are required for the interpretation of current experiments as well as the design of new devices.

core-edge integration↗

DOE OSTI2020

Machine Learning Characterization of Alfvénic and Sub-Alfvénic Chirping and Correlation With Fast-Ion Loss at NSTX

Abrupt large events in the Alfvenic and sub-Alfvenic frequency bands in tokamaks are typically correlated with increased fast ion loss. Here, machine learning is used to speed up the laborious process of characterizing the behaviour of magnetic perturbations from corresponding frequency spectrograms that are typically identified by humans. Analysis allows for comparison between different mode character (such as quiescent, fixedfrequency, chirping, avalanching) and plasma parameters obtained from the TRANSP code such as the ratio of the neutral beam injection (NBI) velocity and the Alfven velocity (v inj. /v A ), the q-profile, and the ratio of the neutral beam beta and the total plasma beta (β beam,i /β). In agreement with previous work by Fredrickson et al., in this work, we find correlation between β beam,i and mode character. In addition, previously unknown correlations are found between moments of the spectrograms and mode character. Character transition from quiescent to non-quiescent behaviour for magnetic fluctuations in the 50 - 200 kHz frequency band is observed along the boundary v Φ ≲ (1/4) (v inj. - 3v A ) where v Φ is the rotation

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Dynamic RNA acetylation revealed by quantitative cross-evolutionary mapping

N 4 -acetylcytidine (ac 4 C) is an ancient and highly conserved RNA modification that is present on tRNA and rRNA and has recently been investigated in eukaryotic mRNA. However, the distribution, dynamics and functions of cytidine acetylation have yet to be fully elucidated. In this work, we report ac 4 C-seq, a chemical genomic method for the transcriptome-wide quantitative mapping of ac 4 C at single-nucleotide resolution. In human and yeast mRNAs, ac 4 C sites are not detected but can be induced—at a conserved sequence motif—via the ectopic overexpression of eukaryotic acetyltransferase complexes. By contrast, cross-evolutionary profiling revealed unprecedented levels of ac 4 C across hundreds of residues in rRNA, tRNA, non-coding RNA and mRNA from hyperthermophilic archaea. Ac 4 C is markedly induced in response to increases in temperature, and acetyltransferase-deficient archaeal strains exhibit temperature-dependent growth defects. Visualization of wild-type and acetyltransferase-deficient archaeal ribosomes by cryo-electron microscopy provided structural insights into the temperature-dependent distribution of ac 4 C and its potential thermoadaptive role. Our studies quantitatively define the ac 4 C landscape, providing a technical and conceptual foundation for elucidating the role of this modification in biology and disease.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Cross-beam energy transfer in direct-drive ICF. II. Theory and simulation of mitigation through increased laser bandwidth

The response of the cross-beam energy transfer instability (CBET) to laser bandwidth is investigated through a combination of theory and simulation. Existing linear theory is generalized to treat broadband lasers, demonstrating that CBET is most effectively suppressed when the bandwidth exceeds the ion-acoustic wave (IAW) frequency. It is shown that for such bandwidths, reverse (seed to pump) transfer becomes possible, which reduces the net energy transfer rapidly as bandwidth is increased. The CBET gain exponent in this regime scales with bandwidth (Δω) as Δω−3 for Gaussian or Lorentzian laser spectra with different scalings possible for other spectra. Comparison of our theory with linearized fluid and particle-in-cell simulations, performed with the laser-plasma simulation environment (LPSE) and vector particle in cell (VPIC) codes, respectively, finds that the model is accurate in the absence of nonlinear processes. However, linear analysis also finds that the IAW energy density scales as Δω−1, implying that nonlinear effects may be more difficult to control than the CBET scaling would suggest. Indeed, nonlinear effects are found to be present in VPIC simulations with high-intensity lasers, despite minimal apparent CBET. Nonlinear processes in the VPIC cases include particle trapping, the two-ion wave decay, and ion wave self-focusing. In some high intensity VPIC cases, these effects lead to net energy transfer from seed to pump and increases to backscatter stimulated Brillouin scattering reflectivities. Finally, for a given bandwidth, we show that improved control of nonlinear processes can be achieved via smoothing by spectral dispersion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Effects of NBI-induced energetic ions on internal kink stability in HL-2M

A systematic investigation is carried out, studying the effect of the neutral beam injection induced energetic particles (EPs) on the n = 1 (n is the toroidal mode number) internal kink (IK) instability in the HL-2M tokamak, utilizing the MARS-F/K code [Liu et al., Phys. Plasmas 7, 3681 (2000) and 15, 112503 (2008)]. A high-beta sawteething HL-2M scenario, simulated by the TRANSP code [Breslau et al. Computer Software (2018)], is chosen for this study. Compared to the fluid model, non-perturbative magnetohydrodynamic (MHD)-kinetic hybrid computations with MARS-K show a generally stabilization effect on the IK, due to drift kinetic resonances associated with EPs. Further, the bounce resonance of trapped EPs has minor influence on the mode stability. In the absence of the plasma equilibrium flow and with the assumed particle pitch distribution, the transit resonance of co-current (countercurrent) passing EPs destabilizes (stabilizes) the IK. With plasma flow, both co- and countercurrent passing EPs tend to stabilize the mode, but the effect is stronger with the countercurrent particles. These modeling results provide useful guidance for interpreting MHD instabilities in the future high-performance experiments in HL-2M.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

A study of convective core overshooting as a function of stellar mass based on two-dimensional hydrodynamical simulations

ABSTRACT We perform two-dimensional (2D) numerical simulations of core convection for zero-age main-sequence stars covering a mass range from 3 to 20 M⊙. The simulations are performed with the fully compressible time-implicit code music. We study the efficiency of overshooting, which describes the ballistic process of convective flows crossing a convective boundary, as a function of stellar mass and luminosity. We also study the impact of artificially increasing the stellar luminosity for 3 M⊙ models. The simulations cover hundreds to thousands of convective turnover time-scales. Applying the framework of extreme plume events previously developed for convective envelopes, we derive overshooting lengths as a function of stellar masses. We find that the overshooting distance (dov) scales with the stellar luminosity (L) and the convective core radius (rconv). We derive a scaling law $d_{\rm ov} \propto L^{1/3} r_{\rm conv}^{1/2}$, which is implemented in a one-dimensional stellar evolution code and the resulting stellar models are compared to observations. The scaling predicts values for the overshooting distance that significantly increase with stellar mass, in qualitative agreement with observations. Quantitatively, however, the predicted values are underestimated for masses ≳10 M⊙. Our 2D simulations show the formation of a nearly adiabatic layer just above the Schwarzschild boundary of the convective core, as exhibited in recent three-dimensional simulations of convection. The most luminous models show a growth in size with time of the nearly adiabatic layer. This growth seems to slow down as the upper edge of the nearly adiabatic layer gets closer to the maximum overshooting length and as the simulation time exceeds the typical thermal diffusive time-scale in the overshooting layer.

79 ASTRONOMY AND ASTROPHYSICS↗