SEARCH · Engineering Papers

Results for “numeric profiling”

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 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Analysis of thermal radiation of a gas-particulate-cloud in a pressurized oxy-coal combustor

Pressurized oxy-combustion (POC) aims to reduce carbon emissions from coal-fired power plants. In the POC process, pulverized coal is burned at elevated pressure in an O2-CO2 environment. The elevated pressure and high CO2 and H2O concentrations strongly impact thermal radiation, thereby significantly distinguishing POC system from conventional atmospheric-pressure air-fired and oxy-fuel combustion. Consequently, a thorough understanding of thermal radiation in POC is required to promote the development of novel combustors. The present computational work analyzes thermal radiation of the gaseous and particle phases in the pilot-scale POC combustor in Washington University in St. Louis. Particularly, the radiation characteristics of the gaseous and particulate cloud were estimated by employing the statistical narrow band model and the large-particle model. It is found that thermal radiation of a particulate cloud dominates in the combustor under the conditions of a furnace temperature of 1500 K and no substantial particle loss. Another important result is that the gas-and-particulate cloud can be approximately treated as a graybody under the same conditions. To be more specific, the results of spectral radiation intensity of a gas comprised of 40% (vol) H2O and 60% CO2 show that when the pressure is increased to 15 bar, and the radiation pathlength to 100 cm, the spectral radiation profile of the gas phase approaches that of a blackbody at the respective temperature. In addition, the emissivity of the particulate cloud has been evaluated as a function of the particle concentration and diameter by means of the large particle model. It is shown that the emissivity grows with the particle concentration but decreases with the particle size for the same mass of a particulate cloud. The emissivity of a particulate cloud in oxy- combustion exceeds that in air-fired combustion. This work could be used to validate the assumption of the gray gas model adopted in numerical simulations.

Li, Lei↗

DOE OSTI2020

Ultra-high aspect ratio pores milled in diamond via laser, ion and electron beam mediated processes

Microfabrication techniques are critical for the rapid prototyping and development of applications for cutting edge materials. Recently diamond has gained considerable interest for quantum photonic, biosensing, inertial confinement fusion and magnetometer applications. In this article, ultra-high aspect ratio milling of diamond micropores by photon, ion and electron based methods is reported. A multiphoton absorption laser ablation approach is revealed to rapidly produce sub-10 μm diameter micropores in diamond with an aspect ratio of 14:1 and a tapered profile at the surface interface. Chemically active, oxygen focused ion beam milling produces high-aspect ratio pores in diamond with an aspect ratio of 65:1 and minimal tapering over the length of the pore, overcoming the physical interaction volume limitations imposed in conventional gallium based focused ion beam milling and laser ablation methods. Oxygen-mediated electron beam induced etching is revealed to negate the limitations imposed by physical sputtering mechanisms utilized in focused ion beam milling via the direct initiation of chemical reactions at the receding surface, producing aspect ratios on the order of 200:1. Numerical simulations reveal the physical basis for the superior aspect-ratio pore milling of the oxygen focused ion beam milling and electron beam induced etching methods. Our results demonstrate direct-write methods for the fabrication of ultra-high aspect micropores in diamond and provide insight into the underlying mechanisms of these physical processes. The three methods demonstrated here can be interchanged for applications based on the desired characteristic aspect ratio and process throughput.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Revealing the dynamics of ultrarelativistic non-equilibrium many-electron systems with phase space tomography

The description of physical processes with many-particle systems is a key approach to the modeling of numerous physical systems. For example in storage rings, where ultrarelativistic particles are agglomerated in dense bunches, the modeling and measurement of their phase-space distribution is of paramount importance: at any time the phase-space distribution not only determines the complete space-time evolution but also provides fundamental performance characteristics for storage ring operation. Here, we demonstrate a non-destructive tomographic imaging technique for the 2D longitudinal phase-space distribution of ultrarelativistic electron bunches. For this purpose, we utilize a unique setup, which streams turn-by-turn near-field measurements of bunch profiles at MHz repetition rates. To demonstrate the feasibility of our method, we induce a non-equilibrium state and show that the phase-space distribution microstructuring as well as the phase-space distribution dynamics can be observed in great detail. Our approach offers a pathway to control ultrashort bunches and supports, as one example, the development of compact accelerators with low energy footprints.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Cone beam neutron interferometry: From modeling to applications

Phase-grating moiré interferometers (PGMIs) have emerged as promising candidates for the next generation of neutron interferometry, enabling the use of a polychromatic beam and manifesting interference patterns that can be directly imaged by existing neutron cameras. However, the modeling of the various PGMI configurations is limited to cumbersome numerical calculations and backward propagation models which often do not enable one to explore the setup parameters. Here we generalize the Fresnel scaling theorem to introduce a k -space model for PGMI setups illuminated by a cone beam, thus enabling an intuitive forward propagation model for a wide range of parameters and experimental setups. The interference manifested by a PGMI is shown to be a special case of the Talbot effect, and the optimal fringe visibility is shown to occur at the moiré location of the Talbot distances. We derive analytical expressions for the contrast and the propagating intensity profiles in various conditions and provide the first analysis of the PGMI dark-field imaging signal when considering sample characterization. The model's predictions are compared to experimental measurements and good agreement is found between them. Last, we propose and experimentally verify a method to recover contrast at typically inaccessible PGMI autocorrelation lengths. The presented work provides a toolbox for analyzing and understanding existing PGMI setups and their future applications, for example extensions to two-dimensional PGMIs and characterization of samples with nontrivial structures. Published by the American Physical Society 2024

Sarenac, D. (ORCID:0000000185753367)↗

DOE OSTI2023

DLIO: A DATA-CENTRIC BENCHMARK FOR DEEP LEARNING APPLICATIONS

SF-22-136 Deep learning has been shown as a successful method for various tasks, and its popularity results in numerous open-source deep learning software tools. Deep learning has been applied to a broad spectrum of scientific domains such as cosmology, particle physics, computer vision, fusion, and astrophysics. Scientists have performed a great deal of work to optimize the computational performance of deep learning frameworks. However, the same cannot be said for I/O performance. As deep learning algorithms rely on big-data volume and variety to effectively train neural networks accurately, I/O is a significant bottleneck on large-scale distributed deep learning training. DLIO, is a novel representative benchmark suite built based on the I/O profiling of the selected workloads. DLIO can be utilized to accurately emulate the I/O behavior of modern deep learning applications. Using DLIO, application developers and system software solution architects can identify potential I/O bottlenecks in their applications and guide optimizations to boost the I/O performance leading to lower training times. The storage vendor can also use DLIO as a guide for designing and optimize the storage and filesystem targeting at deep learning application.

ZHENG, HUIHUO↗

DOE OSTI2021

Thermal-striping analysis methodology for sodium-cooled reactor design

Computational fluid dynamics simulations have been performed to study the applicability of engineering CFD methods for thermals striping analysis. Thermal striping is the fluctuating temperature profile in a solid caused by fluctuating fluid temperature, and the resulting fluctuating thermal stresses cause high cycle fatigue and eventual material failure. Here this study presents the methodology for thermal striping analysis that including a transient conjugate heat transfer RANS model and benchmark the results against a sodium triple jet experiment. The multiphysics calculations are implemented to analyze the thermal stresses in the solid domain caused by the coupled heat transfer between the fluid and solid domains. Validation data from a liquid sodium triple jet experiment include time-averaged temperature measurements and power spectra of the temperature signal. The numerical results agree well with these experimental measures, demonstrating key features such as the dominant frequency of temperature fluctuations. The applicability of a low-cost wall treatment method is demonstrated, enabling key computational savings. Finally, the performance of two finite element stress analysis software packages is compared, and the validity of the lower-cost method is confirmed. These results demonstrate the applicability of engineering methods for computational thermal striping calculations, enabling thermal striping estimations in large fluid systems such as the core of a Sodium-cooled nuclear reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Numerical Simulation of a Natural Convection–Driven Air-Cooled Reactor Cavity Cooling System Experiment

Ensuring the efficient removal of decay heat from the reactor vessel is essential for the safety of advanced reactor technologies. Several Generation-IV concepts incorporate variations in the reactor vessel cooling systems to achieve this objective. High-temperature gas-cooled reactors utilize a reactor cavity cooling system (RCCS), a passive ex-vessel system designed to operate without active components or external power during accident conditions. The RCCS removes decay heat primarily through radiative and convective heat transfer mechanisms. Here, this study presents a comprehensive validation of a computational fluid dynamics Reynolds-averaged Navier-Stokes model for the University of Wisconsin-Madison air-cooled RCCS facility. Validation was conducted for both high- and low-power natural convection cases under a uniform heating profile. Near-wall resolution was found to be critical for accurately modeling natural convection in the RCCS; employing an all-𝑦 + wall treatment resulted in wall temperature discrepancies exceeding 50 °⁢𝐶 compared to a wall-resolved mesh. Thermal-hydraulic behaviors under natural and forced convection conditions were compared within the heated cavity and RCCS. A turbulence model sensitivity analysis indicated that low-Reynolds number k-ɛ, k-ω shear stress transport (SST), and Reynolds stress transport models produce similar wall temperature predictions. A buoyancy modeling sensitivity study revealed that the Boussinesq approximation significantly underpredicted thermal-hydraulic behavior in the RCCS. Based on these findings, modeling recommendations are provided. The validated data set along with identified sensitivities refine the modeling of natural convection in the RCCS. The information produced by this study supports RCCS design, optimization, and safety evaluations, enabling the calibration and verification of reduced-order thermal-hydraulic models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Residual Stress Modeling with Phase Transformation for Wire Arc Additive Manufacturing of B91 Steel

Wire arc additive manufacturing (WAAM) is an energy-efficient manufacturing technique used for near-net-shape production of functional industrial components. However, heat accumulation during deposition and the associated mechanical and metallurgical changes result in complex residual stress profiles across the cross section of the fabricated components. These residual stresses are detrimental to the service life of the components. In this study, sequentially coupled thermomechanical analysis of WAAM B91 steel is conducted to quantify the residual stress variation across the component. The thermomechanical analysis includes a transient heat transfer model and a static stress model that incorporates the transformation-induced plasticity due to martensitic phase transformation. Here, the experimentally calibrated heat transfer model mirrors the temperature variation of the system during the deposition. The results from the stress model are validated via x-ray diffraction measurements, and the numerical results are in good agreement with the experimental data.

20 FOSSIL-FUELED POWER PLANTS↗

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

MHD-induced beta limits in the Large Helical Device

Using the extended-magnetohydrodynamics code, M3D-C1, we perform a systematic numerical study of the effect of externally applied heating on the achievable plasma beta in a ten field-period heliotron. Heat sources of varying intensity are applied to a vacuum magnetic field that is representative of the standard configuration of the Large Helical Device, with R 0 = 3.66 m, where R 0 is the radial position of the magnetic axis in vacuum. As the system is driven to a state that is unstable to low-n magnetohydrodynamic (MHD) modes, nonlinear mode interactions lead to the formation of chaotic magnetic fields. With sufficiently strong heating, a collapse of the electron temperature profile is observed. This demonstrates the necessity of simulating the self-consistent evolution of plasma profiles, without imposing assumptions on the structure of the magnetic field, to accurately determine transport properties in stellarator plasmas. It also highlights the value of these advanced simulation capabilities for accelerating the development of high-performance stellarator operating scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Rapid Screening of Liquid Metal Wetting for a Materials Compatibility Library

Wetting behavior of molten metals on solid substrates is a critical phenomenon influencing numerous industrial applications, including welding, anti-corrosion coatings, and metal additive manufacturing (AM). In particular, molten metal jetting (MMJ), an emerging AM technology, requires that the molten metal remain pinned at the nozzle exit. Thus, each new metal requires a specific nozzle material to ensure consistent droplet ejection and deposition, making it important to rapidly identify the appropriate wetting combinations. However, traditional measurements of wetting angles require expensive equipment and only allow one combination of materials to be investigated at a time which can be time consuming. This work introduces a rapid screening method based on sessile droplet experiments to evaluate wetting profiles across multiple metal–substrate combinations simultaneously. This study investigates the wetting interactions of molten Al alloy (Al4008), Cu, and Sn on various ceramic and metal substrates to identify optimal material combinations for MMJ nozzle designs. Results demonstrate that Al4008 achieves wetting on ceramic substrates such as AlN, TiO 2 , and SiC, with varying mechanisms including chemical reactions and weak surface interactions. Additionally, theoretical predictions regarding miscibility gaps and melting point differences were verified for Cu and Sn on refractory metals like Mo and W. Findings from this study contribute to the establishment of a materials compatibility library, enabling the selection of wetting/non-wetting combinations for stable MMJ operation. This resource not only advances MMJ technologies but also provides valuable insights for broader applications such as welding, coating, and printed electronics.

Materials science↗

DOE OSTI2021

Formation of edge pressure pedestal and reversed magnetic shear due to toroidal rotation in a tokamak equilibrium

Toroidal rotation is well known to play significant roles in the edge transport and L-H transition dynamics of tokamaks. Our recent calculation finds that a sufficiently strong localized toroidal rotation can directly bring out the formation of edge pressure pedestal with reversed magnetic shear that is reminiscent of an H-mode plasma, purely through the effects of toroidal rotation on the tokamak MHD equilibrium itself. In particular, the enhanced edge toroidal rotation enables a substantial peaking of the parallel current profile near edge in higher β regimes, which leads to the flattening or reversal of the local q (safety factor) profile. Here the formation of pressure pedestal along with the reversed magnetic shear region is shown to be the natural outcome of the MHD tokamak equilibrium in a self-consistent response to the presence of a localized toroidal rotation typically observed in H-mode or QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

How accurately can we detect the splashback radius of dark matter haloes and its correlation with accretion rate?

ABSTRACT The splashback radius (Rsp) of dark matter haloes has recently been detected using weak gravitational lensing and cross-correlations with galaxies. However, different methods have been used to measure Rsp and to assess the significance of its detection. In this paper, we use simulations to study the precision and accuracy to which we can detect the splashback radius with 3D density, 3D subhalo, and weak lensing profiles. We study how well various methods and tracers recover Rsp by comparing it with the value measured directly from particle dynamics. We show that estimates of Rsp from density and subhalo profiles correspond to different percentiles of the underlying Rsp distribution of particle orbits. At low accretion rates, a second caustic appears and can bias results. Finally, we show that upcoming lensing surveys may be able to constrain the Rsp–Γdyn relation directly.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Resolving Mesoscale Convective Systems: Grid Spacing Sensitivity in the Tropics and Midlatitudes

Abstract Mesoscale convective systems (MCSs) are a critical global water cycle component and drive extreme precipitation events in tropical and midlatitude regions. However, simulating deep convection remains challenging for modern numerical weather and climate models due to the complex interactions of processes from microscales to synoptic scales. Recent models with kilometer‐scale horizontal grid spacings offer notable improvements in simulating deep convection compared to coarser‐resolution models. Still, deficiencies in representing key physical processes, such as entrainment, lead to systematic biases. Additionally, evaluating model outputs using process‐oriented observational data remain difficult. This study presents an ensemble of MCS simulations with spanning the deep convective gray zone ( from 12 km to 125 m) in the Southern Great Plains of the U.S. and the Amazon Basin. Comparing these simulations with Atmospheric Radiation Measurement (ARM) wind profiler observations, we find greater sensitivity in the Amazon Basin compared to the Great Plains. Convective drafts converge structurally at sub‐kilometer scales, but some deficiencies remain. In both regions, simulated up and downdrafts are too deep and extreme downdrafts are not strong enough. Furthermore, Amazonian updrafts are too strong. Overall, we observe higher sensitivity in the tropics, including an artificial buildup in vertical kinetic energy at scales of , suggesting a need for 250 m in this region. Nevertheless, bulk convergence—agreement of storm‐average statistics—is achievable with kilometer‐scale simulations within a 10% error margin with 1 km providing a good balance between accuracy and computational cost.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Convective Updraft and Downdraft Characteristics of Continental Mesoscale Convective Systems in the Model Gray Zone

Here, the "gray zone" of convective modeling is defined as the range of horizontal grid spacings (Δx) at which turbulent transport processes are only partially resolved by the dynamics of the numerical model. This zone typically covers Δx from a few kilometers to several hundred meters, wherein the realistic representation of convective cloud processes can be challenging. This study characterizes the convective draft behaviors at multiple Δx across the gray zone and determines the appropriate Δx that can reliably capture these salient convective properties. We perform an ensemble of idealized simulations of mesoscale convective systems (MCS) using the Weather Research and Forecasting model at various Δx from 4 km to 250 m over the central U.S. An evaluation of key MCS kinematic properties is constrained using unique, long-term vertical velocity estimates obtained by radar wind profilers deployed by the Department of Energy Atmospheric Radiation Measurement user facility. MCS simulations for all Δx tested overestimate (underestimate) the probabilities of convective updrafts (downdrafts) compared to the observations. In terms of the convective draft intensity, finer-Δx models overestimate the updraft intensity, while the opposite is found for downdrafts. Moving from Δx = 4 km to 250 m, downdrafts become stronger and more frequent especially at middle and upper levels, attributed to additional drag from an increasing graupel frequency and compensation of enhanced updraft velocity. Simulated draft characteristics, including core size, intensity, and probability of occurrence, exhibit pronounced changes at Δx greater than 500 m, but suggest less sensitivity when Δx is reduced below 500 m.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Measurement of laser absorption in underdense plasmas using near-field imaging of the incident and transmitted beams

Measurements of laser absorption in high-temperature, underdense plasmas produced at the Omega Laser Facility are made using two near-field imaging detectors that diagnose the spatial profile and energy of the port P9 beam before and after it transmits through the plasma. The incident beam is sampled using a partial reflection from a full-aperture, (30 cm-diam) uncoated wedge pickoff located before the target chamber vacuum window and final focus lens assembly. A concave mirror reduces the reflected beam size, allowing it to be recorded directly using a charged-coupled device (CCD) camera. The P9 transmitted beam diagnostic (P9TBD) characterizes the transmitted light by terminating the expanded beam on a semi-transparent diffuser and imaging the illuminated surface using a lens and CCD camera. The P9TBD samples a numerical aperture twice as large as the input beam, allowing the energy of transmitted beams with moderate levels of beam spray to be measured accurately. Calibration shots with no plasma provide a path to infer absorption without absolute photometric calibration of either detector. Lastly, the cross-calibration between the two detectors was measured to remain stable at ±200 ppm, enabling measurements of total beam absorption below 1% with ±0.07% error.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗

DOE OSTI2026

Final report: Insights from ARM observations into aerosol processing and transport by extratropical cyclones and aerosol effects on cyclone clouds

This project advanced understanding of aerosol–cloud interactions in extratropical cyclones using observations from the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) observatory and ACE-ENA field campaigns, along with satellite datasets and numerical modeling. We have quantified the susceptibility of frontal clouds to aerosols, showing that this susceptibility is similar to that in non-frontal clouds despite substantial (factor 2) differences in droplet concentrations between these cloud types, and also significant differences in cloud albedo. A key finding was that Aitken-mode (here, 60-100nm diameter) aerosols play a disproportionately important role in droplet activation in frontal clouds, and concentrations of these particles are frequently strongly biased in climate models. Activation of aerosols depends strongly on updraft speeds, and via collaboration with Argonne National Laboratory we contributed to quantifying these better with the ARM radar wind profiler at the ENA site. Our studies in the North Atlantic have analogs in the Southern Ocean, and to draw these out we examined data from the MARCUS and SOCRATES field campaigns. We identified new particle formation within extratropical cyclones as a potentially important source of cloud condensation nuclei in these pristine marine environments. By combining ARM observations with perturbed parameter ensemble simulations, we demonstrated that surface observations can effectively constrain aerosol–cloud adjustments and reduce uncertainty in simulated cloud liquid water path responses by approximately 15%, yielding stronger constraints on historical aerosol cooling effects. The work further established precipitation processes as a dominant control on aerosol–cloud interactions and Earth system predictability.

Gordon, Hamish [Department of Chemical Engineering↗