SEARCH · Search NASA

Results for “traversibility”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Mercury linked to Deccan Traps volcanism, climate change and the end-Cretaceous mass extinction

Mercury (Hg) anomalies linked to Large Igneous Provinces (LIP) volcanism have been identified in sediments across all five major mass extinctions in Earth's history. This study tests whether Hg in marine sediments is a reliable proxy linking Deccan Traps volcanic eruptions to late Maastrichtian global climate warming and the mass extinction at the Cretaceous-Paleogene boundary (KPB). Our primary test site is the Elles section in Tunisia, the auxiliary Global Stratotype Section and Point (GSSP) to El Kef. Elles has the most complete marine sedimentary record and a high average sedimentation rate of ~4.7 cm/ky. Here, we chose the Hor Hahar section in Israel to corroborate the geographic distribution of Hg fallout from Deccan volcanism. Reliability of the Hg proxy over the last 550 ky of the Maastrichtian to early Danian was evaluated based on high-resolution age control (orbital cyclostratigraphy), stable isotope climate record, Hg concentrations, biotic turnover and mass extinction. These results were correlated with the pulsed Deccan eruptive history constrained previously by U-Pb zircon geochronology. Our results support Hg as robust proxy for Deccan volcanism with large Hg spikes marking “extreme event” (EE) pulsed eruptions correlative with climate warming peaks separated by steady, less intense eruptions. Long-term global climate warming began near ~350 ky pre-KPB, reached maximum warming (3–4 °C) between 285 and 200 ky pre-KPB, followed by gradual cooling and rapid temperature drop between 45 and 25 ky pre-KPB. During the last 25 ky before the KPB, multiple Hg EE eruptions correlate with hyperthermal warming that culminated in the rapid mass extinction at Elles during ≤1000 years of the Cretaceous. These latest Cretaceous Hg peaks may correlate with massive, distal, Deccan-sourced lava flows (>1000 km long) that traversed the Indian subcontinent and flowed into the Bay of Bengal, bracketing the mass extinction. These results support Deccan volcanism as a primary driver of the end-Cretaceous mass extinction.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A multi-dimensional dataset for two-phase instability in low pressure natural circulation based on direct transient local measurement

Thermal hydraulic instability can challenge the control and the safety of nuclear reactors. This experimental research captures and investigates the instability in a low-pressure (1–10 atmospheres) natural circulation loop with an adiabatic chimney installed above the heated section. Direct transient local measurement of two-phase flow is performed at multiple axial locations along the flow channel with thermocouples and conductivity probes traversed by customized high-precision mechanisms. Based on the periodicity of oscillations, a post-processing method is proposed to align the data captured non-simultaneously to reconstruct the instability behavior in one representative period. The averaging operation over multiple oscillation cycles is performed to reduce the random errors. Cyclic phenomena are clearly demonstrated by the post-processed data, and rapid void fraction change induced by the flashing is observed to be significant in the driving mechanism for low pressure oscillations. Parametric effects of inlet subcooling, heat flux, inlet flow restriction and system pressure are studied within the experiment range, and qualitative results are explained and compared with existing literature. The ratio between the oscillation period and the travelling time of the fluid particle through the test section is compared with the theoretical value for the density wave oscillation. The new dataset contains directly measured transient multi-dimensional two-phase data along the flow channel for natural circulation under density wave oscillations. The generated dataset is valuable for the validation of system analysis codes and CFD codes in the simulation of natural circulation phenomena and dynamic two-phase flow.

42 ENGINEERING↗

DOE OSTI2022

Numerical simulation of falling film sensible heat transfer over round horizontal tubes

Numerical simulations are performed to explore the heat transfer characteristics of falling films over horizontal round tubes with uniform heat flux imposed for a range of Reynolds numbers spanning the droplet, jet, and sheet regimes. Simulations results agree well with available experimental measurements. The study analyzes the local as well as average heat transfer behavior under the different flow modes. The numerical results show that the local Nusselt number (Nu) distribution depends on the flow features in each mode and varies substantially in all directions for the respective mode. In the droplet mode, the Nu value varies significantly as the droplet impinges and the remnant liquid-bridge retracts (peak instantaneous Nu near 6), followed by wave propagation over the tube surface with peak Nu around 0.25. For the jet modes, the local maximum in the heat transfer occurs off-center to the impingement location with magnitudes of peak Nu = 3.1 for the inline jet mode and Nu = 2.7 for the staggered jet mode, while on the rest of the tube surface, it has an inverse relation with the liquid film thickness. Substantial variations in the heat transfer value are also recorded in the middle of the two impinging jets with Nu = 0.95 in the inline jet mode where the neighboring jets do not interact, and Nu = 0.60 in the crest region of the staggered jet mode where the neighboring jets interact with each other. In the sheet mode, the Nu was seen to depend on the thickness of the liquid waves traversing over the tube surface. Lower Nu values were recorded beneath the crest location of the liquid waves, which increases (1.4–11.6% depending on circumferential location) abruptly in magnitude at the advancing fronts of the waves. The temperature distribution in the liquid film in each of the modes was examined to evaluate the mechanism of heat transfer process. This study also compares the local heat transfer coefficient distribution with the analytical heat transfer models derived to predict heat transfer performance over horizontal tube surfaces.

42 ENGINEERING↗

DOE OSTI2022

Estimating the energy requirement for hydrogen production in proton exchange membrane electrolysis cells using rapid operando hydrogen crossover analysis

Hydrogen (H 2 ) crossover in proton exchange membrane water electrolyzers refers to the process by which hydrogen produced at the cathode traverses the membrane and mixes with the oxygen produced at the anode. This phenomenon reduces efficiency and may pose flammability hazards. In this work we present a method for quantifying the H 2 content of the anode exhaust gas using a gas chromatograph that is capable of sampling data every 2 min. Subsequent theory is presented to calculate the crossover flux, overall H 2 efficiency, and H 2 energy requirements. Results the effects of membrane thickness using Nafion TM N117 (178 um) and Nafion TM NR212 (51 um) membranes. Furthermore, it was found that thinner membranes lead to improved VI performance but exhibit higher crossover rates. Despite their increased crossover, leading to decreased hydrogen efficiency, the calculated required energy for NR212 membrane-electrode assemblies (MEAs) was significantly lower than that of N117 MEAs.

08 HYDROGEN↗

DOE OSTI2022

Dual-scale folding in cutting of commercially pure aluminum alloys

We examine a, hitherto, little-studied and curious machining chip morphology, with tell-tale signs of folding at two different length scales, that is common in cutting of certain ductile and highly strain-hardening metals like soft aluminum alloys, tantalum and niobium. This chip morphology does not appear in the usual catalogues of common chip types. The mechanics of formation of the “dual-scale folded chip” is studied in model material systems of commercially pure aluminum alloys (AA 1100 and AA 8040), that prominently exhibit this chip morphology. The flow, folding and associated plastic instabilities are investigated using micro/macro structure observations of the chip in a plane-strain cutting framework, with high-speed in situ imaging and image analysis of material flow; and force measurements. The smaller-scale folding is shown to develop in the primary deformation zone while the larger-scale folding occurs as the chip traverses the rake face of the tool. The resulting chip is composed of irregularly-spaced large folds, superimposed onto which are the quasi periodic small folds. The representative wavelengths of the two folds differ on average by an order of magnitude, 0.1 mm vs. 2 mm. The observations reveal a direct coupling between the material flow and chip morphology, and how specific attributes of the dual-scale folded chip arise from the flow mechanism. Plastic buckling is found to play a key role in the folding at both length scales. The small-scale folds are characteristic of a sinuous plastic flow mode, while the large-scale folding is characterized by buckling and stick-slip along the tool rake face, triggered by adhesive pinning of the chip to the tool. Important consequences of the dualscale folding are very large cutting forces, and force oscillations of large amplitude, despite the alloys being very soft, only ~ 25 HV. Here, the dual-scale folding is why many of these alloys are classified as “gummy” to machine. Since the dual-scale folded chip is associated with large cutting forces and poor surface quality, there is much to be gained by disrupting this flow type in practical machining applications. Methods for controlling the folding to improve machining performance with the gummy alloys are briefly discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Dynamic crushing of metal lattice metamaterials: Shock mode diagrams and transition to topology-independent compaction regime

Additively manufactured lattice metamaterials offer design versatility in strength and energy absorption and provide an additional degree of freedom through the selection of the lattice topology. Under quasistatic loading, the unit cell structure can strongly affect the stiffness, yield, and post-yield behavior, but whether and to what degree the effect of lattice topology persists into dynamic loading scenarios, up to the compaction shock regime, has not been established. LLNL’ s ALE3D hydrocode was used to perform a computational investigation of dynamic loading in multiple lattice types, including the gyroid, octet, Schwarz D, and rhombic dodecahedron, under impact velocities from 0.25 to 2.25 km/s. Shock Hugoniots for each lattice topology are generated and compared, suggesting that above a critical velocity, distinctions between architectures may not persevere and compacted lattices behave similarly. Here, to investigate the transition between topology-dependent quasistatic compression and the topology-independent regime above the critical velocity, a one-dimensional elastic-linear hardening plasticity-densified solid (E-LHP-DS) shock model for lattice materials was developed that relies upon confined compression to link the quasistatic and shock mechanics. Unlike similar works, the model does not assume rigid behavior prior to yield or locking behavior at densification, allowing a richer exploration of lattice mechanics. With only six parameters, the analytical model simultaneously fit quasistatic confined compression simulations for relative densities 0.1 $≤ \bar{ρ} ≤$ 0.9 and predicted dynamic compaction behavior to traverse several distinct shock modes, each defined by a critical impact speed (equivalently, critical stresses). Comparing the numerical results to the one-dimensional E-LHP-DS shock model predictions suggests that the topology-independence under strong shocks is linked to the onset of densification, which can be predicted based on quasistatic confined compression results.

Cellular material↗

DOE OSTI2024

Exploring constituent redistribution in irradiated U-19Pu-14Zr fuel via electron probe microanalysis

Here, the phenomena of constituent redistribution, wherein a previously homogeneous metallic fuel forms discrete, radially concentric compositional zones upon irradiation was investigated by examining an irradiated U-19Pu-14Zr fuel (where numbers represent wt. %) with a burnup of 11.5 at.% with electron probe microanalysis (EPMA) and quadruple inductively coupled plasma mass spectroscopy (Q-ICP-MS). EPMA-generated U, Pu, and Zr compositional data obtained from a diameter traverse of the sample was converted to mass and was used to: 1) compare the overall fuel element analysis results between the two methods, 2) determine the number of compositionally distinct zones forming as a result of constituent redistribution; and 3) quantify the post-irradiation loss or gain of U, Pu, and Zr atoms in each distinct compositional zone. Weight percent concentrations of U, Pu, and Zr for the overall cross section compare favorably between the two analytical methods, suggesting that the spatially resolved EPMA analysis complements bulk chemical analysis. Among the four identified compositional zones, post-irradiation quantification of U, Pu, and Zr elemental atom content changes shows that the quantity of U atoms lost from the innermost zone is slightly less than the quantity of U atoms gained by the middle two zones, and the quantity of Zr atoms lost from the high-U third zone is slightly less than is gained by the two innermost zones. Pu is lost from all four zones, although the innermost zone and the high-U third zone lose a significantly higher percentage (> 22 %) of their initial Pu atoms than the other two zones. For all three elements, EPMA cannot distinguish between atoms lost due to transport to a different zone from atoms lost due to nuclear processes; however, the insight gained from using this process can be used to experiment with new modeling techniques to predict constituent redistribution in U-Pu-Zr fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Microstructure and texture evolution of friction stir lap welded dissimilar multi-aluminum stack

Here this work investigates the microstructural evolution of a dissimilar multiple-aluminum alloy stack (AA7055, AA7055, AA6022) produced by friction stir lap welding (FSLW). A through-thickness analyses of grain orientation map, grain size and distribution, strain distribution, and texture evolution by electron backscattered diffraction-based characterization technique were performed. Two different welding parameters (combination of tool rotation and traverse speed) were used to study the effects of the welding parameters on these microstructural characteristics. The study shows that higher welding speed and higher rotation speed of the FSLW tool yield finer grains. FSLW resulted in shear texture formation in the weld regions due to high welding speeds. This texture development in the nugget zone is consistent with shear deformation, and (111) pole figures show ideal shear texture components. The Stir Zone (SZ) exhibits shear texture (A/A ̅ and B/B ̅ components), consistent with other shear-assisted processes. Recrystallization mechanisms observed are CDRX in SZ and DDRX in TMAZ (Thermo-Mechanically Affected Zone). Interface analysis reveals the presence of oxide layers between different aluminum sheets, affecting joint mechanical performance. The mechanical property of the lap welded joints was presented in terms of microhardness distribution across the weld cross section. Slightly higher hardness in higher speed welding setup was observed throughout the SZ and AA6022 layer, indicating greater grain refinement and uniform grain structure distribution.

36 MATERIALS SCIENCE↗

DOE OSTI2023

High-speed friction stir butt welding of 25.4 mm thick 7175-T79 aluminum alloy

This presents the first experimental demonstration of friction stir butt welding of 25.4 mm thick high-strength, precipitate-strengthened 7XXX aluminum alloy at a welding speed above 500 mm/min and penetration depth greater than 10 mm. Tool pin threads that terminate away from the shoulder region, helped mitigate nearby stress concentration during tool traversing, enabling a higher welding speed while avoiding pin fracture near shoulder. Two kinds of welds—one-sided with full penetration of the plate thickness and double-sided with partial penetration of the plate thickness—were conducted. Significant grain refinement (grain size 1.3 ± 0.8 µm) was evident at the junction of double-sided friction stir welding (FSW) compared to other regions of weld-nugget or base material. Cross-weld tensile testing with digital image correlation revealed high local strains at the junction of double-sided FSW before macro-yield. This strain concentration impeded the typical early necking on advancing/retreating sides of FSW joints, improving the yield strength of as-welded double-sided FSW by 20–24% over slow- to medium-speed one-sided FSW. In conclusion, the joint efficiency of as-welded, double-sided high-speed FSW was about 76% with respect to the base material’s ultimate tensile strength.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Friction surface layer deposition of triple-phase Al 10 Cr 12 Fe 35 Mn 23 Ni 20 high entropy alloy: Process optimization and microstructural evolution

A high-strength Co-free triple-phase Al 10 Cr 12 Fe 35 Mn 23 Ni 20 high-entropy alloy (HEA) was successfully fabricated using Friction Surface Layer Deposition (FSLD), a bulk manufacturing method. Multiple single-layer deposits were produced by varying forging force (F) and traverse speeds (T r ) to optimize the process parameters. The optimized conditions (F = 40 kN & T r = 200 mm/min) were then applied to manufacture a scaled-up multi-layer specimen. The initial microstructure of the HEA consisted of coarse grains of the soft FCC-phase, long columnar dendrites of the hard BCC-phase, and small precipitates of the harder B2-phase within the BCC-dendrites. During FSLD, the FCC-matrix underwent continuous dynamic recrystallization due to high-temperature severe plastic deformation, forming finer equiaxed grains. Simultaneously, the BCC-dendrites fractured into smaller fragments, some of which experienced partial growth and coarsening under applied stress, resulting in an hourglass morphology. In contrast, the small B2-precipitates within the BCC-fragments dissolved during the elevated temperatures of FSLD and reprecipitated as substantially finer precipitates during continuous cooling post-FSLD. Additionally, the orientation relationships between the FCC and BCC/B2 phases were completely destroyed by the severe thermoplastic deformation inherent to FSLD. The microstructural refinements led to a substantial improvement in hardness from 177 HV to 283 HV, driven by Hall-Petch strengthening. The increased number of interfaces, including coherent BCC-B2 interfaces, potentially enhances the sink strength and radiation tolerance of the HEA, making it a promising candidate for nuclear applications. In conclusion, this study also highlights FSLD as a versatile technique for achieving tunable properties in HEAs, with detailed schematics illustrating the complex mechanisms of phase transformations during processing.

Additive Manufacturing↗

DOE OSTI2020

Horizons of modern molecular dynamics simulation in digitalized solid freeform fabrication with advanced materials

Our ability to shape and finish a component by combined methods of fabrication including (but not limited to) subtractive, additive and/or no theoretical mass-loss/addition during the fabrication is now popularly known as solid free form fabrication. Fabrication of a telescope mirror is a typical example where grinding and polishing processes are first applied to shape the mirror and thereafter an optical coating is usually applied to enhance its optical performance. The area of nanomanufacturing cannot grow without a deep knowledge of the fundamentals of materials and consequently, the use of computer simulations is becoming ubiquitous. This article is intended to introduce the most recent advances in the computation benefit specific to the area of solid free form fabrication as these systems are traversing through the journey of digitalisation and Industry-4.0. Specifically, this article demonstrates that the application of the latest materials modelling approaches, based on techniques such as molecular dynamics, are enabling breakthroughs in applied precision manufacturing techniques.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Differences in the response of two light guide technologies and two readout technologies after an exchange of liquid argon in the dewar

In this investigation the response to the scintillation light generated by through-going cosmic muons in liquid argon (LAr) was measured by two light guide technologies and two readout technologies after five weeks of running in the TallBo dewar at Fermilab. The response was remeasured after the dewar was drained of LAr, refilled, and then run again for an additional four weeks. After the dewar was refilled, there was clear evidence that the scintillation signal had dropped significantly. The two light guide technologies were developed at Indiana University and MIT/Fermilab. The two readout technologies were boards that passively or actively ganged 12 Hamamatsu MPPCs. Two possible explanations were identified for the degraded signal: the response of the two light guide technologies degraded due to damage caused by thermal cycling, and/or unknown differences in the trace residual Xe contamination in the fills of LAr led to the observed drop in scintillation light. Neither absorption nor quenching by N 2 , O 2 , and H 2 O contamination can account for the degradation. Neither the individual Hamamatsu MPPCs nor the passive/active ganging boards appear to have been affected by the thermal cycling. The path length distributions of the cosmics traversing the dewar appear quite similar in both event samples.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

STOPGAP—A time-of-flight extension for the Belle II TOP barrel PID system

The Belle II barrel region is instrumented with the Time of Propagation (TOP) particle identification system. Due to its mechanical design, the individual TOP modules do not overlap, leaving a gap of around 2 cm between them in the azimuthal direction. This leads to a 6%–9% drop in acceptance, depending on the track’s momentum. We propose a solution to remedy these gaps by instrumenting them with fast silicon detectors to directly measure the time-of-flight of traversing particles. We present here a simulation study discussing the performance requirements and the possible sensor technologies, and we demonstrate that such a project could be realized with novel, fast monolithic CMOS sensors, or alternatively AC-LGAD sensors, both of which are expected to reach MIP timing resolutions of down to 50 ps or better.

Hartbrich, Oskar↗

DOE OSTI2021

The Problem of overlapping formation times: In-medium virtual corrections for QCD

High energy particles traversing through medium primarily lose energy by showering through hard bremsstrahlung and pair production. These splitting processes are coherent over large distances in the very high energy limit, leading to suppression from the Landau-Pomeranchuk-Migdal (LPM) effect. Avoiding soft-emission approximations, we study the cases where the coherence lengths of two consecutive splittings overlap (which is important for calculating corrections to LPM effect in QCD). Here, we will show how to compute in-medium virtual corrections to the leading order LPM emission rates for QCD. These loop corrections will be necessary for calculating properties of in-medium high energy parton showers. To simplify the calculations, we will focus on the all-gluon case and work in the large-Nc limit of QCD, where Nc is the number of colors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Modification of charged-particle jets in event-shape engineered Pb–Pb collisions at $\sqrt{S_{NN}}$ = 5.02 TeV

Charged-particle jet yields have been measured in semicentral Pb–Pb collisions at center-of-mass energy per nucleon–nucleon collision $\sqrt{S_{NN}}$ = 5.02 TeV with the ALICE detector at the LHC. These yields are reported as a function of the jet transverse momentum, and further classified by their angle with respect to the event plane and the event shape, characterized by ellipticity, in an effort to study the path-length dependence of jet quenching. Jets were reconstructed at midrapidity from charged-particle tracks using the anti-k T algorithm with resolution parameters R = 0.2 and 0.4, with event-plane angle and event-shape values determined using information from forward scintillating detectors. The results presented in this letter show that, in semicentral Pb–Pb collisions, there is no significant difference between jet yields in predominantly isotropic and elliptical events. However, out-of-plane jets are observed to be more suppressed than in-plane jets. Further, this relative suppression is greater for low transverse momentum (< 50 GeV/c) R = 0.2 jets produced in elliptical events, with out-of-plane to in-plane jet-yield ratios varying up to 5.2σ between different event-shape classes. These results agree with previous studies indicating that jets experience azimuthally anisotropic suppression when traversing the QGP medium, and can provide additional constraints on the path-length dependence of jet energy loss.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Exploration of mass splitting and muon/tau mixing parameters for an eV-scale sterile neutrino with IceCube

We present the first three-parameter fit to a 3+1 sterile neutrino model using 7.634 years of data from the IceCube Neutrino Observatory on $ν$ μ + $\overline{ν}$ μ charged-current interactions in the energy range 500–9976 GeV. Our analysis is sensitive to the mass-squared splitting between the heaviest and lightest mass state (Δm$^{2}_{41}$), the mixing matrix element connecting muon flavor to the fourth mass state (|U μ4 | 2 ), and the element connecting tau flavor to the fourth mass state (|U $\tau$4 | 2 ). Predicted propagation effects in matter enhance the signature through a resonance as atmospheric neutrinos from the Northern Hemisphere traverse the Earth to the IceCube detector at the South Pole. The remaining sterile neutrino matrix elements are left fixed, with |U e4 | 2 = 0 and δ 14 = 0, as they have a negligible effect, and δ 24 = π is set to give the most conservative limits. The result is consistent with the no-sterile neutrino hypothesis with a probability of 4.3%. Profiling the likelihood of each parameter yields the 90% confidence levels: 2.4eV 2 < Δm$^{2}_{41}$) < 9.6eV 2 , 0.0081< |U μ4 | 2 < 0.10, and |U $\tau$4 | 2 < 0.035, which narrows the allowed parameter-space for |U $\tau$4 | 2 . However, the primary result of this analysis is the first map of the 3+1 parameter space exploring the interdependence of Δm$^{2}_{41}$, |U μ4 | 2 , and |U $\tau$4 | 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Structure of the Extracellular Region of the Bacterial Type VIIb Secretion System Subunit EsaA

Gram-positive bacteria use type VII secretion systems (T7SSs) to export effector proteins that manipulate the physiology of nearby prokaryotic and eukaryotic cells. Several mycobacterial T7SSs have established roles in virulence. By contrast, the genetically distinct T7SSb pathway found in Firmicutes bacteria more often functions to mediate bacterial competition. A lack of structural information on the T7SSb has limited the understanding of effector export by this protein secretion apparatus. Here, we present the 2.4 angstrom crystal structure of the extracellular region of the T7SSb subunit EsaA from Streptococcus gallolyticus. Our structure reveals that homodimeric EsaA is an elongated, arrow-shaped protein with a surface-accessible "tip", which in some species of bacteria serves as a receptor for lytic bacteriophages. Because it is the only T7SSb subunit large enough to traverse the peptidoglycan layer of Firmicutes, we propose that EsaA plays a critical role in transporting effectors across the entirety of the Gram-positive cell envelope.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Toward data assimilation of ship-induced aerosol–cloud interactions

Satellite imagery can detect temporary cloud trails or ship tracks formed from aerosols emitted from large ships traversing our oceans, a phenomenon that global climate models cannot directly reproduce. Ship tracks are observable examples of marine cloud brightening, a potential solar climate intervention that shows promise in helping combat climate change. In this paper, we demonstrate a simulation-based approach in learning the behavior of ship tracks based upon a novel stochastic emulation mechanism. Our method uses wind fields to determine the movement of aerosol–cloud tracks and uses a stochastic partial differential equation (SPDE) to model their persistence behavior. This SPDE incorporates both a drift and diffusion term which describes the movement of aerosol particles via wind and their diffusivity through the atmosphere, respectively. We first present our proposed approach with examples using simulated wind fields and ship paths. We then successfully demonstrate our tool by applying the approximate Bayesian computation method-sequential Monte Carlo for data assimilation.

54 ENVIRONMENTAL SCIENCES↗