Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Singularity”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Boundary-consistent B-spline filtering schemes and application to high-fidelity simulations of turbulence

A filtering operation, based on B-spline discretizations, is introduced to target weakly growing mesh-scale oscillations that can arise in high-fidelity turbulence simulations. This is a spectral regularization that can be described using the singular values of a banded matrix operator, with the filtering strength set by a scalar- or vector-valued penalty parameter. The penalty parameter can be specified though it can also be advantageously selected to minimize the generalized cross validation (GCV) measure of distance between the pre- and post-filtered solutions. Efficient algorithms are developed to compute both the scalar and vector penalty parameters. The B-spline filter has a sharper localization to high-wavenumber than compact or explicit filters of the same stencil width and is demonstrated for solutions of the Burgers' equation, decaying Burgers' turbulence, and compressible Navier–Stokes turbulent channel flow. Furthermore, these simulations confirm the scheme's numerical stability and ability to narrowly target the high wavenumber components of numerical solutions. An advantage over finite-difference filters is that these B-spline filters are stable on bounded domains and even preserve formal order of accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Accurate numerical, integral methods for computing drift-kinetic Trubnikov-Rosenbluth potentials

A novel numerical method is employed to compute the integral form of the axi-symmetric Trubnikov-Rosenbluth potentials. Two methods for quadrature in pitch-angle are described and their convergence properties are studied. Careful attention is given to quadrature over a singular Green's function. Here it is shown that an infinite series representation of the Green's function can be used more efficiently than its closed form involving complete elliptic integrals. Then a collocation method in speed, with its associated quadrature scheme, is laid out and its convergence properties are studied. Using the proposed scheme, accurate low-order moments of the field collision operator are obtained using relatively few velocity space degrees of freedom. The scheme is showcased by solving for the equilibrium, axi-symmetric bootstrap current in tokamaks. A C 0 Gauss-Lobatto-Legendre finite element pitch-angle basis with vertex nodes at the trapped/passing boundary is shown, in the context of the integral methods used, to be much more efficient than the more common Legendre polynomial expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Parametric dynamic mode decomposition for reduced order modeling

Dynamic Mode Decomposition (DMD) is a model-order reduction approach, whereby spatial modes of fixed temporal frequencies are extracted from numerical or experimental data sets. The DMD low-rank or reduced operator is typically obtained by singular value decomposition of the temporal data sets. For parameter-dependent models, as found in many multi-query applications such as uncertainty quantification or design optimization, the only parametric DMD technique developed was a stacked approach, with data sets at multiple parameter values were aggregated together, increasing the computational work needed to devise low-rank dynamical reduced-order models. Here in this paper, we present two novel approach to carry out parametric DMD: one based on the interpolation of the reduced-order DMD eigen-pair and the other based on the interpolation of the reduced DMD (Koopman) operator. Numerical results are presented for diffusion-dominated nonlinear dynamical problems, including a multiphysics radiative transfer example. All three parametric DMD approaches are compared.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Compression of tokamak boundary plasma simulation data using a maximum volume algorithm for matrix skeleton decomposition

This report demonstrates satisfactory data compression of SOLPS-ITER simulation output ranging from 2D fields, 1D profiles, and 0D scalar variables with a novel matrix decomposition approach. The singular value decomposition (SVD) scales poorly for large matrix sizes and is unsuited to the application on high dimensional data common to fusion plasma physics simulation. In this work, we employ the columns-submatrix-rows (CUR) matrix factorization technique in order to compute a low-rank approximation up to two orders of magnitude faster than the SVD, but within a nominal L2-norm relative error of ε = 10 –2 . In addition, the CUR approach maintains the original format of the data, in its extracted columns and rows, allowing for interpretable data storage at the original resolution of the simulation. We utilize an iterative algorithm to compute the CUR decomposition of simulation output by maximizing the volume, or linearly independent information content, of a low-rank submatrix contained within the data. Experiments over $\textit{n} × \textit{n}$ randomized test matrices with embedded rank-deficient features show that this maximum volume implementation of CUR matrix approximation has reduced asymptotic computational complexity on the order of n compared to the SVD, which scales approximately as $n^3$. These results show that the CUR technique can be used to effectively select time step snapshots (columns) of over 140 SOLPS-ITER output variables and the associated discretized coordinate timeseries (rows) allowing for reconstruction of the complete simulation dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Permutation-adapted complete and independent basis for atomic cluster expansion descriptors

Atomic cluster expansion (ACE) methods provide a systematic way to describe particle local environments of arbitrary body order. For practical applications it is often required that the basis of cluster functions be symmetrized with respect to rotations and permutations. Existing methodologies yield sets of symmetrized functions that are over-complete. These methodologies thus require an additional numerical procedure, such as singular value decomposition (SVD), to eliminate redundant functions. In this work, it is shown that analytical linear relationships for subsets of cluster functions may be derived using recursion and permutation properties of generalized Wigner symbols. From these relationships, subsets (blocks) of cluster functions can be selected such that, within each block, functions are guaranteed to be linearly independent. It is conjectured that this block-wise independent set of permutation-adapted rotation and permutation invariant (PA-RPI) functions forms a complete, independent basis for ACE. Along with the first analytical proofs of block-wise linear dependence of ACE cluster functions and other theoretical arguments, numerical results are offered to demonstrate this. The utility of the method is demonstrated in the development of an ACE interatomic potential for tantalum. Using the new basis functions in combination with Bayesian compressive sensing sparse regression, some high degree descriptors are observed to persist and help achieve high-accuracy models.

Angular momentum↗
DOE OSTI2025

Estimating QSVT angles for matrix inversion with large condition numbers

Quantum Singular Value Transformation (QSVT) is a state-of-the-art, near-optimal quantum algorithm that can be used for matrix inversion. The QSVT circuit is parameterized by a sequence of angles that must be pre-calculated classically, with the number of angles increasing as the matrix condition number grows. Computing QSVT angles for ill-conditioned problems is a numerically challenging task. Here, we propose a numerical technique for estimating QSVT angles for large condition numbers. This technique allows one to avoid expensive numerical computations of QSVT angles and to emulate QSVT circuits for solving ill-conditioned problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Goal-oriented real-time Bayesian inference for linear autonomous dynamical systems with application to digital twins for tsunami early warning

We present a goal-oriented framework for constructing digital twins with the following properties: (1) they employ discretizations of high-fidelity partial differential equation (PDE) models governed by autonomous dynamical systems, leading to large-scale forward problems; (2) they solve a linear inverse problem to assimilate observational data to infer uncertain model components followed by a forward prediction of the evolving dynamics; and (3) the entire end-to-end, data-to-inference-to-prediction computation is carried out without approximation and in real time through a Bayesian framework that rigorously accounts for uncertainties. Several challenges must be overcome to realize this framework, including the large scale of the forward problem, the high dimensionality of the parameter space, and for a class of problems including those we target, the slow decay of the singular values of the parameter-to-observable map. Here we introduce a methodology to overcome these challenges by exploiting the autonomous structure of the forward model to decompose the solution of the inverse problem into a one-time-only offline phase in which the PDE model is solved a limited number of times (equal to the number of sensors), and an online phase that maps well onto GPUs and computes the parameter inference and prediction of quantities of interest in real time, given observational data. Our ultimate goal is to apply this framework to construct digital twins for subduction zones, including Cascadia, to provide early warning for tsunamis generated by megathrust earthquakes. To this end, we demonstrate how our methodology can be used to employ seafloor pressure observations, along with the coupled acoustic–gravity wave equations, to infer the earthquake-induced spatiotemporal seafloor motion (discretized with $\mathscr{O}$ (10 9 ) parameters) and forward predict the tsunami propagation. We present results of an end-to-end inference, prediction, and uncertainty quantification for a representative test problem with $\mathscr{O}$ (10 8 ) inversion parameters for which goal-oriented Bayesian inference is accomplished exactly and in real time, that is, in a matter of seconds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

What do physics-informed DeepONets learn? Understanding and improving training for scientific computing applications

Physics-informed deep operator networks (DeepONets) have emerged as a promising approach toward numerically approximating the solution of partial differential equations (PDEs). In this work, we aim to develop further understanding of what is being learned by physics-informed DeepONets by assessing the universality of the extracted basis functions and demonstrating their potential toward model reduction with spectral methods. Results provide clarity about measuring the performance of a physics-informed DeepONet through the decays of singular values and expansion coefficients. In addition, we propose a transfer learning approach for improving training for physics-informed DeepONets between parameters of the same PDE as well as across different, but related, PDEs where these models struggle to train well. This approach results in significant error reduction and learned basis functions that are more effective in representing the solution of a PDE.

Deep operator networks↗
DOE OSTI2025

An innovative radial gradient material design using hot isostatic pressing for applications in extreme environments

Functionally graded materials (FGMs) are highly advanced continuous or discontinuous structures whose structural and material properties vary along a singular geometric dimension either in the axial or radial direction. Here, the radial gradient FGM design makes for an optimal structural design to incorporate a bi-metallic structure with a copper-based high entropy alloy (Cu-HEA) with good mechanical properties and high irradiation resistance, and Chromium (Cr) with great corrosion resistance. This study focuses on the experimental design of a metal powder loading mechanism to fabricate a bi-metallic radial gradient structure using Cu-HEA and 99.9 % pure Cr metal powders. The powder loading strategy uses custom-designed concentric cylindrical dividers to separate the individual compositions. Two benchtop trial runs were performed for design optimization. The optimized design was then implemented to eventually load the HEA and Cr powders for consolidation via powder metallurgy hot isostatic pressing (PM-HIP). The electron microscopy analysis reveals the successful fabrication of the radial gradient structure with the chemical mapping analysis, demonstrating the gradual composition shift from the HEA at the center to the pure-Cr at the periphery via a three-step gradient.

High Entropy Alloys (HEAs)↗
DOE OSTI2020

A predictive strain-gradient model with no undetermined constants or length scales

A general meso–scale (GM) crystal plasticity (CP) model was developed that accounts for lower-order (strain hardening) and higher-order (internal stress) effects of geometrically necessary dislocations (GNDs). It is predictive: no arbitrary parameters or length scales were invoked and no ad hoc numerical techniques were employed. It uses general stress field equations for GND content and a novel harmonization technique to enforce consistency of elastic long-range singular defect fields with applied elastic-plastic fields. The model facilitates implementation in commercial finite element programs without requiring special elements, special boundary conditions, or access to element shape functions. GM simulations confirmed, with improved accuracy, previously published predictions of the Hall-Petch effect, Bauschinger effect, and anelasticity. Previously unpredicted phenomena were also predicted: anelasticity and hysteresis for single Ta crystals and strain-hardening stagnation. The internal stresses (higher-order effect) dominate at large length scales, while at small length scales, the GND density hardening (lower-order effect) dominates. As a result, GM predicts that strain heterogeneity and consequent GND internal stresses are important factors in anelasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A general solution for accelerating screw dislocations in arbitrary slip systems with reflection symmetry

Solutions to the differential equations of linear elasticity in the continuum limit in arbitrary crystal symmetry are known only for steady-state dislocations of arbitrary character, i.e. line defects moving at constant velocity. Troubled by singularities at certain ‘critical’ velocities (typically close to certain sound speeds), these dislocation fields are thought to be too idealized, and divergences are usually attributed to neglecting the finite size of the core and to the restriction to constant velocity. In the isotropic limit, accelerating pure screw and edge dislocations were studied some time ago. A generalization to anisotropic crystals has been attempted for pure screw and edge dislocations only for some special cases. This work aims to fill the gap of deriving a general anisotropic solution for pure screw dislocations applicable to slip systems featuring a reflection symmetry, a prerequisite to studying pure screw dislocations without mixing with edge dislocations. Finally, further generalizations to arbitrary mixed dislocations as well as regularizations of the dislocation core are beyond the scope of this paper and are left for future work.

42 ENGINEERING↗
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

DESPERATE: A Python Library for Processing and Denoising NMR Spectra

NMR spectroscopy is an inherently insensitive technique with respect to the amount of observable signal. A common element in all NMR spectra is random thermal noise that is often characterized by a signal-to-noise ratio (SNR). SNR can be generically improved experimentally with repetitive signal averaging or during post-processing with apodization; the former of which often results in long experimental times and the latter results in the loss of spectral resolution. Denoising techniques can instead be used during post-processing to enhance SNR without compromising resolution. The most common approach relies on the singular-value decomposition (SVD) to discard noisy components of NMR data. SVD-based approaches work well, such as Cadzow and PCA, but are computationally expensive when used for large datasets that are often encountered in NMR (e.g., Carr-Purcell/Meiboom-Gill and nD datasets). Herein, we describe the implementation of a new wavelet transform (WT) routine for the fast and robust denoising of 1D and 2D NMR spectra. Several simulated and experimental datasets are denoised with both SVD-based Cadzow or PCA and WT’s, and the resulting SNR enhancements and spectral uniformity are compared. WT denoising offers similar and improved denoising compared with SVD and operates faster by several orders-of-magnitude in some cases. Further, all denoising and processing routines used in this work are included in a free and open-source Python library called DESPERATE.

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

Thermal performance and energy consumption validation of an occupied local government office building outfitted with ceiling tile phase change materials

Buildings present an opportunity for energy conservation and the modulation of peak energy demand through controlled Heating, Ventilation, and Air Conditioning (HVAC) energy use. The administrative and office building stock in the United States holds potential to achieve energy and demand savings through retrofits such as insulation, weatherization, and thermal energy storage. Specifically, there is a need to validate passive phase change material (PCM) applications in full scale in aging administrative buildings in the US to evaluate the energy benefits. Aim of this study was to conduct a whole building level thermal and energy validation of an operational building and explore an alternative method for evaluating energy efficiency. To accomplish this, the study employed PCMs in the drop ceiling and carry out an energy audit and on-site measurement of HVAC systems' energy demand and consumption. A full-scale EnergyPlus energy model, modeled by the authors, served as a baseline for evaluation. The results show that calibrated model's envelope temperature measures fall within the accepted errors. HVAC energy simulation results also fall within the accepted errors for monthly and hourly pre- and post- PCM retrofit electricity and natural gas data. The novelty of this study is that it employees energy scales per Heating Degree Hour and Cooling Degree Hour, in contrast to the commonly used Heating Degree Days and Cooling Degree Days as reported in the literature to analyze energy savings. These findings underscore the pivotal role of a calibrated model in assessing the efficacy of a singular energy measure, like a PCM-retrofitted ceiling, in an occupied office building.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Loss of biological control of enamel mineralization in amelogenin-phosphorylation-deficient mice

Amelogenin, the most abundant enamel matrix protein, plays several critical roles in enamel formation. Importantly, we previously found that the singular phosphorylation site at Ser16 in amelogenin plays an essential role in amelogenesis. Studies of genetically knock-in (KI) modified mice in which Ser16 in amelogenin is substituted with Ala that prevents amelogenin phosphorylation, and in vitro mineralization experiments, have shown that phosphorylated amelogenin transiently stabilizes amorphous calcium phosphate (ACP), the initial mineral phase in forming enamel. Furthermore, KI mice exhibit dramatic differences in the enamel structure compared with wild type (WT) mice, including thinner enamel lacking enamel rods and ectopic surface calcifications. Here, we now demonstrate that amelogenin phosphorylation also affects the organization and composition of mature enamel mineral. In this work, we compared WT, KI, and heterozygous (HET) enamel and found that in the WT elongated crystals are co-oriented within each rod, however, their c-axes are not aligned with the rods' axes. In contrast, in rod-less KI enamel, crystalline c-axes are less co-oriented, with misorientation progressively increasing toward the enamel surface, which contains spherulites, with a morphology consistent with abiotic formation. Furthermore, we found significant differences in enamel hardness and carbonate content between the genotypes. ACP was also observed in the interrod of WT and HET enamel, and throughout aprismatic KI enamel. In conclusion, amelogenin phosphorylation plays crucial roles in controlling structural, crystallographic, mechanical, and compositional characteristics of dental enamel. Thus, loss of amelogenin phosphorylation leads to a reduction in the biological control over the enamel mineralization process.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A numerical and experimental approach to compare the effect of sample thickness in small in-situ SEM and large ex-situ tensile testing in Alloy 709

ASTM standards for tensile tests define specific sample size requirements regardless of grain size of the material. However, sample size requirements for testing should be considered in conjunction with the number of grains within its reduced cross-sectional area. This is particularly important for in-situ Scanning Electron Microscope (SEM) tensile tests, as they have to be conducted on smaller samples. In this study, a comprehensive experimental and numerical evaluation of the effect of specimen thickness (and the resulting number of grains within the cross-section) were conducted using in-situ SEM tensile test (on submillimeter thick samples) and ex-situ tensile tests (on samples of 0.68–5.9 mm thickness) and the results are compared with FEM simulations outcome. All tests were conducted at room temperature and the results are correlated to the number of grains within the thickness of specimens. The tensile test results indicated that even though the 0.2% proof stress and the tensile strength do not vary for the current range of thicknesses (with number of grains in the cross-section varying from 13 to 118), a difference in necking mechanism exists. Post tensile strength, thinner samples undergo more shear failure and diagonal localized necking whereas thicker samples experience more diffuse necking indicated by decreased area of shear failure at the edges and increased area of dimpled tensile failure at the center of the sample. Further, FEM results complement the experimental findings by showing the formation of conjugated localized shear bands on the upper and lower surface in the 0.68 mm thick sample and shear bands combining to form singular shear bands in thicker samples. These results also confirm the validity of the in-situ SEM tensile tests conducted on thinner samples as long as the required minimum number of grains exist within the cross-section.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Tunable Second Harmonic Generation in Twisted Bilayer Graphene

Twisted stacking of van der Waals (vdW) materials introduces a new way in bandstructure engineering and has given rise to numerous extraordinary physical phenomena. Despite the absence of second harmonic generation (SHG) in non-gated monolayer graphene, artificially twisted bilayer graphene (tBLG) possesses more possible point-group symmetries including those with broken inversion-symmetry. Here, we report twist-angle-dependent SHG from tBLG, which is the first demonstration of an elemental material with intrinsically tunable nonlinearity. We show that depending on the twist angle, the susceptibility of the dominant chiral tensor component of tBLGs can vary from 0 to 28×10 4 pm 2 /V, which is at the same order of magnitude as on-resonance susceptibility of monolayer MoS 2 . These results shed light on the underlying symmetry of tBLG systems and electronic band-structure near van Hove singularities (vHs). More importantly, they introduce a new degree of freedom, the twisting, in creating efficient second-order nonlinear material form centrosymmetric constituents

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of Fermi surface topology change on the Kagome superconductor CeRu 2 under pressure

The cubic Laves phase compound CeRu 2 with a Kagome substructure of Ru has been investigated to understand myriad fascinating phenomena resulting from competition among its various physical and geometric features. Such phenomena include flat bands, van Hove singularities, Dirac cones, reentrant superconductivity, magnetism, the Fulde–Ferrell–Larkin–Ovchinnikov state, valence fluctuations, time-irreversible anisotropic s-state superconductivity, etc. Extensive studies have thus been carried out since 1958 when the highly unusual coexistence of superconductivity and ferromagnetism was proposed for the mixed compounds (Ce,Gd)Ru 2 . Activity has accelerated in recent years due to increasing interest in topological states in superconductors. However, there has been little investigation of the mutual influence of these fascinating states. Therefore, we systematically investigated the superconductivity and possible Fermi surface topological change in CeRu 2 via magnetic, resistivity, and structural measurements under pressure up to ~168 GPa. An unusual phase diagram that suggests an intriguing interplay between the compound’s superconducting order and Fermi surface topological order has been constructed. A resurgence in its superconducting transition temperature was observed above 28 GPa. Our experiments have identified a structural transition above 76 GPa and a few tantalizing phase transitions driven by high pressure. Our high-pressure results further suggest that superconductivity and Fermi surface topology in CeRu 2 are strongly intertwined.

36 MATERIALS SCIENCE↗