Search NASA⌕ Search

SEARCH · Search NASA

Results for “mesoscale simulations”

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 487 records · Page 27

Soft matter roadmap

Soft materials are usually defined as materials made of mesoscopic entities, often self-organised, sensitive to thermal fluctuations and to weak perturbations. Archetypal examples are colloids, polymers, amphiphiles, liquid crystals, foams. The importance of soft materials in everyday commodity products, as well as in technological applications, is enormous, and controlling or improving their properties is the focus of many efforts. From a fundamental perspective, the possibility of manipulating soft material properties, by tuning interactions between constituents and by applying external perturbations, gives rise to an almost unlimited variety in physical properties. Together with the relative ease to observe and characterise them, this renders soft matter systems powerful model systems to investigate statistical physics phenomena, many of them relevant as well to hard condensed matter systems. Understanding the emerging properties from mesoscale constituents still poses enormous challenges, which have stimulated a wealth of new experimental approaches, including the synthesis of new systems with, e.g. tailored self-assembling properties, or novel experimental techniques in imaging, scattering or rheology. Theoretical and numerical methods, and coarse-grained models, have become central to predict physical properties of soft materials, while computational approaches that also use machine learning tools are playing a progressively major role in many investigations. This Roadmap intends to give a broad overview of recent and possible future activities in the field of soft materials, with experts covering various developments and challenges in material synthesis and characterisation, instrumental, simulation and theoretical methods as well as general concepts.

36 MATERIALS SCIENCE↗

Resonant soft X-ray scattering in polymer science

Resonant soft X-ray scattering (RSoXS) is an emerging, powerful technique to probe the nano-to-mesoscale structure of polymers and other molecules. It joins together small-angle X-ray scattering (a statistical nanoprobe) with X-ray spectroscopy that brings with it unique chemical and bond-orientation sensitivity. Through over a decade of discovery and development, RSoXS is moving from a niche technique applied to organic electronic thin films to a mature tool applicable to a plethora of polymeric and molecular systems, encompassing new modalities, analyses, and simulation methods. This development promises to deliver increasingly quantitative answers to challenging questions in polymer science as well as expand its usefulness to complementary fields. Here, this review presents a full synopsis of the technique, including background on the theoretical underpinnings, measurement best practices, and examples of recent RSoXS applications and discoveries provided here to accelerate the transition to a broader range of soft matter and polymeric fields.

36 MATERIALS SCIENCE↗

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↗

Particle interspacing effects on the mechanical behavior of a Fe–TiB 2 metal matrix composite using FFT-based mesoscopic field dislocation mechanics

This paper presents an application to metal matrix composites (MMCs) of an enhanced elasto-viscoplastic Fast Fourier Transform (EVP-FFT) formulation coupled with a phenomenological continuum Mesoscale Field Dislocation Mechanics (MFDM) theory. Contrary to conventional crystal plasticity, which only accounts for plastic flow and hardening induced by statistically stored dislocations (SSDs), MFDM-EVP-FFT also describes the evolution of polarized geometrically necessary dislocation (GND) density and its effect on both plastic flow and hardening. Numerical results for a Fe–TiB 2 MMC made of a ferrite matrix (α-Fe) and elastic ceramic particles (TiB 2 ) are presented. Full-field simulations are performed using synthetic periodic unit cells representative of the MMC, with single-crystalline and polycrystalline matrix, for different particle interspacing distances. A strong dependence of the predicted equivalent stress, cumulated plastic strain and GND density fields with particle interspacing distance is observed, in contrast with conventional crystal plasticity. Correlations between these mechanical fields and microstructural features, and their influence on local and global mechanical behavior are examined for the different MMC microstructures.

42 ENGINEERING↗

Pair distribution function analysis of nano-object assemblies

The atomic pair distribution function (aPDF) analysis technique, also known as the total scattering method, which considers both Bragg and diffuse scattering, has been used extensively to probe local atomic arrangements in crystalline and disordered materials. In contrast, there have been limited applications of the PDF in self-assembled nanomaterials, which represent a class of materials built from nanoscale objects, such as nano-colloids, micelles and proteins. As distinguished from atoms, nano-objects have polydispersity in size and shape, and such form-factor effects complicate the application of PDF analysis to nano-systems. Here, in this paper, the application of the PDF is extended to spherical nano-object assemblies and the formulae for the nano-PDF (nPDF) are derived, showing some differences from the aPDF. By numerical simulations, the properties of the nPDF (peak broadening and pattern profile) are studied systematically as a function of structural features, such as nano-object parameters (size and size polydispersity) and assembly structural features (size, shape, structure type and lattice disorder), and of data processing parameters ( q cut-off and `missing' data in ultra-small-angle regions). The nPDF analysis method is found to provide an effective route to revealing not only nanoscale but also mesoscale structural properties, for example the morphology of a nano-assembly.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Alumina – Stabilized SEI and CEI in Potassium Metal Batteries

Aluminum oxide (Al 2 O 3 ) nanopowder is spin-coated onto both sides of commercial polypropene separator to create artificial solid-electrolyte interphase (SEI) and artificial cathode electrolyte interface (CEI) in potassium metal batteries (KMBs). Here, this significantly enhances the stability, including of KMBs with Prussian Blue (PB) cathodes. For example, symmetric cells are stable after 1,000 cycles at 0.5 mA/cm 2 –0.5 mAh/cm 2 and 3.0 mA/cm 2 –0.5 mAh/cm 2 . Alumina modified separators promote electrolyte wetting and increase ionic conductivity (0.59 vs. 0.2 mS/cm) and transference number (0.81 vs. 0.23). Cryo-stage focused ion beam (cryo-FIB) analysis of cycled modified anode demonstrates dense and planar electrodeposits, versus unmodified baseline consisting of metal filaments (dendrites) interspersed with pores and SEI. Alumina-modified CEI also suppresses elemental Fe crossover and reduces cathode cracking. Mesoscale modeling of metal – SEI interactions captures crucial role of intrinsic heterogeneities, illustrating how artificial SEI affects reaction current distribution, conductivity and morphological stability.

25 ENERGY STORAGE↗

Mechanisms of heat flux across the Southern Greenland continental shelf in 1/10° and 1/12° ocean/sea ice simulations

The increased presence of warm Atlantic water on the Greenland continental shelf has been connected to the accelerated melting of the Greenland Ice Sheet, particularly in the southwest and southeast shelf regions. Results from two high-resolution coupled ocean-sea ice simulations that utilized either the 1/10-degree Parallel Ocean Program (POP) or the 1/12-degree HYbrid Coordinate Ocean Model (HYCOM) are used to understand the flux of heat on and off the southern Greenland shelf. The analysis reveals that the region of greatest heat flux onto the shelf is southeast Greenland. On the southwestern shelf, heat is mainly exported from the shelf to the interior basins. We identify differences in the shelf break current structure and on-shelf heat content between the two simulations. Just south of the Denmark strait, there is a seasonally persistent pattern of multi-day variability in the cross-shelf heat flux in both simulations. In the POP simulation, this high-frequency signal results in net on-shore heat flux. In the HYCOM simulation, the signal is weaker and results in net off-shelf heat flux. This variability is consistent with Denmark Strait Overflow eddies traveling along the shelf break.

58 GEOSCIENCES↗

Multilayer interface tracking model of pure tungsten oxidation

Here, we present a numerical model to predict oxide scale growth on tungsten surfaces under exposure to oxygen at high temperatures. The model captures the formation of four thermodynamically-compatible oxide sublayers, WO 2 , WO 2.72 , WO 2.9 , and WO 3 , on top of the metal substrate. Oxide layer growth is simulated by tracking the oxide/oxide and oxide/metal interfaces using a sharp-interface Stefan model coupled to diffusion kinetics. The model is parameterized using selected experimental measurements and electronic structure calculations of the diffusivities of all the oxide subphases involved. We simulate oxide growth at temperatures of 600°C and above, extracting the power law growth exponents in each case, which we find to deviate from classical parabolic growth in several cases. We conduct a comparison of the model predictions with an extensive experimental data set, with reasonable agreement at most temperatures. While many gaps in our understanding still exist, this work is a first attempt at embedding the thermodynamic and kinetic complexity of tungsten oxide growth into a comprehensive mesoscale kinetic model that attempts to capture the essential features of tungsten oxidation to fill existing knowledge gaps and guide and enhance future tungsten oxidation models.

36 MATERIALS SCIENCE↗

Spatiotemporal Variations in the Characteristics of Mesoscale Convective Systems Over Indian Monsoon Zone

Here, this study examines cloud and precipitation properties of mesoscale convective systems (MCSs) in the Indian monsoon zone using 6 years of high-resolution satellite-based tracking data. MCSs exhibit strong regional and seasonal variations, peaking during the summer monsoon (June–September). Concentrated along a southeast-northwest band of the monsoon trough, MCS contribute 40%–70% of total precipitation and predominantly move westward at an average speed of 6.3 m s -1 . Long-lived MCSs (> 40 hr), although 10 times less frequent than short-lived MCSs (5–20 hr), produce 2–3 times more rainfall. Longer MCS lifespans and higher rainfall are associated with stronger deep-layer wind shear and enhanced lower free-tropospheric moisture at initiation. Land and ocean MCSs exhibit distinct diurnal cycles. Most land-based MCSs initiate around 16 LT and mature at 18 LT, whereas Bay of Bengal (BoB) MCSs initiate between 04 and 06 LT and peak by 14 LT. Convective precipitations from MCSs significantly surpasses that from non-MCSs systems, emphasizing their critical role in intense monsoon rainfall. Short- and long-lived MCSs showed distinct lifecycle evolution. Long-lived MCSs over the BoB show 37% larger cold-cloud shields, 44% larger precipitation feature areas, and 33% higher rain rates than their land counterparts, yielding 67% more total rainfall. These BoB MCSs also grow faster and experience more frequent cloud mergers early in their lifecycle. These findings offer critical insights into monsoonal MCSs lifecycle evolution, providing observational constraints for improving model simulations of monsoonal convection.

54 ENVIRONMENTAL SCIENCES↗

Considering computational speed vs. accuracy: Choosing appropriate mesoscale RVE boundary conditions

Modeling a material’s microstructure using continuum theories allows for inspection of the relationship between coarse scale and fine scale behaviors. Computational limits generally require selection of a sub-volume from a bulk sample in order to directly model the microstructure. Boundary conditions are applied to the sub-volume to mimic the excluded bulk material. Appropriate selection of boundary conditions helps effectively determine the appropriate spatial scale required of the sub-volume. Applicable boundary conditions include direct displacement, periodic, and uniform traction. While direct displacement and periodic boundary conditions are commonly used, uniform traction boundary conditions have seen limited use due to rigid body stability issues in simulations of compression or shear deformation. A new application of uniform traction boundary conditions was developed through linear constraint equations, similar to approaches employed by direct displacement and periodic boundary conditions, to quench rigid body motions with minimal interference of the relative deformation of the model. These boundary conditions were tested by compressing several synthetically generated periodic microstructures using the finite element method. Evaluating the effective stiffness along the compression axis, the direct displacement boundary condition produced the stiffest response, whereas the uniform traction boundary condition produced the most compliant. Periodic boundary conditions produced the same response for all volumes analyzed and both the direct displacement and uniform traction boundary conditions trended toward the periodic response as the domain volume increased. Computational performance was also evaluated for each boundary condition using implicit and explicit solvers. Direct displacement boundary conditions presented the lowest computational cost of all of the boundary conditions followed by periodic then uniform traction. The computational expense of periodic and uniform traction boundary conditions limited the viable spatial scale and mesh resolutions able to be simulated. Selection of appropriate boundary conditions for specific uses need to be a balance between allowable computational expense and accuracy of the method. Techniques for evaluating which boundary conditions to use are discussed.

42 ENGINEERING↗

Monte-Carlo modeling of phonon thermal transport using DFT-based anisotropic dispersion relations over the full Brillouin zone

In this work, we present a Monte Carlo (MC) approach to solve the phonon Boltzmann transport equation (BTE) in which the anisotropic phonon dispersion relations over the full Brillouin zone (BZ) are used. In this approach, the discretization of the BZ used to compute the phonon relaxation time places constraints on the direction of scattered phonons in the real-space simulation domain. The phonon dispersion and phonon relaxation times are calculated using the density functional theory (DFT) approach. The modified MC approach is validated by a close examination of its ability to simulate phonon transport in both the ballistic and diffusive regimes for multiple materials including GaAs, InAs, ThO 2 , and α-U. In doing so, the phonon thermal conductivities from 100 K to 1000 K are calculated and compared with traditional non-transport solution of the phonon BTE. It is found that the phonon thermal conductivities of α-U and ThO 2 obtained from MC simulations using isotropic dispersion are larger than the values obtained using anisotropic phonon dispersion relations over the full BZ. The effect of phonon-defect scattering on the thermal conductivity of ThO 2 is also studied as an application of the current MC approach and found to agree with previously computed values in the literature. The MC solver developed here has been parallelized as a step to demonstrate its potential to solving computationally intensive phonon thermal transport problems at the mesoscale.

36 MATERIALS SCIENCE↗

Reconciliation of asynchronous satellite-based $\mathrm{NO_2}$ and $\mathrm{XCO_2}$ enhancements with mesoscale modeling over two urban landscapes

Fossil fuel carbon dioxide (CO 2ff ), the main driver of global warming and climate change, is often co-emitted with nitrogen oxides (NO x ) and precursors to ground-level ozone from anthropogenic sources like power plants or vehicles. In urban and suburban areas, satellite-based NO 2 can be used as a proxy to track the emissions of CO 2ff . Because of NO 2 ’s shorter lifetime, urban NO 2 plumes are more distinguishable from backgrounds and more sensitive to variations in emissions. However, the combination of these two gases is limited by the asynchrony among NO 2 and CO 2 monitoring satellites. We used CO 2ff simulated by the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) model to reconcile the tropospheric NO 2 vertical column density (VCD) from the Tropospheric Monitoring Instrument (TROPOMI) and column-averaged dry-air mole fractions of carbon dioxide enhancements (ΔXCO 2 ) from Orbiting Carbon Observatory 3 (OCO-3) Snapshot Area Maps (SAMs) over a multicity area, Washington D.C.-Baltimore (DC-Balt), and a basin city, Mexico City. NO 2 /CO 2ff ratios over DC-Balt are smaller than Mexico City, indicative of stricter emission restrictions, a more combustion-efficient vehicle fleet, and higher combustion efficiency due to lower altitude in DC-Balt. For single-track cases, the spatial correlations between NO 2 and ΔXCO 2 over Mexico City are stronger than DC-Balt because the NO 2 and CO 2 are mostly trapped in the valley of Mexico City, while DC-Balt is severely affected by distant sources (i.e., US East Coast cities). Using multi-track averaging, spatial correlation coefficients increase with the number of days used for averaging. The correlations reached a maximum when averaging >12 continuous images for DC-Balt and > 10 continuous images for Mexico City. This finding indicates that multi-track averaging using modeled CO 2ff as a proxy is helpful to filter the noise in single-track images, to cancel the interference from distant sources, and to magnify correlations between NO 2 and CO 2ff . Mexico City showed stronger spatial correlations but weaker temporal correlations than DC-Balt due to biomass burning hot spots and large transport errors caused by the trapping effects of the surrounding mountains. Tracking the 20-day moving average of CO 2ff emissions using TROPOMI NO 2 seems technically feasible, considering the relationship between correlation coefficients and the number of available satellite images.

54 ENVIRONMENTAL SCIENCES↗

Future Changes in Midwest Extreme Precipitation Depend on Storm Type

Midwestern U.S. extreme precipitation is associated with multiple storm types including mesoscale convective systems (MCSs) and/or training thunderstorms, tropical cyclone (TC) remnants, and winter storms. Anthropogenic warming is expected to increase climatological precipitation globally, however, there may be little correspondence with regional storm-based changes. Furthermore, uncertainty remains in precipitation-temperature scaling due to use of convective parameterization in most global models. In this study, we investigated historically impactful extreme precipitation events from multiple types of Midwest storms using the Weather Research and Forecasting model at convection-permitting resolution. We simulated five-member ensembles of historical hindcasts and experiments representing the storms in the future using the pseudo-global warming method. We found that future precipitation changes depend on storm type, with increases near Clausius-Clapeyron (CC) for winter storms, no consensus for MCSs and/or training thunderstorms, and sub-CC increases for TC remnants. This research highlights the importance of considering storm type in future extreme precipitation projections.

54 ENVIRONMENTAL SCIENCES↗

Examining Ohmic Losses in Fuel-Cell Catalyst Layers with Different Pt/C Ratios

A model of a fuel-cell cathode with distinct proton conductivities at the agglomerate and electrode scales was used to simulate polarization curves under conditions with facile oxygen transport. Here, the model was applied to electrodes containing agglomerates with different ratios of platinum to carbon in order to examine the hypothesis that intra-agglomerate proton conductivity is substantially lower than ionomer conductivity. The model predicts decreasing performance and increasing sensitivity to Pt/C ratio as intra-agglomerate conductivity declines relative to a constant through-electrode conductivity. Experiments on catalyst layers with Vulcan supports do not show the predicted trend indicating that mesoscale intra-agglomerate conductivity is not orders of magnitude lower than bulk ionomer conductivity.

08 HYDROGEN↗

Resolved gravity waves in the tropical stratosphere: Impact of horizontal resolution and deep convection parametrization

Convectively generated gravity waves (CGGWs) are important for numerical weather prediction due to their effect on the quasi-biennial oscillation (QBO) in the stratosphere. Here, using global ECMWF IFS simulations at TCo7999 (or 1.25 km), TCo2559 (or 3.9 km) and TCo1279 (or 7.8 km) horizontal resolutions, sensitivity of resolved CGGWs to the horizontal resolution and to the explicit versus parametrized representation of deep convection is elucidated during the westerly shear phase of the QBO. Parametrized deep convection is found to inhibit CGGWs, resulting in a twofold reduction in CGGW forcing. When deep convection is explicitly resolved, the total CGGW forcing is almost unchanged across the horizontal resolutions. However, the contribution of long and mesoscale CGGWs (with horizontal wavelengths 100 km ≤ $λ_h$ < $1,900$ km) to the total CGGW forcing decreases and the contribution of smaller-scale CGGWs (with $λ_h$ < $100$ km) increases as the horizontal resolution increases. At the maximum CGGW forcing altitude, at TCo7999 resolution 43% of the total CGGW forcing is due to long and mesoscale waves, whereas at TCo2559 and TCo1279 resolutions their contribution is 70% and 90%, respectively. While CGGW forcing from long and mesoscale waves is similar at TCo7999 resolution with explicit deep convection and at TCo1279 resolution with parametrized deep convection, CGGW forcing from these waves is artificially enhanced at TCo1279 and TCo2559 resolutions with explicit deep convection. This is due to the explicit deep convection being too strong and having too much variance for 100 km ≤ $λ_h$ < $1,900$ km. Therefore, parametrizations of deep convection and CGGWs (to account for forcing from waves with $λ_h$ < $100$ km) are required even at TCo2559 resolution. Additionally, resolved CGGW forcing at TCo7999 resolution is examined for the easterly shear phase of the QBO; similar to the westerly shear phase, the smaller-scale waves contribute 55% to the total CGGW forcing at the maximum CGGW forcing altitude.

54 ENVIRONMENTAL SCIENCES↗

Uranium recovery from seawater using amidoxime-based braided polymers synthesized from acrylic fibers

Global demand for nuclear energy is expected to rise in the coming decades. To meet these growing needs, new uranium resources must be explored. One of the potential alternatives to traditional uranium mining is oceanic uranium. The capture and recovery of uranium from the ocean has been under investigation for some time, with many recent studies focused on amidoxime based adsorbents. These adsorbents while able to achieve high uranium recovery capacities are, nevertheless, expensive to produce and adsorb a significant amount of hard to remove vanadium. The purpose of this study is to evaluate the adsorption performance of amidoxime based polymer braids synthesized from acrylic fibers which are designed to significantly cut polymer synthesis and conditioning costs. Adsorption experiments were performed in a recirculating raceway flume system at environmental conditions, approximately 10.8 oC, with 40-micron prefiltered seawater over 28 days for small and large sized LCW polymer braids. For both braid sizes, the adsorption of vanadium was far lower on the LCW materials than on Oak Ridge National Laboratory’s amidoxime based AF1 and AI8 small braids, with the LCW small braids outperforming both materials with respect to uranium adsorption under similar conditions. Adsorption modeling was used to simulate the performance at higher temperatures based on 20oC experiments previously performed with similar materials. Simulation results indicated that the small LCW braids would have a less significant advantage with respect to uranium adsorption over both AF1 and AI8 at 20oC and adsorb slightly less uranium than AF1 at 31oC. Vanadium adsorption by LCW small braids was less than one third of vanadium adsorption by either AF1 or AI8 for all temperatures. This behavior indicates that the LCW materials are able to achieve comparable uranium adsorption to other amidoxime adsorbents previously developed, while being cheaper to produce and adsorbing significantly less vanadium.

Wiechert, Alexander I.↗

Segmentation of RDX and TNT in X‐Ray Computed Tomography Reconstructions of Melt‐Cast Explosives

ABSTRACT Three‐dimensional mesoscale characterization of heterogeneous melt‐cast high explosives is challenging because of the difficulty differentiating binder from explosive crystals: two functionally different materials which are typically similar in density by design. Here, we report an algorithm which can differentiate hexahydro‐1,3,5‐trinitro‐1,3,5‐triazine (RDX) from 2,4,6‐trinitrotoluene (TNT) in x‐ray computed tomography (CT) volumes with tens of microns resolution. This method allows us to quantify RDX/TNT content, porosity, and RDX domain size. We calibrated the segmentation algorithm using simulated x‐ray CT volumes containing object models of RDX crystals within a TNT matrix. We then segmented and analyzed CT data for Composition B (Comp B), a 60/40 RDX/TNT mixture, and Cyclotol, a 75/25 RDX/TNT mixture. We examined melt‐cast samples fabricated with 100% theoretical maximum density (TMD) and 85% TMD. For the 100% TMD Comp B and Cyclotol samples, the RDX content values calculated by segmentation were 3% and 9% lower, respectively, than the values measured by high‐performance liquid chromatography on material from the same synthesis lots. This result is consistent with the expected underreporting of RDX content resulting from x‐ray CT resolution limits on RDX particles with diameters smaller than 25 µm. The 85% TMD samples were less accurately segmented with our algorithm due to the confounding presence of voids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atmospheric Rivers in the Eastern and Midwestern United States Associated With Baroclinic Waves

Atmospheric rivers (ARs) significantly impact the hydrological cycle and associated extremes in western continental regions. Recent studies suggest ARs also influence water resources and extremes in continental interiors. AR detection tools indicate that AR conditions are relatively frequent in areas east of the Rocky Mountains. The origin of these ARs, whether from synoptic-scale waves or mesoscale processes, is unclear. This study uses meteorological composite maps and transects of AR conditions during the four seasons. The analysis reveals that ARs east of the Rockies are associated with long-wave, baroclinic Rossby waves. This result demonstrates that eastern North American ARs are dynamically similar to their western coastal counterparts, though mechanisms for vertical moisture flux differ between the two. These findings provide a foundation for understanding future climate change and ARs in this region and offer new methods for evaluating climate model simulations.

58 GEOSCIENCES↗