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,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Exploring a possible origin of the QCD critical endpoint

We develop a new model of the QCD critical endpoint by matching the deconfinement phase transition line of quark-gluon bags with the similar line at which the bag surface tension coefficient vanishes. Unlike all previous studies of such models the deconfined phase in our approach is defined not by an essential singularity of the isobaric partition function but its simple pole. As an unexpected result we find out that the first-order phase transition which is usually defined by a discontinuity of the first derivative of the bag system pressure results from a discontinuity of the derivative of surface tension coefficient of quark-gluon bags.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Variational perturbation theory and nonperturbative calculations in QCD

A nonperturbative approach based on the variational perturbation theory in quantum chromodynamics is developed. The variational series is different from the conventional perturbative expansion and can be used to go beyond the weak-coupling regime. The approach suggested takes into account the summation of threshold singularities and the involvement of nonperturbative light quark masses. Phenomenological applications of this approach to describe physical quantities connected with the hadronic {tau}-decay data: the R{sub {tau}} ratio, the light-quark Adler function, and the smeared R{sub {Delta}} function are presented. The description of examined quantities includes an infrared region and, therefore, they cannot be directly calculated within the standard perturbation theory. It is shown that in spite of this fact the approach suggested gives a rather good result for these quantities down to the lowest energy scale.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

The early universe and cosmogenesis

We extrapolate the Cosmological Standard Model to the past, determine initial geometrical conditions in the early universe, and consider a new cosmogenesis paradigm based on the concept of black-and-white holes with integrable singularities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Noncompact Riemannian spaces with G{sub 2}, Spin(7) and SU(2m) holonomies

The survey contains a description of the author’s results on constructing Riemannian metrics with exceptional holonomies on noncompact manifolds. We give a detailed description of how to reduce a condition of the holonomy to be exceptional to a system of ODEs. We describe the behavior of the trajectories and the ways to resolve singularities of such systems, and also study concrete examples. Introduction and first chapter contain necessary definitions and well-known facts about holonomies. Conclusion contains a description of other results proven with the help of the same technique and a minimal survey on Lorentzian holonomies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Iterative algorithms for processing experimental data

The need to solve linear and nonlinear integral equations arise, e.g., in recovering plasma parameters from the data of multichannel diagnostics. The paper presents an iterative method for solving integral equations with a singularity at the upper limit of integration. The method consists in constructing successive approximations and calculating the integral by quadrature formulas in each integration interval. An example of application of the iterative algorithm to numerically solve an integral equation similar to those arising in recovering the plasma density profile from reflectometry data is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Density of states of narrow superconducting channels in the regime of quantum fluctuations of the order parameter

The current–voltage characteristics of superconductor–insulator–semiconductor (S1–I–S2) tunnel junctions, where superconducting electrode S2 is a thin nanowire, are studied experimentally. The observed blurring of the gap singularities is interpreted as a manifestation of the order parameter quantum fluctuations. We propose a model taking into account the broadening of the density of states due to the interaction of electrons with the Mooij–Schön plasmon mode emerging in a quasi-one-dimensional superconducting channel in the regime of quantum fluctuations of the order parameter. The model gives results that are in a reasonable qualitative agreement with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Features of the Percolation Scheme of Vibrational Spectrum Reconstruction in the Ga{sub 1–} {sub x}Al{sub x}P Alloy

Specific features of the properties of Ga–P lattice vibrations have been investigated using the percolation model of a mixed Ga{sub 1–} {sub x}Al{sub x}P crystal (alloy) with zero lattice mismatch between binary components of the alloy. In contrast to other two-mode alloy systems, in Ga{sub 1–} {sub x}Al{sub x}P a percolation splitting of δ ~ 13 cm{sup –1} is observed for the low-frequency mode of GaP-like vibrations. An additional GaP mode (one of the percolation doublet components) split from the fundamental mode is observed for the GaP-rich alloy, which coincides in frequency with the gap corresponding to the zero density of one-phonon states of the GaP crystal. The vibrational spectrum of impurity Al in the GaP crystal has been calculated using the theory of crystal lattice dynamics. Upon substitution of lighter Al for the Ga atom, the calculated spectrum includes, along with the local mode, a singularity near the gap with the zero density of phonon states of the GaP crystal, which coincides with the mode observed experimentally at a frequency of 378 cm{sup –1} in the Ga{sub 1–} {sub x}Al{sub x}P (x < 0.4) alloy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Modified Eddington-inspired-Born-Infeld gravity with a trace term

In this study, a modified Eddington-inspired-Born-Infeld (EiBI) theory with a pure trace term g μν R being added to the determinantal action is analysed from a cosmological point of view. It corresponds to the most general action constructed from a rank two tensor that contains up to first order terms in curvature. This term can equally be seen as a conformal factor multiplying the metric g μν . This very interesting type of amendment has not been considered within the Palatini formalism despite the large amount of works on the Born-Infeld-inspired theory of gravity. This model can provide smooth bouncing solutions which were not allowed in the EiBI model for the same EiBI coupling. Most interestingly, for a radiation filled universe there are some regions of the parameter space that can naturally lead to a de Sitter inflationary stage without the need of any exotic matter field. Finally, in this model we discover a new type of cosmic “quasi-sudden” singularity, where the cosmic time derivative of the Hubble rate becomes very large but finite at a finite cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Continuous Integration and Delivery for HPC

Continuous integration, delivery, and deployment (CICD) is widely used in DevOps communities, as it allows for teams of all sizes to deploy rapidly-changing hardware and software resources quickly and confidently. In this paper, we will describe how University of Colorado Boulder Research Computing has adopted these practices on the RMACC Summit supercomputer to allow system engineers and researchers alike to capitalize on the benefits of CICD-centric development workflows. We will introduce the topic of CICD at a high level and describe how such practices can ease common software management challenges for High-Performance Computing (HPC) resources. We will then document the infrastructure deployed for Summit, and explain how software such as Jenkins and Singularity enabled adaptation for an HPC environment. We will conclude with two case studies discussing the use of our CICD infrastructure: one case study from the perspective of a system engineer maintaining user-facing resources, and the other case study from the perspective of a researcher developing, maintaining, and using the MFiX-Exa codebase.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A parallel hierarchical blocked adaptive cross approximation algorithm

This article presents a low-rank decomposition algorithm based on subsampling of matrix entries. The proposed algorithm first computes rank-revealing decompositions of submatrices with a blocked adaptive cross approximation (BACA) algorithm, and then applies a hierarchical merge operation via truncated singular value decompositions (H-BACA). The proposed algorithm significantly improves the convergence of the baseline ACA algorithm and achieves reduced computational complexity compared to the traditional decompositions such as rank-revealing QR. Numerical results demonstrate the efficiency, accuracy, and parallel scalability of the proposed algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Decadal Prediction and Stochastic Simulation of Hydroclimate Over Monsoonal Asia

We developed further our advanced methods of time series analysis and empirical model reduction (EMR) and applied them to climatic time series relevant to hydroclimate over Monsoonal Asia. The EMR methodology was both generalized further and laid on a rigorous mathematical basis via multilayered stochastic models (MSMs). We identified easily testable conditions that imply the existence of a global random attractor for MSMs and allow for non-polynomial predictors. This existence, in turn, guarantees the numerical stability of the MSMs so obtained. We showed that, in the presence of low-frequency variability (LFV), EMR prediction can be improved further by including information from selected times in the system’s past. This prediction method, dubbed Past-Noise Forecasting (PNF), was successfully applied to the Madden-Julian Oscillation (MJO). Our time series analysis and forecasting methods, based on singular-spectrum analysis (SSA) and its enhancements, were applied to several multi-centennial proxy records provided by the Lamont team. These included the Palmer Drought Severity Index (PDSI) for 1300–2005 from the Monsoonal Asia Drought Atlas (MADA), and a 300-member ensemble of pseudo-reconstructions of Indus River discharge for 1702–2005. The latter was shown to exhibit a robust 27-yr low-frequency mode, which helped multi-decadal retroactive forecasts with no look-ahead over this 300-year interval.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

What do we mean by the word “Shock”?

From one vantage point, a shock is a continuous but drastic change in state variables that occurs over very small time and length scales. These scales and associated changes in state variables can be measured experimentally. From another vantage point, a shock is a mathematical singularity consisting of instantaneous changes in state variables. This more mathematical view gives rise to analytical solutions to idealized problems. And from a third vantage point, a shock is a structure in a hydrocode prediction. Its width depends on the simulation’s grid resolution and artificial viscosity. These three vantage points can be in conflict when ideas from the associated fields are combined, and yet combining them is an important goal of an integrated modeling program. This presentation explores an example of how models for real materials in the presence of real shocks react to a hydrocode’s numerical shocks of finite width. The presentation will include an introduction to plasticity for the novice, an historical view of plasticity algorithms, a demonstration of how pursuing the meaning of “shock” has resulted in hydrocode improvements, and will conclude by answering some of the questions that arise from that pursuit. After the technical part of the presentation, a few slides advertising LANL’s Computational Physics Student Summer Workshop will be shown.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Limited-memory adaptive snapshot selection for proper orthogonal decomposition

Reduced order models are useful for accelerating simulations in many-query contexts, such as optimization, uncertainty quantification, and sensitivity analysis. However, offline training of reduced order models can have prohibitively expensive memory and floating-point operation costs in high-performance computing applications, where memory per core is limited. To overcome this limitation for proper orthogonal decomposition, we propose a novel adaptive selection method for snapshots in time that limits offline training costs by selecting snapshots according an error control mechanism similar to that found in adaptive time-stepping ordinary differential equation solvers. The error estimator used in this work is related to theory bounding the approximation error in time of proper orthogonal decomposition-based reduced order models, and memory usage is minimized by computing the singular value decomposition using a single-pass incremental algorithm. Results for a viscous Burgers’ test problem demonstrate convergence in the limit as the algorithm error tolerances go to zero; in this limit, the full order model is recovered to within discretization error. The resulting method can be used on supercomputers to generate proper orthogonal decomposition-based reduced order models, or as a subroutine within hyperreduction algorithms that require taking snapshots in time, or within greedy algorithms for sampling parameter space.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Fracture Capabilities in Grizzly with the extended Finite Element Method (X-FEM)

Efforts are underway to develop fracture mechanics capabilities in the Grizzly code to enable it to be used to perform deterministic fracture assessments of degraded reactor pressure vessels (RPVs). A capability was previously developed to calculate three-dimensional interaction- integrals to extract mixed-mode stress-intensity factors. This capability requires the use of a finite element mesh that conforms to the crack geometry. The eXtended Finite Element Method (X-FEM) provides a means to represent a crack geometry without explicitly fitting the finite element mesh to it. This is effected by enhancing the element kinematics to represent jump discontinuities at arbitrary locations inside of the element, as well as the incorporation of asymptotic near-tip fields to better capture crack singularities. In this work, use of only the discontinuous enrichment functions was examined to see how accurate stress intensity factors could still be calculated. This report documents the following work to enhance Grizzly’s engineering fracture capabilities by introducing arbitrary jump discontinuities for prescribed crack geometries; X-FEM Mesh Cutting in 3D: to enhance the kinematics of elements that are intersected by arbitrary crack geometries, a mesh cutting algorithm was implemented in Grizzly. The algorithm introduces new virtual nodes and creates partial elements, and then creates a new mesh connectivity; Interaction Integral Modifications: the existing code for evaluating the interaction integral in Grizzly was based on the assumption of a mesh that was fitted to the crack geometry. Modifications were made to allow for the possibility of a crack front that passes arbitrarily through the mesh; and Benchmarking for 3D Fracture: the new capabilities were benchmarked against mixed-mode three-dimensional fracture problems with known analytical solutions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Structural Aspects of Neutron Survival Probabilities

The neutron survival probability (and related quantities including probabilities of extinction and initiation) is a central element of the broader stochastic theory of neutron populations and finds application in fields including reactor start-up, analysis of reactor power bursts and criticality accidents, and safeguards. In a full neutron transport formulation, the equation governing the single-neutron survival probability is a backward or adjoint-like integro-partial differential equation with the added complexity of being highly nonlinear. Analogous formulations of this equation exist in the context of many approximate theories of neutron transport, with the point kinetics formulation having received significant theoretical attention since the 1940s. This work continues this tradition by providing a novel analysis of the single-neutron survival probability equation using the tools of boundary layer theory. The analysis reveals that the “fully dynamic” solution of the single-neutron survival probability equation—and some key probability distributions derived from it—may be cast as a singular perturbation around the underlying quasi-static single-neutron probability of initiation. In this perturbation solution, the expansion parameter is the ratio of the neutron generation time to a macroscopic time scale characterizing the overall system evolution; this interpretation illuminates some of the fundamental structural aspects of neutron survival phenomena.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Machine Learning-Based Model Predictive Control of Two-Time-Scale Systems

In this study, we present a general form of nonlinear two-time-scale systems, where singular perturbation analysis is used to separate the dynamics of the slow and fast subsystems. Machine learning techniques are utilized to approximate the dynamics of both subsystems. Specifically, a recurrent neural network (RNN) and a feedforward neural network (FNN) are used to predict the slow and fast state vectors, respectively. Moreover, we investigate the generalization error bounds for these machine learning models approximating the dynamics of two-time-scale systems. Next, under the assumption that the fast states are asymptotically stable, our focus shifts toward designing a Lyapunov-based model predictive control (LMPC) scheme that exclusively employs the RNN to predict the dynamics of the slow states. Additionally, we derive sufficient conditions to guarantee the closed-loop stability of the system under the sample-and-hold implementation of the controller. A nonlinear chemical process example is used to demonstrate the theory. In particular, two RNN models are constructed: one to model the full two-time-scale system and the other to predict solely the slow state vector. Both models are integrated within the LMPC scheme, and we compare their closed-loop performance while assessing the computational time required to execute the LMPC optimization problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A SECOND-ORDER DIVERGENCE-CONSTRAINED MULTIDIMENSIONAL NUMERICAL SCHEME FOR RELATIVISTIC TWO-FLUID ELECTRODYNAMICS

A new multidimensional simulation code for relativistic two-fluid electrodynamics (RTFED) is described. The basic equations consist of the full set of Maxwell’s equations coupled with relativistic hydrodynamic equations for separate two charged fluids, representing the dynamics of either an electron–positron or an electron–proton plasma. It can be recognized as an extension of conventional relativistic magnetohydrodynamics (RMHD). Finite resistivity may be introduced as a friction between the two species, which reduces to resistive RMHD in the long wavelength limit without suffering from a singularity at infinite conductivity. A numerical scheme based on HLL (Harten–Lax–Van Leer) Riemann solver is proposed that exactly preserves the two divergence constraints for Maxwell’s equations simultaneously. Several benchmark problems demonstrate that it is capable of describing RMHD shocks/discontinuities at long wavelength limit, as well as dispersive characteristics due to the two-fluid effect appearing at small scales. This shows that the RTFED model is a promising tool for high energy astrophysics application.

79 ASTRONOMY AND ASTROPHYSICS↗