Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mesoscale dynamics”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Flow reduction of hydrocarbon liquid in silica nanochannel: Insight from many-body dissipative particle dynamics simulations

A modified many-body dissipative particle dynamics (mDPD) model recently developed for realistic mesoscale multiphase flow simulations is rigorously parameterized, calibrated, and applied for elucidating the flow mechanisms of hydrocarbon liquids (i.e., heptane in this work) in amorphous silica cylindrical nanochannels with inner diameters ranging from 4.5 to 22.5 nm. The simulation results suggest the presence of a strong threshold of pressure gradient under which heptane cannot be driven to flow. The threshold for the 4.5 nm diameter pore is 10 to 100 times as high as for the 9–22.5 nm diameter pore, highlighting a remarkable nanoconfinement effect. Fluid viscosity is found to exhibit a shear-thinning phenomenon with intensity to weaken with increasing channel diameter — a phenomenon not observed in nanochannel flow of liquid water and gas in literature. Most remarkably, the radial profiles of average longitudinal flow velocity fitted by the modified Hagen-Poiseuille equation showed a negative slip length (–2.5% to –0.5% relative to the diameter) and a reduction of apparent permeability by 16% to 23%. Furthermore, this finding suggests silica nanochannels tend to deter hydrocarbon flow, a phenomenon that is opposed to the flow enhancement reported in most of the prior nanochannel flow studies in literature.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

A Novel Method for Diagnosing Land–Atmosphere Coupling Sensitivity in a Single-Column Model

The response of boundary layer properties and cloudiness to changes in surface evaporative fraction (EF) is investigated in a single-column model to quantify the locally coupled impact of subgrid surface variations on the atmosphere during summer. Sensitive coupling days are defined when the model atmosphere exhibits large variations across a range of EFs centered on the analyzed value. Coupling sensitivity exists as both positive feedback (cloudiness increases with EF) and negative feedback (clouds increase with decreasing EF) regimes. The positive regime manifests in shallow convection situations, which are capped by a strengthened inversion and subsidence, restricting the vertical extent of convection to just above the boundary layer. Surfaces with larger EF (greater surface latent heat flux) can inject more moisture into the vertically confined system, lowering the cloud base and an increasing cloud liquid water path (LWP). Negative feedback regimes tend to manifest when large-scale deep convection, such as from mesoscale convective systems and fronts, is advected through the domain, where convection strengthens over surfaces with a lower EF (greater surface sensible heat flux). The invigoration of these systems by the land surface leads to an increase in LWP through strengthened updrafts and stronger coupling between the boundary layer and the free atmosphere. These results apply in the absence of heterogeneity-induced mesoscale circulations, providing a one-dimensional dynamical perspective on the effect of surface heterogeneity. Finally, this study provides a framework of intermediate complexity, lying between parcel theory and high-resolution coupled land–atmosphere modeling, and therefore isolates the relevant first-order processes in land–atmosphere interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Perpendicular crossing chains enable high mobility in a noncrystalline conjugated polymer

The nature of interchain π-system contacts, and their relationship to hole transport, are elucidated for the high-mobility, noncrystalline conjugated polymer C16-IDTBT by the application of scanning tunneling microscopy, molecular dynamics, and quantum chemical calculations. The microstructure is shown to favor an unusual packing motif in which paired chains cross-over one another at near-perpendicular angles. By linking to mesoscale microstructural features, revealed by coarse-grained molecular dynamics and previous studies, and performing simulations of charge transport, it is demonstrated that the high mobility of C16-IDTBT can be explained by the promotion of a highly interconnected transport network, stemming from the adoption of perpendicular contacts at the nanoscale, in combination with fast intrachain transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2025

Active-noise-induced dynamic clustering of passive colloidal particles

Active fluids generate spontaneous, often chaotic mesoscale flows. Harnessing these flows to drive soft materials embedded within an active fluid into structures with controlled length scales and lifetimes is a key challenge at the interface between the fields of active matter and nonequilibrium self-assembly. Here, we present a simple and efficient computational approach to model soft materials advected by active fluids, by simulating particles moving in a spatiotemporally correlated noise field. To illustrate our approach, we simulate the dynamical self-organization of repulsive colloids within such an active noise field. The colloids form structures whose sizes and dynamics can be tuned by the correlation time and length of the active fluid, and range from small rotating droplets to clusters with internal flows and system-spanning sizes that vastly exceed the active correlation length. Our results explain how the interplay between active fluid time and length scales and emergent driven assembly can be used to rationally design functional assemblies. More broadly, our approach can be used to efficiently simulate diverse active fluids and other systems with spatiotemporally correlated noise.

Brownian dynamics↗
DOE OSTI2023

Mesoscale Convective Systems in DYAMOND Global Convection‐Permitting Simulations

Abstract This study examines the deep convection populations and mesoscale convective systems (MCSs) simulated in the DYAMOND (DYnamics of the atmospheric general circulation modeled on non‐hydrostatic domains) winter project. A storm tracking algorithm is applied to six DYAMOND simulations and a global high‐resolution satellite cloud and precipitation data set for comparison. The simulated frequencies of tropical deep convection and organized convective systems vary widely among models and regions, although robust MCSs are generally underestimated. The diurnal cycles of MCS initiation and mature stages are well simulated, but the amplitudes are exaggerated over land. Most models capture the observed MCS lifetime, cloud shield area, rainfall volume and movement speed. However, cloud‐top height and convective rainfall intensity are consistently overestimated, and stratiform rainfall area and amount are consistently underestimated. Possible causes for the model differences compared to observations and implications for future model developments are discussed.

58 GEOSCIENCES↗
DOE OSTI2024

Ab initio generalized Langevin equation

We introduce a machine learning–based approach called ab initio generalized Langevin equation (AIGLE) to model the dynamics of slow collective variables (CVs) in materials and molecules. In this scheme, the parameters are learned from atomistic simulations based on ab initio quantum mechanical models. Force field, memory kernel, and noise generator are constructed in the context of the Mori–Zwanzig formalism, under the constraint of the fluctuation–dissipation theorem. Combined with deep potential molecular dynamics and electronic density functional theory, this approach opens the way to multiscale modeling in a variety of situations. Here, we demonstrate this capability with a study of two mesoscale processes in crystalline lead titanate, namely the field-driven dynamics of a planar ferroelectric domain wall, and the dynamics of an extensive lattice of coarse-grained electric dipoles. In the first case, AIGLE extends the reach of ab initio simulations to a regime of noise-driven motions not accessible to molecular dynamics. In the second case, AIGLE deals with an extensive set of CVs by adopting a local approximation for the memory kernel and retaining only short-range noise correlations. The scheme is computationally more efficient than molecular dynamics by several orders of magnitude and mimics the microscopic dynamics at low frequencies where it reproduces accurately the dominant far-infrared absorption frequency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Mesoscale Convective Systems Tracking Method Intercomparison (MCSMIP): Application to DYAMOND Global km‐Scale Simulations

Abstract Global kilometer‐scale models represent the future of Earth system modeling, enabling explicit simulation of organized convective storms and their associated extreme weather. Here, we comprehensively evaluate tropical mesoscale convective system (MCS) characteristics in the DYAMOND (DYnamics of the atmospheric general circulation modeled on non‐hydrostatic domains) simulations for both summer and winter phases. Using 10 different feature trackers applied to simulations and satellite observations, we assess MCS frequency, precipitation, and other key characteristics. Substantial differences (a factor of 2–3) arise among trackers in observed MCS frequency and their precipitation contribution, but model‐observation differences in MCS statistics are more consistent across trackers. DYAMOND models are generally skillful in simulating tropical mean MCS frequency, with multi‐model mean biases ranging from −2%–8% over land and −8%–8% over ocean (summer vs. winter). However, most DYAMOND models underestimate MCS precipitation amount (23%) and their contribution to total precipitation (17%). Biases in precipitation contributions are generally smaller over land (13%) than over ocean (21%), with moderate inter‐model variability. While models better simulate MCS diurnal cycles and cloud shield characteristics, they overestimate MCS precipitation intensity and underestimate stratiform rain contributions (up to a factor of 2), particularly over land, albeit observational uncertainties exist. Additionally, models exhibit a wide range of precipitable water in the tropics compared to reanalysis and satellite observations, with many models showing exaggerated sensitivity of MCS precipitation intensity to precipitable water. The MCS metrics developed here provide process‐oriented diagnostics to guide future model development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Micro-beam bending of FCC bicrystals: A comparison between defect dynamics simulations and experiments

To understand the role of the grain boundary (GB) in plasticity at small scale, a concurrently coupled mesoscale plasticity model was developed to simulate micro-bending of bicrystalline micron-sized beams. By coupling dislocation dynamics (DD) with a finite element model (FEM), a novel defect dynamics model provides the means to investigate intricate interactions between dislocations and GBs under various loading conditions. Our simulations of micro-bending agree well with corresponding micro-bending experiments, and they show that mechanical response of bicrystals could have not only hardening but also softening depending on the characters of the GB. In addition, changing the location of the GB in the microbeams results in different mechanical responses; GBs located at the neutral plane show softening compared to single crystals, while inclined GBs located halfway along the length of the beam show little effect. In conclusion, simulation results could provide a clear picture on detailed dislocation-GB interactions, and quantitative resolved shear stress analysis supplemented by dislocation density distribution is used to analyze the mechanical response of bicrystalline samples.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Structure of Complex Liquid–Liquid Extraction Organic Phases for Rare Earth Separations

Complex, multicomponent liquids with hierarchical structure and phase transitions are encountered in many natural and industrial processes, including in chemical separations. One notable example is aggregation and organic phase splitting in liquid–liquid extraction (LLE) of metal ions. While these two phenomena that have long been closely associated, a mechanistic link between mesoscale structure and the capacity-limiting organic phase splitting remains elusive due to complexity of these systems. Here, in this study, we combine small-angle X-ray scattering (SAXS), X-ray photon correlation spectroscopy (XPCS), and molecular dynamics simulation to reveal a comprehensive picture of structure at the nano- and mesoscale in these complex solutions. For the representative case of rare earth extraction from an acidic aqueous phase by a malonamide extractant in dodecane, we investigate a wide range of process-relevant extractant and acid concentrations to provide a complete picture of how aggregation depends on composition. We decompose organic phase structure from SAXS into two contributions, which together can capture the scattering at all compositions: composition fluctuations described by the Ornstein–Zernike equation at low wavenumber Q, and nanostructure modeled by a “pre-peak” at intermediate Q. The former contains information about the thermodynamics of demixing, while the latter reflects nanoscopic self-assembly of the extractant and extracted solutes. While fluctuations have typically not been considered in the literature, we find they in fact dominate the total structure for nearly all practical conditions. As only the fluctuations have a strong temperature response, we confirm this attribution with temperature-dependent SAXS measurements, including for extracted europium nitrate complexes. SAXS and XPCS measurements near the critical point find static and dynamic scaling consistent with theory. Overall, this new paradigm for understanding LLE organic phases connects composition, nanoscale, and mesoscale structuring to phase behavior, providing both a comprehensive picture of solution structure and a quantitative link between aggregation and third phase formation.

Peroutka, Allison A. [Argonne National Laboratory ↗
DOE OSTI2020

An in situ shearing x-ray measurement system for exploring structures and dynamics at the solid–liquid interface

Revealing interfacial structure and dynamics has been one of the essential thematic topics in material science and condensed matter physics. Synchrotron-based x-ray scattering techniques can deliver unique and insightful probing of interfacial structures and dynamics, in particular, in reflection geometries with higher surface and interfacial sensitivity than transmission geometries. We demonstrate the design and implementation of an in situ shearing x-ray measurement system, equipped with both inline parallel-plate and cone-and-plate shearing setups and operated at the advanced photon source at Argonne National Laboratory, to investigate the structures and dynamics of end-tethered polymers at the solid-liquid interface. With a precise lifting motor, a micrometer-scale gap can be produced by aligning two surfaces of a rotating upper shaft and a lower sample substrate. A torsional shear flow forms in the gap and applies tangential shear forces on the sample surface. The technical combination with nanoscale rheology and the utilization of in situ x-ray scattering allow us to gain fundamental insights into the complex dynamics in soft interfaces under shearing. In this work, we demonstrate the technical scope and experimental capability of the in situ shearing x-ray system through the measurements of charged polymers at both flat and curved interfaces upon shearing. Furthermore, through the in situ shearing x-ray scattering experiments integrated with theoretical simulations, we aim to develop a detailed understanding of the short-range molecular structure and mesoscale ionic aggregate morphology, as well as ion transport and dynamics in soft interfaces, thereby providing fundamental insight into a long-standing challenge in ionic polymer brushes with a significant technological impact.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid–liquid mixtures

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Multiscale computational understanding and growth of 2D materials: a review

Abstract The successful discovery and isolation of graphene in 2004, and the subsequent synthesis of layered semiconductors and heterostructures beyond graphene have led to the exploding field of two-dimensional (2D) materials that explore their growth, new atomic-scale physics, and potential device applications. This review aims to provide an overview of theoretical, computational, and machine learning methods and tools at multiple length and time scales, and discuss how they can be utilized to assist/guide the design and synthesis of 2D materials beyond graphene. We focus on three methods at different length and time scales as follows: (i) nanoscale atomistic simulations including density functional theory (DFT) calculations and molecular dynamics simulations employing empirical and reactive interatomic potentials; (ii) mesoscale methods such as phase-field method; and (iii) macroscale continuum approaches by coupling thermal and chemical transport equations. We discuss how machine learning can be combined with computation and experiments to understand the correlations between structures and properties of 2D materials, and to guide the discovery of new 2D materials. We will also provide an outlook for the applications of computational approaches to 2D materials synthesis and growth in general.

Momeni, Kasra↗
DOE OSTI2020

Simulation and Emulation of X-Ray Diffraction from Dynamic Compression Experiments

Many important aspects of the dynamic thermo-mechanical response of materials occur at the mesoscale, i.e. a physical scale of interactions smaller than what can be adequately described by homogenous behaviors, yet larger than the scale of the atomic lattice. Concurrent advancements in computational power, continuum theory, and experimental diagnostics are enabling unprecedented understanding of such interactions. However, we cannot develop a sufficient level of confidence in such mesoscale capability until the constitutive description of the underlying constituents is reliably representative of their actual physical behavior. Therefore, there is a strong need to combine experimental, modeling, and data-science techniques to validate models of the thermomechanical response of individual single crystals. One experimental diagnostic with high potential impact to shock physics and materials science is in-situ x-ray diffraction. This paper is primarily focused on simulation of x-ray diffraction in shock physics, but with an aim toward quantifying parametric uncertainty of simulation models. Here, we develop and demonstrate a data-science and model-driven approach to constrain the parameterization of continuum models of crystal lattice deformation associated with the shock response of crystalline materials. The framework is built around the connection between continuum hydrodynamic simulations of lattice deformation and a new Bragg diffraction simulation code, BarberShop. The dynamic deformation of a crystal lattice is modeled using the DiscoFlux model within an arbitrary Lagrangian-Eulerian hydrodynamic code, FLAG. These detailed continuum simulations of lattice deformation can be computationally slow, thus a statistical model is used to emulate the evolution of lattice deformation fields in time and across the considered model parameter space. Emulated lattice deformation fields can then be generated rapidly for any combination of physics model parameters. In turn, these fields can be fed into BarberShop to realize a rapid prediction of Bragg diffraction patterns associated with particular values of physics model parameters. The framework enables parameterization of the single crystal model to obtain Bragg diffraction patterns that most closely resemble a corresponding measurement. Furthermore, the framework naturally provides sensitivities of the lattice deformation to the physics parameters. We highlight the utility of this framework through the application to a synthetic closed-loop inverse problem leading to the parameterization of a single crystal material model. As a model problem, we consider the dynamic response of the energetic molecular crystal, cyclotrimethylenetrinitramine (or RDX), under dynamic compression induced by simulated flyer plate impact experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electronic 1/f noise as a probe of dimensional effects on spin-glass dynamics

Over the past decade, spin-glass simulations have improved to the point that they now access time- and length-scales comparable to experiments at the mesoscale. A recent series of thin-film field-cooled/zero-field-cooled magnetization (FC/ZFC) experiments demonstrated activated spin dynamics, with a temperature-independent activation energy proportional to the logarithm of the film thickness and with coefficients in remarkable agreement with the simulation. These measurements require the application of small magnetic fields, which has been shown to affect the spin-glass energy landscape. Measurements of the 1/f noise in metallic spin-glasses have been previously shown to be a sensitive probe of the spin dynamics, and the measurements can be made without applying a magnetic field. In this mini-review, we review these techniques and discuss how transport measurements can fit into the current landscape of spin-glass measurements. We compare previous measurements to more recent measurements on similar films, made with ostensibly different cooling protocols, and compare both the previous and recent measurements to the magnetometry. The transport measurements—taken over a wider range of temperature than magnetometry—suggest that the maximum spin-glass energy barrier height is temperature-dependent, not fixed, possibly due to two-dimensional dynamics. We discuss this possibility, along with future measurements, which may be able to resolve this mystery.

1/f noise↗
DOE OSTI2026

Nonlocal Effect of Percolated Particle Networks on Viscoelasticity of Polymer–Filler Nanocomposites: A Mesoscale Simulation Study

With nanoparticles (NPs) as fillers, polymer nanocomposites (PNCs) usually exhibit enhanced mechanical properties. However, a direct connection between the microscopic structural relaxation and macroscopic mechanical properties of PNCs remains to be established. To investigate the micro-to-macro connection, we develop a mesoscale model, in which the NPbridging polymer chains are represented by a dynamic bonded interaction between NPs, and the bulk polymer matrix is implicitly modeled by overdamped Langevin dynamics. Extensive equilibrium simulations are performed to quantify the microscopic dynamics of model PNCs. Systematic analyses of modified Rouse dynamics, dynamic structure factor, and relaxation modulus uncover that the microscopic relaxation dynamics of PNCs are significantly decelerated across different length scales because of nonlocal effects of percolated particle networks a phenomenon that has not been adequately captured in prior simulation studies. We find that NPvolume- fraction and NP-bonding-energy barrier are the two critical variables that affect bulk viscoelasticity the most. The proposed mesoscale model is versatile and provides a powerful framework for studying structure−property relations of different PNCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Comparing Synoptic Pattern Evolution for Flash‐Flood‐Producing and Non‐Flash‐Flood‐Producing Mesoscale Convective Systems in the United States

Understanding how the short-term evolution of synoptic weather patterns influence Mesoscale Convective Systems (MCSs) is essential, as these systems are responsible for over half of central U.S. flash floods, leading to substantial socioeconomic and water resource management impacts. This study analyzes long-term MCS data, flash flood reports, and atmospheric reanalyses from 2007 to 2017 using a machine learning clustering algorithm to examine how the synoptic weather patterns evolve prior to MCS initiation. While the clusters reflect seasonal and regional differences in MCS occurrence, they do not consistently distinguish between MCSs that do and do not produce flash floods. Systems in the southern Great Plains are more flood-prone when a synoptic-scale forcing, located near the system, drives strong water vapor transport from the nearby moisture source. More generally under different synoptic weather patterns, a broader precipitating area is the most dominant factor governing MCS flash flood potential.

atmospheric dynamics↗
DOE OSTI2021

Ocean fronts and eddies force atmospheric rivers and heavy precipitation in western North America

Atmospheric rivers (ARs) are responsible for over 90% of poleward water vapor transport in the mid-latitudes and can produce extreme precipitation when making landfall. However, weather and climate models still have difficulty simulating and predicting landfalling ARs and associated extreme precipitation, highlighting the need to better understand AR dynamics. Here, using high-resolution climate models and observations, we demonstrate that mesoscale sea-surface temperature (SST) anomalies along the Kuroshio Extension can exert a remote influence on landfalling ARs and related heavy precipitation along the west coast of North America. Inclusion of mesoscale SST forcing in the simulations results in approximately a 40% increase in landfalling ARs and up to a 30% increase in heavy precipitation in mountainous regions and this remote impact occurs on two-week time scales. The asymmetrical response of the atmosphere to warm vs. cold mesoscale SSTs over the eddy-rich Kuroshio Extension region is proposed as a forcing mechanism that results in a net increase of moisture flux above the planetary boundary layer, prompting AR genesis via enhancing moisture transport into extratropical cyclones in the presence of mesoscale SST forcing.

54 ENVIRONMENTAL SCIENCES↗