Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex 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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simulated plant-mediated oxygen input has strong impacts on fine-scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands

Methane (CH 4 ) emissions from wetland ecosystems are controlled by redox conditions in the soil, which are currently underrepresented in Earth system models. Plant-mediated radial oxygen loss (ROL) can increase soil O 2 availability, affect local redox conditions, and cause heterogeneous distribution of redox-sensitive chemical species at the root scale, which would affect CH 4 emissions integrated over larger scales. In this study, we used a subsurface geochemical simulator (PFLOTRAN) to quantify the effects of incorporating either spatially homogeneous ROL or more complex heterogeneous ROL on model predictions of porewater solute concentration depth profiles (dissolved organic carbon, methane, sulfate, sulfide) and column integrated CH 4 fluxes for a tidal coastal wetland. From the heterogeneous ROL simulation, we obtained 18% higher column averaged CH 4 concentration at the rooting zone but 5% lower total CH 4 flux compared to simulations of the homogeneous ROL or without ROL. This difference is because lower CH 4 concentrations occurred in the same rhizosphere volume that was directly connected with plant-mediated transport of CH 4 from the rooting zone to the atmosphere. Sensitivity analysis indicated that the impacts of heterogeneous ROL on model predictions of porewater oxygen and sulfide concentrations will be more important under conditions of higher ROL fluxes or more heterogeneous root distribution (lower root densities). Despite the small impact on predicted CH 4 emissions, the simulated ROL drastically reduced porewater concentrations of sulfide, an effective phytotoxin, indicating that incorporating ROL combined with sulfur cycling into ecosystem models could potentially improve predictions of plant productivity in coastal wetland ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Exascale granular microstructure reconstruction in 3D volumes of arbitrary geometries with generative learning

Reconstructing 3D granular microstructures within volumes of arbitrary geometries from limited 2D image data is crucial for predicting the material properties, as well as performances of structural components accounting for material microstructural effects. We present a novel generative learning framework that enables exascale reconstruction of granular microstructures within complex 3D geometric volumes. Building upon existing transfer learning techniques using pre-trained convolutional neural networks (CNN), we introduce several key innovations to overcome the difficulties inherent in arbitrary geometries. Our framework incorporates periodic boundary conditions using circular padding techniques, ensuring continuity and representativeness of the reconstructed microstructures. We also introduce a novel seamless transition reconstruction (STR) method that creates statistically equivalent transition zones to integrate multiple pre-existing 3D microstructure volumes. Based on STR, we propose a cost-effective strategy for reconstructing microstructures within complex geometric volumes, minimizing computational waste. Validation through numerical experiments using kinetic Monte Carlo simulations demonstrates accurate reproduction of grain statistics, including grain size distributions and morphology. A case study involving the reconstruction of a 4-blade propeller microstructure illustrates the method’s capability to efficiently handle complex geometries. In conclusion, the proposed framework significantly reduces computational demands while maintaining high reconstruction quality, paving the way for scalable microstructure reconstruction in materials design and analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electrostatic Interactions Control the Nanostructure of Conjugated Polyelectrolyte–Polymeric Ionic Liquid Blends

Polyelectrolyte complexation offers unique opportunities to compatibilize polymers with very different backbone chemistries and to control the morphology of the resulting blend via electrostatic manipulation. In this study, we demonstrate the ability to formulate homogeneous complexes of a conjugated polyelectrolyte with a polymeric ionic liquid, utilizing the electrostatic attraction among their oppositely charged side chains. Variation of electrostatic parameters, such as counterion concentration or polymer charge fraction, tunes the morphology of these polymer complexes from homogeneously disordered blend to weakly structured microemulsion where the local ordering arises from backbone-immiscibility-induced microphase segregation. Our experimental observations are in qualitative agreement with both field-theoretic simulation and random-phase approximation calculations. Simulated morphology snapshots suggest and experimental evidence also indicates that the microphase-segregated complex likely takes on a cocontinuous microemulsion structure. Our findings show that ionic interactions are an effective pathway to compatibilize polymers at macroscopic length scales while achieving controlled nanostructures in these ionic blends. Such systems have great potential for engineering the nanostructure of polymers to tailor applications such as nanofiltration, catalysis, and energy storage, where local ordering can enhance the physical properties of an otherwise macroscopically homogeneous structure.

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

Examining metastatic behavior within 3D bioprinted vasculature for the validation of a 3D computational flow model

Understanding the dynamics of circulating tumor cell (CTC) behavior within the vasculature has remained an elusive goal in cancer biology. To elucidate the contribution of hydrodynamics in determining sites of CTC vascular colonization, the physical forces affecting these cells must be evaluated in a highly controlled manner. To this end, we have bioprinted endothelialized vascular beds and perfused these constructs with metastatic mammary gland cells under physiological flow rates. By pairing these in vitro devices with an advanced computational flow model, we found that the bioprinted analog was readily capable of evaluating the accuracy and integrated complexity of a computational flow model, while also highlighting the discrete contribution of hydrodynamics in vascular colonization. This intersection of these two technologies, bioprinting and computational simulation, is a key demonstration in the establishment of an experimentation pipeline for the understanding of complex biophysical events.

42 ENGINEERING↗
DOE OSTI2022

Experimental Procedure and Results for Contact Thermal Conductance Measurements performed during FY 2021

Weapons Test Engineering groups at LANL are responsible for testing and simulating weapon assemblies and subassemblies in support of qualification and certification of nuclear weapon systems. The complex nature of the assemblies themselves coupled with the loading environments, gives rise to low confidence and high uncertainties in the interpretation of both experimental and simulated results. In an effort to reduce these uncertainties, a project was started at LANL to perform focused and systematic experiments combined with simulations to increase confidence in both arenas. A key goal of the work performed under the Delivery Environments (DE) Testbeds to Reduce Uncertainties in Simulations and Tests TRUST program, focuses on identifying and minimizing sources of uncertainty associated with experimental and computational techniques. Specifically, with the above as motivation, a measurement procedure was developed by MST-8 to perform thermal conductivity measurements on cylinders of both similar and dissimilar metallic systems as a function of loading conditions and varying material interface morphologies between these components.

36 MATERIALS SCIENCE↗
DOE OSTI2023

QES-Plume v1.0: a Lagrangian dispersion model

Low-cost simulations providing accurate predictions of transport of airborne material in urban areas, vegetative canopies, and complex terrain are demanding because of the small-scale heterogeneity of the features influencing the mean flow and turbulence fields. Common models used to predict turbulent transport of passive scalars are based on the Lagrangian stochastic dispersion model. The Quick Environmental Simulation (QES) tool is a low-computational-cost framework developed to provide high-resolution wind and concentration fields in a variety of complex atmospheric-boundary-layer environments. Part of the framework, QES-Plume, is a Lagrangian dispersion code that uses a time-implicit integration scheme to solve the generalized Langevin equations which require mean flow and turbulence fields. Here, QES-Plume is driven by QES-Winds, a 3D fast-response model that computes mass-consistent wind fields around buildings, vegetation, and hills using empirical parameterizations, and QES-Turb, a local-mixing-length turbulence model. In this paper, the particle dispersion model is presented and validated against analytical solutions to examine QES-Plume’s performance under idealized conditions. In particular, QES-Plume is evaluated against a classical Gaussian plume model for an elevated continuous point-source release in uniform flow, the Lagrangian scaling of dispersion in isotropic turbulence, and a non-Gaussian plume model for an elevated continuous point-source release in a power-law boundary-layer flow. In these cases, QES-Plume yields a maximum relative error below 6 % when compared with analytical solutions. In addition, the model is tested against wind-tunnel data for a uniform array of cubical buildings. QES-Plume exhibits good agreement with the experiment with 99 % of matched zeros and 59 % of the predicted concentrations falling within a factor of 2 of the experimental concentrations. Furthermore, results also emphasize the importance of using high-quality turbulence models for particle dispersion in complex environments. Finally, QES-Plume demonstrates excellent computational performance.

58 GEOSCIENCES↗
DOE OSTI2023

KinCat v.1.0

SAND2024-02099O The software is designed to allow researchers to perform kinetic Monte Carlo (KMC) simulations of catalytic reactions on a 2D lattice. The code is written efficiently to run on a variety of shared memory computing architectures (e.g. GPU, multi-core) and to natively express the full complexity of lateral interactions on reaction rates. The software allows researchers to perform KMC simulations of catalytic reactions on a 2D lattice. It uses parallel shared-memory computing architectures to reduce run-times and allows for simultaneous simulation of multiple independent runs. 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.

Najm, Habib↗
DOE OSTI2024

Quantum Simulation of SU(3) Lattice Yang-Mills Theory at Leading Order in Large- N c Expansion

Quantum simulations of the dynamics of QCD have been limited by the complexities of mapping the continuous gauge fields onto quantum computers. By parametrizing the gauge invariant Hilbert space in terms of plaquette degrees of freedom, we show how the Hilbert space and interactions can be expanded in inverse powers of N c . At leading order in this expansion, the Hamiltonian simplifies dramatically, both in the required size of the Hilbert space as well as the type of interactions involved. Adding a truncation of the resulting Hilbert space in terms of local energy states we give explicit constructions that allow simple representations of SU(3) gauge fields on qubits and qutrits. This formulation allows a simulation of the real time dynamics of a SU(3) lattice gauge theory on a 5 × 5 and 8 × 8 lattice on ibm_torino with a CNOT depth of 113.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Multiscale simulations reveal TDP-43 molecular-level interactions driving condensation

The RNA-binding protein TDP-43 is associated with mRNA processing and transport from the nucleus to the cytoplasm. TDP-43 localizes in the nucleus as well as accumulating in cytoplasmic condensates such as stress granules. Aggregation and formation of amyloid-like fibrils of cytoplasmic TDP-43 are hallmarks of numerous neurodegenerative diseases, most strikingly present in >90% of amyotrophic lateral sclerosis (ALS) patients. If excessive accumulation of cytoplasmic TDP-43 causes, or is caused by, neurodegeneration is presently not known. In this work, we use molecular dynamics simulations at multiple resolutions to explore TDP-43 self- and cross-interaction dynamics. A full-length molecular model of TDP-43, all 414 amino acids, was constructed from select structures of the protein functional domains (N-terminal domain, and two RNA recognition motifs, RRM1 and RRM2) and modeling of disordered connecting loops and the low complexity glycine-rich C-terminus domain. All-atom CHARMM36m simulations of single TDP-43 proteins served as guides to construct a coarse-grained Martini 3 model of TDP-43. The Martini model and a coarser implicit solvent Cα model, optimized for disordered proteins, were subsequently used to probe TDP-43 interactions; self-interactions from single-chain full-length TDP-43 simulations, cross-interactions from simulations with two proteins and simulations with assemblies of dozens to hundreds of proteins. Our findings illustrate the utility of different modeling scales for accessing TDP-43 molecular-level interactions and suggest that TDP-43 has numerous interaction preferences or patterns, exhibiting an overall strong, but dynamic, association and driving the formation of biomolecular condensates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Update on the Simulation of Commercial Drying of Spent Nuclear Fuel

The purpose of this report is to document improvements in the simulation of commercial vacuum drying procedures at the Nuclear Energy Work Complex at Sandia National Laboratories. Validation of the extent of water removal in a dry spent nuclear fuel storage system based on drying procedures used at nuclear power plants is needed to close existing technical gaps. Operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates additional, well-designed investigations of drying process efficacy and water retention. Scaled tests that incorporate relevant physics and well-controlled boundary conditions are essential to provide insight and guidance to the simulation of prototypic systems undergoing drying processes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Response of a Pressurized Water Reactor Dashpot Region to Commercial Drying Cycles

The purpose of this report is to document updates to the simulation of commercial vacuum drying procedures at the Nuclear Energy Work Complex at Sandia National Laboratories. Validation of the extent of water removal in a dry spent nuclear fuel storage system based on drying procedures used at nuclear power plants is needed to close existing technical gaps. Operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates additional, well-designed investigations of drying process efficacy and water retention. Scaled tests that incorporate relevant physics and well-controlled boundary conditions are essential to provide insight and guidance to the simulation of prototypic systems undergoing drying processes. This report documents testing updates for the Dashpot Drying Apparatus (DDA), an apparatus constructed at a reduced scale with multiple Pressurized Water Reactor (PWR) fuel rod surrogates and a single guide tube dashpot. This apparatus is fashioned from a truncated 5×5 section of a prototypic 17×17 PWR fuel skeleton and includes the lowest segment of a single guide tube, often referred to as the dashpot region. The guide tube in this assembly is open and allows for insertion of a poison rod (neutron absorber) surrogate.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Multibody for Everybody (M4E): A Symbolic Dynamics Modeling Tool with Applications in Simulation, Control, and Optimization

Developing the analytical model of a multibody system is often the initial step in control and optimization. The analytical model (equations of motion) describes a system’s time evolution under specified forcing conditions. Although developing these equations is easy for simple systems, this process becomes more complex for systems composed of multiple bodies. Deriving equations of motion for complex multibody systems requires specialized expertise in multibody dynamics, is time-consuming, and is susceptible to error. To address this issue, this paper presents an open-source, easy-to-use, systematic framework to derive symbolic equations of motion in both Python and MATLAB using the joint coordinate formulation. This formulation results in a set of ordinary differential equations that use the minimum set of coordinates needed to model a system. The symbolic representation provides better insight into the influence of design parameters on system performance, facilitates sensitivity analysis and parameter studies, and supports direct implementation of control and optimization routines. The tool enables numerical simulation for specified parameter sets, is modular for straightforward integration with other tools and libraries, and allows incorporation of hydrodynamics, mooring, and other external forces. The result is a reproducible, extensible pipeline for modeling, simulation, and design of complex multibody systems. The proposed tool is versatile and can be applied to domains such as robotics, control, and design. In addition, we integrated external libraries that provide capabilities for modeling offshore systems such as underwater robots and marine energy converters.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Mechanistic Insights into the Demethylation of Lignin‐Derived Structures Using Protic Ionic Liquids: A Density Functional Theory Study

Lignin valorization is restricted by the stability of its methoxy groups, creating a critical need for efficient demethylation strategies. Here, in this study, density functional theory (DFT) is employed to dissect the mechanistic pathways of demethylation in lignin model compounds, guaiacol and syringol, using protic ionic liquids (PILs) that act as both solvent and catalyst. Conductor‐like Screening Model for Real Solvents (COSMO‐RS) analysis identifies monoethanolammonium acetate ([MEOA][Ace]) as the most promising medium, attributed to its strong hydrogen bonding network and solvation ability. By integrating implicit and explicit solvation models, it is revealed that an acid‐catalyzed hydrolytic mechanism governs demethylation, with PILs stabilizing crucial transition states and intermediates. Complementary electronic structure evaluations, including highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO‐LUMO) gap analysis, charge distribution, and electrostatic potential mapping, demonstrate how PILs lower energetic barriers and enhance reactivity. To simulate realistic environments, this study is extended to lignin dimer complexes with varying water content, uncovering how water‐bridged solvation reverses demethylation preference from guaiacyl to syringyl units. This mechanistic shift aligns with experimental observations showing faster S‐unit reactivity in hydrated systems. Together, these findings provide atomic‐level insight into lignin demethylation dynamics and highlight how tuning acid concentration in PILs can accelerate kinetics, enabling a rational pathway toward next‐generation biomass conversion technologies.

COSMO-RS↗
DOE OSTI2025

Unraveling the Hsp70-ROS-autophagy axis in pentachlorophenol-challenged lung and liver epithelial cells

Pentachlorophenol (PCP) was extensively utilized as an organochlorine pesticide and wood preservative in the United States from the 1930s until the Environmental Protection Agency (EPA) imposed restrictions due to concerns about its toxicity and potential carcinogenic properties. Although it is no longer widely used, PCP remains a concern due to its environmental persistence and potential for long-term health effects. Significant occupational and environmental exposures have likely occurred, with the health and economic costs of PCP exposure potentially being substantial given its known toxicity. Notably, PCP exhibits rapid absorption through both the skin and respiratory system and has been shown to cause hepatotoxicity, developmental toxicity, immunotoxicity, irritation, and carcinogenicity in laboratory animal studies. PCP exposure induces oxidative stress, a key mechanism underlying its inflammatory and toxic effects, which can activate cellular stress responses including upregulation of heat shock protein 70 (Hsp70). Previous studies in lung and liver epithelial cells have shown that Hsp70 and oxidative stress play pivotal roles in triggering autophagy. This study establishes the critical role of the Hsp70-reactive oxygen species (ROS)-autophagy axis in regulating cellular responses to PCP exposure in human alveolar (A549) and liver carcinoma (HepG2) epithelial cells. Our research elucidated the molecular mechanisms underlying PCP's cellular effects, demonstrating that its exposure resulted in increased expression of autophagy-related proteins (Beclin-1, LC3B, ATG12, and ATG16), subunits of NADPH oxidase (NCF-1, NCF-2, NOX2, and Rac), and antioxidant proteins (SOD and GPx) in both lung and liver cell types. Notably, PCP augmented the interaction between Hsp70 and the autophagy regulator Beclin-1. Pretreatment with the ROS inhibitor N-acetylcysteine or Hsp70 knockdown markedly reversed PCP-induced responses. Our in-silico protein–protein docking analysis and molecular dynamics simulation studies revealed enhanced interactions and/or stable confirmations maintained throughout the simulations for TLR4-Hsp70 and Hsp70-Beclin-1 complexes in the presence of PCP. These findings provide a strong foundation for future studies, employing in vivo experimental models and human populations to identify promising targets for PCP-induced toxicity and cellular injury. As a result, these findings may have far-reaching implications for public health and environmental policy, ultimately leading to the identification of biomarkers and the development of more effective interventions for environmentally induced toxicity and diseases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Topography-Induced TKE Budget Behavior Over an Amazon Forest

Data from three different heights (35, 50 and 81 m) of one of the ATTO Project towers in the Amazon forest were used to calculate the TKE (turbulence kinetic energy) budget and some other statistics within the RSL (roughness sublayer). The statistical analyses were carried out for unstable and stable cases. The vertical transport and the vertical advection terms do not explain the imbalances found in the TKE budget, highlighting the impor- tance of horizontal terms in complex terrain. Indeed, patterns independent of the time of day were found for vertical transport, horizontal and vertical velocity skewness, the mean vertical component velocity, and momentum and vertical TKE flux. Here, the highest values of mean vertical velocity were observed from the direction where there is a valley near the tower, which may indicate a topographical effect. Monin-Obukhov Similarity Theory (MOST) predictions for the mean velocity gradient, strictly applicable only in the inertial sublayer (ISL), gave reasonable predictions at 35 and 50 m, but dimensionless dissipation rates displayed large scatter. At 81 m, there were clear deviations from MOST for the two functions, disclosing the effect of topography. In addition, a neutral LES (large-eddy simu- lation) study showed good agreement with tower data for the same general wind direction.

ATTO project↗
DOE OSTI2022

A level-of-details framework for representing occupant behavior in agent-based models

We report agent-based modeling is an advanced computational technique capable of representing complex and dynamic processes of human behavior in building performance simulation. Though the agent-based approach supports diverse applications concerning human behavior modeling within the built environment, there is no consensus on the optimal amount of information or level of granularity needed for occupant information representation. This paper attempts to formalize the level of details (LoD) needed for occupant behavior representation in agent-based environments. A novel framework, grounded on the concept of LoD, is proposed to select the required details in representing occupants in agent-based models. Ten attributes related to occupants' presence, movement, behavioral processes, and repertoire are considered to define the LoD. The framework identifies use case parameters as the guiding principle and allows a hybrid approach for selecting varying degrees of occupant attributes to serve the purpose of simulation. A discussion on the pertinence of different occupant behavior LoDs in relation to the desired objective and simulation context is also presented. The study intends to support the occupant behavior research by advancing agent-based occupant modeling in building performance simulation.

42 ENGINEERING↗
DOE OSTI2021

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗
DOE OSTI2024

Elucidate the molecular basis of ampholytic chitosan as a high-performance cryoprotectant to myosin denaturation: The importance of saccharide charges

The uses of charged poly/oligo-saccharides enable the retarding of protein denaturation against various environmental stresses during food storage and manufacturing. However, at subzero temperatures, the molecular basis of such stabilization behaviors, i.e., cryoprotections, remain less explored. Here, in this study, we introduced an ampholytic saccharide, carboxymethyl chitooligosaccharide (CMCO) that effectively inhibited the freezing-induced myosin denaturation. The in-depth cryoprotective mechanism was systematically investigated by using molecular dynamic simulation and multispectral characterizations. Results showed that CMCO may interact with myosin through hydrogen bonding and electrostatic interactions, which caused the expelling of water at protein surfaces and the reduced conformational flexibility of myosin molecules. Due to this water replacement event, both secondary and ternary structures of myosin became freezing-resistant, leading to the inhibited protein aggregations and retained functionalities, such as solubility, Ca 2+ -ATPase activity, and gelling properties. Moreover, cryoprotective behaviors of CMCO were charge-dependent. CMCO with a higher degree of carboxymethyl substitution (DS: 1.2) was inclined to bind and stabilize myosin molecules better than the low-DS (DS: 0.8) one, even though both outperformed other cryoprotective saccharides. Therefore, this investigation not only introducing a high-performance myosin cryoprotectant, but also elaborated the cryoprotective mechanism of ampholytic saccharides.

59 BASIC BIOLOGICAL SCIENCES↗