Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NATURAL CONVECTION”

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 1,027 records · Page 57

Imaging Plasma Density Structures in the Soft X-Rays Generated by Solar Wind Charge Exchange with Neutrals

Both heliophysics and planetary physics seek to understand the complex nature of the solar wind's interaction with solar system obstacles like Earth's magnetosphere, the ionospheres of Venus and Mars, and comets. Studies with this objective are frequently conducted with the help of single or multipoint in situ electromagnetic field and particle observations, guided by the predictions of both local and global numerical simulations, and placed in context by observations from far and extreme ultraviolet (FUV, EUV), hard X-ray, and energetic neutral atom imagers (ENA). Each proposed interaction mechanism (e.g., steady or transient magnetic reconnection, local or global magnetic reconnection, ion pick-up, or the Kelvin- Helmholtz instability) generates diagnostic plasma density structures. The significance of each mechanism to the overall interaction (as measured in terms of atmospheric/ionospheric loss at comets, Venus, and Mars or global magnetospheric/ionospheric convection at Earth) remains to be determined but can be evaluated on the basis of how often the density signatures that it generates are observed as a function of solar wind conditions. This paper reviews efforts to image the diagnostic plasma density structures in the soft (low energy, 0.1-2.0 keV) X-rays produced when high charge state solar wind ions exchange electrons with the exospheric neutrals surrounding solar system obstacles.

Sibeck, David G.↗

Some aspects of core formation in Mercury

Some questions dealing with the nature and history of a large metallic core within Mercury are considered. These include the existence of a core, its size, whether it is fluid or solid, the timescale for core formation, the geological consequences of core formation, and whether such consequences are consistent with the surface geology. Several indirect lines of evidence are discussed which suggest the presence of a large iron-rich core. A core-formation model is examined in which core infall is accompanied by an increase of 17 km in planetary radius, an increase of 700 K in mean internal temperature, and substantial melting of the mantle. It is argued that if the core differentiated from an originally homogeneous planet, that event must have predated the oldest geological units comprising most of the planetary surface. A convective dynamo model for the source of Mercury's magnetic field is shown to conflict with cosmochemical models that do not predict a substantial radiogenic heat source in the core.

Solomon, S. C.↗

Exploring the Inorganic Composition of the Asian Tropopause Aerosol Layer Using Medium-Duration Balloon Flights

Satellite observations have revealed an enhanced aerosol layer near the tropopause over Asia during the summer monsoon, called the “Asian Tropopause Aerosol Layer” (ATAL). In this work, aerosol particles in the ATAL were collected with a balloon-borne impactor near the tropopause region over India, using extended-duration balloon flights, in summer 2017 and winter 2018. The chemical composition of these particles was further investigated by quantitative analysis using offline ion chromatography. Nitrate (NO_3 ^1) and nitrite (NO_2 ^-) were found to be the dominant ions in the collected aerosols with values ranging between 87 and 343 ng m^(−3) at STP (standard temperature and pressure) during the summer campaign. In contrast, sulfate (SO_4 ^(2-)) levels were found to be above the detection limit (>10 ng m^(−3) at STP) only in winter. In addition, we determined the origin of the air masses sampled during the flights using the analysis of back trajectories as well as a convective proxy from cloud-top temperature fields derived from a geostationary satellite. The results obtained from this analysis were put into the context of large-scale transport and aerosol distribution using GEOS-Chem chemical transport model simulations. The first flight in summer 2017 which sampled an air mass within the Asian monsoon anticyclone (AMA), influenced by convection over Western China, was associated with particle size diameters from 0.05 to 0.15 µm. In contrast, the second flight sampled air masses at the edge of the AMA associated with a larger particle size radius (>2 µm) with a higher NO_2 ^- concentration. The sampled air masses in winter 2018 were likely affected by smoke from the Pacific Northwest fire event in Canada, which occurred 7 months before our campaign, associated with concentration enhancements of SO_4 (^2-) and Ca^(2+). Overall, our results suggest that nitrogen-containing particles represent a large fraction of cloud-free and in-cloud aerosols populating the ATAL, which is partially in agreement with the results from aircraft measurements during the StratoClim (Stratospheric and upper tropospheric processes for better climate predictions) campaign. The exact nature of those particles is still unknown, but their coincidences with subvisible cirrus clouds and their sizes suggest nitric acid trihydrate (NAT) as a possible candidate, as NAT has already been observed in the tropical upper troposphere and lower stratosphere in other studies. Furthermore, GEOS-Chem model simulations indicate that lightning NOx emissions could have significantly impacted the production of nitrate aerosols sampled during the summer of 2017.

Hazel Vernier↗

The role of radiation-dynamics interaction in regional numerical weather prediction

The role of radiation-dynamics interaction in regional numerical weather prediction of severe storm environment and mesoscale convective systems over the United States is researched. Based upon the earlier numerical model simulation experiments, it is believed that such interaction can have a profound impact on the dynamics and thermodynamics of regional weather systems. The research will be carried out using real-data model forecast experiments performed on the Cray-X/MP computer. The forecasting system to be used is a comprehensive mesoscale prediction system which includes analysis and initialization, the dynamic model, and the post-forecast diagnosis codes. The model physics are currently undergoing many improvements in parameterizing radiation processes in the model atmosphere. The forecast experiments in conjunction with in-depth model verification and diagnosis are aimed at a quantitative understanding of the interaction between atmospheric radiation and regional dynamical processes in mesoscale models as well as in nature. Thus, significant advances in regional numerical weather prediction can be made. Results shall also provide valuable information for observational designs in the area of remote sensing techniques to study the characteristics of air-land thermal interaction and moist processes under various atmospheric conditions.

Chang, Chia-Bo↗

Coupling the Solar-Wind/IMF to the Ionosphere through the High Latitude Cusps

Magnetic merging is a primary means for coupling energy from the solar wind into the magnetosphere-ionosphere system. The location and nature of the process remain as open questions. By correlating measurements form diverse locations and using large-scale MHD models to put the measurements in context, it is possible to constrain out interpretations of the global and meso-scale dynamics of magnetic merging. Recent evidence demonstrates that merging often occurs at high latitudes in the vicinity of the cusps. The location is in part controlled by the clock angle in the interplanetary magnetic field (IMF) Y-Z plane. In fact, B(sub Y) bifurcated the cusp relative to source regions. The newly opened field lines may couple to the ionosphere at MLT locations of as much as 3 hr away from local noon. On the other side of noon the cusp may be connected to merging sites in the opposite hemisphere. In face, the small convection cell is generally driven by opposite hemisphere merging. B(sub X) controls the timing of the interaction and merging sites in each hemisphere, which may respond to planar features in the IMF at different times. Correlation times are variable and are controlled by the dynamics of the tilt of the interplanetary electric field phase plane. The orientation of the phase plane may change significantly on time scales of tens of minutes. Merging is temporally variable and may be occurring at multiple sites simultaneously. Accelerated electrons from the merging process excite optical signatures at the foot of the newly opened field lines. All-sky photometer observations of 557.7 nm emissions in the cusp region provide a "television picture" of the merging process and may be used to infer the temporal and spatial variability of merging, tied to variations in the IMF.

Maynard, Nelson C.↗

Estimation of the viscosities of silicate liquids at high pressure from measurements of oxygen diffusivity

The dynamics and evolution of a magma ocean depend significantly on the structure and physical properties of silicate liquids at high pressure. Viscosity is particularly important because of its effect on the dynamics of convection and crystal setting. Traditionally, the viscosities of silicate liquids have been measured at high pressure by falling sphere viscometry but, due to some limitations of this technique, few measurements have been made at pressures exceeding 2.5 GPa. An alternative approach, which has previously been proven up to 2 GPa, is to use high pressure measurements of oxygen self-diffusivity (D) to estimate viscosity (eta) from the Eyring relationship eta = kT/D lambda (where k is the Boltzmann constant, T is absolute temperature, and lambda is the diffusive jump distance). We have performed preliminary experiments on Na2Si4O9 liquid in order to explore the feasibility of this method up to pressures in excess of 10 GPa. Diffusion couples were prepared from glass starting materials with one half of the sample enriched in O-18 and the other half containing the natural abundance (0.2 wt percent). Diffusion experiments were performed at 1600 to 1800 C at pressures in the range of 2.5 to 10 GPa for times up to 6 minutes using a 1200 ton multianvil apparatus. Oxygen diffusivities were derived from the resulting O-18 concentration profiles, which were analyzed using an ion microprobe. At 1800 C, oxygen diffusivities increase continuously from 1 x 10(exp -10) m(exp -2) s(exp -1) at 2.5 Gpa to 5 x 10(exp -10) m(exp -2) s(exp -1) at 10 GPa. From the Eyring relationship, the viscosity of Na2SiO9 liquid is predicted to decrease by approx. 0.7 log units over this pressure range at 1800 C. These trends, which can be related to changes in the structure of the liquid at high pressure, are in agreement with results of molecular dynamics calculations. These results suggest that it should be possible to estimate the viscosities of silicate liquids up to at least 15 GPa and 2200 C by this method.

Rubie, D. C.↗

Role of the ocean in climate changes

The present program aimed at the study of ocean climate change is prepared by a group of scientists from State Oceanographic Institute, Academy of Science of Russia, Academy of Science of Ukraine and Moscow State University. It appears to be a natural evolution of ideas and achievements that have been developed under national and international ocean research projects such as SECTIONS, WOCE, TOGA, JGOFS and others. The two primary goals are set in the program ROCC. (1) Quantitative description of the global interoceanic 'conveyor' and it's role in formation of the large scale anomalies in the North Atlantic. The objectives on the way to this goal are: to get the reliable estimates of year-to-year variations of heat and water exchange between the Atlantic Ocean and the atmosphere; to establish and understand the physics of long period variations in meridianal heat and fresh water transport (MHT and MFWT) in the Atlantic Ocean; to analyze the general mechanisms, that form the MHT and MFWT in low latitudes (Ekman flux), middle latitudes (western boundary currents) and high latitudes (deep convection) of the North Atlantic; to establish and to give quantitative description of the realization of global changes in SST, surface salinity, sea level and sea ice data. (2) Development of the observational system pointed at tracing the climate changes in the North Atlantic. This goal merges the following objectives: to find the proper sites that form the inter annual variations of MHT; to study the deep circulation in the 'key' points; to develop the circulation models reflecting the principle features of interoceanic circulation; and to define global and local response of the atmosphere circulation to large scale processes in the Atlantic Ocean.

Gulev, Sergey K.↗

New Empirical Formulation for the Sublimational Breakup of Graupel and Dendritic Snow

Abstract Ice fragments are generated by sublimation of ice particles in subsaturated conditions in natural clouds. Conceivably, such sublimational breakup would be expected to cause ice multiplication in natural clouds. Any fragment that survives will grow to become ice precipitation that may sublimate and fragment further. As a first step toward assessing this overlooked process, a formulation is proposed for the number of ice fragments from sublimation of ice particles for an atmospheric model. This is done by amalgamating laboratory observations from previously published studies. The concept of a “sublimated mass activity spectrum” for the breakup is applied to the dataset. The number of ice fragments is determined by the relative humidity over ice and the initial size of the parent ice particles. The new formulation applies to dendritic crystals and heavily rimed particles only. Finally, a thought experiment is performed for an idealized scenario of subsaturation with in-cloud descent. Scaling analysis yields an estimate of an ice enhancement ratio of about 5 (10) within a weak deep convective downdraft of about 2 m s −1 , for an initial monodisperse population of dendritic snow (graupel) particles of 3 L −1 and 2 mm. During descent, there is a dynamic equilibrium between continual emission of fragments and their depletion by sublimation. A simplified bin microphysics parcel model exhibits this dynamical quasi equilibrium, consistent with the thought experiment. The fragments have average lifetimes of around 90 and 70 s for dendrites and graupel, respectively. Sublimational breakup is predicted to cause significant secondary ice production.

54 ENVIRONMENTAL SCIENCES↗

Interaction between Convection and Heat Transfer in Crystal Growth

Crystals are integral components in some of our most sophisticated and rapidly developing industries. Single crystals are solids with the most uniform structures that can be obtained on an atomic scale. Because of their structural uniformity, crystals can transmit acoustic and electromagnetic waves and charged particles with essentially no scattering or interferences. This transparency, which can be selectively modified by controlled additions of impurities known as dopants, is the foundation of modern electronic industry. It has brought about widespread application of crystals in transistors, lasers, microwave devices, infrared detectors, magnetic memory devices, and many other magnets and electro-optic components. The performance of a crystal depends strongly on its compositional homogeneity. For instance, in modern microcircuitry, compositional variations of a few percent (down to a submicron length scale) can seriously jeopardize predicted yields. Since crystals are grown by carefully controlled phase transformations, the compositional adjustment in the solid is often made during growth from the nutrient. Hence, a detailed understanding of mass transfer in the nutrient is essential. Moreover, since mass transfer is often the slowest process during growth, it is usually the rate limiting mechanism. Crystal growth processes are usually classified according to the nature of the parent phase. Nevertheless, whether the growth occurs by solidification from a melt (melt growth), nucleation from a solution (solution growth), condensation from a vapor (physical vapor transport) or chemical reaction of gases (chemical vapor deposition), the parent phase is a fluid. As is with most non-equilibrium processes involving fluids, liquid or vapor, fluid motion plays an important role, affecting both the concentration and temperature gradients at the soli-liquid interface.

Source record↗

Enabling residential heating decarbonization through hydronic low-temperature thermal distribution using forced-air assistive devices

Space heating represents approximately one-tenth of the United States’ energy use and has a breadth of potential for emission reduction. An element of space heating, hydronic heat distribution methods, use water supply temperatures up to 82.2 °C (180 °F). However, this operating temperature can be incompatible with high-efficiency heat generation systems, which typically provide heating temperatures of up to 60 °C (140 °F). In this work, a low-cost retrofit solution is developed using experimentally validated computational tools by incorporating an airflow distributor that preferentially directs the airflow from a fan to enhance heat transfer over the finned-tube heat exchanger found in conventional baseboards. The same heat output of traditional baseboard heat distribution systems operating at higher water supply temperatures (71.1 – 82.2 °C) can also be achieved at lower temperatures (≤60 °C). Specifically, results indicate that the technology can produce more than a 46.7% improvement in the heat transfer output at temperatures as low as 60 °C (140 °F), effectively matching the same output range achieved by 71.1 – 82.2 °C (160–180 °F) water supply temperatures. An important benefit of this solution is that it utilizes the existing infrastructure as opposed to requiring the replacement of the entire distribution system. Additionally, this technology enables existing building infrastructure to be coupled with newer, high-efficiency heating systems such as condensing boilers, solar-thermal systems, and geothermal/air-to-water heat pumps that may produce lower water supply temperatures (60 °C). A geospatially resolved techno-economic and environmental analysis is completed and presented to further understand the equivalent carbon footprint of enabling higher efficiency heat generation systems with the improved efficiency heat distribution system disclosed. Using 2021 grid-average emission factors, an emission reduction of up to a 67.5% decrease (2658 kgCO 2 /yr) for a single-family home, depending on state and climate region, could be realized by replacing a traditional natural gas-fired boiler with an air-to-water heat pump coupled to the same high-temperature heat distribution system along with the low-cost retrofit solution. A complete CO 2 emission reduction of 6688 kgCO 2 /yr, depending on state and climate region, could be realized if all electricity is further renewably sourced. Thus, this study provides a possible pathway towards enabling the reduction of operational emissions in space heating.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Evaluating 3 decades of precipitation in the Upper Colorado River basin from a high-resolution regional climate model

Abstract. Convection-permitting regional climate models (RCMs) have recently become tractable for applications at multi-decadal timescales. These types of models have tremendous utility for water resource studies, but better characterization of precipitation biases is needed, particularly for water-resource-critical mountain regions, where precipitation is highly variable in space, observations are sparse, and the societal water need is great. This study examines 34 years (1987–2020) of RCM precipitation from the Weather Research and Forecasting model (WRF; v3.8.1), using the Climate Forecast System Reanalysis (CFS; CFSv2) initial and lateral boundary conditions and a 1 km × 1 km innermost grid spacing. The RCM is centered over the Upper Colorado River basin, with a focus on the high-elevation, 750 km2 East River watershed (ERW), where a variety of high-impact scientific activities are currently ongoing. Precipitation is compared against point observations (Natural Resources Conservation Service Snow Telemetry or SNOTEL), gridded climate datasets (Newman, Livneh, and PRISM), and Bayesian reconstructions of watershed mean precipitation conditioned on streamflow and high-resolution snow remote-sensing products. We find that the cool-season precipitation percent error between WRF and 23 SNOTEL gauges has a low overall bias (x^ = 0.25 %, s = 13.63 %) and that WRF has a higher percent error during the warm season (x^ = 10.37 %, s = 12.79 %). Warm-season bias manifests as a high number of low-precipitation days, though the low-resolution or SNOTEL gauges limit some of the conclusions that can be drawn. Regional comparisons between WRF precipitation accumulation and three different gridded datasets show differences on the order of ± 20 %, particularly at the highest elevations and in keeping with findings from other studies. We find that WRF agrees slightly better with the Bayesian reconstruction of precipitation in the ERW compared to the gridded precipitation datasets, particularly when changing SNOTEL densities are taken into account. The conclusions are that the RCM reasonably captures orographic precipitation in this region and demonstrates that leveraging additional hydrologic information (streamflow and snow remote-sensing data) improves the ability to characterize biases in RCM precipitation fields. Error characteristics reported in this study are essential for leveraging the RCM model outputs for studies of past and future climates and water resource applications. The methods developed in this study can be applied to other watersheds and model configurations. Hourly 1 km × 1 km precipitation and other meteorological outputs from this dataset are publicly available and suitable for a wide variety of applications.

Geology↗

Effect of volcanic emissions on clouds during the 2008 and 2018 Kilauea degassing events

Volcanic eruptions in otherwise clean environments are “natural experiments” wherein the effects of aerosol emissions on clouds and climate can be partitioned from meteorological variability and anthropogenic activities. In this work, we combined satellite retrievals, reanalysis products, and atmospheric modeling to analyze the mechanisms of aerosol–cloud interactions during two degassing events at the Kilauea volcano in 2008 and 2018. The eruptive nature of the 2008 and 2018 degassing events was distinct from long-term volcanic activity for Kilauea. Although previous studies assessed the modulation of cloud properties from the 2008 event, this is the first time such an analysis has been reported for the 2018 event and that multiple degassing events have been analyzed and compared at this location. Both events resulted in significant changes in cloud effective radius and cloud droplet number concentration that were decoupled from local meteorology and in line with an enhanced cloud albedo. However, it is likely that the effects of volcanic emissions on liquid water path and cloud fraction were largely offset by meteorological variability. Comparison of cloud anomalies between the two events suggested a threshold response of aerosol–cloud interactions to overcome meteorological effects, largely controlled by aerosol loading. In both events, the ingestion of aerosols within convective parcels enhanced the detrainment of condensate in the upper troposphere, resulting in deeper clouds than observed under pristine conditions. Accounting for ice nucleation on ash particles led to enhanced ice crystal concentrations at cirrus levels and a slight decrease in ice water content, improving the correlation of the model results with the satellite retrievals. Overall, aerosol loading, plume characteristics, and meteorology contributed to changes in cloud properties during the Kilauea degassing events.

Katherine H. Breen↗

Sooting Limits Of Diffusion Flames With Oxygen-Enriched Air And Diluted Fuel

Oxygen-enhanced combustion permits certain benefits and flexibility that are not otherwise available in the design of practical combustors, as discussed by Baukal. The cost of pure and enriched oxygen has declined to the point that oxygen-enhanced combustion is preferable to combustion in air for many applications. Carbon sequestration is greatly facilitated by oxygen enrichment because nitrogen can be eliminated from the product stream. For example, when natural gas (or natural gas diluted with CO2) is burned in pure oxygen, the only significant products are water and CO2. Oxygen-enhanced combustion also has important implications for soot formation, as explored in this work. We propose that soot inception in nonpremixed flames requires a region where C/O ratio, temperature, and residence time are above certain critical values. Soot does not form at low temperatures, with the threshold in nonpremixed flames ranging from about 1250-1650 K, a temperature referred to here as the critical temperature for soot inception, Tc. Soot inception also can be suppressed when residence time is short (equivalently, when the strain rate in counterflow flames is high). Soot induction times of 0.8-15 ms were reported by Tesner and Shurupov for acetylene/nitrogen mixtures at 1473 K. Burner stabilized spherical microgravity flames are employed in this work for two main reasons. First, this configuration offers unrestricted control over convection direction. Second, in steady state these flames are strain-free and thus can yield intrinsic sooting limits in diffusion flames, similar to the way past work in premixed flames has provided intrinsic values of C/O ratio associated with soot inception limits.

Sunderland, P. B.↗

Predicting the Evolution of Shallow Cumulus Clouds With a Lotka‐Volterra Like Model

Abstract In numerical weather prediction and climate models, boundary‐layer clouds are controlled by a wide range of subgrid‐scale processes. However, understanding the nature of these processes and their role in the evolution of the cloud size distribution as a whole has been elusive. To address this issue, we adopt a novel empirical framework from the field of population dynamics to model the evolution of cloud size statistics by using the shallow cumulus properties obtained from a large‐eddy simulation (LES). Our approach involves representing the cloud size distribution and the total cloud area using a revised Lotka‐Volterra model and ridge linear model, respectively. The physical interpretation of the total cloud area and coefficients obtained from the optimization of the models reveals three stages probably interpreted by dominant processes: the formation of new clouds, the growth of single clouds, and a steady state with organized transitions involving the growth and decay of multiple clouds. Furthermore, we showcase the potential of this framework to serve as a component of scale‐aware parameterizations of shallow‐convective clouds in atmospheric models.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of the 7-km GEOS-5 Nature Run

This report documents an evaluation by the Global Modeling and Assimilation Office (GMAO) of a two-year 7-km-resolution non-hydrostatic global mesoscale simulation produced with the Goddard Earth Observing System (GEOS-5) atmospheric general circulation model. The simulation was produced as a Nature Run for conducting observing system simulation experiments (OSSEs). Generation of the GEOS-5 Nature Run (G5NR) was motivated in part by the desire of the OSSE community for an improved high-resolution sequel to an existing Nature Run produced by the European Centre for Medium-Range Weather Forecasts (ECMWF), which has served the community for several years. The intended use of the G5NR in this context is for generating simulated observations to test proposed observing system designs regarding new instruments and their deployments. Because NASA's interest in OSSEs extends beyond traditional weather forecasting applications, the G5NR includes, in addition to standard meteorological components, a suite of aerosol types and several trace gas concentrations, with emissions downscaled to 10 km using ancillary information such as power plant location, population density and night-light information. The evaluation exercise described here involved more than twenty-five GMAO scientists investigating various aspects of the G5NR performance, including time mean temperature and wind fields, energy spectra, precipitation and the hydrological cycle, the representation of waves, tropical cyclones and midlatitude storms, land and ocean surface characteristics, the representation and forcing effects of clouds and radiation, dynamics of the stratosphere and mesosphere, and the representation of aerosols and trace gases. Comparisons are made with observational data sets when possible, as well as with reanalyses and other long model simulations. The evaluation is broad in scope, as it is meant to assess the overall realism of basic aspects of the G5NR deemed relevant to the conduct of OSSEs. However, because of the relatively short record and other practical considerations, these comparisons cannot provide a definitive, statistically sound assessment of all model deficiencies, or guarantee the G5NR's suitability for all OSSE applications. Differences between the observed and simulated behavior also must be judged in the context of basic internal atmospheric variability which can introduce variations that are not necessarily controlled by the prescribed sea surface temperatures used in generating the G5NR. The results show that the G5NR performs well as measured by the majority of metrics applied in this evaluation. Particular benefits derived from the 7-km resolution of G5NR include realistic representations of extreme weather events in both the tropics and extratropics including tropical cyclones, Nor'easters and mesoscale convective complexes; improved representation of the diurnal cycle of precipitation over land; well-resolved surface-atmosphere interactions such as katabatic wind flows over Antarctica and Greenland; and resolution of orographically generated gravity waves that propagate into the upper atmosphere and influence the large scale circulation. Obvious deficiencies in the G5NR include a "splitting" of the inter-tropical convergence zone, which leads to a weaker-than-observed Hadley circulation and related deficiencies in the depiction of stationary wave patterns. Also, while the G5NR captures global cloud features and radiative effects well in general, close comparison with observations reveals higher-than-observed cloud brightness, likely due to an overabundance of cloud condensate; less distinct cloud minima in subtropical subsidence zones, consistent with a weak Hadley circualtion; and too few near-coastal marine stratocumulus clouds.

GEOS-5↗

Review of Recent Results in Global MHD Modeling: ISTP Project Scientist for Theory and Ground-Based Observations

Global MHD (magnetohydrodynamic) simulations have shown a remarkable ability to describe the global dynamics of geospace. The limitations of the physical approximations underlying MHD would seem to limit the effectiveness of these codes, since kinetic and hybrid effects should manifest themselves by cross-scale coupling from microscales to mesoscales to global scales. However three effects appear to allow the codes to operate much more successfully than one would at first believe. They are:(l) the globally self-consistent nature of the codes with very well defined exterior boundary conditions (the solar wind) which allows the proper intercommunication between magnetospheric regions on MHD scales, (2) the control by global dynamics of the boundary layer locations where micro and meso scale processes operate, and (3) the critical role of numerical diffusion and with a sufficiently high resolution grid, the use of an empirical resistivity term, which if set at a level where the major magnetosphere boundaries properly calibrate against their observed locations, appear to well represent the effects of kinetic and hybrid processes on the global dynamics. The effectiveness of the global MHD codes, which have been developed under the ISTP mission, in describing Wind, Polar and Geotail observations, as well as ground-based observations are described. Particular emphasis is placed upon the Polar imaging data which when combined with ground-based data and global MHD-based synthetic aurora and convection patterns provide a powerful tool in understanding the final link in the solar-terrestrial chain: coupling into the atmosphere and ionosphere.

Curtis, Steven↗

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗

ONERA-NASA Cooperative Effort on Liner Impedance Eduction

As part of a cooperation between ONERA and NASA, the liner impedance eduction methods developed by the two research centers are compared. The NASA technique relies on an objective function built on acoustic pressure measurements located on the wall opposite the test liner, and the propagation code solves the convected Helmholtz equation in uniform ow using a finite element method that implements a continuous Galerkin discretization. The ONERA method uses an objective function based either on wall acoustic pressure or on acoustic velocity acquired above the liner by Laser Doppler Anemometry, and the propagation code solves the linearized Euler equations by a discontinuous Galerkin discretization. Two acoustic liners are tested in both ONERA and NASA ow ducts and the measured data are treated with the corresponding impedance eduction method. The first liner is a wire mesh facesheet mounted onto a honeycomb core, designed to be linear with respect to incident sound pressure level and to grazing ow velocity. The second one is a conventional, nonlinear, perforate-over-honeycomb single layer liner. Configurations without and with ow are considered. For the nonlinear liner, the comparison of liner impedance educed by NASA and ONERA shows a sensitivity to the experimental conditions, namely to the nature of the source and to the sample width.

Primus, Julien↗