Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational framework”

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

A molecular dynamics study of domain switching in BiFeO 3 nanofilm under DC electric field

The performance of ferroelectric materials is deeply rooted in the domain switching process triggered by an applied electric field. To understand the effect of controlling factors on the domain switching characteristics, nanoscale studies are essential. In this work, an atomistic-scale computational framework is developed to study the domain switching process of bismuth ferrite (BiFeO 3 ) nanofilms. Through a series of molecular dynamics simulations, the effects of duration and magnitude of an electric field, domain size and nanofilm thickness on the domain switching behavior are investigated. The evolution of domains under different magnitudes and durations of the applied electric field is studied, and the results of unswitched and switched cases are summarized in a switching map. Analysis of various nanofilms reveals that by increasing the film thickness, the switching time is increased. Additionally, the nanofilms with different domain sizes and configurations are examined, and the minimum domain size for complete reversal switching is obtained. Our results are helpful to predict and tune the switching process in ferroelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Screening green solvents for multilayer plastic film recycling processes

Multilayer (ML) plastic films are essential packaging materials that help protect products from diverse external factors; however, only 5% of all ML films are recycled in the United States. Solvent-based technologies are a promising alternative for recycling ML films because they enable recovery of constituent polymer resins. For example, the Solvent Targeted Recovery and Precipitation (STRAPTM) process sequentially dissolves and separates polymer components using a series of targeted solvent washes. A crucial design aspect of this process is the impact of selected solvents on human health and on the environment. Here, this work introduces a computational framework that integrates molecular modeling, process modeling, techno-economic analysis (TEA), and life-cycle analysis (LCA) to quickly screen green solvents for solvent-based ML recycling processes. Initial screening for solvents based on selectivity is performed by estimating temperature-dependent solubilities using molecular-scale models. Subsequent screening uses basic estimates of energy use and octanol-water partition coefficients (logP) as key measures of health, safety, and environmental hazards. Detailed process modeling, TEA, and LCA are used on a reduced set of promising solvents identified in early screening steps to more accurately determine how solvent selection and associated operating conditions impact overall economics and environmental impacts. The framework is used for the identification of green solvents (from a database of 1,000 solvents) that separate an industrial ML film composed of polyethylene (PE), ethylene vinyl alcohol (EVOH), and polyethylene terephthalate (PET). Our analysis shows the effectiveness of the framework and reveals fundamental trade-offs between solvent greenness, solubility, and economics. Our work emphasizes the importance of taking a holistic systems view during solvent design and aims to inform the development of new processes for ML film recycling and the identification of new ML films that are easier to recycle.

economics↗
DOE OSTI2025

A deep learning and finite element approach for exploration of inverse structure–property designs of lightweight hybrid composites

Hybrid composites have important applications, such as high-performance and lightweight materials in aerospace and automotive industries. Hybrid composites utilize the synergy of diverse fillers to achieve desired material properties, but usually have more complicated microstructures. While topology optimization can optimize a particular property, designing hybrid composites for customized mechanical performances, e.g. full-range stress–strain curve, remains challenging. Here, a computational framework that integrated finite element analysis (FEA) and artificial intelligence (AI) methods of Conditional Generative Adversarial Networks (cGAN) deep learning and transfer learning was developed to establish inverse structure–property relationships and design tailor-made hybrid composites. Based on FEA-generated datasets of hybrid fiber-particle–matrix microstructures and their corresponding full-range stress–strain curves, a cGAN architecture was trained to generate tailored microstructures and establish structure–property relationships. Similarity in microstructural features and well-matched stress–strain curves based on the AI-generated composites were achieved. In conclusion, transfer learning was used to expand the pre-trained model for designing different materials systems.

Hybrid composites↗
DOE OSTI2025

Resolving Lonsdaleite's decade-long controversy: Atomistic insights into a metastable diamond polymorph

Lonsdaleite, a theoretically proposed hexagonal diamond polymorph, has remained at the center of a five-decade scientific controversy since its 1967 identification. While some studies claim it exhibits superior hardness through compression-induced structural changes, others contend it is merely a stacking-faulted cubic diamond. Meteoritic samples and synthetic preparations have yielded conflicting evidence, with even advanced characterisation techniques like XRD and TEM failing to provide definitive proof. In this work, we employ first-principles density functional theory (DFT) and molecular dynamics (MD) simulations to generate unambiguous theoretical fingerprints through XRD, Raman, and SAED patterns that distinguish true Lonsdaleite from cubic diamond and its defective variants. Our atomistic approach quantifies the thermodynamic metastability of Lonsdaleite under realistic pressure-temperature conditions, reveals distinct spectral signatures through simulated Raman and resolves the structural ambiguity through generalised stacking fault energy analysis. By establishing clear criteria for definitive identification, this study provides long-awaited clarity to the Lonsdaleite debate while offering a robust computational framework for characterising metastable carbon phases in meteoritic, synthetic and industrial materials.

DFT↗
DOE OSTI2022

Optimal high-throughput virtual screening pipeline for efficient selection of redox-active organic materials

As global interest in renewable energy continues to increase, there has been a pressing need for developing novel energy storage devices based on organic electrode materials that can overcome the shortcomings of the current lithium-ion batteries. One critical challenge for this quest is to find materials whose redox potential (RP) meets specific design targets. In this study, we propose a computational framework for addressing this challenge through the effective design and optimal operation of a high-throughput virtual screening (HTVS) pipeline that enables rapid screening of organic materials that satisfy the desired criteria. Starting from a high-fidelity model for estimating the RP of a given material, we show how a set of surrogate models with different accuracy and complexity may be designed to construct a highly accurate and efficient HTVS pipeline. We demonstrate that the proposed HTVS pipeline construction and operation strategies substantially enhance the overall screening throughput.

36 MATERIALS SCIENCE↗
DOE OSTI2026

scPlantAnnotate: an accurate and robust transformer-based model for plant cell type annotation

Accurate cell type annotation remains a major bottleneck in plant single-cell RNA sequencing (scRNA-seq), where existing tools are often adapted from animal studies and perform sub-optimally on plant data. The lack of plant-specific computational frameworks limits the construction of plant cell atlases and downstream biological discovery. We develop and evaluate scPlantAnnotate, a Transformer-based reference annotation framework tailored for plant scRNA-seq data, and benchmark it against state-of-the-art deep learning and conventional methods across multiple plant species. Species-specific scPlantAnnotate models were trained using curated datasets from Arabidopsis thaliana, Zea mays, Oryza sativa, and Glycine max. We compared scPlantAnnotate with leading baselines under both standard random-split evaluation and a more stringent leave-one-dataset-out setting, which tests robustness to completely unseen datasets and tissue types. scPlantAnnotate consistently outperforms existing approaches across all four species under random-split evaluation. In the leave-one-dataset-out setting for A. thaliana, where performance drops markedly for all methods due to strong batch effects and dataset heterogeneity, scPlantAnnotate nonetheless achieves the highest Accuracy, Macro-F1, Balanced Accuracy, and Macro-AUROC on average and ranks first on most held-out datasets. These results demonstrate improved robustness to dataset shifts, a critical yet underexplored challenge in plant scRNA-seq analysis. A freely accessible web server enables users to annotate their own datasets using pretrained models. scPlantAnnotate provides a plant-specific, Transformer-based framework for single-cell annotation that delivers state-of-the-art performance and enhanced robustness to unseen datasets. By addressing limitations of existing tools and enabling scalable reference-based annotation, scPlantAnnotate supports the development of comprehensive plant cell atlases and facilitates broader use of single-cell genomics in plant biology.

Bioinformatics↗
DOE OSTI2024

Computational modeling of coupled mechanical damage and electrochemistry in ternary oxide composite electrodes

Performance degradation of ternary layered oxide cathodes largely originates from their loss of structural integrity in cyclic usage. Mechanical damage, such as intergranular fracture of the active particles, is not only a mechanical cleavage process but also interferes with electrochemical kinetics such as infiltration of liquid electrolyte, surface corrosion of the constituent primary particles, and may eventually isolate the primary grains from the electron conducting network. Here, in this work, we develop a computational framework that integrates electrochemistry of a LiNi x Mn y Co 1−x−y O 2 (NMC) composite cathode with mechanical damage of the active particles. To fully examine the intricate chemomechanical behavior of the electrode, we evaluate the effects of the anisotropic material properties, the influence of mechanical potential on Li transport, and the concurrent intergranular fracture and electrolyte penetration along the grain boundaries upon multiple cycles. Electrolyte infiltration benefits capacity retention but aggravates further mechanical damage by corrosion. Structural failure mostly occurs in the first charging due to the anisotropic mechanical strain between the primary grains, while the resulting damage remains stable in the later few cycles. The results are consistent with experimental observations and the integration of electrochemistry and mechanical failure enables a step further understanding of the complex mechanism of battery degradation.

Battery degradation↗
DOE OSTI2025

Seismic response of vertical dry storage casks under three-dimensional earthquake motions

Ensuring the long-term seismic safety of dry storage casks (DSCs) is becoming increasingly critical as these systems evolve from temporary to de facto permanent repositories for spent nuclear fuels. Traditional seismic soil–structure interaction (SSI) assessment methods use one-dimensional deconvolution or simplified boundary conditions to model incident waves. Although computationally appealing, simplifying assumptions may alter the seismic risk by neglecting the full complexity of three-dimensional (3D) wave propagation effects. To address this challenge, this paper introduces a novel high-fidelity computational framework that leverages the Domain Reduction Method (DRM) with perfectly matched layers (PML) to accurately transfer complex, 3D seismic wavefields from regional-scale fault-rupture simulations into local-scale finite element models of DSCs. Using broadband, physics-based ground motions from a generic M w 7.0 strike-slip event, both single-cask and multi-cask configurations were investigated under near- and far-field conditions. Emphasis is placed on capturing complex SSI, spatial variability in the ground motion, and nonlinear phenomena such as cask rocking and sliding. Numerical results demonstrate that near-field conditions, where forward directivity and fling-step effects dominate, lead to significantly higher DSC rocking and sliding. Far-field cases, by contrast, generally exhibit modest responses. Incorporating SSI tends to amplify or alter DSC response spectra and introduce response variability, which underscores the need for site-specific evaluations and robust modeling approaches to ensure the seismic integrity of DSCs in interim spent fuel storage installations.

Das, Tonmoy↗
DOE OSTI2023

Stochastic multiscale modeling for quantifying statistical and model errors with application to composite materials

This paper provides a coherent and efficient computational framework for stochastic multiscale analysis of material systems in the presence of parametric uncertainties and modeling errors. Uncertainty in those model parameters that are not deduced as upscaled quantities is attributed to an uncertainty “germ”. While such parameters can appear at any scale, they are predominant at the finest analysis scale. Additional uncertainties stemming from statistical estimation, attributed to lack of data and model error, are associated with each submodel contributing to the multiscale system. Here, a robust and efficient framework based on a generalized extended polynomial chaos expansion (gEPCE) is proposed to simultaneously propagate all these uncertainties in order to provide a probabilistic representation of specific quantities of interest (QoI). We characterize the full probability distribution of the QoI and the uncertainty in the failure probability pertaining to its tails. By combining gEPCE with kernel density estimation (KDE) and directional derivatives, we construct sensitivity measures that connect these statistical metrics of QoI to the various sources of uncertainty to assess their individual and combined impacts. An illustrative problem featuring three-point bending of a composite beam is investigated to demonstrate the presented approach.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Thermomechanical analysis and modeling of a high-temperature light trapping planar cavity solar receiver

The development of durable particle-based high-temperature solar receivers is critical for advancing concentrating solar-thermal (CST) technologies to enable high-efficiency power generation and industrial process heat. Here, this study presents a computational framework to evaluate the thermomechanical performance of a proposed enclosed light-trapping planar cavity receiver designed for particle-based thermal energy systems. The receiver incorporates absorptive cavities and fluidized particle-bed channels to enhance heat capture and reduce thermal losses. Finite element analysis (FEA) is employed to assess stress, strain, and creep-fatigue behavior under concentrated solar flux using realistic thermal boundary conditions derived from coupled system models and experimental assembly parameters. The analysis investigates the influence of particle-to-wall heat transfer coefficients (HTC) ranging from 800 to 1800 W/m 2 .K on the thermomechanical response of six candidate high-temperature alloys: Alloy 740H, Alloy 282, Alloy 617, 316H, Alloy 230, and 800H. Results show that increasing HTC reduces thermal gradients, leading to lower stresses and strains and extended minimum predicted creep life. While all materials satisfy fatigue life requirements under the investigated conditions, significant differences in creep resistance are observed. Alloy 740H consistently exhibits the longest minimum predicted creep life and the most favorable durability margins, followed by Alloy 282, with the remaining materials showing reduced creep resistance under identical loading. The reported creep lives are conservative lower-bound estimates intended for comparative material evaluation. This framework highlights the critical roles of material selection and geometry optimization in improving mechanical durability and reliability of solar-thermal receivers, forming a foundation for future experimental validation and design optimization.

14 SOLAR ENERGY↗
DOE OSTI2020

Logic-based analysis of gene expression data predicts association between TNF, TGFB1 and EGF pathways in basal-like breast cancer

For breast cancer, clinically important subtypes are well characterized at the molecular level in terms of gene expression profiles. In addition, signaling pathways in breast cancer have been extensively studied as therapeutic targets due to their roles in tumor growth and metastasis. However, it is challenging to put signaling pathways and gene expression profiles together to characterize biological mechanisms of breast cancer subtypes since many signaling events result from post-translational modifications, rather than gene expression differences. We designed a logic-based computational framework to explain the differences in gene expression profiles among breast cancer subtypes using Pathway Logic and transcriptional network information. Pathway Logic is a rewriting-logic-based formal system for modeling biological pathways including post-translational modifications. Our method demonstrated its utility by constructing subtype-specific path from key receptors (TNFR, TGFBR1 and EGFR) to key transcription factor (TF) regulators (RELA, ATF2, SMAD3 and ELK1) and identifying potential association between pathways via TFs in basal-specific paths, which could provide a novel insight on aggressive breast cancer subtypes. Lastly, codes and results are available at http://epigenomics.snu.ac.kr/PL/.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Charge Transport in Solvated Donor–Acceptor Functionalized Peptoids: Molecular Dynamics and Rate Theory

Scalable solar-energy conversion requires photoactive materials that combine the efficiency of natural photosynthetic systems with the stability and processability needed for practical applications. Achieving reliable charge transport in soft, self-assembled organic materials remains challenging, as structural fluctuations and environmental effects strongly influence charge-transfer (CT) rates. Here, we present a broadly applicable computational framework for evaluating CT rates in the condensed phase, combining Fermi’s golden rule rate theory with inputs from all-atom molecular dynamics (MD) simulations and first-principles electronic-structure calculations. The approach does not rely on system-specific parametrization and is applicable to a wide range of soft and disordered materials. We demonstrate the applicability and usefulness of the framework on redox-active peptoids functionalized with iron–porphyrin (Fe–P) complexes, a bioinspired platform with programmable donor–acceptor units and tunable three-dimensional organization. The calculated CT rates exhibit strong sensitivity to molecular conformation, with variations spanning several orders of magnitude. This dependence is shown to arise from the pronounced variation in diabatic electronic coupling with the relative orientations and separations of the Fe–P complexes across the conformational ensemble. The framework provides a consistent route for connecting atomistic structure to CT kinetics in the condensed phase and enables analysis of structure–rate relations in organic semiconducting systems.

Charge transfer↗
DOE OSTI2025

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗
DOE OSTI2021

Comparative Studies of the Structural and Transport Properties of Molten Salt FLiNaK Using the Machine-Learned Neural Network and Reparametrized Classical Forcefields

Despite surging interest in molten salt reactors and thermal storage systems, knowledge of the physicochemical properties of molten salts are still inadequate due to demanding experiments that require high temperature, impurity control, and corrosion mitigation. Therefore, the ability to predict these properties for molten salts from first-principles computations is urgently needed. Herein, we developed and compared a machine-learned neural network force field (NNFF) and a reparametrized rigid ion model (RIM) for a prototypical molten salt LiF–NaF–KF (FLiNaK). We found that NNFF was able to reproduce both the structural and transport properties of the molten salt with first-principles accuracy and classical-MD computational efficiency. Furthermore, the correlation between the local atomic structures and the dynamics was identified by comparing with RIMs, suggesting the significance of polarization of anions implicitly embedded in the NNFF. Furthermore, this work demonstrated a computational framework that can facilitate the screening of molten salts with different chemical compositions, impurities, and additives, and at different thermodynamic conditions suitable for the next-generation nuclear reactors and thermal energy storage facilities.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Deciphering the Cofilin Oligomers via Intermolecular Disulfide Bond Formation: A Coarse-Grained Molecular Dynamics Approach to Understanding Cofilin’s Regulation on Actin Filaments

Cofilin, a key actin-binding protein, orchestrates the dynamics of the actomyosin network through its actin-severing activity and by promoting the recycling of actin monomers. Recent experimental work suggests that cofilin also forms functionally distinct oligomers through thiol post-translational modification (PTM) that encourages actin nucleation and assembly. Despite these advances, the structural conformations of cofilin oligomers that modulate actin activity remain elusive because there are combinatorial ways to oxidize thiols in cysteines to form disulfide bonds rapidly. This study employs molecular dynamics simulations to investigate human cofilin 1 as a case study for exploring cofilin dimers via disulfide bond formation. Using the free energy profiling, our simulations unveil a range of probable cofilin dimer structures not represented in current Protein Data Bank entries. These candidate dimers are characterized by their distinct population distributions and relative free energies. Of particular note is a dimer featuring an interface between cysteines 139 and 147 residues, which demonstrates stable free energy characteristics and intriguingly symmetrical geometry. In contrast, the experimentally proposed dimer structure exhibits a less stable free energy profile. Here, we also evaluate frustration quantification based on the energy landscape theory in the protein-protein interactions at the dimer interfaces. Notably, the 39-39 dimer configuration emerges as a promising candidate for forming cofilin tetramers, as substantiated by frustration analysis. Additionally, docking simulations with actin filaments further evaluate the stability of these cofilin dimer-actin complexes. Our findings thus offer a computational framework for understanding the role of thiol post-translational modification cofilin proteins in regulating oligomerization, and the subsequent cofilin-mediated actin dynamics in the actomyosin network.

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

Wetting of a Hydrophobic Surface: Far-IR Action Spectroscopy and Dynamics of Microhydrated Naphthalene

The interaction of water and polycyclic aromatic hydrocarbons is of fundamental importance in areas as diverse as materials science and atmospheric and interstellar chemistry. The interplay between hydrogen bonding and dipole–π interactions results in subtle dynamics that are challenging to describe from first principles. Here, we employ far-IR action vibrational spectroscopy with the infrared free-electron laser FELIX to investigate naphthalene with one to three water molecules. We observe diffuse bands associated with intermolecular vibrational modes that serve as direct probes of the loose binding of water to the naphthalene surface. These signatures are poorly reproduced by static DFT or Møller–Plesset computations. Instead, a rationalization is achieved through Born–Oppenheimer Molecular Dynamics simulations, revealing the active mobility of water over the surface, even at low temperatures. Therefore, our work provides direct insights into the wetting interactions associated with shallow potential energy surfaces while simultaneously demonstrating a solid experimental–computational framework for their investigation.

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

Molecular Dynamics Investigation of Oil Wetting on Synthetic Polymer Substrates

Understanding the interaction of polymer surfaces with nonpolar, low surface tension liquids, or whether a substrate is oleophobic versus oleophilic, is critical for applications ranging from antifingerprint coatings to oil–water separation membranes and oil spill remediation. Despite its technological importance, the molecular mechanisms underpinning polymer–oil wetting are not well characterized. Here, we employ molecular dynamics simulations to investigate the behavior of n-hexadecane in contact with chemically distinct polymer surfaces, spanning eight constitutional unit chemistries as well as amorphous and crystalline morphologies. This permits a critical examination of both thermodynamic and dynamic descriptors of oleophobicity, including oil contact angle, dewetting free energy, interfacial diffusivity, and a proposed “ghost probe energy,” which does not require explicit simulation of oil–polymer interactions. We find that the oil contact angle does not reliably distinguish oleophobic behavior across polymer chemistries, whereas other metrics provide clearer and more consistent differentiation. Analysis of the results reveals that polymer–oil wetting behavior is primarily governed by interfacial van der Waals interactions and modulated by surface flexibility and morphology. Collectively, this work establishes a computational framework for characterizing oil wetting and provides additional insight into what molecular-level factors dictate trends in oleophobicity.

lipids↗
DOE OSTI2022

Te/CdTe and Al/CdTe Interfacial Energy Band Alignment by Atomistic Modeling

A synergistic approach that incorporates first-principles atomistic modeling with numerical device simulations is used to systematically evaluate the role of heterointerfaces within metal-chalcogenide-based photovoltaic technologies. Two interfaces involving either a tellurium back contact or aluminum back electrode combined with a cadmium telluride absorber layer within cadmium-telluride-based thin-film solar cells are investigated on an atomic scale to determine the mechanisms contributing to variations in device performance. Electronic structures and predicted charge transport behavior with respect to cadmium and tellurium termination of the absorber layer are studied along the polar oriented CdTe{111} facets. The computational methodology reveals a noticeable contrast between the Schottky barrier forming Al/CdTe interface versus the Type I Te/CdTe heterojunction. Greater band bending features are exhibited by the cadmium termination as opposed to the tellurium termination for each interface case. Subsequent device modeling suggests that 3.6% higher photovoltaic conversion efficiency is achievable for the cadmium termination relative to the tellurium termination of the Te/CdTe interface. Based strictly on an idealistic representation, both interface models show the importance of atomic-scale interfacial properties for cadmium telluride solar cell device performance with their bulk properties being validated in comparison to published experimental data. Here, the synergistic approach offers a suitable method to analyze solar cell interfaces through a predictive computational framework for the engineering and optimization of metal-chalcogenide-based thin-film photovoltaic technologies.

Green’s function↗