Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Physics”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effect of Ambient Organic Acids on the Water Structure at ${\rm TiO}_{2}$ Interfaces

A molecular-level understanding of the effects of ambient organic compounds on the wettability of titanium dioxide ( ${\rm TiO}_{2}$ ) surfaces is relevant to many of its energy-related and environmental applications. Herein, we focus on two common atmospheric carboxylic acids, formic and acetic acid, and characterize their adsorption/ desorption at the aqueous interfaces of anatase and rutile ${\rm TiO}_{2}$ using molecular dynamics with an ab initio deep neural network potential. Our simulations show that these acids prefer to be localized in the interfacial water layers close to the ${\rm TiO}_{2}$ surface where they are stabilized by the interaction/exchange of their acid proton with a surface oxygen, rather than chemisorb at the surface Ti sites by displacing the adsorbed water. Notably, these acids make the surface of anatase hydrophobic, whereas the larger fraction of adsorbed water dissociation can offset their effect on rutile. Furthermore, these results provide a picture where carboxylic acids control the wettability of ${\rm TiO}_{2}$ largely through acid-base chemistry at the interface rather than chemisorption on the oxide surface, a finding that can help improve the design of self-cleaning surfaces and photocatalytic devices.

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

Practical choices for space-filling designs

Space-filling designs are now commonly used as a flexible model-free strategy for providing good coverage throughout an input space of interest for a variety of computer and physical experiment design scenarios. Some of the preliminary choices about how to frame the problem and which type of design to use can have a substantial impact on the success or failure of the experiment, and yet how to make these critical choices is often under-emphasized in the literature. In this paper, we explore several of the practical choices required by the experimenter and describe a sequence of steps to help create an ideal design that matches the goals and constraints of the experiment. These choices include the specification of the input space, the scaling of the variables, the degree of uniformity of the design points across the input space, and the space-filling characteristics. In addition, some new tools for defining the weights to implement non-uniform space-filling designs are provided. Finally, the methods are demonstrated with several illustrative examples and a real-world chemical engineering experiment for carbon capture.

42 ENGINEERING↗
DOE OSTI2020

Quantum Spectral Methods for Differential Equations

Recently developed quantum algorithms address computational challenges in numerical analysis by performing linear algebra in Hilbert space. Such algorithms can produce a quantum state proportional to the solution of a d-dimensional system of linear equations or linear differential equations with complexity poly(logd). While several of these algorithms approximate the solution to within ϵ with complexity poly(log(1/ϵ)), no such algorithm was previously known for differential equations with time-dependent coefficients. In this work, we develop a quantum algorithm for linear ordinary differential equations based on so-called spectral methods, an alternative to finite difference methods that approximates the solution globally. Using this approach, we give a quantum algorithm for time-dependent initial and boundary value problems with complexity poly(logd, log(1/ϵ)).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Convex Programs for Minimal-Area Problems

The minimal-area problem that defines string diagrams in closed string field theory asks for the metric of least area on a Riemann surface with the condition that all non-contractible closed curves have length at least 2π. This is an extremal length problem in conformal geometry as well as a problem in systolic geometry. Here, we consider the analogous minimal-area problem for homology classes of curves and, with the aid of calibrations and the max flow-min cut theorem, formulate it as a local convex program. We derive an equivalent dual program involving maximization of a concave functional. These two programs give new insights into the form of the minimal-area metric and are amenable to numerical solution. We explain how the homology problem can be modified to provide the solution to the original homotopy problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Minimal-Area Metrics on the Swiss Cross and Punctured Torus

The closed string theory minimal-area problem asks for the conformal metric of least area on a Riemann surface with the condition that all non-contractible closed curves have length greater than or equal to an arbitrary constant that can be set to one. Through every point in such a metric there is a geodesic that saturates the length condition, and the saturating geodesics in a given homotopy class form a band. The extremal metric is unknown when bands of geodesics cross, as it happens for surfaces of non-zero genus. Here, we use recently proposed convex programs to numerically find the minimal-area metric on the square torus with a square boundary, for various sizes of the boundary. For large enough boundary the problem is equivalent to the “Swiss cross” challenge posed by Strebel. We find that the metric is positively curved in the two-band region and flat in the single-band regions. For small boundary the metric develops a third band of geodesics wrapping around it, and has both regions of positive and of negative curvature. This surface can be completed to provide the minimal-area metric on a once-punctured torus, representing a closed-string tadpole diagram.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

MFNets: data efficient all-at-once learning of multifidelity surrogates as directed networks of information sources

We present an approach for constructing a surrogate from ensembles of information sources of varying cost and accuracy. The multifidelity surrogate encodes connections between information sources as a directed acyclic graph, and is trained via gradient-based minimization of a nonlinear least squares objective. While the vast majority of state-of-the-art assumes hierarchical connections between information sources, our approach works with flexibly structured information sources that may not admit a strict hierarchy. The formulation has two advantages: (1) increased data efficiency due to parsimonious multifidelity networks that can be tailored to the application; and (2) no constraints on the training data—we can combine noisy, non-nested evaluations of the information sources. Finally, numerical examples ranging from synthetic to physics-based computational mechanics simulations indicate the error in our approach can be orders-of-magnitude smaller, particularly in the low-data regime, than single-fidelity and hierarchical multifidelity approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Improved Architectures and Training Algorithms for Deep Operator Networks

Operator learning techniques have recently emerged as a powerful tool for learning maps between infinite-dimensional Banach spaces. Trained under appropriate constraints, they can also be effective in learning the solution operator of partial differential equations (PDEs) in an entirely self-supervised manner. In this work we analyze the training dynamics of deep operator networks (DeepONets) through the lens of Neural Tangent Kernel theory, and reveal a bias that favors the approximation of functions with larger magnitudes. To correct this bias we propose to adaptively re-weight the importance of each training example, and demonstrate how this procedure can effectively balance the magnitude of back-propagated gradients during training via gradient descent. We also propose a novel network architecture that is more resilient to vanishing gradient pathologies. Taken together, our developments provide new insights into the training of DeepONets and consistently improve their predictive accuracy by a factor of 10-50x, demonstrated in the challenging setting of learning PDE solution operators in the absence of paired input-output observations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Green’s Function Perspective on the Nonlinear Density Response of Quantum Many-Body Systems

We derive equations of motion for higher order density response functions using the theory of thermodynamic Green’s functions. We also derive expressions for the higher order generalized dielectric functions and polarization functions. Moreover, we relate higher order response functions and higher order collision integrals within the Martin–Schwinger hierarchy. We expect our results to be highly relevant to the study of a variety of quantum many-body systems such as matter under extreme temperatures, densities, and pressures.

Density Functional Theory↗
DOE OSTI2024

Ghost states and surface structures of the charge density wave kagome metal ScV 6 Sn 6

In this work, we investigate the high-temperature phase of the kagome metal ScV 6 Sn 6 using scanning tunneling microscopy/spectroscopy (STM/S) and density functional theory calculations. STM topographic images of the cleaved sample reveal two distinct surface terminations: flat islands with Sn termination and trenches terminated by kagome layers with Sn as the outermost atomic layer. STS measurements on the Sn-terminated and kagome-terminated surfaces show significant differences, in particular the presence of large density of states near the Fermi level in the former case. Our first-principles calculations reveal that the charge density on the kagome-terminated surface gives rise to “ghost states” which show intensity away from surface atoms, arising due to hybridization of orbitals above the surface. These states can obscure the intrinsic properties of the surface, potentially leading to misattribution of the surface termination. This underscores the need for careful interpretation in STM studies, especially when discerning surface states of localized states. Understanding the surface structure of this versatile quantum material provides essential information for interpreting surface-sensitive experiments, tailoring material properties, engineering interfaces, and controlling stability and reactivity. This knowledge paves the way for further exploration and potential applications of kagome lattice materials in various fields, including quantum computing, topological physics, and advanced electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantifying the atomistic free-volume morphology of materials with graph theory

Here, we introduce a new computational methodology for the identification and characterization of free volume within/around atomistic configurations. This scheme employs a three-stage workflow, by which spheres are iteratively grown inside of voxels, and ultimately converted to planar graphs, which are then characterized via a graph-based order parameter. Our approach is computationally efficient, physically intuitive, and universally transferable to any material system. Validation of our methodology is performed on several sets of materials problems: (1) classification of unique free volumes in various crystal phases, (2) autonomous detection and classification of complex surface defects during epitaxial growth simulations, (3) characterization of free volume defects in metals/alloys, and (4) quantification of the spatio-temporal behavior of nano-scale free volume morphologies as a function of both temperature and free-volume size. Our method accurately identifies and characterizes unique free volumes over a multitude of systems and length scales, indicating its potential for future use in understanding the relationship between free volume morphology and material properties under both static and dynamic conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Calibrating uncertain parameters in melt pool simulations of additive manufacturing

Melt pool scale numerical modeling of additive manufacturing (AM) processes can provide predictive capabilities and theoretical insight into the process-property-structure-performance relationships for AM parts. Despite capabilities of numerical models to solve complex multi-physics problems, it is often important to consider a tradeoff between detailed physics and computational cost. Therefore, sources of uncertainty in both experimental conditions and the parameters needed for modeling require models to be validated against empirical evidence. Here, a method is proposed to calibrate uncertain parameters used in continuum-scale melt pool models for powder bed fusion (PBF) AM. Both a simplified heat transfer model and a heat transfer and fluid flow model were investigated. A surrogate model and Markov chain-based optimization algorithm calibrated melt pool geometry for models within experimental variation of the target melt pool width and depth from the NIST AM-Bench 2018-02 dataset. The melt pool temperature distributions, solidification parameters, and simulated multi-layer solidification microstructures were compared between the two models. Similar results from both models indicate that calibrated, lower fidelity numerical models may be used in place of higher fidelity models to generate melt pool solidification data. Finally, these calibrated models therefore enable lower computational cost melt pool simulations without a noticeable decrease in simulation accuracy for grain-scale microstructure simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Benchmarking collective effects of electron interactions in a wiggler with OPAL-FEL

OPAL-FEL is a recently developed tool for the modeling of particle accelerators containing wigglers or undulators. It extends the well established 3D electrostatic particle-tracking code OPAL, by merging it with the finite-difference time-domain electromagnetic solver MITHRA. We present results of two benchmark cases where OPAL-FEL simulations are compared to experimental results. Both experiments concern electron beamlines where the longitudinal phase space is modulated with a short magnetic wiggler. Good agreement was found in both the space charge and radiation dominated regimes.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Structural and Dynamic Heterogeneity of Deep Eutectic Solvents Composed of Choline Chloride and Ortho-Phenol Derivatives

Structural, thermal, and dynamic properties of four deep eutectic solvents comprising choline chloride paired with ortho-phenolic derivative hydrogen-bond donors were probed using experiments and molecular simulations. The hydrogen-bond donors include phenol, catechol, o-chlorophenol, and o-cresol, in a 3:1 mixture with the hydrogen-bond acceptor choline chloride. Density, viscosity, and pulsed-field gradient NMR diffusivity measurements were conducted over a range of temperatures. Classical and ab initio molecular dynamics simulation results match experimental data reasonably well. Furthermore, the simulation results were then used to perform a more detailed analysis of the local structure and dynamics of these systems.

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

Nuclear–Electronic Orbital Quantum Mechanical/Molecular Mechanical Real-Time Dynamics

Simulating the nuclear–electronic quantum dynamics of large-scale molecular systems in the condensed phase is key for studying biologically and chemically important processes such as proton transfer and proton-coupled electron transfer reactions. Herein, the real-time nuclear–electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) approach is combined with a hybrid quantum mechanical/molecular mechanical (QM/MM) strategy to enable the accurate description of coupled nuclear–electronic quantum dynamics in the presence of heterogeneous environments such as solvent or proteins. The densities of the electrons and quantum protons are propagated in real time, while the other nuclei are propagated classically on the instantaneous electron–proton vibronic surface. This approach is applied to phenol bound to lysozyme, intramolecular proton transfer in malonaldehyde, and nonequilibrium excited-state intramolecular proton transfer in o-hydroxybenzaldehyde. Furthermore, these examples illustrate that the RT-NEO-TDDFT framework, coupled with an atomistic representation of the environment, allows the simulation of condensed-phase systems that exhibit significant nuclear quantum effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Direct Ab Initio Simulation of the Synthesis of BaZrO 3 and the Microstructure Impacts on Proton Transport

Controlling and predicting the processing-structure-performance relationship in functional materials is a grand challenge in materials science, with important implications for a wide range of emerging applications; a high fidelity understanding of the performance impact of microstructures formed under synthesis conditions is required to develop advanced materials, such as solid-state fuel cells and electrolyzers. Using the ceramic BaZrO 3 as a case study, we directly simulate the synthesis and investigate how proton transport is dictated by microstructures. We develop a framework that couples density functional theory (DFT), machine-learning interatomic potential (MLIP) driven molecular dynamics, and grand canonical Monte Carlo to perform large-scale, microstructure-resolved, atomistic simulations of proton transport in experimentally representative polycrystalline structures. Our fully ab initio approach, using a MLIP as a proxy for DFT, allows us to quantify the competition between two distinct diffusion mechanisms: one associated with grain-boundary regions and another within grains. When the impacts of grain boundaries are taken into account, proton transport exhibits substantial deviation from the bulk oxide limit. This addresses long-standing discrepancies between theory and experiments. Our integrated approach provides atomistic insight into microstructure-dependent proton pathways in BaZrO 3 and establishes a general protocol for predicting processing-structure-performance relationships.

organic↗