Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SIMULATION METHODS”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Dynamic energy absorption characteristics of hollow microlattice structures

Hollow microlattice structures are promising candidates for advanced energy absorption and their characteristics under dynamic crushing are explored. The energy absorption can be significantly enhanced by inertial stabilization, shock wave effect and strain rate hardening effect. In this paper we combine theoretical analysis and comprehensive finite element method simulation to decouple the three effects, and then obtain a simple model to predict the overall dynamic effects of hollow microlattice structures. Inertial stabilization originates from the suppression of sudden crushing of the microlattice and its contribution scales with the crushing speed, v. Shock wave effect comes from the discontinuity across the plastic shock wave front during dynamic loading and its contribution scales with e. The strain rate effect increases the effective yield strength upon dynamic deformation and increases the energy absorption density. A mechanism map is established that illustrates the dominance of these three dynamic effects at a range of crushing speeds. Compared with quasi-static loading, the energy absorption capacity a dynamic loading of 250 m/s can be enhanced by an order of magnitude. The study may shed useful insight on designing and optimizing the energy absorption performance of hollow microlattice structures under various dynamic loads. (C) 2014 Elsevier Ltd. All rights reserved.

Liu, YL↗
DOE OSTI2020

A large-scale, agent-based simulation of metropolitan freight movements with passenger and freight market interactions

This study summarizes the first stage in the implementation of an agent-based freight modeling system that has a global representation of agents and detailed modeling of a large-scale transportation network. The model is used to evaluate the transportation and energy impacts of goods movement across urban and national scales. The framework is implemented within POLARIS, a C++-based Planning and Operations Language for Agent-based Regional Integrated Simulation, which consists of an activity-based modeling (ABM) and dynamic traffic assignment (DTA) system that has robust features for passenger travel. This platform provides a tool to model interactions among consumers, producers, and the transportation system. The main objective of this initial implementation is to implement a freight model within POLARIS following an agent-based paradigm with behavioral and simulation methods. This paper presents the initial framework and illustrates the application of the model. Building upon earlier works, a parcel location assignment algorithm for business establishments in the population is documented, along with a method for estimating establishment production and consumption volumes. In addition to population generation, other features of the model include push-pull supply chains, multimodal path choice, choice of transportation logistics node, and dynamic traffic assignment. A module with e-commerce supply and demand was also developed to analyze the effects of e-commerce delivery on last-mile energy use and congestion.

agent-based freight↗
DOE OSTI2021

Interpretable artificial intelligence and exascale molecular dynamics simulations to reveal kinetics: Applications to Alzheimer's disease

The rapid increase in computing power, especially with the integration of graphics processing units, has dramatically increased the capabilities of molecular dynamics simulations. To date, these capabilities extend from running very long simulations (tens to hundreds of microseconds) to thousands of short simulations. However, the expansive data generated in these simulations must be made interpretable not only by the investigator who performs them but also by others as well. Here, we demonstrate how integrating learning techniques, such as artificial intelligence, machine learning, and neural networks, into analysis pipelines can reveal the kinetics of Alzheimer's disease (AD) protein aggregation. Finally, we review select AD targets, describe current simulation methods, and introduce learning concepts and their application in AD, highlighting limitations and potential solutions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Investigation of inter-subject variation in ultrafine particle deposition across human nasal airways: A study involving children, adults, and the elderly

Occupational and environmental exposure to toxic nanoparticles, driven by the rapid expansion of nanotechnology, raises significant respiratory health concern. Numerous studies have explored airflow and particle dynamics in adult nasal airways, but understanding the impact of age-related anatomical changes in children and the elderly remains limited. This study systematically investigates age-related anatomical variations and associated influence on nasal airflow dynamics and ultrafine particle deposition characteristics. Using Computational Fluid-Particle Dynamics (CFPD) method, simulation was conducted under diverse inhalation conditions spanning a wide age range, including: two children (5 years old), two young adults (in their twenties), and two elderly (over 77 years old). Our results reveal distinctive variations across age groups in anatomical dimensions, which affect distribution of wall shear stress where the elderly and children display unique patterns distinct from the young adults. While total deposition efficiency differs significantly between children and adults, filtration efficiency in the subregion with most deposition, main respiratory, remains consistent. However, inter-subject differences are observed in the vestibular and olfactory regions,emphasizing nuanced impact of age-related anatomical variations. Overall and subregional empirical equations for deposition efficiency were developed by incorporating the combined diffusion parameter, Sc a Δ b , corroborating the use of geometrical characteristic parameters for each specific subject in predicting nasal deposition efficiency across age groups. Our findings contribute to predictive nanoparticle exposure analysis in nasal airways across different age groups, thereby enhancing respiratory healthcare for individuals across the life span.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Exploring strain rate effects upon 3D materials using high speed in situ X-ray tomoscopy

Cellular materials are ubiquitous in our modern society. They may be stochastic gas-blown foams (e.g., polyurethane), foamed starches (e.g., cereals), or, in this case, 3D printed microlattices. Failure in these materials is often driven by surface or sub-surface defects, which may be nucleated at a surface roughness, an interior void, or inclusion interfaces that may not be typically observable. Obfuscating our understanding further, bulk materials are known to exhibit strain-rate-dependent mechanical response, making a subsurface understanding of damage even more critical. For the first time, an in situ uniaxial mechanical loading stage that simultaneously rotates specimens up to 18 Hz was fielded at a synchrotron for 3D tomographic imaging. This capability opens a plethora of materials science opportunities to explore strain rate effects in materials and examining deformation, fracture, and delamination’s (in composites) for a complete 3D picture (movie) of material response. We demonstrate the deformation of 3D printed polymer lattice structures, of three different material types, at 0.25, 1.1, and 2.2 s −1 strain rates. We successfully imaged the 3D deformation of these materials and can directly compare the same printed structure to the three material types at three strain rates, all in 3D. Material point method simulations were applied to one of the materials to better understand the role of voids on the 3D printed structure’s performance.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Economic consequences of aviation system disruptions: A reduced-form computable general equilibrium analysis

The state of the art approach to economic consequence analysis (ECA) is computable general equilibrium (CGE) modeling. However, such models contain thousands of equations and cannot readily be incorporated into computerized systems used by policy analysts to yield estimates of economic impacts of various types of transportation system failures due to natural hazards, human related attacks or technological accidents. This paper presents a reduced-form approach to simplify the analytical content of CGE models to make them more transparent and enhance their utilization potential. The reduced-form CGE analysis is conducted by first running simulations one hundred times, varying key parameters, such as magnitude of the initial shock, duration, location, remediation, and resilience, according to a Latin Hypercube sampling procedure. Statistical analysis is then applied to the “synthetic data” results in the form of both ordinary least squares and quantile regression. The analysis yields linear equations that are incorporated into a computerized system and utilized along with Monte Carlo simulation methods for propagating uncertainties in economic consequences. Although our demonstration and discussion focuses on aviation system disruptions caused by terrorist attacks, the approach can be applied to a broad range of threat scenarios.

Aviation System Disruptions↗
DOE OSTI2024

Langmuir turbulence in suspended kelp farms

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

Bo, Tong (ORCID:0000000260300561)↗
DOE OSTI2026

Transition from Vehicular to Structural Ionic Transport in Electrified Alkali Aqueous Solutions

A molecular understanding of the solvation and dynamics of ions under static electric fields is crucial for modelling a wide range of natural and technological processes. Yet, traditional simulation methods suffer from a trade-off that has to be made between accuracy and statistical convergence. To bridge this gap, herein we extend our recently introduced Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) approach to investigate the solvation structures and ionic transport mechanisms of electrified alkali cationic solutions. We obtain ionic conductivities for Li+, Na+ and Cs+ from the field dependence of the ionic current density in good agreement with experiment. Surprisingly, the migration mechanism is found to be strikingly different for the three ions despite their similar ionic conductivities. While Li+ conducts predominantly through vehicular migration of a stable 4-fold coordinated ion at all field strengths, Cs+ conducts strictly through a structural diffusion mechanism, where 9-12 transient first shell water coordination bonds are broken and renewed. Notably, aqueous Na+ emerges as a “Goldilocks” ion: its ion-water interactions are strong enough to maintain distinct 5-6 fold coordination shells at zero field (unlike Cs+) yet labile enough to be strongly perturbed by electric fields (unlike Li+). As a consequence, we observe an electric field-induced transition from vehicular to structural ionic transport for Na+ that is accompanied by a marked increase in the ionic current density. Our results imply that the conductance mechanism of ions with moderate ion-solvent interactions can be effectively tuned by external electric fields.

Joll, Kit↗
DOE OSTI2020

Utilizing Site Disorder in the Development of New Energy-Relevant Semiconductors

We report controlling site disorder in ternary and multinary compounds enables tuning optical and electronic properties at fixed lattice constants and stoichiometries, moving beyond many of the challenges facing binary alloy systems. Here, we consider possible enhancements to energy-related applications through the integration of disorder-tunable materials in devices such as light-emitting diodes, photonics, photovoltaics, photocatalytic materials, batteries, and thermoelectrics. However, challenges remain in controlling and characterizing disorder. Focusing primarily on II–IV–V 2 materials, we identify three metrics for experimentally characterizing cation site disorder. Complementary to these experiments, we discuss simulation methods to understand disordered materials. Nonidealities, such as off-stoichiometry and oxygen incorporation, can occur while synthesizing metastable disordered materials. While nonidealities may seem undesirable, we describe how if harnessed they could provide another knob for tuning disorder and subsequently properties. To illustrate the effects of disorder on device-relevant properties, we provide case examples of disordered materials and their potential in device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Plasmon 3D Electron Tomography and Local Electric-Field Enhancement of Engineered Plasmonic Nanoantennas

Plasmonic nano-antennas are pushing the limits of optical imaging resolution capabilities in near-field scanning optical microscopy (NSOM). Accordingly, these techniques are driving the basic understanding of photonic and optoelectronic nanoscale devices with applications in sensing, energy conversion, solid-state lighting and information technology. Imaging the localized surface plasmon resonance (LSPR) at the nanoscale is a key to understanding the optical responses of a given tip geometry in order to engineer better plasmonic nano-antennas for near-field experiments. In recent years the advancement of focused ion beam technology provides the ability to directly modify plasmonic structures with nanometer resolution. Also, scanning transmission electron microscopy (STEM) with electron energy loss spectroscopy (EELS) is an established technique allowing imaging of LSPR. Specifically, the combination of these two techniques provides spectrally sensitive two-dimensional (2D) imaging information to better visualize and understand LSPR on the nanometer scale. This can be combined with electron tomography to provide the three-dimensional LSPR distribution. Here in this paper we demonstrate the fabrication of Au nano-pyramids using helium ion microscopy, and analyze the LSPR in 3D reconstructions produced by total variation (TV)-norm minimization of a set of 2D STEM-EELS maps. Additionally, a boundary element simulation method was used to verify the experimentally observed nanopyramid LSPR modes. Finally, we show that the point-spread-functions (PSF) of LSPR mode hot spots in nanopyramids differ to local electric-field enhancement under optical excitation making direct comparison to NSOM experimental resolution difficult. Furthermore, the STEM-EELS results show how LSPR modes are influenced by the tip characteristics, which can inform the development of new nano-antenna designs.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Evolution of Metastable Structures at Bimetallic Surfaces from Microscopy and Machine-Learning Molecular Dynamics

The restructuring of interfaces plays a crucial role in materials science and heterogeneous catalysis. Bimetallic systems, in particular, often adopt very different compositions and morphologies at surfaces compared to the bulk. For the first time, we reveal a detailed atomistic picture of long-time scale restructuring of Pd deposited on Ag using microscopy, spectroscopy, and novel simulation methods. By developing and performing accelerated machine-learning molecular dynamics followed by an automated analysis method, we discover and characterize previously unidentified surface restructuring mechanisms in an unbiased fashion, including Pd-Ag place exchange and Ag pop-out as well as step ascent and descent. Remarkably, layer-by-layer dissolution of Pd into Ag is always preceded by an encapsulation of Pd islands by Ag, resulting in a significant migration of Ag out of the surface and a formation of extensive vacancy pits within a period of microseconds. These metastable structures are of vital catalytic importance, as Ag-encapsulated Pd remains much more accessible to reactants than bulk-dissolved Pd. We report our approach is broadly applicable to complex multimetallic systems and enables the previously intractable mechanistic investigation of restructuring dynamics at atomic resolution.

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

In Situ Insights into Enhanced Cooperative Ligand Exchange Kinetics via Solvent-Induced Restacking in a 2D Metal–Organic Framework

Understanding the reaction kinetics at catalytically active sites is crucial for integrating catalytic two-dimensional (2D) materials into industrial processes. This study focuses on in situ observation of ligand exchange kinetics and solvent-assisted structural restacking transition in the 2D paddle wheel-based MOF [Cu 2 (dttc) 2 ] n (DUT-134(Cu), dttc = dithieno[3,2-b:2′,3′- d]thiophene-2,6-dicarboxylate). The ligand exchange process, involving the replacement of dimethylformamide (DMF) with nitriles such as acetonitrile (ACN), pentanenitrile, and heptanenitrile, was investigated using advanced in situ characterization techniques with high temporal resolution, including powder X-ray diffraction and Raman spectroscopy. The larger analytes exhibited reduced exchange rates, consistent with enhanced steric hindrance and greater diffusion constraints. Interestingly, the study revealed that the exchange of DMF with ACN induces a structural transition to higher symmetry within few seconds, a transition from AB to AA stacking mode of the layers, and a widening of the interlayer distance. Crucially, this structural transition dramatically accelerates the solvent exchange process through cooperative effects, offering critical advantages for catalytic applications. Notably, the reverse exchange from ACN to DMF proceeds more slowly and does not reverse the structural changes, but a new phase is formed with preserved AA stacking. By isotope labeling of linker molecules in combination with two complementary theoretical vibrational simulation methods, the precise assignment of Raman bands and the vibrational modes associated with the ligand exchange process could be achieved. These pioneering insights into the dynamic behavior of 2D MOFs, coupled with ligand exchange, establish a highly promising and transformative approach to achieving enhanced tunability and responsiveness in future catalytic applications.

Layers↗
DOE OSTI2023

Urban Flood Modeling: Uncertainty Quantification and Physics‐Informed Gaussian Processes Regression Forecasting

Abstract Estimating uncertainty in flood model predictions is important for many applications, including risk assessment and flood forecasting. We focus on uncertainty in physics‐based urban flooding models. We consider the effects of the model's complexity and uncertainty in key input parameters. The effect of rainfall intensity on the uncertainty in water depth predictions is also studied. As a test study, we choose the Interconnected Channel and Pond Routing (ICPR) model of a part of the city of Minneapolis. The uncertainty in the ICPR model's predictions of the floodwater depth is quantified in terms of the ensemble variance using the multilevel Monte Carlo (MC) simulation method. Our results show that uncertainties in the studied domain are highly localized. Model simplifications, such as disregarding the groundwater flow, lead to overly confident predictions, that is, predictions that are both less accurate and uncertain than those of the more complex model. We find that for the same number of uncertain parameters, increasing the model resolution reduces uncertainty in the model predictions (and increases the MC method's computational cost). We employ the multilevel MC method to reduce the cost of estimating uncertainty in a high‐resolution ICPR model. Finally, we use the ensemble estimates of the mean and covariance of the flood depth for real‐time flood depth forecasting using the physics‐informed Gaussian process regression method. We show that even with few measurements, the proposed framework results in a more accurate forecast than that provided by the mean prediction of the ICPR model.

Kohanpur, Amir H.↗
DOE OSTI2018

Fragile aspects of topological transition in lossy and parity-time symmetric quantum walks

Quantum walks often provide telling insights about the structure of the system on which they are performed. Here, in PT-symmetric and lossy dimer lattices, the topological properties of the band structure manifest themselves in the quantization of the mean displacement of such a walker. We investigate the fragile aspects of a topological transition in these two dimer models. We find that the transition is sensitive to the initial state of the walker on the Bloch sphere, and the resultant mean displacement has a robust topological component and a quasiclassical component. In PT symmetric dimer lattices, we also show that the transition is smeared by nonlinear effects that become important in the PT-symmetry broken region. By carrying out consistency checks via analytical calculations, tight-binding results, and beam-propagation-method simulations, we show that our predictions are easily testable in today’s experimental systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

The relationship between grain boundary structure, defect mobility, and grain boundary sink efficiency

Nanocrystalline materials have received great attention due to their potential for improved functionality and have been proposed for extreme environments where the interfaces are expected to promote radiation tolerance. However, the precise role of the interfaces in modifying defect behavior is unclear. Using long-time simulations methods, we determine the mobility of defects and defect clusters at grain boundaries in Cu. We find that mobilities vary significantly with boundary structure and cluster size, with larger clusters exhibiting reduced mobility, and that interface sink efficiency depends on the kinetics of defects within the interface via the in-boundary annihilation rate of defects. Thus, sink efficiency is a strong function of defect mobility, which depends on boundary structure, a property that evolves with time. Further, defect mobility at boundaries can be slower than in the bulk, which has general implications for the properties of polycrystalline materials. Finally, we correlate defect energetics with the volumes of atomic sites at the boundary.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2018

Manipulating molecular order in nematic liquid crystal capillary bridges via surfactant adsorption: guiding principles from dissipative particle dynamics simulations

The ability of liquid crystals (LCs) to change orientational order is used in applications, ranging from sensors to displays. The aim of this work is to computationally investigate how surfactant adsorption on cylindrical LC bridges can be used to control such orientational order. Building from classical fundamental lessons, understanding the ordering of mesogens along a preferred axis with the help of molecular modelling contributes to investigations of systems that could be a platform for LC-based sensing applications. The coarse-grained dissipative particle dynamics (DPD) simulation method is implemented here, because it allows us to quantify the effect of molecular features on the properties of meso-scopic systems containing LC bridges, an aqueous solvent, and surfactants at various concentrations. Three surfactant types are modelled with short, medium, and long tail lengths, respectively. All surfactants adsorb at the LC–water interface. It is found that the length of the surfactant hydrophobic tail determines the effectiveness by which the LC order is affected. Short tails are not as effective as long ones. Surfactants with long tails affect the LC order, but, in agreement with experiments, predominantly only within a short distance from the LC–water interface. For these surfactants, the surface density at the LC–water interface is an important knob that can be used to control the order of the LCs. As the effective LC–surfactant interactions change, so does the distribution of the surfactants at the interface. Consistent with theoretical expectations, the results presented here elucidate the effect of molecular features on the anchoring mechanism between surfactants and mesogens within cylindrical bridges dispersed in aqueous systems and could be helpful for designing novel surface-active compounds in the development of advanced sensing devices based on LCs.

74 ATOMIC AND MOLECULAR PHYSICS↗