Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble prediction”

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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Amplified drought trends in Nepal increase the potential for Himalayan wildfires

Here we report in the spring 2021, Nepal underwent a record wildfire season in which active fires were detected at a rate 10 times greater than the 2002–2020 average. Prior to these major wildfire events, the country experienced a prolonged precipitation deficit and extreme drought during the post-monsoon period (starting in October 2020). An analysis using observational, reanalysis, and climate model ensemble data indicates that both climate variability and climate change-induced severe drought conditions were at play. Further analysis of climate model outputs suggests the likely reoccurrence of drought conditions, thus favoring active wildfire seasons in Nepal throughout the twenty-first century. While the inter-model uncertainty is large and direct modeling of wildfire spread and suppression has not been completed, the demonstrated relationship between a drought index (the standardized precipitation and evapotranspiration index) and subsequent fire activity may offer actionable opportunities for forest managers to employ the monitoring and projection of climate anomalies at sub-seasonal to decadal timescales to inform their management strategies for Nepal’s wildlands.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Phenological shifts in lake stratification under climate change

One of the most important physical characteristics driving lifecycle events in lakes is stratification. Already subtle variations in the timing of stratification onset and break-up (phenology) are known to have major ecological effects, mainly by determining the availability of light, nutrients, carbon and oxygen to organisms. Despite its ecological importance, historic and future global changes in stratification phenology are unknown. Here, we used a lake-climate model ensemble and long-term observational data, to investigate changes in lake stratification phenology across the Northern Hemisphere from 1901 to 2099. Under the high-greenhouse-gas-emission scenario, stratification will begin 22.0 ± 7.0 days earlier and end 11.3 ± 4.7 days later by the end of this century. It is very likely that this 33.3 ± 11.7 day prolongation in stratification will accelerate lake deoxygenation with subsequent effects on nutrient mineralization and phosphorus release from lake sediments. Further misalignment of lifecycle events, with possible irreversible changes for lake ecosystems, is also likely.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Hydrodynamic computations of high-power laser drives generating metal ejecta jets from surface grooves

Understanding dynamic fragmentation in shock-loaded metals and predicting properties of the resulting ejecta are of considerable importance for both basic and applied science. The nature of material ejection has been shown to change drastically when the free surface melts on compression or release. In this work, we present hydrodynamic simulations of laser-driven microjetting from micron-scale grooves on a tin surface. We study microjet formation across a range of shock strengths from drives that leave the target solid after release to drives that induce shock melting in the target. The shock-state particle velocity (Up) varies from 0.3 to 3 km/s and the shock breakout pressure is 3–120 GPa. The microjet tip velocity is 1–8 km/s and the free-surface velocity varies from 0.1 to 5 km/s. Two tin equations of state are examined: a “soft” model (LEOS 501) where the target melts for U p > 1 km/s and a more detailed multiphase model (SESAME 2161) that melts for U p > 1.4 km/s. We use these two models to examine the influence of phase change and the choice of the material model on microjet formation and evolution. We observe in our computational results that jet formation can be classified into three regimes: a low-energy regime where material strength affects jet formation, a moderate-energy regime dominated by the changing phase of tin material, and a high-energy regime where results are insensitive to the material model and jet formation is described by an idealized steady-jet theory. Using an ensemble of 2D simulations, we show that these trends hold across a wide range of drive energies and groove angles.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mid-century climate change impacts on tornado-producing tropical cyclones

Tornadoes are a co-occurring extreme that can be produced by landfalling tropical cyclones (TCs). These tornadoes can exacerbate the loss of life and property damage caused by the TC from which they were spawned. It is uncertain how the severe weather environments of landfalling TCs may change in a future climate and how this could impact tornado activity from TCs. In this study, we investigated four TCs that made landfall in the U.S. and produced large tornado outbreaks. We performed four-member ensembles of convective-allowing (4-km resolution) regional climate model simulations representing each TC in the historical climate and a mid-twenty-first century future climate. To identify potentially tornadic storms, or TC-tornado (TCT) surrogates, we used thresholds for three-hourly maximum updraft helicity and radar reflectivity, as tornadoes are not resolved in the model. We found that the ensemble-mean number of TCT-surrogates increased substantially (56–299%) in the future, supported by increases in most-unstable convective available potential energy, surface-to-700-hPa bulk wind shear, and 0–1-km storm-relative helicity in the tornado-producing region of the TCs. On the other hand, future changes in most-unstable convective inhibition had minimal influence on future TCT-surrogates. This provides robust evidence that tornado activity from TCs may increase in the future. Furthermore, TCT-surrogate frequency between 00Z and 09Z increased for three of the four cases, suggesting enhanced tornado activity at night, when people are asleep and more likely to miss warnings. All of these factors indicate that TC-tornadoes may become more frequent and a greater hazard in the future, compounding impacts from future increases in TC winds and precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Semileptonic form factors for $$B\rightarrow D^*\ell \nu $$ at nonzero recoil from $$2+1$$-flavor lattice QCD: Fermilab Lattice and MILC Collaborations

Abstract We present the first unquenched lattice-QCD calculation of the form factors for the decay $$B\rightarrow D^*\ell \nu $$ B → D ∗ ℓ ν at nonzero recoil. Our analysis includes 15 MILC ensembles with $$N_f=2+1$$ N f = 2 + 1 flavors of asqtad sea quarks, with a strange quark mass close to its physical mass. The lattice spacings range from $$a\approx 0.15$$ a ≈ 0.15 fm down to 0.045 fm, while the ratio between the light- and the strange-quark masses ranges from 0.05 to 0.4. The valence b and c quarks are treated using the Wilson-clover action with the Fermilab interpretation, whereas the light sector employs asqtad staggered fermions. We extrapolate our results to the physical point in the continuum limit using rooted staggered heavy-light meson chiral perturbation theory. Then we apply a model-independent parametrization to extend the form factors to the full kinematic range. With this parametrization we perform a joint lattice-QCD/experiment fit using several experimental datasets to determine the CKM matrix element $$|V_{cb}|$$ | V cb | . We obtain $$\left| V_{cb}\right| = (38.40 \pm 0.68_{\text {th}} \pm 0.34_{\text {exp}} \pm 0.18_{\text {EM}})\times 10^{-3}$$ V cb = ( 38.40 ± 0 . 68 th ± 0 . 34 exp ± 0 . 18 EM ) × 10 - 3 . The first error is theoretical, the second comes from experiment and the last one includes electromagnetic and electroweak uncertainties, with an overall $$\chi ^2\text {/dof} = 126/84$$ χ 2 /dof = 126 / 84 , which illustrates the tensions between the experimental data sets, and between theory and experiment. This result is in agreement with previous exclusive determinations, but the tension with the inclusive determination remains. Finally, we integrate the differential decay rate obtained solely from lattice data to predict $$R(D^*) = 0.265 \pm 0.013$$ R ( D ∗ ) = 0.265 ± 0.013 , which confirms the current tension between theory and experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A multi-dimensional parametric study of variability in multi-phase flow dynamics during geologic CO 2 sequestration accelerated with machine learning

Successful geologic CO 2 storage projects depend on numerical simulations to predict reservoir performance during site selection, injection verification, and post-injection monitoring phases of the project. These numerical simulations solve non-linear sets of coupled partial differential equations, while accounting for multi-phase fluid dynamics on the basis of constitutive equations that are embedded into the solution scheme. As a consequence, individual simulations often require tens to hundreds of hours to complete on high-performance computing clusters. Moreover, laboratory experiments reveal that parametric functions for capillary pressure and relative permeability exhibit substantial variability, even within the same rock type. This combination of computational expense and wide-ranging parametric variability means that there remains substantial uncertainty in the behavior of multi-phase CO 2 -water systems, particularly in the context of feedbacks between relative permeability and capillary pressure. To bridge this knowledge gap, here we develop a novel workflow that utilizes physics-based numerical simulation to train an artificial neural network (ANN) emulator for interrogating the multivariate parameter space that governs both capillary pressure and relative permeability. With this approach, the ANN is trained to emulate both fluid pressure distribution and CO 2 saturation, which are then interrogated quantitatively to generate parametric response surface mappings with high-fidelity resolution. Results from this study initially show that capillary entry pressure is the dominant control on both CO 2 plume geometry and fluid pressure propagation when considering the combined effects of capillary pressure and relative permeability, particularly when phase interference is low and residual CO 2 saturation is high. Moreover, the ANN emulator provides tremendous computational speed-up by computing 2691 individual simulations in several minutes; whereas, the same simulation ensemble would have required ~3 years of simulation time using only physics-based simulation methods (25,000 times speed up).

58 GEOSCIENCES↗
DOE OSTI2025

High-Resolution Model Intercomparison Project phase 2 (HighResMIP2) towards CMIP7

Abstract. Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. Consequently, the role of enhanced horizontal resolution in improved process representation in all components of the climate system continues to be of great interest. Recent simulations suggest the possibility of significant changes in both large-scale aspects of the ocean and atmospheric circulations and in the regional responses to climate change, as well as improvements in representations of small-scale processes and extremes, when resolution is enhanced. The first phase of the High-Resolution Model Intercomparison Project (HighResMIP1) was successful at producing a baseline multi-model assessment of global simulations with model grid spacings of 25–50 km in the atmosphere and 10–25 km in the ocean, a significant increase when compared to models with standard resolutions on the order of 1° that are typically used as part of the Coupled Model Intercomparison Project (CMIP) experiments. In addition to over 250 peer-reviewed manuscripts using the published HighResMIP1 datasets, the results were widely cited in the Intergovernmental Panel on Climate Change report and were the basis of a variety of derived datasets, including tracked cyclones (both tropical and extratropical), river discharge, storm surge, and impact studies. There were also suggestions from the few ocean eddy-rich coupled simulations that aspects of climate variability and change might be significantly influenced by improved process representation in such models. The compromises that HighResMIP1 made should now be revisited, given the recent major advances in modelling and computing resources. Aspects that will be reconsidered include experimental design and simulation length, complexity, and resolution. In addition, larger ensemble sizes and a wider range of future scenarios would enhance the applicability of HighResMIP. Therefore, we propose the High-Resolution Model Intercomparison Project phase 2 (HighResMIP2) to improve and extend the previous work, to address new science questions, and to further advance our understanding of the role of horizontal resolution (and hence process representation) in state-of-the-art climate simulations. With further increases in high-performance computing resources and modelling advances, along with the ability to take full advantage of these computational resources, an enhanced investigation of the drivers and consequences of variability and change in both large- and synoptic-scale weather and climate is now possible. With the arrival of global cloud-resolving models (currently run for relatively short timescales), there is also an opportunity to improve links between such models and more traditional CMIP models, with HighResMIP providing a bridge to link understanding between these domains. HighResMIP also aims to link to other CMIP projects and international efforts such as the World Climate Research Program lighthouse activities and various digital twin initiatives. It also has the potential to be used as training and validation data for the fast-evolving machine learning climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Determination of confinement regime boundaries via separatrix parameters on Alcator C-Mod based on a model for interchange-drift-Alfvén turbulence

The separatrix operational space (SepOS) model (Eich and Manz 2021 Nucl. Fusion 61 086017) is shown to predict the L–H transition, the L-mode density limit, and the ideal magnetohydrodynamic ballooning limit in terms of separatrix parameters for a wide range of Alcator C-Mod plasmas. The model is tested using Thomson scattering measurements across a wide range of operating conditions on C-Mod, spanning $\overline{n}_{\mathrm{e}}$ = 0.3–5.5 $\times 10^{20}\,$m −3 , $B_{\mathrm{t}} = 2.5$–8.0 T, and $B_{\mathrm{p}}$= 0.1–1.2 T. An empirical regression for the electron pressure gradient scale length, $\lambda_{{p}_{\mathrm{e}}}$, against a turbulence control parameter, $\alpha_{\mathrm{t}}$, and the poloidal fluid gyroradius, $\rho_{\mathrm{s,p}}$, is constructed for H-modes and found to require positive exponents for both regression parameters, indicating turbulence widening of near-scrape-off layer widths at high $\alpha_{\mathrm{t}}$ and an inverse scaling with $B_{\mathrm{p}}$, consistent with results on ASDEX Upgrade. The SepOS model is also tested in the unfavorable drift direction and found to apply well to all three boundaries, including the L–H transition as long as a correction to the Reynolds energy transfer term, $\alpha_\mathrm{RS} \lt 1$ is applied. I-modes typically exist in the unfavorable drift direction for values of $\alpha_{\mathrm{t}} \lesssim 0.3$. Finally, an experiment studying the transition between the Type-I ELMy and EDA H-mode is analyzed using the same framework. It is found that a recently identified boundary $\alpha_{\mathrm{t}} = 0.55$ at excludes most EDA H-modes but that the balance of wavenumbers responsible for the L-mode density limit, namely $k_\mathrm{EM} = k_\mathrm{RBM}$, may better describe the transition on C-Mod. The ensemble of boundaries validated and explored is then applied to project regime access and limit avoidance for the SPARC primary reference discharge parameters.

ELM suppression↗
NASA NTRS2016

The Mexican Drought Atlas: Tree-Ring Reconstructions of the Soil Moisture Balance During the Late Pre-Hispanic, Colonial, and Modern Eras

Mexico has suffered a long history and prehistory of severe sustained drought. Drought over Mexico is modulated by ocean-atmospheric variability in the Atlantic and Pacific, raising the possibility for long-range seasonal climate forecasting, which could help mediate the economic and social impacts of future dry spells. The instrumental record of Mexican climate is very limited before 1920, but tree-ring chronologies developed from old-growth forests in Mexico can provide an excellent proxy representation of the spatial pattern and intensity of past moisture regimes useful for the analysis of climate dynamics and climate impacts. The Mexican Drought Atlas (MXDA) has been developed from an extensive network of 252 climate sensitive tree-ring chronologies in and near Mexico. The MXDA reconstructions extend from 1400 CE-2012 and were calibrated with the instrumental summer (JJA) self-calibrating Palmer Drought Severity Index (scPDSI) on a 0.5deg latitude/longitude grid extending over land areas from 14 to 34degN and 75-120degW using Ensemble Point-by-Point Regression (EPPR) for the 1944-1984 period. The grid point reconstructions were validated for the period 1920-1943 against instrumental gridded scPDSI values based on the fewer weather station observations available during that interval. The MXDA provides a new spatial perspective on the historical impacts of moisture extremes over Mexico during the past 600-years, including the Aztec Drought of One Rabbit in 1454, the drought of El Ano de Hambre in 1785-1786, and the drought that preceded the Mexican Revolution of 1909-1910. The El Nino/Southern Oscillation (ENSO) is the most important ocean-atmospheric forcing of moisture variability detected with the MXDA. In fact, the reconstructions suggest that the strongest central equatorial Pacific sea surface temperature (SST) teleconnection to the soil moisture balance over North America may reside in northern Mexico. This ENSO signal has stronger and more time-stable correlations than computed for either the Atlantic Multidecadal Oscillation or Pacific Decadal Oscillation. The extended Multivariate ENSO Index is most highly correlated with reconstructed scPDSI over northern Mexico, where warm events favor moist conditions during the winter, spring, and early summer. This ENSO teleconnection to northern Mexico has been strong over the past 150 years, but it has been comparatively weak and non-stationary in the MXDA over central and southern Mexico where eastern tropical Pacific and Caribbean/tropical Atlantic SSTs seem to be more important. The ENSO teleconnection to northern Mexico is weaker in the available instrumental PDSI, but analyses based on the millennium climate simulations with the Community Earth System Model suggest that the moisture balance during the winter, spring, and early summer over northern Mexico may indeed be particularly sensitive to ENSO forcing. Nationwide drought is predicted to become more common with anthropogenic climate change, but the MXDA reconstructions indicate that intense "All Mexico" droughts have been rare over the past 600 years and their frequency does not appear to have increased substantially in recent decades.

El Nino↗
DOE OSTI2024

A comparison of model validation approaches for echo state networks using climate model replicates

As global temperatures continue to rise, climate mitigation strategies such as stratospheric aerosol injections (SAI) are increasingly discussed, but the downstream effects of these strategies are not well understood. As such, there is interest in developing statistical methods to quantify the evolution of climate variable relationships during the time period surrounding an SAI. Feature importance applied to echo state network (ESN) models has been proposed as a way to understand the effects of SAI using a data-driven model. This approach depends on the ESN fitting the data well. If not, the feature importance may place importance on features that are not representative of the underlying relationships. Typically, time series prediction models such as ESNs are assessed using out-of-sample performance metrics that divide the times series into separate training and testing sets. However, this model assessment approach is geared towards forecasting applications and not scenarios such as the motivating SAI example where the objective is using a data driven model to capture variable relationships. Here, in this paper, we demonstrate a novel use of climate model replicates to investigate the applicability of the commonly used repeated hold-out model assessment approach for the SAI application. Simulations of an SAI are generated using a simplified climate model, and different initialization conditions are used to provide independent training and testing sets containing the same SAI event. The climate model replicates enable out-of-sample measures of model performance, which are compared to the single time series hold-out validation approach. For our case study, it is found that the repeated hold-out sample performance is comparable, but conservative, to the replicate out-of-sample performance when the training set contains enough time after the aerosol injection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

flat10MIP: an emissions-driven experiment to diagnose the climate response to positive, zero and negative CO2 emissions

Abstract. The proportionality between global mean temperature and cumulative emissions of CO2 predicted in Earth system models (ESMs) is the foundation of carbon budgeting frameworks. Deviations from this behavior could impact estimates of required net-zero timings and negative emissions requirements to meet the Paris Agreement climate targets. However, existing ESM diagnostic experiments do not allow for direct estimation of these deviations as a function of defined emissions pathways. Here, we perform a set of climate model diagnostic experiments for the assessment of transient climate response to cumulative CO2 emissions (TCRE), the Zero Emissions Commitment (ZEC), and climate reversibility metrics in an emissions-driven framework. The emissions-driven experiments provide consistent independent variables simplifying simulation, analysis and interpretation, with emissions rates more comparable to recent levels than existing protocols using model-specific compatible emissions from the CMIP DECK 1pctCO2 experiment, where emissions rates tend to increase during the experiment, such that at the time of CO2 doubling in year 70, emissions are much greater than present-day values. A base experiment, “esm-flat10”, has constant emissions of CO2 of 10 GtC per year (near-present-day values), and initial results show that the TCRE estimated in this experiment is about 0.1 K less than that obtained using 1pctCO2. A subset of ESMs exhibit land carbon sinks that saturate during this experiment. A branch experiment, esm-flat10-zec, illustrates that both positive and negative ZEC effects are less pronounced under esm-flat10 than under 1pctCO2 – the magnitude of ZEC50 in ESMs is, on average, reduced by 30 % compared with 1pctCO2 branch experiments. A final experiment, esm-flat10-cdr, assesses climate reversibility under negative emissions, where we find that peak warming may occur before or after net zero and that the asymmetry in temperature at a given level of cumulative emissions between the positive and negative emissions phases is well described by ZEC in most models. Further, we find that existing probabilistic simple climate model (SCM) ensembles tend to overestimate temperature reversibility compared with ESMs, highlighting the need for additional constraints. We propose a set of climate diagnostic indicators to quantify various aspects of climate reversibility. These experiments were suggested as potential candidates in CMIP7 and have since been adopted as “fast track” simulations.

Sanderson, Benjamin M↗
DOE OSTI2023

Improving Tropical Cyclogenesis Forecasts of Hurricane Irma (2017) through the Assimilation of All-Sky Infrared Brightness Temperatures

The assimilation of satellite all-sky infrared (IR) brightness temperatures (BTs) has been shown in previous studies to improve intensity forecasts of tropical cyclones. Here, in this study, we examine whether assimilating all-sky IR BTs can also potentially improve tropical cyclogenesis forecasts by improving the pregenesis cloud and moisture fields. By using an ensemble-based data assimilation system, we show that the assimilation of upper-tropospheric water vapor channel BTs observed by the Meteosat-10 SEVIRI instrument two days before the formation of a tropical depression improves the genesis forecast of Hurricane Irma (2017), a classic Cape Verde storm, by up to 24 h while also capturing its later rapid intensification in deterministic forecasts. In an experiment that withholds the assimilation of all-sky IR BTs, the assimilation of conventional observations from the Global Telecommunications System (GTS) leads to the premature genesis of Hurricane Irma by at least 24 h. This premature genesis is shown to result from an overestimation of the spatial coverage of deep convection within the African easterly wave (AEW) from which Irma eventually forms. The gross overestimation of deep convection without all-sky IR BTs is accompanied by higher column saturation fraction, stronger low-level convergence, and the earlier spinup of a low-level meso- β -scale vortex within the AEW that ultimately becomes Hurricane Irma. Through its adjustment to the initial moisture and cloud conditions, the assimilation of all-sky IR BTs leads to a more realistic convective evolution in forecasts and ultimately a more realistic timing of genesis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2004

Starting Vortex Identified as Key to Unsteady Ejector Performance

Unsteady ejectors are currently under investigation for use in some pulse-detonation-engine-based propulsion systems. Experimental measurements made in the past, and recently at the NASA Glenn Research Center, have demonstrated that thrust augmentation can be enhanced considerably when the driver is unsteady. In ejector systems, thrust augmentation is defined as = T(sup Total)/T(sup j), where T(sup Total) is the total thrust of the combined ejector and driving jet and T(sup j) is the thrust due to the driving jet alone. There are three images in this figure, one for each of the named thrust sources. The images are color contours of measured instantaneous vorticity. Each image is an ensemble average of at least 150 phase-locked measurements. The flow is from right to left, and the shape and location of each driver is shown on the far right of each image. The emitted vortex is a clearly defined "doughnut" of highly vortical (spinning) flow. In these planar images, the vortex appears as two distorted circles, one above, and one below the axis of symmetry. Because they are spinning in the opposite direction, the two circles have vorticity of opposite sign and thus are different colors. There is also a rectangle shown in each image. Its width represents the ejector diameter that was found experimentally to yield the highest thrust augmentation. It is apparent that the optimal ejector diameter is that which just "captures" the vortex: that is, the diameter bounding the outermost edge of the vortex structure. The exact mechanism behind the enhanced performance is unclear; however, it is believed to be related to the powerful vortex emitted with each pulse of the unsteady driver. As such, particle imaging velocimetry (PIV) measurements were obtained for three unsteady drivers: a pulsejet, a resonance tube, and a speaker-driven jet. All the drivers were tested with ejectors, and all exhibited performance enhancement over similarly sized steady drivers. The characteristic starting vortices of each driver are shown in these images. The images are color contours of measured instantaneous vorticity. Each image is an ensemble average of at least 150 phase-locked measurements. The flow is from right to left. The shape and location of each driver is shown on the far right of each image. The rectangle shown in each image represents the ejector diameter that was found experimentally to yield the highest thrust augmentation. It is apparent that the optimal ejector diameter is that which just "captures" the vortex: that is, the diameter bounding the outermost edge of the vortex structure. Although not shown, it was observed that the emitted vortex spread as it traveled downstream. The spreading rate for the pulsejet is shown as the dashed lines in the top image. A tapered ejector was fabricated that matched this shape. When tested, the ejector demonstrated superior performance to all those previously tested at Glenn (which were essentially of straight, cylindrical form), achieving a remarkable thrust augmentation of 2. The measured thrust augmentation is shown as a function of ejector length. Also shown are the thrust augmentation values achieved with the straight, cylindrical ejectors of varying diameters. Here, thrust augmentation is plotted as a function of ejector length for several families of ejector diameters. It can be seen that large thrust augmentation values are indeed obtained and that they are sensitive to both ejector length and diameter, particularly the latter. Five curves are shown. Four correspond to straight ejector diameters of 2.2, 3.0, 4.0, and 6.0 in. The fifth curve corresponds to the tapered ejector contoured to bound the emitted vortex. For each curve, there are several data points corresponding to different lengths. The largest value of thrust augmentation is 2.0 for the tapered ejector and 1.81 for the straight ejectors. Regardless of their diameters, all the ejectors trend toward peak performance at a particular leng. That the cross-sectional dimensions of optimal ejectors scaled precisely with the vortex dimensions on three separate pulsed thrust sources demonstrates that the action of the vortex is responsible for the enhanced ejector performance. The result also suggests that, in the absence of a complete understanding of the entrainment and augmentation mechanisms, methods of characterizing starting vortices may be useful for correlating and predicting unsteady ejector performance.

Paxson, Daniel E.↗
NASA NTRS1994

A comparison between general circulation model simulations using two sea surface temperature datasets for January 1979

Simulations with the UCLA atmospheric general circulation model (AGCM) using two different global sea surface temperature (SST) datasets for January 1979 are compared. One of these datasets is based on Comprehensive Ocean-Atmosphere Data Set (COADS) (SSTs) at locations where there are ship reports, and climatology elsewhere; the other is derived from measurements by instruments onboard NOAA satellites. In the former dataset (COADS SST), data are concentrated along shipping routes in the Northern Hemisphere; in the latter dataset High Resolution Infrared Sounder (HIRS SST), data cover the global domain. Ensembles of five 30-day mean fields are obtained from integrations performed in the perpetual-January mode. The results are presented as anomalies, that is, departures of each ensemble mean from that produced in a control simulation with climatological SSTs. Large differences are found between the anomalies obtained using COADS and HIRS SSTs, even in the Northern Hemisphere where the datasets are most similar to each other. The internal variability of the circulation in the control simulation and the simulated atmospheric response to anomalous forcings appear to be linked in that the pattern of geopotential height anomalies obtained using COADS SSTs resembles the first empirical orthogonal function (EOF 1) in the control simulation. The corresponding pattern obtained using HIRS SSTs is substantially different and somewhat resembles EOF 2 in the sector from central North America to central Asia. To gain insight into the reasons for these results, three additional simulations are carried out with SST anomalies confined to regions where COADS SSTs are substantially warmer than HIRS SSTs. The regions correspond to warm pools in the northwest and northeast Pacific, and the northwest Atlantic. These warm pools tend to produce positive geopotential height anomalies in the northeastern part of the corresponding oceans. Both warm pools in the Pacific produce large-scale circulation anomalies with a pattern that resembles that obtained using COADS SSTs as well as EOF 1 of the control simulation; the warm pool in the Atlantic does not. These results suggest that the differences obtained with COADS SSTs and HIRS SSTs are mostly due to the differences in the datasets over the northern Pacific. There was a blocking episode near Greenland in late January 1979. Both simulations with warm SST anomalies over the northwest and northeast Pacific show a tendency toward increased incidence of North Atlantic blocking; the simulation with warm SST anomalies over the northwest Atlantic shows a tendency toward decreased incidence. These results suggest that features in both SST datasets that do not have a counterpart in the other dataset contribute signficantly to the differences between the simulated and observed fields. The results of this study imply that uncertainties in current SST distributions for the world oceans can be as important as the SST anomalies themselves in terms of their impact on the atmospheric circulation. Caution should be exercised, therefore, when linking anomalous circulation and SST patterns, especially in long-range prediction.

Ose, Tomoaki↗
DOE OSTI2025

The Influence of African Easterly Waves on Atlantic Tropical Cyclone Tracks and Landfall in Large Ensembles

African easterly waves (AEWs) are an important precursor or “seed” for Atlantic tropical cyclones (TCs), with 60%–80% of major hurricanes observed to originate from AEWs. However, climate model simulations indicate that AEWs are not necessary to maintain annual Atlantic TC frequency. Furthermore, small ensembles suggest that AEWs may impact the spatial distribution and landfall of Atlantic TCs. Here, we investigated the influence of AEWs on the spatial distribution of Atlantic TC tracks and landfall using 50-member ensembles of TC-permitting regional model simulations for five hurricane seasons characterized by different levels of TC activity. The control simulations are seasonal hindcasts in which AEWs were prescribed through the eastern lateral boundary condition using reanalysis. In the experiments, we suppressed AEWs by applying a 2–10 day filter to the eastern lateral boundary condition. In response to AEW suppression, we discovered statistically significant increases in Atlantic TC frequency (10%–26%) and landfall (16%–44%), a westward shift in TC genesis location and tracks with increased landfall over the Caribbean Islands, southwestern Gulf Coast, and southeastern US coast, and increases in mid-tropospheric relative humidity in the main development region. In addition, we evaluated TC genesis mechanisms in the absence of AEWs and found evidence that the intertropical convergence zone (ITCZ) intensified and extended northward, resulting in ITCZ wave-breaking that shed vortices which served as TC seeds. By uncovering the connections between TC seed types and the likelihood of TC landfall, this research can provide Atlantic coastal and island communities with useful information to prepare for TC impacts.

African easterly waves↗
NASA NTRS2010

Tropical Cyclones, Hurricanes, and Climate: NASA's Global Cloud-Scale Simulations and New Observations that Characterize the Lifecycle of Hurricanes

One of the primary interests of Global Change research is the impact of climate changes and climate variability on extreme weather events, such as intense tropical storms and hurricanes. Atmospheric climate models run at resolutions of global weather models have been used to study the impact of climate variability, as seen in sea surface temperatures, on the frequency and intensity of tropical cyclones. NASA's Goddard Earth Observing System Model, version 5 (GEOS-5) in ensembles run at 50 km resolution has been able to reproduce the interannual variations of tropical cyclone frequency seen in nature. This, and other global models, have found it much more difficult to reproduce the interannual changes in intensity, a result that reflects the inability of the models to simulate the intensities of the most extreme storms. Better representation of the structures of cyclones requires much higher resolution models. Such improved representation is also fundamental to making best use of satellite observations. In collaboration with NOAA's Geophysical Fluid Dynamics Laboratory, GEOS-5 now has the capability of running at much higher resolution to better represent cloud-scale resolutions. Global simulations at cloud-permitting resolutions (10- to 3.5-km) allows for the development of realistic tropical cyclones from tropical storm 119 km/hr winds) to category 5 (>249km1hr winds) intensities. GEOS-5 has produced realistic rain-band and eye-wall structures in tropical cyclones that can be directly analyzed against satellite observations. For the first time a global climate model is capable of representing realistic intensity and track variability on a seasonal scale across basins. GEOS-5 is also used in assimilation mode to test the impact of NASA's observations on tropical cyclone forecasts. One such test, for tropical cyclone Nargis in the Indian Ocean in May 2008, showed that observations from Atmospheric Infrared Sounder (AIRS) and the Advanced Microwave Sounding Unit (AMSU-A) on Aqua substantially reduced forecast track errors. Tropical cyclones in the northern Indian Ocean pose serious challenges to operational weather forecasting systems, partly due to their shorter lifespan and more erratic track, compared to those in the Atlantic and the Pacific. SA is also bringing several state of the art instruments in recent field campaigns to peer under the clouds and study the inner workings of the tropical storms. With the Genesis and Rapid Intensification Processes (GRIP) experiment, a NASA Earth science field experiment in 2010 that includes the Global Hawk Unmanned Airborne System (UAS) configured with a suite of in situ and remote sensing instruments that are observing and characterizing the lifecycle of hurricanes, we expect significant improvement in our understanding of the track and intensification processes with the assimilation of the satellite and field campaign observations of meteorological parameters in the numerical prediction models.

Putman, William M.↗