Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesoscale modeling”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Comparison of Energetic Electron Fluxes Measured by GPS and THEMIS Spacecraft in the Inner Magnetosphere

Key elements of space weather models are energetic electron fluxes in the inner magnetosphere and the outer radiation belt. Flux depletion is driven by various loss processes: scattering into atmosphere, magnetopause shadowing. Flux enhancement is driven by various acceleration processes: local wave-particle interactions, radial transport, plasma sheet injections. Many of these processes operate on ~ hour timescales. Such mesoscale flux variations are not well traced by equatorial spacecraft with much longer orbits. Energetic electron detectors onboard the Global Positioning System (GPS) constellation provide a unique opportunity for probing such ~ hour-scale flux variations. Measurements from up to 23 identically instrumented GPS satellites cover a wide energy and L-shell range with a subhour time resolution. However, their orbits are inclined and thus all measurements at L-shell >4.3 are off-equatorial. Here, in this report, we present a comparison of equatorial THEMIS and nonequatorial GPS measurements of omnidirectional ≤600 keV electron fluxes. Such a comparison allows us to derive coefficients for using off-equatorial GPS fluxes to infer the equatorial values. These coefficients depend on particle energy and L-shell. We demonstrate a new data set derived from GPS measurements and discuss how it can be used to investigate mesoscale dynamics of energetic electron fluxes in the inner magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Integration of seismic-pressure-petrophysics inversion of continuous active-seismic monitoring data for monitoring and quantifying CO 2 plume (Final Report)

The overall objective of this project is to develop and validate an integrated package of joint seismic-pressure-petrophysics inversion (jSPPI) of continuous active-source seismic monitoring dataset capable of providing real-time monitoring of CO 2 plume during geologic carbon sequestration (GCS). The three specific developments include: (a) the methodologies for fast seismic full waveform inversion of continuous active source seismic monitoring, (CASSM) datasets for simultaneously estimating velocity and attenuation, and with data assimilation; (b) joint Bayesian petrophysical inversion of seismic models and pressure data for providing and updating CO 2 saturation models; (c) the methods using multiple datasets including (Crainfield and Frio-II borehole) synthetic, laboratory, and field CASSM datasets. The outcomes of jSPPI include (a) a workflow for processing CASSM data, (b) Bayesian inversion algorithms using CASSM data and pressure response data, and (c) integration with data assimilation algorithms for continuously updating site-specific models used for prediction and reservoir management. The validation of joint FWI will be conducted using synthetic models based on the Cranfield and Frio experiments as well as field CASSM datasets collected as part of the Frio-II pilot injection. To quantify and map the mass and distribution of CO 2 (saturation), we will jointly invert velocity and attenuation measurements from the FWI with a Bayesian approach using a rock physics model for attenuation (e.g., White’s attenuation model with two selected patch sizes (White, 1976; Dutta and Seriff, 1979)). The Bayesian inversion will be applied to each time step in the CASSM survey in an updating scheme, which integrates with an ensemble of reservoir simulations at each step. A more complete experimental validation dataset will be collected as part of a mesoscale (2-3 m) gas-CO 2 injection experiment utilizing a higher frequency version of the CASSM system developed for laboratory studies; the integrated inversion will be demonstrated using this dataset which will provide both a dense geometry as well as more precise secondary confirmation measurements (e.g. saturation) typically not available in the field. The resulting real-time map of CO 2 saturation is able to provide a deeper scientific understanding of the complex, time-varying dynamics of subsurface fluid flow migration path as well as the rapid detection of CO 2 leakage hazards.

25 ENERGY STORAGE↗
DOE OSTI2025

The Importance of an Ensemble Approach for Modeling Aerosol‐Convection Interactions

Aerosol-convection interactions modulate cloud microphysics, thermodynamics, and updraft intensity, contributing to climate-scale aerosol-radiative forcing. However, quantifying aerosol indirect effects in mixed-phase deep convection remains challenging due to uncertainties in parameterized physics and initial conditions driving nonlinear evolution of convective processes. This study investigates the convective updraft sensitivity to random initial temperature perturbations using an idealized ensemble modeling framework informed by in situ thermodynamic and aerosol observations from the DOE TRACER field campaign. We analyze the impact of small-scale initial perturbations on updraft velocity and supersaturation, and determine the ensemble size required to minimize stochastic internal variability. Results show that minor thermodynamic perturbations can produce updraft variability comparable to aerosol-induced changes reported in prior work. An ensemble of 10 members sufficiently reduces variability, enabling robust investigation of aerosol-related updraft invigoration and informing the design of future ensemble-based aerosol-convection interaction studies in environments with significant spatiotemporal mesoscale thermodynamic and aerosol heterogeneity.

aerosol-convection interactions↗
DOE OSTI2023

MEMPHIS (Mesoscle Multiphysics Phase Field Simulator v.0

SAND2024-01427O Mesoscale Multiphysics Phase Field Simulator (MEMPHIS) is used for solving interfacial problems and can be applied to many evolutionary problems. Some examples are solidification dynamics, viscous fingering, fracture mechanics, hydrogen embrittlement, vesicle dynamics, electro-mechanical degradation of batteries, aging of microstructures, and segregation. MEMPHIS is written in Fortran 90 in a modular fashion, where model and numerical solvers are decoupled from one another. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Dingreville, Remi↗
DOE OSTI2020

Grand potential sintering simulations of doped UO 2 accident-tolerant fuel concepts

Chromium-doped UO is a widely-studied near-term deployable accident-tolerant fuel concept because it results in a dense, large-grain structure that increases the fuel resistance to densification, swelling, and fission gas release. A new charged-interstitial mechanism was recently proposed to describe the behavior of dopants like chromium in sintered UO 2 . Based on that mechanism, manganese was suggested as an even stronger dopant than chromium. In the current work we use mesoscale sintering simulations in an effort to validate the new mechanism. We compare the relative behavior of Cr-doped and undoped UO 2 against experimental data in the literature. We also make predictions of the relative behavior of Mn-doped UO 2 . Simulations are done using the phase field-based grand potential sintering model. Dopants have two effects on sintered UO 2 . They marginally increase the densification rate and greatly increase the average grain size. Both of these effects are individually tested in small-scale simulations. Then large 3D sintering simulations are used to test the combined behavior of both effects. The results for the densification rate simulations are consistent with experiments for Cr-doped fuels. However, we note the grain growth rates are lower than what is found in the literature for Cr-doped fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Improved Simulation of East Asian Summer Monsoon in the High-resolution CESM1 and Its Causes

Based on the high- and low-resolution Community Earth System Model, version 1 (CESM1), and corresponding simulations from phase 6 of the Coupled Model Intercomparison Project (CMIP6), we compare the interannual variability of the East Asian summer monsoon (EASM). The EASM interannual variability is characterized by the anomalous western North Pacific anticyclone (WNPAC) circulation and the dipole rainfall pattern with a negative southern lobe over the western North Pacific and a positive northern lobe along the Meiyu–Baiu region, which is better reproduced by the high-resolution models. Here, the reason for the improvement in the high-resolution models has been attributed to the better simulation of the warm temperature advection from the wind anomalies on the climatological temperature gradient. Positive sea surface temperature (SST) anomalies over the tropical Indian Ocean are the key to the improved wind anomalies featuring a WNPAC in the high-resolution models. The warm SST anomalies over the tropical Indian Ocean strengthen the WNPAC by triggering a Kelvin-wave response to the enhanced heat release induced by the increased precipitation. Based on the mixed-layer heat budget analysis, the warm SST anomalies over the western Indian Ocean in the high-resolution CESM1 are tied to the anomalous easterly wind along the equator, which reduces surface evaporation and upwelling. Therefore, the better simulations of air-sea feedback and the oceanic mesoscale eddy over the western Indian Ocean are the key for the improved simulation of the EASM interannual variations in the high-resolution CESM1.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Persistent meanders and eddies lead to quasi-steady Lagrangian transport patterns in a weak western boundary current

The Brazil Current (BC) is a weak western boundary current flowing along the Southwestern Atlantic Ocean. It is frequently described as a flow with intense mesoscale activity and relatively low volume transport between 5.0 and 10.0 Sv. We use a 13-year eddy-resolving primitive-equation simulation to show that the presence of persistent meanders and eddies leads to characteristic quasi-steady Lagrangian transport patterns, aptly extracted through climatological Lagrangian Coherent Structures (cLCS). The cLCS position the surface expression of the BC core along the 2000-m isobath, in excellent visual agreement with high resolution satellite sea-surface temperature and the model Eulerian mean velocity. The cLCS deformation pattern also responds to zonally persistent cross-shelf SSH transition from positive (high) values near coastline to low (negative) values between 200- and 2000-m and back to positive (high) offshore from the 2000-m isobath. Zonally-paired cyclonic and anticyclonic structures are embedded in this transition, also causing the cLCS to deform into chevrons. An efficient transport barrier is identified close to the 200-m isobath confirmed by limited inshore movement of drogued drifters and accurately indicated by an along slope maxima of climatological strength of attraction. We also show that the persistent cyclonic and anticyclonic structures may induce localized cross-shelf transport. Regions of low climatological strength of attraction coincide with large shelves and with stagnant synthetic trajectories. We also show that cLCS accurately depict trajectories initiated at the location of Chevron’s spill (November 2011) as compared to synthetic and satellite-tracked trajectories, and the outline of the oil from that accident. There is also an agreement between the large-scale oil slicks reaching the Brazilian beaches (from August 2019 to February 2020) and the strength of climatological attraction at the coast. Our work also clarifies the influence of persistent mesoscale structures on the regional circulation. The identification and quantitative description of climatological Lagrangian coherent structures is expected to improve the effectiveness of future emergency response to oil spills, contingency planning, rescue operations, larval and fish connectivity assessment, drifter launch strategies, waste pollutant and marine debris dispersion and destination.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Coastal-Urban-Rural Atmospheric Gradient Experiment (CoURAGE) Science Plan

Understanding the mechanisms governing the urban atmospheric environment is critical for informing urban populations regarding the impacts of climate change and associated mitigation and adaptation measures. Earth system (climate and weather) models have not yet been adapted to provide accurate predictions of climate and weather variability within cities, nor do they provide well-tested representations of the impacts of urban systems on the atmospheric environment. These limitations are largely due to limited field data available for testing and development of these models. We will deploy the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility’s first Mobile Facility (AMF1) to the mid-Atlantic region surrounding the city of Baltimore for the Coast-Urban-Rural Atmospheric Gradient Experiment (CoURAGE). This deployment will create a four-node regional atmospheric observatory network including Baltimore and its three primary surrounding environments – rural, urban, and bay. CoURAGE investigators will study the interactions among the Earth’s surface, the atmospheric boundary layer, aerosols and atmospheric composition, clouds, radiation, and precipitation at each site, and examine how the spatial gradients across the region interact to create the climate conditions in Baltimore. This study will determine the degree to which Baltimore’s atmospheric environment depends on interactive feedbacks in the atmospheric system and the degree to which conditions in Baltimore depend on the surrounding environment. Some topics of interest include how urban land management exacerbates heat waves, the impact of regional mesoscale winds (nocturnal jet, bay breeze) on urban air pollution and cloud cover, and the impact of the urban heat island and aerosol production on heavy precipitation events. Understanding this integrated coast-urban-rural system quantitatively and with good accuracy and precision is critical to informing climate adaptation and mitigation efforts in the city of Baltimore. The understanding gained should be applicable to many similar coastal, mid-latitude urban centers. Another important objective of CoURAGE is to improve the representation of the climate of coastal cities in Earth systems models (ESMs). CoURAGE investigators will use the observations to test current ESMs, identify weaknesses and work towards improved simulations of this complex environment. The ARM core facility will be deployed in the city of Baltimore, complementing the Baltimore Social-Environmental Collaborative (BSEC), a DOE urban integrated field laboratory (UIFL). Ancillary sites will be deployed to rural Maryland northwest of Baltimore, and to the southern end of Kent Island within Chesapeake Bay. The fourth node will be a long-term atmospheric observatory operated in Beltsville, Maryland by Howard University and the Maryland Department of the Environment. Measurements will be conducted for one year, starting in December of 2024. There will be two intensive operational periods (IOPs), one in summer and one in winter, when the ancillary sites will be enhanced with additional balloon launches, tethered balloon system (TBS) operation, and added atmospheric composition measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Global High‐Resolution Mesoscale Convective System Database Using Satellite‐Derived Cloud Tops, Surface Precipitation, and Tracking

Abstract A new methodology is developed to construct a global (60°S–60°N) long‐term (2000–2019) high‐resolution (∼10‐km h) mesoscale convective system (MCS) database by tracking MCS jointly using geostationary satellite infrared brightness temperature ( T b ) and precipitation feature (PF) characteristics from the Integrated Multi‐satellitE Retrievals for GPM precipitation data sets. Independent validation shows that the satellite‐based MCS data set is able to reproduce important MCS statistics derived from ground‐based radar network observations in the United States and China. We show that by carefully considering key PF characteristics in addition to T b signatures, the new method significantly improves upon previous T b ‐only methods in detecting MCSs in the midlatitudes for all seasons. Results show that MCSs account for over 50% of annual total rainfall across most of the tropical belt and in selected regions of the midlatitudes, with a strong seasonality over many regions of the globe. The tracking database allows Lagrangian aspects such as MCS lifetime and translational speed and direction to be analyzed. The longest‐lived MCSs preferentially occur over the subtropical oceans. The land MCSs have higher cloud‐tops associated with more intense convection, and oceanic MCSs have much higher rainfall production. While MCSs are observed in many regions of the globe, there are fundamental differences in their dynamic and thermodynamic structures that warrant a better understanding of processes that control their evolution. This global database provides significant opportunities for observational and modeling studies of MCSs, their characteristics, and roles in regional and global water and energy cycles, as well as their hydrologic and other impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Two-dimensional patterning of mesoscale fibers using acoustophoresis

The performance of functional composites can rely critically on the arrangement of secondary phases; for example, patterned networks of conductive particles can impart anisotropic thermal, electric or ionic conductivity while preserving flexibility in the matrix. We demonstrate the use of standing acoustic waves to generate periodic patterns of short fibers. We extend the range of possible patterns with the first demonstration of both rectangular grids and arrays of octagons interspersed with rectangles. These newly demonstrated patterns are rationalized using theoretical models of acoustic forces and torques on fibers that account for two-dimensional spatial variations arising from applied acoustic fields. The models enable simulations of fiber motion, which are used to (i) map out final fiber positions as a function of initial position and orientation, and (ii) corroborate experiments visualizing fiber motion and final patterns. This approach provides a fast and accurate way to predict emergent fiber patterns as a function of excitation modality and fiber length. The theory and experiments clearly indicate strong coupling between the length of the fibers and the spacing of the acoustic nodes. Furthermore, this coupling is used to estimate reductions in percolation thresholds associated with the ratio of fiber length and acoustic wavelength.

25 ENERGY STORAGE↗
DOE OSTI2023

Accelerating discrete dislocation dynamics simulations with graph neural networks

Discrete dislocation dynamics (DDD) is a widely employed computational method to study plasticity at the mesoscale that connects the motion of dislocation lines to the macroscopic response of crystalline materials. However, the computational cost of DDD simulations remains a bottleneck that limits its range of applicability. Here, we introduce a new DDD-GNN framework in which the expensive time-integration of dislocation motion is entirely substituted by a graph neural network (GNN) model trained on DDD trajectories. As a first application, we demonstrate the feasibility and potential of our method on a simple yet relevant model of a dislocation line gliding through an array of obstacles. We show that the DDD-GNN model is stable and reproduces very well unseen ground-truth DDD simulation responses for a range of straining rates and obstacle densities, without the need to explicitly compute nodal forces or dislocation mobilities during time-integration. Our approach opens new promising avenues to accelerate DDD simulations and to incorporate more complex dislocation motion behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Two-Day Westward-Propagating Inertio-Gravity Waves during GoAmazon

The dynamical and thermodynamical features of Amazonian 2-day westward-propagating inertia–gravity (WIG) waves are examined. On the basis of a linear regression analysis of satellite brightness temperature and data from the 2014–15 Observations and Modeling of the Green Ocean Amazon (GoAmazon) field campaign, it is shown that Amazonian WIG waves exhibit structure and propagation characteristics consistent with the n = 1 WIG waves from shallow water theory. These WIG waves exhibit a pronounced seasonality, with peak activity occurring from March to May and a minimum occurring from June to September. Evidence is shown that mesoscale convective systems over the Amazon are frequently organized in 2-day WIG waves. Results suggest that many of the Amazonian WIG waves come from preexisting 2-day waves over the Atlantic, which slow down when coupled with the deeper, more intense convection over tropical South America. In contrast to WIG waves that occur over the ocean, Amazonian 2-day WIG waves exhibit a pronounced signature in surface temperature, moisture, and heat fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Understanding the distinct impacts of MCS and non-MCS rainfall on the surface water balance in the central US using a numerical water-tagging technique

Warm-season rainfall associated with mesoscale convective systems (MCSs) in the central US is characterized by higher intensity and nocturnal timing compared to rainfall from non-MCS systems, suggesting their potentially different footprints on the land surface. To differentiate the impacts of MCS and non-MCS rainfall on the surface water balance, a water tracer tool embedded in the Noah land surface model with multi-parameterization options (WT-Noah-MP) is used to numerically “tag” water from MCS and non-MCS rainfall separately during April to August (1997-2018) and track their transit in the terrestrial system. From the water-tagging results, over 50% of warm-season rainfall leaves the surface-subsurface system through evapotranspiration by the end of August, but non-MCS rainfall contributes a larger fraction. However, MCS rainfall plays a more important role in generating surface runoff. These differences are mostly attributed to the rainfall intensity differences. The higher intensity MCS rainfall tends to produce more surface runoff through infiltration excess flow and drives a deeper penetration of the rainwater into the soil. Over 70% of the top 10 percentile runoff is contributed by MCS rainfall, demonstrating its important contribution to local flooding. In contrast, lighter intensity non-MCS rainfall resides mostly in the top layer and contributes more to evapotranspiration through soil evaporation. Diurnal timing of rainfall has negligible effects on the flux partitioning for both MCS and non-MCS rainfall. Differences in soil moisture profiles for MCS and non-MCS rainfall and the resultant evapotranspiration suggest differences in their roles in soil moisture-precipitation feedbacks and ecohydrology.

Hu, Huancui↗
DOE OSTI2020

Characterizing Dislocation Defects in Dark Field X-ray Microscopy Images of Bulk Aluminum

Material defects play a large role in material response under shock loading, yet the process of how these defects initiate, propagate, and annihilate is not well understood. Using a newly developed diagnostic, dark field X-ray microscopy (DFXM), we can now visualize the behavior of dislocation defects in materials at the mesoscale under varying conditions. Using data from the European Synchrotron Radiation Facility, we apply a variety of image processing techniques to capture relevant features to locate and characterize size and orientation of dislocation defects in bulk aluminum DFXM images. Beyond simple visualization, this analysis drives statistical characterization of the defects and their dynamic response to further improve the relevant physics models.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Differences in Lightning—Storm Microphysical Relationships in Subtropical South America Versus the Continental U.S. and Their Impact on Model Lightning Parameterizations

Satellite studies have shown that some of the deepest convection with high lightning flash rates occurs downstream of the Andes in subtropical South America (SSA). Using high resolution ground‐based observations in SSA, we explore the relationship between lightning and radar‐inferred storm microphysics. Results show these storms have significantly larger 35‐dBZ volumes and more graupel per lightning flash compared to storms in the U.S. We then test microphysics‐based lightning parameterizations for storms in SSA on a 1‐km convection‐permitting simulation of a mesoscale convective system. Results from the simulation match well with observed flash rates. If parameterizations based on U.S. storms had been applied, flash rates would have been overestimated by a factor of eight. These results suggest that storms in SSA are fundamentally different compared to storms in the U.S., and caution should be exercised when applying microphysics‐based lightning parameterizations in different regions of the world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A modified many-body dissipative particle dynamics model for mesoscopic fluid simulation: methodology, calibration, and application for hydrocarbon and water

The many-body dissipative particle dynamics (mDPD) is a prominent mesoscopic multiphase model for fluid transport in mesoconfinement. However, it has been a long-standing challenge for mDPD (and other multiphase-enabled DPD models) to accurately predict real-fluid static and dynamic properties simultaneously. We have developed a modified mDPD model to overcome the issue and a rigorous calibration approach that uses reference data, including experimental and/or molecular dynamics (MD) simulations to parameterise the modified mDPD for real fluids. We choose heptane as a representative example of hydrocarbon in source rocks to demonstrate the model's capability to accurately predict the equation of state (EOS), free surface tension, diffusivity, and viscosity. Our timing test shows that the modified mDPD is 400–500 times faster than its MD counterpart for simulating bulk heptane in equivalent volumes. To further demonstrate the robustness of the model, we revisited the benchmark problem of mesoscopic modelling of water, in which all the previous works on DPD reported only a limited portion of the water properties. Here, we show that the modified mDPD can provide accurate modelling of water static and dynamic properties and an EOS that matches the experimental data to a large range of confinement pressure.

42 ENGINEERING↗