SEARCH · Engineering Papers

Results for “Radiative properties”

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 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Tethered balloon-borne profile measurements of atmospheric properties in the cloudy atmospheric boundary layer over the Arctic sea ice during MOSAiC: Overview and first results

The tethered balloon-borne measurement system BELUGA (Balloon-bornE moduLar Utility for profilinG the lower Atmosphere) was deployed over the Arctic sea ice for 4 weeks in summer 2020 as part of the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition. Using BELUGA, vertical profiles of dynamic, thermodynamic, aerosol particle, cloud, radiation, and turbulence properties were measured from the ground up to a height of 1,500 m. BELUGA was operated during an anomalously warm period with frequent liquid water clouds and variable sea ice conditions. Three case studies of liquid water phase, single-layer clouds observed on 3 days (July 13, 23, and 24, 2020) are discussed to show the potential of the collected data set to comprehensively investigate cloud properties determining cloud evolution in the inner Arctic over sea ice. Simulated back-trajectories show that the observed clouds have evolved within 3 different air masses (“aged Arctic,” “advected over sea ice,” and “advected over open ocean”), which left distinct fingerprints in the cloud properties. Strong cloud top radiative cooling rates agree with simulated results of previous studies. The weak warming at cloud base is mostly driven by the vertical temperature profile between the surface and cloud base. In-cloud turbulence induced by the cloud top cooling was similar in strength compared to former studies. From the extent of the mixing layer, it is speculated that the overall cloud cooling is stronger and thus faster in the warm oceanic air mass. Larger aerosol particle number concentrations and larger sizes were observed in the air mass advected over the sea ice and in the air mass advected over the open ocean.

Arctic↗

DOE OSTI2020

Long-range transported North American wildfire aerosols observed in marine boundary layer of eastern North Atlantic

Wildfire is a major source of biomass burning aerosols, which greatly impact Earth climate. Tree species in North America (NA) boreal forests can support high-intensity crown fires, resulting in elevated injection height and longer lifetime (on the order of months) of the wildfire aerosols. Given the long lifetime, the properties of aged NA wildfire aerosols are required to understand and quantify their effects on radiation and climate. Here we present comprehensive characterization of climatically relevant properties, including optical properties and cloud condensation nuclei (CCN) activities of aged NA wildfire aerosols, emitted from the record-breaking Canadian wildfires in August 2017. Despite the extreme injection height of ~12 km, some of the wildfire plumes descended into the marine boundary layer in the eastern North Atlantic over a period of ~2 weeks, owing to the dry intrusions behind mid-latitude cyclones. The aged wildfire aerosols have high single scattering albedos at 529 nm (ω 529 ; 0.92–0.95) while low absorption Ångström exponents (Å abs ) at 464 nm/648 nm (0.7–0.9). In comparison, Å abs of fresh/slightly aged ones are typically 1.4–3.5. This low Å abs 529 and low Å abs ~110 nm diameter) with thick non-absorbing coatings. The accelerated descent of the wildfire plume also led to strong increase of CCN concentration at the supersaturation levels representative of marine low clouds. The hygroscopicity parameter, κ CCN , of the aged wildfire aerosols varies from 0.2 to 0.4, substantially lower than that of background marine boundary layer aerosols. However, the high fraction of particles with large diameter (i.e., within accumulation size ranges, ~100–250 nm) compensates for the low values of κ, and as a result, the aged NA wildfire aerosols contribute more efficiently to CCN population. These results provide direct evidence that the long-range transported NA wildfires can strongly influence CCN concentration in remote marine boundary layer, therefore the radiative properties of marine low clouds. Given the expected increases of NA wildfire intensity and frequency and regular occurrence of dry intrusion following mid-latitude cyclones, the influence of NA wildfire aerosols on CCN and clouds in remote marine environment need to be further examined.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Uncertainty of SW Cloud Radiative Effect in Atmospheric Models Due to the Parameterization of Liquid Cloud Optical Properties

Clouds are largely responsible for the spread of climate models predictions. Here we focus on the uncertainties in cloud shortwave radiative effect due to the parameterization of liquid cloud single scattering properties (SSPs) from liquid water content (LWC) and droplet number concentration (N), named parameterization of cloud optical properties. Uncertainties arise from not accounting for the droplet size distribution (DSD)—which affects the estimation of the effective radius (r eff ) and modulates the r eff -dependency of the SSPs—and from averaging SSPs over wide spectral bands. To assess these uncertainties a series of r eff -dependent SSPs parameterizations corresponding to various DSDs and spectral averaging methods are derived and implemented in a radiative code. Combined with the DSD-dependent estimation of r eff they are used to compute the bulk radiative properties (reflectance, transmittance, absorptance) of various clouds (defined in terms of LWC and N), including a homogeneous cloud, more realistic case studies, and outputs of a climate model. The results show that the cloud radiative forcing can vary up to 20% depending on the assumed DSD. Likewise, differences up to 20% are obtained for heating rates. The estimation of r eff is the main source of uncertainty, while the SSPs parameterization contributes to around 20% of the total uncertainty. Spectral averaging is less an issue, except for atmospheric absorption. Overall, global shortwave cloud radiative effect can vary by 6 W m –2 depending on the assumed DSD shape, which is about 13% of the best observational estimate.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Vibration Studies for the CLARA Interferometer

The purpose of the CLARA experiment is to study the nature of undulator radiation emitted by single electrons circulating in the IOTA storage ring. The classical and quantum properties of the radiation are investigated by measuring its coherence length, intensity fluctuations and time correlations. The experiment took data in IOTA Run 4 (2022-2023). The key component of the apparatus is a Mach-Zehnder interferometer (MZI), in which the optical length of one of the arms can be precisely controlled. For some measurements, the fine regulation of the arm length must be smaller than the radiation wavelength. For this reason, the apparatus is particularly sensitive to mechanical vibrations. In this note, we model and measure the effect of vibrations on the performance of the MZI under various conditions. Several improvements of the setup were implemented to minimize systematic distortions of the observed interference patterns. In this report, we present a mathematical model of the effect of mechanical vibrations on observed detector and coincidence rates. We also describe the measurements that were made to estimate the magnitude and spectra of rate fluctuations and their sources in the CLARA MZI under various conditions. Finally, we estimate the magnitude of arm length fluctuations and we deduce the sensitivity to coincidence-rate variations in our apparatus.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Reconstruction of beam parameters and betatron radiation spectra measured with a Compton spectrometer

The photon flux resulting from high-energy electron beam interactions with high-field systems, such as those found in the upcoming FACET-II experiments at the SLAC National Accelerator Laboratory, yields deep insight into the electron beam’s underlying dynamics during the interaction. However, extracting this information is an intricate process. To demonstrate how to approach this challenge using modern methods, this paper utilizes simulated data that models plasma wakefield acceleration-derived betatron radiation in experiments to determine reliable methods of reconstructing key beam and beam-plasma interaction properties. For betatron radiation measurements, translating the observed 200⁢ keV to 30⁢ MeV photon double-differential energy-angle spectra obtained from an advanced Compton spectrometer requires testing multiple methods to optimize the pipeline from its response to incident electron beam information. The paper compares maximum likelihood estimation and machine learning to refine the translation of photon spectra into precise electron beam metrics, such as spot size, energy, and emittance, enhancing the understanding of beam behavior within these dense, high-field environments. We also introduce machine learning and the expected maximization algorithm to reconstruct the primary photon spectrum, employing a multilayer neural network for regression analysis of the energy and angle spectra. With appropriate modifications, the advanced methods reproduce relevant incident beam parameters with high accuracy, even for beam sizes in the <10 μ⁢m range. This capacity is critical to understanding intense beam propagation and its optimization in plasma.

Beam code development & simulation techniques↗

DOE OSTI2022

Radiation hard gallium oxide scintillators for high count rate radiation detection

Despite decades of materials research, the availability of appropriate crystalline radiation-hard scintillators with ultrafast decay times (~10ns) and high light yields (>2000 ph/MeV) are still limited. In this study, we demonstrate the potential of gallium oxide (β-Ga 2 O 3 ) based scintillators for high count rate applications. The low-cost rugged β-Ga 2 O 3 scintillators were grown using the optical float zone (FZ) technique. Several dopants and growth atmospheres were used to demonstrate the balance between ultrafast primary decay time and the light yield of the scintillators. Light yields as high as 6446 ± 716 ph/MeV were obtained with 18.3 ns decay time for 662 keV gamma rays. Decay times as fast as 12ns were obtained with a 3212 ± 357 ph/MeV light yield. An excellent 662 keV gamma energy resolution of 7% was obtained using the β-Ga 2 O 3 :Ce crystals. Even with these excellent scintillation properties, the gamma radiation hardness of β-Ga 2 O 3 crystals was better than the leading radiation-hard lead tungstate (PbWO 4 ) crystals. Twenty-six scintillators fabricated from the grown β-Ga 2 O 3 crystals were tested for radiation hardness with a 60 Co gamma irradiation rate of 14 krad/h with a total radiation dose of 2 Mrad. The scintillation properties of these crystals remained unaffected. To date, no other scintillator has demonstrated such high radiation resistance. Furthermore, these FZ-grown low-cost, non-hygroscopic, and high-density β-Ga 2 O 3 scintillators can replace the existing scintillators in many applications ranging from nuclear and high energy physics experiments to nuclear security.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Beryllium (Be) Handbook

Beryllium (Be), atomic number 4, is a silver gray metal of low density (1.85g/cm 3 ), moderately high melting point (1289°C), and quite good stability in the atmosphere. Favorable mechanical properties, particularly specific stiffness (elastic modulus/density), have resulted in a number of weight-critical structural applications. Other notable properties of beryllium include good dimensional stability, high specific heat and extremely high transparency to X-radiation. The nuclear properties of beryllium include a high neutron scattering cross section (6 barns for thermal neutrons) and a low neutron absorption cross section (0.009 barns for thermal neutrons) making it a good choice for a neutron reflector. The (n, 2n) reaction in beryllium also makes it suitable for a neutron multiplier.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Evaporation due to infrared heating and natural convection

The results obtained from a drying experiment under infrared heating are compared to those obtained from correlations for the mass transfer coefficients, which are commonly used to obtain evaporation rates for drying applications. It was shown that the mass transfer coefficients evaluated by natural convection correlations are approximately ten times lower than those measured under infrared (IR) exposure. In this work, an attempt to explain this discrepancy, which has been also observed in other studies and explained on the basis of thermal radiation effects, optical properties of NMP are presented and the implications of the IR and near infrared (NIR) spectrum on the thermal radiation effects and ensuing evaporation rates are discussed. The results show that expected raise in the film temperature, which would be due to the thermal radiation effects, cannot be used to describe the measured evaporation rates for this solvent. Therefore, this study provides important evaporation mass flux data, which cannot be predicted using traditional correlations, for the design of IR assisted drying systems used in membrane fabrication.

42 ENGINEERING↗

DOE OSTI2021

Defect chemistry of Cr-B binary and Cr-Al-B MAB phases: Effects of covalently bonded B networks

Transition metal borides, which are three-dimensional (3D) layered materials containing covalently bonded B networks, have shown a number of excellent properties, such as radiation resistance and the ability to act as a diffusion barrier in integrated circuits. However, defect behavior, which controls many of the materials’ properties, has remained unknown in these materials. In this work, we investigate the effects of the B networks on the defect chemistry in both binary borides (CrB, Cr 3 B 4 , Cr 2 B 3 ) and ternary MAB phases (Cr 2 AlB 2 , Cr 3 AlB 4 , Cr 4 AlB 6 ) using first-principles calculations. We find that increasing the number of B rings in the structure leads to lower formation energies and higher concentrations of Frenkel pairs. The results can be explained by the fact that the strongest Cr-B bond is weakened when borides have more B rings, leading to a reduction in the formation energy of Cr and B vacancies. Also, the bonds associated with Cr atoms bonded within B rings are softer in structures containing more B rings, which allows Cr interstitials to form with a lower energy cost and contributes to an increase in the concentration of Cr interstitials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Mixed Material Scintillator Systems Particle ID Modelling Report

AMixed-Material Scintillator System (MMSS) is a radiation detector using a scintillator made of a heterogeneous, structured mixture of two or more materials. MMSSs are designed so that the structured mixing of materials encodes properties of the radiation detected in the scintillator. This allows for new radiation detectors with advantages over detectors using traditional homogeneous scintillators. This report contains results from this project’s second task, to model and quantify the promise of the MMSSs invented earlier in the project. Those inventions have been divided into two categories, and this report focuses on the first category: MMSS inventions the use zones of differently-colored scintillators to detect and characterize neutron sources. We call this category particle-ID, or PID, detectors. Inventions using gradients designed to encode the position of radiation interactions using gradients will be discussed in a future report.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Mixed Material Scintillator Systems Position Reconstruction Modelling Report

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Improving Predictive Capability in REHEDS Simulations with Fast, Accurate, and Consistent Non-Equilibrium Material Properties

Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Role of chemical disorder on radiation-induced defect production and damage evolution in NiFeCoCr

Understanding chemical disorder in many concentrated solid solution alloys (CSAs) at the levels of electrons and atoms has attracted increasing attention as a path forward to reveal and identify underlying mechanisms for extraordinary mechanical properties and improved radiation tolerance. Single-phase NiFeCoCr CSA is a common base for many high-entropy alloys (HEAs) that have shown improved mechanical strength and radiation tolerance. In this study, defect production and damage evolution in NiFeCoCr under ion irradiation at room temperature to dose over 20 dpa are determined using ion channeling technique along both <100> and <110> directions utilizing multiple probing beam energies. The results obtained from the multi-axial and multi-energy channeling analysis are compared with those previously obtained for Ni crystals irradiated under similar conditions. The influence of chemical complexity on defect production and clustering at early-stage under room temperature irradiation up to dose of 1 dpa is discussed based on positron annihilation spectroscopy results. Defect structure evaluation in Ni and NiFeCoCr is also discussed based on transmission electron microscopy results over a prolonged irradiation at both room and elevated temperatures. Compared with chemically complex NiFeCoCr, larger dislocation loops thus less lattice strain are expected to form in pure Ni. Moreover, the role of chemical disorder in this CSA is also investigated based on ab initio calculations using large supercells. Finally, to understand the impact of chemical complexity effect on defect structure evolution, this integrated research effort attempts to link the relatively large charge redistribution due to difference in valence electron counts resulting from alloying different 3d transition metal elements, moderate lattice distortion arising from similar adaptable atomic size, and notable suppressed or delayed damage evolution in NiFeCoCr.

36 MATERIALS SCIENCE↗

DOE OSTI2021

3D radiation, density, and MHD structures following neon shattered pellet injection into stable DIII-D Super H-mode discharges

Six nominally repeat neon shattered pellet injection (SPI) shutdowns of stable DIII-D Super H-modes are studied to understand the 3D properties of the radiation and impurity transport. The radiation efficiency and radiation peaking determine whether first wall melting is expected following disruption mitigation in ITER. Previous studies make use of axisymmetric approximations to infer radiation efficiencies, but validating the high efficiency required by ITER necessitates improved accuracy, and this work contributes by exploring the 3D radiation and density structures that will inform forward modeling. When the neon shatter plume produced by the SPI reaches the plasma edge, m/n = 3/1 and 2/1 island O-points are observed to align with the injection trajectory in five out of six cases, suggesting that the injected material seeds the island O-points. Field aligned neon structures emitting Ne-I line radiation drift at 1 km/s in the ion diamagnetic drift direction during the pre-TQ, tracking the motion of the m/n = 2/1 island O-point. Neon fragments penetrate to the q = 2 surface by the time of the TQ. Techniques to constrain the 3D emissivity are explored, and one method constrains a 3D flux tube that is consistent with the radiation data, and when mapped to the interferometers, intersects the lasers that measure the highest density. The resulting structure derived from the radiation measurements exists near the 2/1 island X-point. In five repeatable discharges, the peak of the radiation in the toroidal direction exists in a 120° toroidal sector where the injection occurs, in contrast with the outlier discharge where the toroidal peak exists in the complementary 240° toroidal sector far from the injector, and where a 50% lower density rise is observed. In conclusion, the n = 1 phase behavior is markedly different in the outlier discharge, suggesting a possible dependence of the radiation structure and the assimilation efficiency on MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

A machine learning framework for accurate and robust analysis of radiation detector pulses

The microscopic properties of atomic nuclei are used to study various scientific questions. They are essential for understanding the fundamental forces of nature and the chemical evolution of the universe. Detecting decay radiation from radioactive nuclei makes it possible to probe these fundamental nuclear properties. Detector waveform traces may contain additional information about the radiation. Generally, advanced signal processing techniques are needed to extract this additional information, often involving fitting the waveform with model response functions using non-linear least-squares optimization with second-order gradient methods. While this is a powerful technique, it is also computationally expensive, leading to slow processing time, which scales with the volume of data. To address this problem, we have developed a machine learning (ML) approach that infers the characteristics of traces from a model detector response function. In particular, we are interested in classifying whether a single recorded trace consists of one or two pulse constituents and estimating the pulse parameters. Furthermore, our proposed ML method can precisely extract the pulses’ parameters, such as energy and timing information, and accurately classify the pulse multiplicity of a trace. Unlike non-learning-based approaches, our ML approach uses neural networks that are significantly faster at inference, as they do not require any optimization during this stage.

Curve fitting↗

DOE OSTI2020

Radiative heat transfer in FLiBe molten salt participating medium in a vertical heated tube under forced and mixed convection laminar flows

The contribution of radiative heat transfer (RHT) to convective heat transfer in a heated vertical pipe with laminar flow and constant wall temperature, for FLiBe molten salt, is investigated computationally. The combined effects of conduction in the fluid, forced convection, buoyancy, temperature-dependent physical properties, and thermal radiation are investigated. The P1 approximation is employed in the discretization of the Radiative Transfer Equation (RTE). The COMSOL Multiphysics software is used to generate the numerical solutions. The theoretical analysis developed here demonstrates that only for intermediate values of the optical thickness, i.e. τ D ~ O 1 , the participating media effects are expected to be important. This analysis is confirmed computationally and a value of τ D ≈ 4 is shown to lead to the highest increase in overall heat transfer behavior. Under forced convection, the Nusselt ratio N u total / N u no - rad reaches a peak value of 1.76 at τ D = 4 and z/D = 200 , and is less than 1.1 for τ D < 0 . 1 and τ D > 60 . Under aiding-flow mixed convection, N u total / N u no - rad reaches a peak of 1.34 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 4 and τ D > 36 . Under opposing-flow mixed convection, Nu total / Nu no - rad reaches a peak of 1.81 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 2 and τ D > 50 . RHT effects on the overall heat transfer are most pronounced in opposing mixed convection, where N u total / N u no - rad is 2.26 observed at z/D = 50 . The sensitivity to wall emissivity is evaluated; in forced convection, the peak N u total / N u no - rad is 1.66 at ε = 0 (reflective wall) and 2.04 at ε = 1 (absorptive wall). Finally, the sensitivity to pipe diameter is also discussed. Overall, for a vertical heated tube, radiative heat transfer effects lead to enhancement of heat transfer by as high as a factor of two, and they depend on the optical thickness of the flow, mixed convection environment ( Gr / R e 2 and direction of flow relative to the gravitational force), surface emissivity, and entrance effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

The impact of cloud microphysics and ice nucleation on Southern Ocean clouds assessed with single-column modeling and instrument simulators

Abstract. Supercooled liquid clouds are common at higher latitudes (especially over the Southern Ocean) and are critical for constraining climate projections. We take advantage of the Macquarie Island Cloud and Radiation Experiment (MICRE) to perform an analysis of observed and simulated cloud processes over the Southern Ocean in a region and season dominated by supercooled liquid clouds. Using a single-column version of the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS), we compare two different cloud microphysical schemes to ground-based observations of cloud, precipitation, and radiation over a 2.5-month period (1 January–17 March 2017). Both schemes are able to reproduce aspects of the cloud and radiation observations during MICRE to within the uncertainty of the data when the thermodynamic profile is prescribed with relaxation. There are differences in water mass and representation of reflectivity between the schemes. A sensitivity study of the cloud microphysics schemes, one a bulk one-moment scheme and the other a two-moment scheme with prediction of mass and number, indicates that several key processes create differences between the schemes. Surface radiative fluxes and total water path are highly sensitive to the formation and fall speed of precipitation. The prediction of hydrometeor number with the two-moment scheme yields a better comparison with observed reflectivity and radiative fluxes, despite predicting higher liquid water contents than observed. With the two-moment scheme, we are also able to test the sensitivity of the results to the input of liquid cloud condensation nuclei (CCN) and ice nuclei (IN). The cloud properties and resulting radiative effects are found to be sensitive to the CCN and IN concentrations. More CCN and IN increase liquid and ice water paths, respectively. Thus, both the dynamic environment and aerosols, integrated through the cloud microphysics, are important for properly representing Southern Ocean cloud radiative effects.

58 GEOSCIENCES↗

DOE OSTI2020

Fully compressible simulations of waves and core convection in main-sequence stars

Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

79 ASTRONOMY AND ASTROPHYSICS↗