SEARCH · Search NASA

Results for “complex system modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Models and Algorithms for Equilibrium Analysis of Mixed-Material Nucleic Acid Systems

Dynamic programming algorithms within the NUPACK software suite enable analysis of equilibrium base-pairing properties for complex and test tube ensembles containing arbitrary numbers of interacting nucleic acid strands. Currently, calculations are limited to single-material systems that are either all-RNA or all-DNA. Here, to enable analysis of mixed-material systems that are critical for modern applications in vitro, in situ, and in vivo, we develop physical models and dynamic programming algorithms that allow the material of the system to be specified at nucleotide resolution. Free energy parameter sets are constructed for both RNA/DNA and RNA/2'OMe-RNA mixed-material systems by combining available empirical mixed-material parameters with single-material parameter sets to enable treatment of the full complex and test tube ensembles. New dynamic programming recursions account for the material of each nucleotide throughout the recursive process. For a complex with N nucleotides, the mixed-material dynamic programming algorithms maintain the O(N 3 ) time complexity of the single-material algorithms, enabling efficient calculation of diverse physical quantities over complex and test tube ensembles (e.g., complex partition function, equilibrium complex concentrations, equilibrium base-pairing probabilities, minimum free energy secondary structure(s), and Boltzmann-sampled secondary structures) at a cost increase of roughly 2.0-3.5×. The results of existing single-material algorithms are exactly reproduced when applying the new mixed-material algorithms to single-material systems. Accuracy is significantly enhanced using mixed-material models and algorithms to predict RNA/DNA and RNA/2'OMe-RNA duplex melting temperatures from the experimental literature as well as RNA/DNA melt profiles from new experiments. In conclusion, mixed-material analyses can be performed online using the NUPACK web app (www.nupack.org) or locally using the NUPACK Python module.

2′OMe-RNA↗

DOE OSTI2025

Neutrons reveal the dynamics of leaf thylakoids in living plants

The study is the first known exploration of photosynthetic membranes dynamics in living plants by high resolution quasielastic neutron scattering spectroscopy. We investigated the mobility and flexibility of thylakoid membranes in common duckweed (Landoltia punctata) and identified dynamics across various length scales corresponding to individual membranes and membranes stack. We employed classical models typically used to study lipid bilayers to characterize the undulation modes and rigidity of the membranes and reveal how structural variations influence the observed complex dynamics. Our findings show that the stacks of thylakoids in duckweed behave as rigid systems, exhibiting an effective bending coefficient in the lower range associated with surfactant membranes. In contrast, the single thylakoid leaflets display greater apparent flexibility and are well situated within the bi-continuous surfactant phase dynamics. While our observations enhance the understanding of the intricate architecture and mobility of photosynthetic cellular machinery, they also highlight the limitations of applying ideal lipid membranes models to describe complex biological systems. This work opens more questions and the need for further investigations across extended length and time scales, as well as the importance of rigorous sample preparation and experimental control.

Applied physics↗

DOE OSTI2024

Mechanisms of Three-Dimensional Solid-Phase Epitaxial Crystallization of Strontium Titanate

Strontium titanate (SrTiO 3 , STO) is a complex metal oxide with a cubic perovskite crystal structure. Due to its easily described and understood crystal structure in the cubic phase, STO is an ideal model system for exploring the mechanistic details of solid-phase epitaxy (SPE) in complex oxides. SPE is a crystallization approach that aims to guide crystal growth at low homologous temperatures to achieve targeted microstructures. Beyond planar thin films, SPE can also exploit the addition of a chemically inert, noncrystallizing, amorphous obstacle in the path of crystallization to generate complex three-dimensional structures. The introduction of this mask fundamentally alters the SPE process, inducing a transition from two- to three-dimensional geometries and from vertical to lateral crystal growth under the influence of the crystal/mask/amorphous boundary. Using a combination of molecular dynamics simulations and experiments, we identify several unique phenomena in the nanoscale growth behaviors in both conventional (unmasked) and masked SPE. Examining conventional SPE of STO, we find that crystallization at the interface is strongly correlated to, and potentially driven by, density fluctuations in the region of the amorphous STO near the crystalline/amorphous interface with a strong facet dependence. In the masked case, we find that the crystalline growth front becomes nonplanar near contact with the mask. We also observe a minimum vertical growth requirement prior to lateral crystallization. Both phenomena depend on the relative bulk and interfacial free energies of the three-phase (crystal/mask/amorphous) system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Efficient mapping between void shapes and stress fields using Deep Convolutional Neural Networks with sparse data

Establishing fast and accurate structure-to-property relationships is an important component in the design and discovery of advanced materials. Physics-based simulation models like the finite element method (FEM) are often used to predict deformation, stress, and strain fields as a function of material microstructure in material and structural systems. Such models may be computationally expensive and time intensive if the underlying physics of the system is complex. This limits their application to solve inverse design problems and identify structures that maximize performance. In such scenarios, surrogate models are employed to make the forward mapping computationally efficient to evaluate. However, the high dimensionality of the input microstructure and the output field of interest often renders such surrogate models inefficient, especially when dealing with sparse data. Deep convolutional neural network (CNN) based surrogate models have shown great promise in handling such high-dimensional problems. In this paper, a single ellipsoidal void structure under a uniaxial tensile load represented by a linear elastic, high-dimensional and expensive-to-query, FEM model. We consider two deep CNN architectures, a modified convolutional autoencoder framework with a fully connected bottleneck and a UNet CNN, and compare their accuracy in predicting the von Mises stress field for any given input void shape in the FEM model. Additionally, a sensitivity analysis study is performed using the two approaches, where the variation in the prediction accuracy on unseen test data is studied through numerical experiments by varying the number of training samples from 20 to 100.

surrogate modeling; convolutional neural networks;↗

DOE OSTI2026

Physics vs structure: A systematic benchmark of learning strategies for multi-zone building thermal dynamics

Recent advances in physics-informed and data-driven machine learning promise improved thermal models for advanced building control, yet there is limited quantitative evidence on when added physics structure and architectural complexity are beneficial. Here, this work presents a systematic benchmark of five representative system identification methods for modeling multi-zone building thermal dynamics: linear state-space models, multi-layer perceptrons, neural state-space models, neural ordinary differential equations, and physically-consistent neural networks. The methods are evaluated across multiple data regimes and zone coupling strategies. Using a high-fidelity multi-zone commercial building emulator, we examine short-term and long-term prediction accuracy, computational efficiency, and ease of development. Our results reveal critical trade-offs between prediction performance, model complexity, and physical consistency. We demonstrate that decoupled, nonlinear black-box models consistently outperform coupled physics-constrained architectures in both predictive accuracy and out-of-distribution robustness in majority of the test cases for the building type considered in the study. Our findings quantify the cost of complexity in building thermal modeling and provide concrete, actionable, scenario-based guidelines for selecting model classes for control-oriented applications.

Building thermal modeling↗

DOE OSTI2024

Analysis of AP1000 Small-Break Loss-of-Coolant Accident Using Reactor Transient Simulator

The Westinghouse Electric Company’s Advanced Passive Reactor (AP1000) is characterized by the incorporation of passive safety systems (PSSs) designed to ensure core cooling during transient events. The assessment of PSSs requires evaluation of their performance through a combination of experiments and simulations employing various thermal-hydraulic codes. In addition, detailed evaluation of PSSs for a specific reactor system transient analysis such as loss-of-coolant-accident analysis supports understanding representative integral effects test facility development and the further evolution model development and assessment process. Developing a reactor system code is a complex and time-consuming process that requires significant engineering expertise and effort. It can take several months to even years to complete in the early stages of reactor system design and analysis. However, this process can be expedited through the use of transient simulator models for similar reactor systems, which can be used for lesson learning and training purposes. This study uses the Personal Computer Transient Analyzer (PCTRAN) code. The main advantage of PCTRAN is its ease of use and ability to run faster than real time. This study presents the results obtained for a small-break loss-of-coolant accident (SBLOCA) for two breaks using the full version (licensed) of PCTRAN. The purpose of this investigation is to evaluate the overall system behavior during the postulated SBLOCA event as well as assess the capability of the PCTRAN code to reproduce the system response during transient events. The obtained results were compared with the Westinghouse NOTRUMP system code. The PCTRAN code proved to be reliable in predicting the qualitative behavior of the system in both transient cases. As for the system response, it was found that it is contingent on the activation time of the PSSs. The differences in reactor coolant system pressure between the two codes were attributed to the critical flow model and simplification of mass and energy balance. Despite PCTRAN’s limitations, it can still provide a reasonable prediction of various reactor parameters such as pressure, mass flow rate, and void fraction during a SBLOCA scenario. It is worth noting that PCTRAN currently employs a bulk approach similar to that of the Modular Accident Analysis Program (MAAP) and MELCOR codes. However, the upcoming version of PCTRAN will include an artificial intelligence–based detection and accident prevention system, as well as different models for different reactor components. Consequently, PCTRAN has the potential to be upgraded to match the system thermal-hydraulic codes of the U.S. Nuclear Regulatory Commission and become more widely used in cybersecurity to safeguard nuclear power plants from cyberattacks.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

Deep potential molecular dynamics simulations of low-temperature plasma-surface interactions

Machine learning approaches to potential generation for molecular dynamics (MD) simulations of low-temperature plasma-surface interactions could greatly extend the range of chemical systems that can be modeled. Empirical potentials are difficult to generalize to complex combinations of multiple elements with interactions that might include covalent, ionic, and metallic bonds. This work demonstrates that a specific machine learning approach, Deep Potential Molecular Dynamics (DeepMD), can generate potentials that provide a good model of plasma etching in the Si-Cl-Ar system. Comparisons are made between MD results using DeepMD models and empirical potentials, as well as experimental measurements. Pure Si properties predicted by the DeepMD model are in reasonable agreement with experimental results. Simulations of Si bombardment by Ar + ions demonstrate the ability of the DeepMD method to predict sputtering yields as well as the depth of the amorphous-crystalline interface. Etch yields as a function of flux ratio and ion energy for simultaneous Cl 2 and Ar + impacts are in good agreement with previous simulation results and experiment. Predictions of etch yields and etch products during plasma-assisted atomic layer etching of Si-Cl 2 -Ar are shown to be in good agreement with MD predictions using empirical potentials and with experiment. Finally, good agreement was also seen with measurements for the spontaneous etching of Si by Cl atoms at 300 K. Further, the demonstration that DeepMD can reproduce results from MD simulations using empirical potentials is a necessary condition to future efforts to extend the method to a much wider range of systems for which empirical potentials may be difficult or impossible to obtain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Spatio-temporal multivariate cluster evolution analysis for detecting and tracking climate impacts

Recent years have seen a growing concern about climate change and its impacts. While Earth System Models (ESMs) can be invaluable tools for studying the impacts of climate change, the complex coupling processes encoded in ESMs and the large amounts of data produced by these models, together with the high internal variability of the Earth system, can obscure important source-to-impact relationships. Here, this paper presents a novel and efficient unsupervised data-driven approach for detecting statistically-significant impacts and tracing spatio-temporal source-impact pathways in the climate through a unique combination of ideas from anomaly detection, clustering and Natural Language Processing (NLP). Using as an exemplar the 1991 eruption of Mount Pinatubo in the Philippines, we demonstrate that the proposed approach is capable of detecting known post-eruption impacts/events. We additionally describe a methodology for extracting meaningful sequences of post-eruption impacts/events by using NLP to efficiently mine frequent multivariate cluster evolutions, which can be used to confirm or discover the chain of physical processes between a climate source and its impact(s).

Anomaly detection↗

DOE OSTI2024

Characterization of Deformational Isomerization Potential and Interconversion Dynamics with Ultrafast X-ray Solution Scattering

Dimeric complexes composed of d 8 square planar metal centers and rigid bridging ligands provide model systems to understand the interplay between attractive dispersion forces and steric strain in order to assist the development of reliable methods to model metal dimer complexes more broadly. [Ir 2 (dimen) 4 ] 2+ (dimen = para-diisocyanomenthane) presents a unique case study for such phenomena, as distortions of the optimal structure of a ligand with limited conformational flexibility counteract the attractive dispersive forces from the metal and ligand to yield a complex with two ground state deformational isomers. Here, in this work, we use ultrafast X-ray solution scattering (XSS) and optical transient absorption spectroscopy (OTAS) to reveal the nature of the equilibrium distribution and the exchange rate between the deformational isomers. The two ground state isomers have spectrally distinct electronic excitations that enable the selective excitation of one isomer or the other using a femtosecond duration pulse of visible light. We then track the dynamics of the nonequilibrium depletion of the electronic ground state population—often termed the ground state hole—with ultrafast XSS and OTAS, revealing a restoration of the ground state equilibrium in 2.3 ps. This combined experimental and theoretical study provides a critical test of various density functional approximations in the description of bridged d 8 –d 8 metal complexes. The results show that density functional theory calculations can reproduce the primary experimental observations if dispersion interactions are added, and a hybrid functional, which includes exact exchange, is used.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Automated Calibration for Rapid Optical Spectroscopy Sensor Development for Online Monitoring

An automated platform has been developed to assist researchers in the rapid development of optical spectroscopy sensors to quantify species from spectral data. This platform performs calibration and validation measurements simultaneously. Real-time, in situ monitoring of complex systems through optical spectroscopy has been shown to be a useful tool; however, building calibration models requires development time, which can be a limiting factor in the case of radiological or otherwise hazardous systems. While calibration time can be reduced through optimized design of experiments, this study approached the challenge differently through automation. The ATLAS (Automated Transient Learning for Applied Sensors) platform used pneumatic control of stock solutions to cycle flow profiles through desired calibration concentrations for multivariate model construction. Additionally, the transients between desired concentrations based on flow calculations were used as validation measurements to understand model predictive capabilities. This automated approach yielded an incredible 76% reduction in model development time and a 60% reduction in sample volume versus estimated manual sample preparation and static measurements. The ATLAS system was demonstrated on two systems: a three-lanthanide system with Pr/Nd/Ho representing a use case with significant overlap or interference between analyte signatures and an alternate system containing Pr/Nd/Ni to demonstrate a use case in which broad-band corrosion species signatures interfered with more distinct lanthanide absorbance profiles. Both systems resulted in strong model prediction performance (RMSEP < 9%). Lastly, ATLAS was demonstrated as a tool to simulate process monitoring scenarios (e.g., column separation) in which models can be further optimized to account for day-to-day changes as necessary (e.g., baseline correction). Ultimately, ATLAS offers a vital tool to rapidly screen monitoring methods, investigate sensor fusion, and explore more complex systems (i.e., larger numbers of species).

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

Unitary Qubit Lattice Algorithms for Plasma Physics

This final technical report summarizes research conducted under DOE Award DE-SC0021653 to develop unitary Quantum Lattice Algorithms for modeling electromagnetic wave propagation and scattering in complex media, including plasmas. The project developed and validated quantum-inspired formulations of Maxwell's equations that preserve unitary evolution and can be evaluated on classical high-performance computing systems while providing a foundation for future quantum-computing implementations. Major accomplishments include the development of two- and three-dimensional algorithms for electromagnetic scattering; scalable, distributed-memory implementations demonstrated on the Perlmutter supercomputer; formulations for nonlinear lossless fluid dynamics and cold, lossless, inhomogeneous magnetized plasmas; and an explicit quantum algorithm for a time-discretized Lorenz model. Simulations reproduced a range of characteristic wave phenomena, including transient effects that are not readily apparent in conventional frequency-domain studies, demonstrating the effectiveness of the proposed approach for modeling complex electromagnetic and plasma systems. The work establishes a unified theoretical and computational framework for quantum and quantum-inspired simulation and provides a foundation for future implementation on fault-tolerant quantum systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Specific iron binding to natural sphingomyelin membrane induced by non-specific co-solutes

Sphingomyelin (SPM), a crucial phospholipid in the myelin sheath, plays a vital role in insulating nerve fibers. We hypothesize that iron ions selectively bind to the phosphatidylcholine (PC) template within the SPM membrane under near-physiological conditions, resulting in disruptions to membrane organization. These interactions could potentially contribute to the degradation of the myelin sheath, thereby playing a role in the development of neurodegenerative diseases. We utilized synchrotron-based X-ray spectroscopy and diffraction techniques to study the interaction of iron ions with a bovine spinal-cord SPM monolayer (ML) at the liquid-vapor interface under physiological conditions. The SPM ML serves as a model system, representing localized patches of lipids within a more complex membrane structure. The experiments assessed iron binding to the SPM membrane both in the presence of salts and with additional evaluation of the effects of various ion species on membrane behavior. Grazing incidence X-ray diffraction was employed to analyze the impact of iron binding on the structural integrity of the SPM membrane. Furthermore, our results demonstrate that iron ions in dilute solution selectively bind to the PC template of the SPM membrane exclusively at near-physiological salt concentrations (e.g., NaCl, KCl, KI, or CaCl 2 ) and are pH-dependent. In-significant binding was detected in the absence of these salts or at near-neutral pH with salts. The surface adsorption of iron ions is correlated with salt concentration, reaching saturation at physiological levels. In contrast, multivalent ions such as La 3+ and Ca 2+ do not bind to SPM under similar conditions. Notably, iron binding to the SPM membrane disrupts its in-plane organization, suggesting that these interactions may compromise membrane integrity and contribute to myelin sheath damage associated with neurological disorders.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Uncovering obscured phonon dynamics from powder inelastic neutron scattering using machine learning

The study of phonon dynamics is pivotal for understanding material properties, yet it faces challenges due to the irreversible information loss inherent in powder inelastic neutron scattering spectra and the limitations of traditional analysis methods. In this study, we present a machine learning framework designed to reveal obscured phonon dynamics from powder spectra. Using a variational autoencoder, we obtain a disentangled latent representation of spectra and successfully extract force constants for reconstructing phonon dispersions. Notably, our model demonstrates effective applicability to experimental data even when trained exclusively on physics-based simulations. The fine-tuning with experimental spectra further mitigates issues arising from domain shift. Analysis of latent space underscores the model’s versatility and generalizability, affirming its suitability for complex system applications. Furthermore, our framework’s two-stage design is promising for developing a universal pre-trained feature extractor. This approach has the potential to revolutionize neutron measurements of phonon dynamics, offering researchers a potent tool to decipher intricate spectra and gain valuable insights into the intrinsic physics of materials.

domain adaptation↗

DOE OSTI2025

DG2DAG: Learning Directed Acyclic Graphs from Functional Priors

Physics-based systems-of-systems models are computationally expensive. Reduced graphical models can decrease computational complexity, but may not proffer an end-to-end model from upstream inputs to downstream outputs. We consequently are interested in reducing models on directed graphs to models on a directed acyclic subgraph such that preserves accurate reconstruction of nodes. The consequence is a model with a topological ordering, providing a one-way flow of computation, and a causal interpr

Voronin, Alexey [Sandia National Laboratories (SNL↗

DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗

DOE OSTI2025

Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST

ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.

17 WIND ENERGY↗

DOE OSTI2024

Resilient water infrastructure partnerships in institutionally complex systems face challenging supply and financial risk tradeoffs

Abstract As regions around the world invest billions in new infrastructure to overcome increasing water scarcity, better guidance is needed to facilitate cooperative planning and investment in institutionally complex and interconnected water supply systems. This work combines detailed water resource system ensemble modeling with multiobjective intelligent search to explore infrastructure investment partnership design in the context of ongoing canal rehabilitation and groundwater banking in California. Here we demonstrate that severe tradeoffs can emerge between conflicting goals related to water supply deliveries, partnership size, and the underlying financial risks associated with cooperative infrastructure investments. We show how hydroclimatic variability and institutional complexity can create significant uncertainty in realized water supply benefits and heterogeneity in partners’ financial risks that threaten infrastructure investment partnership viability. We demonstrate how multiobjective intelligent search can design partnerships with substantially higher water supply benefits and a fraction of the financial risk compared to status quo planning processes. This work has important implications globally for efforts to use cooperative infrastructure investments to enhance the climate resilience and financial stability of water supply systems.

Science & Technology - Other Topics↗

DOE OSTI2025

Computing virtual dark-field X-ray microscopy images of complex discrete dislocation structures from large-scale molecular dynamics simulations

Dark-field X-ray microscopy (DFXM) is a novel diffraction-based imaging technique that non-destructively maps the local deformation from crystalline defects in bulk materials. While studies have demonstrated that DFXM can spatially map 3D defect geometries, it is still challenging to interpret DFXM images of the high-dislocation-density systems relevant to macroscopic crystal plasticity. This work develops a scalable forward model to calculate virtual DFXM images for complex discrete dislocation structure(s) (DDS) obtained from atomistic simulations. Our new DDS-DFXM model integrates a non-singular formulation for calculating the local strain from the DDS and an efficient geometrical optics algorithm for computing the DFXM image from the strain field. We apply the model to complex DDS obtained from a large-scale molecular dynamics simulation of compressive loading on single-crystal silicon. Simulated DFXM images exhibit prominent contrast for dislocation features between the multiple slip systems, demonstrating the potential of DFXM to resolve features from dislocation multiplication. In conclusion, the integrated DDS-DFXM model provides a toolbox for DFXM experimental design and image interpretation in the context of bulk crystal plasticity for a range of measurements across shock plasticity and the broader materials science community.

X-ray imaging↗