Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical mathematics”

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

Inference of phase field fracture models

The phase field approach to modeling fracture uses a diffuse damage field to represent cracks. This representation mollifies singularities that arise in computations with sharp interface models and some of the resultant difficulties in the mathematical and numerical treatment of fracture. Phase field fracture models have proven effective at representing crack propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, lead to an array of possible models. There exists a large body of literature focusing on this class of models and on the impact of model form on the predicted crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here, we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used for the degradation of elastic response due to damage and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). Furthermore, in this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Recent progress in the phase-field dislocation dynamics method

The phase-field dislocation dynamics (PFDD) method, originated in 2002, is a continuum dislocation model that uses order parameters to describe dislocation slips in crystalline materials. In the past two decades, and especially since it was last reviewed in 2016, PFDD was advanced significantly in terms of the mathematical formulation, numerical implementation, and applicability. The main purpose of this short review is to summarize recent progress made to improve the energy functional formulation and numerical techniques of PFDD as well as its recent applications. Additionally, some recommendations for future work to further extend the PFDD method are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stiff neural ordinary differential equations

Neural Ordinary Differential Equations (ODEs) are a promising approach to learn dynamical models from time-series data in science and engineering applications. This work aims at learning neural ODEs for stiff systems, which are usually raised from chemical kinetic modeling in chemical and biological systems. We first show the challenges of learning neural ODEs in the classical stiff ODE systems of Robertson’s problem and propose techniques to mitigate the challenges associated with scale separations in stiff systems. We then present successful demonstrations in stiff systems of Robertson’s problem and an air pollution problem. The demonstrations show that the usage of deep networks with rectified activations, proper scaling of the network outputs as well as loss functions, and stabilized gradient calculations are the key techniques enabling the learning of stiff neural ODEs. The success of learning stiff neural ODEs opens up possibilities of using neural ODEs in applications with widely varying time-scales, such as chemical dynamics in energy conversion, environmental engineering, and life sciences.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Spectral Denoising for Accelerated Analysis of Correlated Ionic Transport

Computation of correlated ionic transport properties from molecular dynamics in the Green-Kubo formalism is expensive, as one cannot rely on the affordable mean square displacement approach. We use spectral decomposition of the short-time ionic displacement covariance to learn a set of diffusion eigenmodes that encode the correlation structure and form a basis for analyzing the ionic trajectories. This allows systematic reduction of the uncertainty and accelerate computations of ionic conductivity in systems with a steady-state correlation structure. We provide mathematical and numerical proofs of the method’s robustness and demonstrate it on realistic electrolyte materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling

In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Incorporating Full Elastodynamic Effects and Dipping Fault Geometries in Community Code Verification Exercises for Simulations of Earthquake Sequences and Aseismic Slip (SEAS)

ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences of earthquakes and aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, and verify and advance the next generation of physics-based earthquake models that reproduce all phases of the seismic cycle. With the goal of advancing SEAS models to robustly incorporate physical and geometrical complexities, here we present code comparison results from two new benchmark problems: BP1-FD considers full elastodynamic effects, and BP3-QD considers dipping fault geometries. Seven and eight modeling groups participated in BP1-FD and BP3-QD, respectively, allowing us to explore these physical ingredients across multiple codes and better understand associated numerical considerations. With new comparison metrics, we find that numerical resolution and computational domain size are critical parameters to obtain matching results. Codes for BP1-FD implement different criteria for switching between quasi-static and dynamic solvers, which require tuning to obtain matching results. In BP3-QD, proper remote boundary conditions consistent with specified rigid body translation are required to obtain matching surface displacements. With these numerical and mathematical issues resolved, we obtain excellent quantitative agreements among codes in earthquake interevent times, event moments, and coseismic slip, with reasonable agreements made in peak slip rates and rupture arrival time. We find that including full inertial effects generates events with larger slip rates and rupture speeds compared to the quasi-dynamic counterpart. For BP3-QD, both dip angle and sense of motion (thrust versus normal faulting) alter ground motion on the hanging and foot walls, and influence event patterns, with some sequences exhibiting similar-size characteristic earthquakes, and others exhibiting different-size events. These findings underscore the importance of considering full elastodynamics and nonvertical dip angles in SEAS models, as both influence short- and long-term earthquake behavior and are relevant to seismic hazard.

Geochemistry & Geophysics↗
DOE OSTI2020

First-principles Modeling and Design of Solid-State Interfaces for the Protection and Use of Lithium Metal Anodes

Li-ion batteries are one of the most advanced energy storage technologies in use today. Li-ion batteries are used in a multitude of applications ranging from consumer electronics, medical devices, sensors and grid storage. However, improving the capacity and energy density delivered by current Li-ion technology requires advanced materials research into novel chemical systems. In this project we have focused particularly in the use of solid-state electrolytes with lithium metal electrodes. Research into all solid-state batteries (ASSB) with Li metal electrodes has significantly expanded in recent years, however most studies reported experimental findings, which left substantial room for theoretical and modeling work as a tool to understand and determine design principles allowing reliable and safe use of ASSBs with Li metal. Among the remaining obstacles preventing reliable use of ASSBs with a Li metal electrode, the stability of the interface between the solid electrolyte and Li metal, and the propagation/dendrite formation of Li metal and resulting mechanical degradation of the electrolyte are key phenomenon that are yet to be fully understood. In the current project we have addressed these two coupled phenomena using first principles calculations and mesoscale continuum modeling. We have obtained chemical and electrochemical stability windows for several solid electrolyte materials. Additionally, from mathematical and numerical modeling of Li protrusion and dendrite initiation during plating and stripping we have determined design criteria in terms of chemical, electrochemical, and mechanical properties and operating conditions for which stable deposition can occur. We also considered the effects of mixed electronic-ionic conduction in solid electrolytes, which has more recently been suggested as another important mechanism involved in ASSB failure. Throughout our work we have successfully addressed important questions necessary for the use of ASSB’s. We have determined guiding principles for materials properties and operating conditions necessary to operate ASSB’s. And have proposed novel solid electrolyte materials with predicted chemical stability an ionic conductivity. Although this represents significant progress in our understanding, open questions remain in order to fully develop reliable and safely operate ASSBs with Li metal. Future work, building on this project will require further experimental, theoretical and simulation efforts to address remaining questions.

25 ENERGY STORAGE↗
DOE OSTI2020

Center for the Geologic Storage of CO 2 (GSCO2) (Final Report)

The Center for Geologic Storage of CO2 (GSCO2) focuses on basic-science approaches to recognized industry and technology challenges for commercial deployment of carbon dioxide (CO2) storage in clastic reservoirs. With this focus, the GSCO2 strives to generate new conceptual, mathematical, and numerical models applicable to geologic storage systems in specific and strategically identified research areas. Over the four-year research period (plus two-year no-cost extension), the GSCO2 developed its focus on one central research question: What is the mechanism of injection-induced microseismicity, and can we control and predict its occurrence?

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Deep learning of free boundary and Stefan problems

Free boundary problems appear naturally in numerous areas of mathematics, science and engineering. These problems present a great computational challenge because they necessitate numerical methods that can yield an accurate approximation of free boundaries and complex dynamic interfaces. In this work, we propose a multi-network model based on physics-informed neural networks to tackle a general class of forward and inverse free boundary problems called Stefan problems. Specifically, we approximate the unknown solution as well as any moving boundaries by two deep neural networks. Besides, we formulate a new type of inverse Stefan problems that aim to reconstruct the solution and free boundaries directly from sparse and noisy measurements. We demonstrate the effectiveness of our approach in a series of benchmarks spanning different types of Stefan problems, and illustrate how the proposed framework can accurately recover solutions of partial differential equations with moving boundaries and dynamic interfaces. All code and data accompanying this manuscript are publicly available at https://github.com/PredictiveIntelligenceLab/DeepStefan.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The DPG Method for the Convection-Reaction Problem Revisited

We study both conforming and non-conforming versions of the practical DPG method for the convection-reaction problem. We determine that the most common approach for DPG stability analysis (construction of a local Fortin operator) is infeasible for the convection-reaction problem. We then develop a line of argument based on the direct construction of a global Fortin operator; we find that employing a polynomial enrichment for the test space does not suffice for this purpose, motivating the introduction of a (two-element) subgrid mesh. The argument combines mathematical analysis with numerical experiments

42 ENGINEERING↗
DOE OSTI2022

The DPG Method for the Convection-Reaction Problem, Revisited

In this work, we study both conforming and non-conforming versions of the practical DPG method for the convection-reaction problem. We determine that the most common approach for DPG stability analysis – construction of a local Fortin operator – is infeasible for the convection-reaction problem. We then develop a line of argument based on a direct proof of discrete stability; we find that employing a polynomial enrichment for the test space does not suffice for this purpose, motivating the introduction of a (two-element) subgrid mesh. The argument combines mathematical analysis with numerical experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Asymptotic preserving methods for fluid electron-fluid models in the large magnetic field limit with mathematically guaranteed properties (Final Report)

The current manuscript is a final report on the activities carried out under the Project LDRD-CIS #226834. In scientific terms, the work reported in this manuscript is a continuation of the efforts started with Project LDRD-express #223796 with final report of activities SAND2021-11481, see [83]. In this section we briefly explain what pre-existing developments motivated the current body of work and provide an overview of the activities developed with the funds provided. The overarching goal of the current project LDRD-CIS #226834 and the previous project LDRD-express #223796 is the development of numerical methods with mathematically guaranteed properties in order to solve the Euler-Maxwell system of plasma physics and generalizations thereof. Even though Project #223796 laid out general foundations of space and time discretization of Euler-Maxwell system, overall, it was focused on the development of numerical schemes for purely electrostatic fluid-plasma models. In particular, the project developed a family of schemes with mathematically guaranteed robustness in order to solve the Euler-Poisson model. This model is an asymptotic limit where only electrostatic response of the plasma is considered. Its primary feature is the presence of a non-local force, the electrostatic force, which introduces effects with infinite speed propagation into the problem. Even though instantaneous propagation of perturbations may be considered nonphysical, there are plenty of physical regimes of technical interest where such an approximation is perfectly valid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Accelerating resonant spectroscopy simulations using multishifted biconjugate gradient

Resonant spectroscopies, which involve intermediate states with finite lifetimes, provide important insights into collective excitations in quantum materials that are otherwise inaccessible. However, theoretical understanding in this area is often limited by the numerical challenges of solving Kramers-Heisenberg-type response functions for large-scale systems. To address this, we introduce a multishifted biconjugate gradient algorithm that exploits the shared structure of Krylov subspaces across spectra with varying incident energies, effectively reducing the computational complexity to that of linear spectroscopies. Both mathematical proofs and numerical benchmarks confirm that this algorithm substantially accelerates spectral simulations, achieving constant complexity independent of the number of incident energies, while ensuring accuracy and stability. This development provides a scalable, versatile framework for simulating advanced spectroscopies in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Diffusioosmotic flow reversals due to ion–ion electrostatic correlations

Existing theories of diffusioosmosis have neglected ion–ion electrostatic correlations, which are important in concentrated electrolytes. Here, we develop a mathematical model to numerically compute the diffusioosmotic mobilities of binary symmetric electrolytes across low to high concentrations in a charged parallel-plate channel. We use the modified Poisson equation to model the ion–ion electrostatic correlations and the Bikerman model to account for the finite size of ions. We report two key findings. First, ion–ion electrostatic correlations can cause a unique reversal in the direction of diffusioosmosis. Such a reversal is not captured by existing theories, occurs at ≈ 0.4 M for a monovalent electrolyte, and at a much lower concentration of ≈ 0.003 M for a divalent electrolyte in a channel with the same surface charge. This highlights that diffusioosmosis of a concentrated electrolyte can be qualitatively different from that of a dilute electrolyte, not just in its magnitude but also its direction. Second, we predict a separate diffusioosmotic flow reversal, which is not due to electrostatic correlations but the competition between the underlying chemiosmosis and electroosmosis. This reversal can be achieved by varying the magnitude of the channel surface charge without changing its sign. However, electrostatic correlations can radically change how this flow reversal depends on the channel surface charge and ion diffusivity between a concentrated and a dilute electrolyte. Furthermore, the mathematical model developed here can be used to design diffusioosmosis of dilute and concentrated electrolytes, which is central to applications such as species mixing and separation, enhanced oil recovery, and reverse electrodialysis.

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

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗
DOE OSTI2025

Quantum Hamiltonian algorithms for maximum independent sets

ABSTRACT We compare two quantum Hamiltonian algorithms that address the maximum independent set problem: one based on the emergent non-Abelian gauge matrix in adiabatic evolution of an energetically isolated manifold of states; the other based on designed application of single-qubit operations. We demonstrate that they are mathematically equivalent in the sense that one is the other’s interaction picture. Despite their mathematical equivalence, our numerical simulations show significant differences between them in performance, which is explained analytically. Intriguingly, this equivalence unveils that the PXP model, recently prominent in quantum dynamics research, can be viewed as quantum diffusion over the median graph of all independent sets governed by the non-Abelian gauge matrix.

Science & Technology - Other Topics↗