SEARCH · Engineering Papers

Results for “Numerical techniques”

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 811 records · Page 45

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Development of a molecularly informed biogeochemical framework for reactive transport modeling of subsurface carbon inventories, transformations and fluxes (Final Report)

The overall objective of the project was to combine new molecular-level characterization strategies with soil carbon flux measurements to develop and evaluate model representations of subsurface carbon cycling. We expanded on existing studies in the East River watershed, Colorado, in collaboration with the Berkeley Lab Watershed Function Scientific Focus Area (SFA), SLAC Groundwater Quality SFA, and Rocky Mountain Biological Laboratory (RMBL) to develop an elevation and vegetation gradient that is now the subject of long-term monitoring by the USGS. To achieve the overall objective, we combined field studies of soil respiration with laboratory analyses, ranging from spectroscopy to incubation studies. This combination of techniques enabled us to develop a new understanding of the drivers of high-elevation soil respiration. In the process, we developed three new modeling approaches that improve our ability to conceptualize and ultimately to represent soil respiration in numerical models. The first approach is a plot-scale transient inverse model that can be used to determine in situ CO 2 production rates from measured concentration profiles and surface fluxes. Application of this method revealed the importance of plant phenology and deep CO 2 production in moderating CO 2 fluxes to the atmosphere. The second modeling approach is a molecular-scale tool that enables spectroscopic and elemental data for carbon speciation to be transformed into functional group abundances, or the ‘SOC-fga model’. This method uniquely enables carbon speciation to be tracked within a reactive transport framework to partition carbon among different pathways and storage zones within the soil. The third modeling tool builds strongly on the previous approaches and captures the microbial processes driving heterotrophic respiration. This ‘dormancy model’ allows the native soil microbial population to respond transiently to the presence or absence of water in order to catalyze carbon respiration. This approach was also compared to the simpler and more widely used first-order model using two experimental datasets with different temporal and spatial resolutions. Our results illustrate that the simpler first order model provides a robust and efficient representation of deep (>1 meter) soil respiration, but that shallow soils, where most respiration occurs, require explicit representation of moisture-dependent activation and dormancy rates. By assessing soil organic carbon turnover at multiple scales, we see a complex array of controls emerge. Atte scale of a hillslope, spatial heterogeneity in soil respiration rates dominates and is uncorrelated with instantaneous soil moisture and plant community. At the profile scale, the balance between plant inputs and water availability is the dominant control. At the microbial to molecular scale, physiological processes associated with carbon use and carbon speciation are important controls. These scale-dependent controls emphasize the need for new modeling approaches that examine their interactions and hierarchies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

A Comprehensive Overview of Postbiotics with a Special Focus on Discovery Techniques and Clinical Applications

The increasing interest in postbiotics, a term gaining recognition alongside probiotics and prebiotics, aligns with a growing number of clinical trials demonstrating positive outcomes for specific conditions. Postbiotics present several advantages, including safety, extended shelf life, ease of administration, absence of risk, and patentability, making them more appealing than probiotics alone. This review covers various aspects, starting with an introduction, terminology, classification of postbiotics, and brief mechanisms of action. It emphasizes microbial metabolomics as the initial step in discovering novel postbiotics. Commonly employed techniques such as NMR, GC-MS, and LC-MS are briefly outlined, along with their application principles and limitations in microbial metabolomics. The review also examines existing research where these techniques were used to identify, isolate, and characterize postbiotics derived from different microbial sources. The discovery section concludes by highlighting challenges and future directions to enhance postbiotic discovery. In the second half of the review, we delve deeper into numerous published postbiotic clinical trials to date. We provide brief overviews of system-specific trial applications, their objectives, the postbiotics tested, and their outcomes. The review concludes by highlighting ongoing applications of postbiotics in extended clinical trials, offering a comprehensive overview of the current landscape in this evolving field.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2026

Effect of macroscopic surface defects on dynamic damage: An experimental and numerical study

This study examines the impact of macroscopic surface defects on the dynamic ductile damage behavior of polycrystalline metals using plate-impact experiments. Defects of various shapes (flat, round, and point) were manufactured on the free or impact surfaces of annealed copper specimens. The experiments were diagnosed with photon Doppler velocimetry measurements and soft recovery techniques. The experimental results revealed that defect shape and location significantly affect velocity–time profiles and void distribution. In order to understand the dynamics of shock propagation and corresponding ductile damage evolution within the target specimens, the impact experiments were modeled with a calibrated Tepla model, a dynamic ductile damage model for polycrystalline metals [Nguyen et al., Int. J. Solids Struct. 329, 113833 (2026)]. Overall, our resulting simulated velocity showed good agreement with measured velocity, and our simulated porosity distributions qualitatively matched experimental data. Based on our simulation results, defects on the free surface were found to distort rarefaction waves and therefore the corresponding spall planes. On the other hand, defects on the impact surface generated a delayed shock when the flyer plate and the defected target area were in contact, leading to distortion of the spall plane. The distortion of the spall plane resulted in a non-uniform distribution of voids within the defected specimens. Lastly, we examined the role of local work hardening due to the defect manufacturing process on the velocity and porosity distribution, using Tepla simulations with a simplified representation of local hardening. This investigation highlights the importance of defect geometry, location, and local hardening associated with defect generation in dynamic ductile damage processes.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Multimodal Defect Imaging of Pure Tungsten Components Fabricated via Electron Beam Powder Bed Fusion

The utilization of additive manufacturing (AM) techniques for refractory materials in high-temperature environments has significantly expanded because of the ability to fabricate geometrically complex components. Electron beam powder bed fusion (EB-PBF), which provides lower residual stress, a cleaner vacuum environment, and better efficiency for high melting point, is one of the best-suited AM methods to produce advanced refractory components. However, the property variation attributed to the heterogeneous microstructure and process-induced defects has hindered the widespread adoption of EB-PBF-produced material like tungsten. While numerous in-situ monitoring and defect detection methods have been demonstrated for EB-PBF, a workflow that compares and evaluates process-induced abnormalities from different imaging perspectives is still limited. This study examines a feature-embedded tungsten component manufactured via the EB-PBF process to demonstrate the defect detection capabilities of a multimodal defect imaging workflow. The predefined and process-induced defects are evaluated by harnessing various imaging techniques, including in-situ electron imaging, layerwise near-infrared (NIR) imaging, post-build high-energy x-ray computed tomography (CT), and conventional destructive metallography. The results highlight the strengths and limitations of distinctive defect imaging techniques concerning specific defect types, sizes, and conditions. It was found that electron imaging can provide more abnormal detection capabilities while maintaining a higher measuring accuracy, against the conventional metallography in this case study, compared with NIR and CT imaging techniques.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Solid‐state nuclear magnetic resonance of spin‐9/2 nuclei 115 In and 209 Bi in functional inorganic complex oxides

Abstract Indium and bismuth are technologically important elements, in particular as oxides for optoelectronic applications. 115 In and 209 Bi are both I = 9/2 nuclei with high natural abundances and moderately high frequencies but large nuclear electric quadrupole moments. Leveraging the quadrupolar interaction as a measure of local symmetry and polyhedral distortions for these nuclei could provide powerful insights on a range of applied materials. However, the absence of reported nuclear magnetic resonance (NMR) parameters on these nuclei, particularly in oxides, hinders their use by the broader materials community. In this contribution, solid‐state 115 In and 209 Bi NMR of three recently discovered quaternary bismuth or indium oxides are reported, supported by density functional theory calculations, numerical simulations, diffraction and additional multinuclear ( 27 Al, 69,71 Ga, and 121 Sb) solid‐state NMR measurements. The compounds LiIn 2 SbO 6 , BiAlTeO 6 , and BiGaTeO 6 are measured without special equipment at 9.4 T, demonstrating that wideline techniques such as the QCPMG pulse sequence and frequency‐stepped acquisition can enable straightforward extraction of quadrupolar tensor information in I = 9/2 115 In and 209 Bi even in sites with large quadrupolar coupling constants. Relationships are described between the NMR observables and local site symmetry. These are amongst the first reports of the NMR parameters of 115 In, 121 Sb, and 209 Bi in oxides.

Griffith, Kent J.↗

DOE OSTI2021

Statistical learning for predicting density–matrix-based electron dynamics

In this study, we consider the problem of learning density-dependent molecular Hamiltonian matrices from time series of electron density matrices, all in the context of Hartree–Fock theory. Prior work developed a solution to this problem for small molecular systems with density and Hamiltonian matrices of size at most 6 × 6. Here, using a battery of techniques, we scale prior methods to larger molecular systems with, for example, 29 × 29 matrices. This includes systems that either have more electrons or are expressed in large basis sets such as 6-311++G**. Scaling the method to larger systems enhances its relevance for realistic applications in chemistry and physics. To achieve this scaling, we apply dimensionality reduction, ridge regression and analytic computation of Hessians. Through the combination of these techniques, we are able to learn Hamiltonians by minimizing an objective function that encodes local propagation error. Importantly, these learned Hamiltonians can then be used to predict electron dynamics for thousands of steps: When we use our learned Hamiltonians to numerically solve the time-dependent Hartree–Fock equation, we obtain predicted dynamics that are in close quantitative agreement with ground truth dynamics. This includes field-off trajectories similar to the training data and field-on trajectories outside of the training data.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Robust second-order approximation of the compressible Euler equations with an arbitrary equation of state

Here, this paper is concerned with the approximation of the compressible Euler equations supplemented with an arbitrary or tabulated equation of state. The proposed approximation technique is robust, formally second-order accurate in space, invariant-domain preserving, and works for every equation of state, tabulated or analytic, provided the pressure is nonnegative. An entropy surrogate functional that grows across shocks is proposed. The numerical method is verified with novel analytical solutions and then validated with several computational benchmarks seen in the literature including problems with composite waves.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Characterization of nonlinear ultrasonic waves behavior while interacting with poor interlayer bonds in large-scale additive manufactured materials

Over the past decades, researchers have developed several nonlinear ultrasonic techniques for quality control of materials commonly used in different applications. Owing to the superior sensitivity of nonlinear ultrasound waves to small defects such as micro-cracks, their applicability in different nondestructive testing (NDT) problems has been investigated in numerous studies. These studies utilize frequency domain analysis to detect the generation of higher harmonics because of the formation of defects in the inspected medium. Frequency domain analysis based on the Fourier transform is a significant approach used in linear systems; however, it may not perform adequately on nonlinear systems. Hence, studies on nonlinear dynamics and physics consider analyzing systems' behavior in the phase-space domain. In contrast to the frequency domain analysis, which can result in information loss, analysis in the phase-space domain retains all the information regarding a system's states. Here, we investigate the nonlinearities induced by poor interlayer bonding in polymer-based additive manufactured parts in the phase-space domain. It is convenient to characterize the nonlinearity in the phase-space domain because it provides a geometrical representation of a system's states. Two types of low quality interlayer bond are considered. The first type is simulated artificially while the second type is manufactured by reducing the bond quality during the printing process. The analysis verified that the received ultrasonic signals exhibit classical nonlinear behavior in the phase-space domain while interacting with simulated poor interlayer bonds. In addition, the results showed that the behavior of ultrasonic waves is amplitude-dependent and evolves into models that have not been previously reported. Furthermore, Largest Lyapunov Exponent (LLE) is used to quantify the behavior of nonlinear ultrasonic waves while interacting with poor interlayer bonds. Using LLE, it was observed that the divergence rate of the phase-space trajectories depends on the amplitude of the excitation. This observation quantitatively proves that nonlinear behavior of ultrasound while interacting with poor interlayer bonds can be amplitude-dependent. The results of both simulated and inherent poor interlayer bond cases showed that LLE can be used as a reliable quantitative damage-sensitive feature to detect and potentially characterize weak bonds, which are difficult to detect using conventional approaches. Additionally, the reported results in the phase-space domain provide a basis for proposing a new mathematical model for ultrasonic waves interacting with poor interlayer bonds.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Direct ab initio calculation of the 4 He nuclear electric dipole polarizability

The calculation of nuclear electromagnetic sum rules by directly diagonalizing the nuclear Hamiltonian in a large basis is numerically challenging and has not been performed for A>2 nuclei. With the significant progress of high performance computing, we show that calculating sum rules using numerous discretized continuum states obtained by directly diagonalizing the ab initio no-core shell model Hamiltonian is achievable numerically. Specifically, we calculate the 4 He electric dipole (E1) polarizability, that is an inverse energy weighted sum rule, employing the Daejeon16 NN interaction. We demonstrate that the calculations are numerically tractable as the dimension of the basis increases and are convergent. Our results for the 4 He electric dipole polarizability are consistent with the most recent experimental data and are compared with those of other theoretical studies employing different techniques and various interactions.

Astronomy & Astrophysics↗

DOE OSTI2025

Least H 2 norm updating of quadratic interpolation models for derivative-free trust-region algorithms

One particular class of derivative-free optimization algorithms is trust-region algorithms based on quadratic models given by the under-determined interpolation. Different techniques in updating the quadratic model from iteration to iteration will give different interpolation models. We propose a new way to update the quadratic model by minimizing the $H^{2}$ norm of the difference between neighboring quadratic models. The motivation for applying the $H^{2}$ norm is given. The theoretical properties of our new updating technique are also presented. We propose the projection in the sense of $H^{2}$ norm and the interpolation error analysis of our model function. We obtain the coefficients of the quadratic model function using the Karush–Kuhn–Tucker (KKT) conditions. Numerical results show the advantages of our model on the test set considered, and the derivative-free algorithms based on our least $H^{2}$ norm updating quadratic model functions can solve test problems with fewer function evaluations than the algorithm based on the least Frobenius norm updating model and the other compared methods.

derivative-free optimization↗

DOE OSTI2020

Accurate Consensus-based Distributed Averaging with Variable Time Delay in Support of Distributed Secondary Control Algorithms

We report that distributed secondary control has been widely used in hierarchical control structures, where multiple distributed generators (DGs) need to coordinate to regulate system voltage and frequency. In these systems, consensus algorithms determine the average of a group of dynamic states (e.g. voltages measured by a group of DGs). To be useful, consensus algorithms must be computationally efficient, stable and accurate. In practice, numerous practical implementation challenges significantly affect the consensus equilibrium. In this paper, we quantify the accuracy deviations of the distributed average observer algorithms proposed in the literature to demonstrate the problems with the state-of-the-art distributed averaging techniques. A novel approach is proposed that achieves accurate average tracking in the presence of time-varying communication delays among agents. In our implementation, time synchronization of all distributed controllers is enabled by a novel software platform, called Resilient Information Architecture Platform for the Smart Grid (RIAPS). The proposed distributed average observer is implemented on hardware controllers and its effectiveness is validated in a controller hardware-in-the-loop testbed.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

Periplasmic biomineralization for semi-artificial photosynthesis

Semiconductor-based biointerfaces are typically established either on the surface of the plasma membrane or within the cytoplasm. In Gram-negative bacteria, the periplasmic space, characterized by its confinement and the presence of numerous enzymes and peptidoglycans, offers additional opportunities for biomineralization, allowing for nongenetic modulation interfaces. We demonstrate semiconductor nanocluster precipitation containing single- and multiple-metal elements within the periplasm, as observed through various electron- and x-ray-based imaging techniques. The periplasmic semiconductors are metastable and display defect-dominant fluorescent properties. Unexpectedly, the defect-rich (i.e., the low-grade) semiconductor nanoclusters produced in situ can still increase adenosine triphosphate levels and malate production when coupled with photosensitization. We expand the sustainability levels of the biohybrid system to include reducing heavy metals at the primary level, building living bioreactors at the secondary level, and creating semi-artificial photosynthesis at the tertiary level. The biomineralization-enabled periplasmic biohybrids have the potential to serve as defect-tolerant platforms for diverse sustainable applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Error Estimates for the Optimal Control of a Parabolic Fractional PDE

In this work, we consider the integral definition of the fractional Laplacian and analyze a linear-quadratic optimal control problem for the so-called fractional heat equation; control constraints are also considered. We derive existence and uniqueness results, first order optimality conditions, and regularity estimates for the optimal variables. To discretize the state equation we propose a fully discrete scheme that relies on an implicit finite difference discretization in time combined with a piecewise linear finite element discretization in space. We derive stability results and a novel $L^2(0,T;L^2(\Omega))$ a priori error estimate. On the basis of the aforementioned solution technique, we propose a fully discrete scheme for our optimal control problem that discretizes the control variable with piecewise constant functions, and we derive a priori error estimates for it. We illustrate the theory with one- and two-dimensional numerical experiments.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Simulating many-engine spacecraft: Exceeding 1 quadrillion degrees of freedom via information geometric regularization

We present an optimized implementation of the recently proposed information geometric regularization (IGR) for unprecedented scale simulation of compressible fluid flows applied to multi-engine spacecraft boosters. We improve upon state-of-the-art computational fluid dynamics (CFD) techniques in terms of computational cost, memory footprint, and energy-to-solution metrics. Unified memory on coupled CPU–GPU or APU platforms increases problem size with negligible overhead. Mixed half/single-precision storage and computation are used on well-conditioned numerics. We simulate flow at 200 trillion grid points and 1 quadrillion degrees of freedom, exceeding the current record by a factor of 20. A factor of 4 wall-time speedup is achieved over optimized baselines. Ideal weak scaling is observed on OLCF Frontier, LLNL El Capitan, and CSCS Alps using the full systems. Strong scaling is near ideal at extreme conditions, including 80% efficiency on CSCS Alps with an 8 node baseline and stretching to the full system.

Wilfong, Benjamin [Georgia Institute of Technology↗

DOE OSTI2025

Generative unfolding with distribution mapping

Machine learning enables unbinned, highly-differential cross section measurements. A recent idea uses generative models to morph a starting simulation into the unfolded data. We show how to extend two morphing techniques, Schrödinger Bridges and Direct Diffusion, in order to ensure that the models learn the correct conditional probabilities. This brings distribution mapping (DM) to a similar level of accuracy as the state-of-the-art conditional generative unfolding methods. Numerical results are presented with a standard benchmark dataset of single jet substructure as well as for a new dataset describing a 22-dimensional phase space of Z+2 -jets.

Butter, Anja↗

DOE OSTI2021

Simulation of Low-Rm physics in complex geometries on GPUs with LGR

Efficient modeling of low magnetic Reynolds number (low-Rm) magnetohydrodynamics is often challenging and requires the implementation of innovative techniques to avoid key barriers experienced with prior approaches. We detail a new paradigm for first-principles simulation of the solution to the low-Rm governing equations in complex geometries. As a result of a number of innovative numerical advances, the next-generation GPU (graphics processing unit) accelerated physics code LGR has been successfully applied to the modeling of exploding wire problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Optimal Equilibrium Selection of Price-Maker Agents in Performance-Based Regulation Market

This paper analyzes the oligopolistic equilibrium of multiple price-maker agents in performance-based regulation (PBR) markets. In these markets, there are price-maker agents representing some frequency regulation (FR) providers and a number of independent price-taker FR providers. An equilibrium problem with equilibrium constraints (EPEC) model is employed in this paper to study the equilibria of a PBR market in the presence of price-maker agents and pricetaker FR providers. Due to the incorporation of the FR providers' dynamics, the proposed model is reformulated as a mixed-integer linear programming (MILP) problem over innovative mathematical techniques. An optimal equilibrium point is also selected for the market, where none of the agents is unique deviator and the dynamic performance of power system is improved simultaneously. The effectiveness of proposed optimal equilibrium point is evaluated in the numerical results section by comparing the outputs with conventional optimal dispatches of the FR providers.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2026

Mixed-precision numerics in scientific applications: survey and perspectives

The explosive demand for artificial intelligence (AI) workloads has led to a significant increase in silicon area dedicated to lower-precision computations on recent high-performance computing hardware designs. However, mixed-precision capabilities, which can achieve performance improvements of up to 8x compared to double-precision in extreme compute-intensive workloads, remain largely untapped in most scientific applications. A growing number of efforts have shown that mixed-precision algorithmic innovations can deliver superior performance without sacrificing accuracy. These developments should prompt computational scientists to seriously consider whether their scientific modeling and simulation applications could benefit from the acceleration offered by new hardware and mixed-precision algorithms. In this survey, we (1) review progress across diverse scientific domains—fluid dynamics, weather and climate, quantum chemistry, and computational genomics—that have begun adopting mixed-precision strategies; (2) examine state-of-the-art algorithmic techniques such as iterative refinement, splitting and emulation schemes, and adaptive precision solvers; (3) assess their implications for accuracy, performance, and resource utilization; and (4) survey the emerging software ecosystem that enables mixed-precision methods at scale. We conclude with perspectives and recommendations on cross-cutting opportunities, domain-specific challenges, and the role of co-design between application scientists, numerical analysts, and computer scientists. Collectively, this survey underscores that mixed-precision numerics can reshape computational science by aligning algorithms with the evolving landscape of hardware capabilities.

Graphics processing units↗