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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2026

Cr plasma-material-interaction in PISCES-RF: D thermal release, retention, and erosion

Pure chromium (Cr) targets were exposed to high-flux deuterium (D) plasmas in the Pisces-RF linear plasma device, and measurements of D retention, release, and erosion were subsequently performed. Post-exposure D retention was quantified using temperature-programmed desorption on Cr targets irradiated by 50 eV ions over a broad range of exposure temperatures (423–873 K) and ion fluences (3 × 10 24 – 3 × 10 26 m −2 ). The retained D inventory was observed to decrease rapidly with increased exposure temperature, from approximately ∼7 × 10 20 m −2 at ∼ 420 K, to then saturate near ∼ 10 20 m −2 for exposure temperatures above ∼550 K. Separately, at fixed exposure temperature (∼450 K), D retention was found to have only a weak dependence on increasing ion fluence. Lastly, the erosion of Cr in D plasma was investigated for ion impact energies in the range 40 ≤ E i ≤ 250 eV. Erosion was inferred using optical emission spectroscopy (OES) from the ratio of emission lines (Cr I (425.4 nm)/D I (656.1 nm)) measured close to the target. Conversion of the OES yield data to net erosion yield was made with singular ion energy target mass-loss measurements. These net erosion yield data were then further corrected to obtain gross erosion yield by accounting for a re-deposition factor, computed using a simple model. The gross erosion yield is found to be 2–4 times lower than predicted by SDTrimSP, consistent with that typically observed for light-ion sputtering under high-flux plasma conditions.

Chromium↗
DOE OSTI2023

Principal Component Analysis of azimuthal flow in intermediate-energy heavy-ion reactions

Principal Component Analysis (PCA) via Singular Value Decomposition (SVD) of large datasets is an adaptive exploratory method to uncover natural patterns underlying the data. Several recent applications of the PCA-SVD to event-by-event single-particle azimuthal angle distribution matrices in ultra-relativistic heavy-ion collisions at RHIC-LHC energies indicate that the sine and cosine functions chosen a priori in the traditional Fourier analysis are naturally the most optimal basis for azimuthal flow studies according to the data itself. We perform PCA-SVD analyses of mid-central Au+Au collisions at $E$ $beam$ / $A$ =1.23 GeV simulated using an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model to address the following two questions: (1) if the principal components of the covariance matrix of nucleon azimuthal angle distributions in heavy-ion reactions around 1 GeV/nucleon are naturally sine and/or cosine functions and (2) what if any advantages the PCA-SVD may have over the traditional flow analysis using the Fourier expansion for studying the EOS of dense nuclear matter. In conclusion, we find that (1) in none of our analyses the principal components come out naturally as sine and/or cosine functions, (2) while both the eigenvectors and eigenvalues of the covariance matrix are appreciably EOS dependent, the PCA-SVD has no apparent advantage over the traditional Fourier analysis for studying the EOS of dense nuclear matter using the azimuthal collective flow in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Conundrum for the free energy of a holonomous gluonic plasma at cubic order

We compute the term ~ g 3 in the free energy for a S U ( N ) gauge theory with nonzero holonomy at nonzero temperature. If the holonomy is generated kinematically by the introduction of gauge invariant sources coupled to Polyakov loops, the contribution of charged (off-diagonal) gluons to the free energy at order g 3 , F ( 3 : c . g . ) , is singular: F ( 3 : c . g . ) ≠ 0 without holonomy, but F ( 3 : c . g . ) = 0 when the holonomy is nonzero, even infinitesimally. We show that the absence of the charged gluon contribution is required by gauge invariance alone and is therefore a universal feature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Physical resurgent extrapolation

Expansions of physical functions are controlled by their singularities, which have special structure because they themselves are physical, corresponding to instantons, caustics or saddle configurations. Resurgent asymptotics formalizes this idea mathematically, and leads to significantly more powerful extrapolation methods to extract physical information from a finite number of terms of an expansion, including precise decoding of non-perturbative effects. We quantify the gain of precision for various extrapolation procedures, showing that significant improvements can be achieved using exactly the same input data, and we illustrate the general method with examples from quantum mechanics and quantum field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Uncertainty estimation of bifurcated solutions in the Rayleigh–Bénard problem for advanced nuclear reactors applications

Multiphysics models of nuclear reactors frequently comprise nonlinear systems of equations. The nonlinear nature of these models could lead to solution bifurcations, where a small change in a certain parameter, e.g., the thermophysical properties of the coolant, can lead to a sudden change in the system’s behavior. At the point in parameter space where this happens, called a critical point, the Jacobian matrix of the model’s nonlinear operator becomes singular potentially permitting multiple solutions to coexist. In this paper, we perform uncertainty estimation (UE) in a parameter range that includes bifurcated solutions within the context of Rayleigh–Bénard problem. We perform this analysis assuming uncertain temperature difference, and tilt angle for the iterative solution algorithm with a unit Prandtl number (Pr = 1). Also, we perform this analysis under uncertain thermophysical properties for both FLiBe molten salt and liquid sodium as working fluid. We deploy two approaches to compute statistical moments for the resulting distributions of selected flow-field variables. The first approach is the blind computation of the mean and the standard deviation without any consideration of solution bifurcation, while the second approach utilizes k-means clustering to cluster each branch’s solutions together and compute separate statistical moments for each branch. The statistical distributions are obtained by perturbing the selected parameters about nominal values that correspond to a solution on one of the valid branches, and that solution is used as initial guess for the iterative solution algorithm. We found that perturbation of any parameter when its nominal value is close to its critical point always leads to branch jumping, i.e., the iterations converge to a solution on a branch different from the branch of the initial guess. This produces a statistical ensemble comprised of fundamentally different solutions leading to wrong mean values and uncertainty estimates, whereas clustering provides an efficient way to deal with this type of computation. This work is important for developing Gen IV nuclear systems because many of these systems rely on natural convection for cooling especially in accident conditions.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Harmonizing solar induced fluorescence across spatial scales, instruments, and extraction methods using proximal and airborne remote sensing: A multi-scale study in a soybean field

Solar-induced chlorophyll fluorescence (SIF) has been widely used to track vegetation photosynthesis at different scales ranging from in-situ measurements to satellite products. Airborne platforms sample SIF data at a spatial scale intermediate between in-situ and satellite, matching that of ground measurement (e.g. flux tower footprints and other field sampling), enabling us to explore causes of SIF variation and validate satellite-based SIF products. However, harmonizing SIF across sensors and platforms (correcting for systematic errors to yield a consistent, comparable SIF product) is challenging because SIF can be retrieved in different absorption windows, with different instruments and methods complicating the comparison between different observational levels (i.e., ground, airborne, satellites) and between sites equipped with different instruments with varying optical properties (spectral resolution and sampling intervals, spatial resolution). Additionally, the spatial and temporal variability of atmospheric properties can influence the retrieval of the weak SIF signal. Because of these complications, direct comparisons of airborne and ground SIF across scales are rarely attempted. Here, in this study, we combined airborne SIF data with simultaneous ‘ground truth’ data collected by stationary and mobile platforms in a soybean field in Nebraska, USA. In this effort, we tested several SIF extraction methods, including Fraunhofer Line Discrimination (FLD), improved Fraunhofer Line Discrimination (iFLD), Spectral Fitting Method (SFM), SpecFit, and a Singular Vector Decomposition (SVD) method. The SpecFit method was sensitive to the 715–740 nm water bands and removing the water bands in the fitting process yielded better agreement between the airborne and ground SIF spectra. Accurate estimation of the ground level downwelling irradiance obtained by ground measurements over a calibration target improved agreement between airborne and ground SIF retrievals at the O 2 A band, and allowed us to derive a SIF dataset with improved agreement across platforms and sampling scales. This experimental approach provided a method for generating comparable SIF signals across instruments, methods and platforms, which is critical to understanding the SIF-GPP relationship at different scales and to cross-validate the diversity of platforms used for satellite products calibration and validation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Odd-viscosity-induced passivation of Moffatt vortices

Odd viscosity is associated with a number of striking physical phenomena: a disk rotating in a viscous liquid experiences a normal compressive or tensile stress depending on its sense of rotation; chiral active matter interfacial deformations level off in a manner akin to surface tension; and unstable thermocapillary films become stabilized. In this article we show that odd viscosity decreases the size and intensity of Moffatt vortices in a corner formed by a solid substrate and a free interface. The effect is present because odd viscosity creates a normal-stress-induced shear stress at the liquid–gas interface. The new shear stress is out of phase and oppositely directed with respect to the interfacial liquid velocity. It squeezes the vortices and damps their rotation. As a consequence, Moffatt vortices can now be observed with greater facility in an experiment. Odd viscosity also increases the range of corner angles over which Moffatt vortices can be observed, and, in general, it shows a tendency to moderate singular behaviour.

Mechanics↗
DOE OSTI2022

Improved accuracy in degenerate variational integrators for guiding centre and magnetic field line flow

First-order-accurate degenerate variational integration (DVI) was introduced in Ellison et al. ( Phys. Plasmas , vol. 25, 2018, 052502) for systems with a degenerate Lagrangian, i.e. one in which the velocity-space Hessian is singular. In this paper we introduce second-order-accurate DVI schemes, both with and without non-uniform time stepping. We show that it is not in general possible to construct a second-order scheme with a preserved two-form by composing a first-order scheme with its adjoint, and discuss the conditions under which such a composition is possible. We build two classes of second-order-accurate DVI schemes. We test these second-order schemes numerically on two systems having non-canonical variables, namely the magnetic field line and guiding centre systems. Variational integration for Hamiltonian systems with non-uniform time steps, in terms of an extended phase space Hamiltonian, is generalized to non-canonical variables. It is shown that preservation of proper degeneracy leads to single-step (one-step) methods without parasitic modes, i.e. to non-uniform time step DVIs. This extension applies to second-order-accurate as well as first-order schemes, and can be applied to adapt the time stepping to an error estimate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A learned closure method applied to phase mixing in a turbulent gradient-driven gyrokinetic system in simple geometry

We present a new method for formulating closures that learn from kinetic simulation data. We apply this method to phase mixing in a simple gyrokinetic turbulent system – temperature-gradient-driven turbulence in an unsheared slab. The closure, called the learned multi-mode (LMM) closure, is constructed by, first, extracting an optimal basis from a nonlinear kinetic simulation using singular value decomposition. Subsequent nonlinear fluid simulations are projected onto this basis and the results are used to formulate the closure. We compare the closure with other closures schemes over a broad range of the relevant two-dimensional parameter space (collisionality and gradient drive). We find that the turbulent kinetic system produces phase-mixing rates much lower than the linear expectations, which the LMM closure is capable of capturing. We also compare radial heat fluxes. A Hammett–Perkins closure, generalized to include collisional effects, is quite successful throughout the parameter space, producing ${\sim }14\,\%$ root-mean-square (r.m.s.) error. The LMM closure is also very effective: when trained at three (two) points (in a 35 point parameter grid), the LMM closure produces $8\,\%$ ( $12\,\%$ ) r.m.s. errors. The LMM procedure can be readily generalized to other closure problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

On reduced modelling of the modulational dynamics in magnetohydrodynamics

This paper explores structure formation in two-dimensional magnetohydrodynamic (MHD) turbulence as a modulational instability (MI) of turbulent fluctuations. We focus on the early stages of structure formation and consider simple backgrounds that allow for a tractable model of the MI while retaining the full chain of modulational harmonics. This approach allows us to systematically examine the validity of popular closures such as the quasilinear approximation and other low-order truncations. We find that, although such simple closures can provide quantitatively accurate approximations of the MI growth rates in some regimes, they can fail to capture the modulational dynamics in adjacent regimes even qualitatively, falsely predicting MI when the system is actually stable. We find that this discrepancy is due to the excitation of propagating spectral waves (PSWs) which can ballistically transport energy along the modulational spectrum, unimpeded until dissipative scales, thereby breaking the feedback loops that would otherwise sustain MIs. The PSWs can be self-maintained as global modes with real frequencies and drain energy from the primary structure at a constant rate until the primary structure is depleted. To describe these waves within a reduced model, we propose an approximate spectral closure that captures them and MIs on the same footing. We also find that introducing corrections to ideal MHD, conservative or dissipative, can suppress PSWs and reinstate the accuracy of the quasilinear approximation. In this sense, ideal MHD is a ‘singular’ system that is particularly sensitive to the accuracy of the closure within mean-field models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗