Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SECOND”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Superconvergence results of Legendre spectral projection methods for Volterra integral equations of second kind

In this paper, Legendre spectral projection methods are applied for the Volterra integral equations of second kind with a smooth kernel. We prove that the approximate solutions of the Legendre Galerkin and Legendre collocation methods converge to the exact solution with the order O(n{sup −r}) in L{sup 2}-norm and order O(n{sup −r+(1/2)}) in infinity norm, and the iterated Legendre Galerkin solution converges with the order O(n{sup −2r}) in both L{sup 2}-norm and infinity norm, whereas the iterated Legendre collocation solution converges with the order O(n{sup −r}) in both L{sup 2}-norm and infinity norm, n being the highest degree of Legendre polynomials employed in the approximation and r being the smoothness of the kernels. We have also considered multi-Galerkin method and its iterated version, and prove that the iterated multi-Galerkin solution converges with the order O(n{sup −3r}) in both infinity and L{sup 2} norm. Numerical examples are given to illustrate the theoretical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Construction of the Nordsieck second derivative methods with RK stability for stiff ODEs

In this paper, we study the construction and implementation of special Nordsieck second derivative general linear methods of order p and stage order q=p in which the number of input and output values is r=p rather than r=p+1. We will construct A- and L-stable methods of orders three and four in this form with Runge–Kutta stability properties. The efficiency of the constructed methods and reliability of the proposed error estimates are shown by implementing of the methods in a variable stepsize environment on some well-known stiff problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Integrated modeling of second phase precipitation in cold-worked 316 stainless steels under irradiation

The current work combines the Cluster Dynamics (CD) technique and CALPHAD-based precipitation modeling to address the second phase precipitation in cold-worked (CW) 316 stainless steels (SS) under irradiation at 300–400 °C. CD provides the radiation enhanced diffusion and dislocation evolution as inputs for the precipitation model. The CALPHAD-based precipitation model treats the nucleation, growth and coarsening of precipitation processes based on classical nucleation theory and evolution equations, and simulates the composition, size and size distribution of precipitate phases. We benchmark the model against available experimental data at fast reactor conditions (9.4 × 10 –7 dpa/s and 390 °C) and then use the model to predict the phase instability of CW 316 SS under light water reactor (LWR) extended life conditions (7 × 10 –8 dpa/s and 275 °C). The model accurately predicts the γ' (Ni 3 Si) precipitation evolution under fast reactor conditions and that the formation of this phase is dominated by radiation enhanced segregation. The model also predicts a carbide volume fraction that agrees well with available experimental data from a PWR reactor but is much higher than the volume fraction observed in fast reactors. We propose that radiation enhanced dissolution and/or carbon depletion at sinks that occurs at high flux could be the main sources of this inconsistency. The integrated model predicts ~1.2% volume fraction for carbide and ~3.0% volume fraction for γ' for typical CW 316 SS (with 0.054 wt% carbon) under LWR extended life conditions. Finally, this work provides valuable insights into the magnitudes and mechanisms of precipitation in irradiated CW 316 SS for nuclear applications.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Energy levels, weighted oscillator strengths, transition rates, lifetimes, hyperfine structures and isotope shifts of He-like-Sn deduced from the relativistic multiconfiguration Dirac–Hartree–Fock and second-order many-body perturbation theory calculations

We disclose relativistic multiconfiguration Dirac–Hartree–Fock (MCDHF) spectrum calculations for SnXLIX. Energy levels, weighted oscillator strengths, isotope shifts, hyperfine structure and Landé g{sub J} factors are calculated for 127 odd- and even-parity states as well as lifetimes and transition rates between these states. To scrutinize the accuracy of our results, we have implemented parallel calculations using a Flexible Atomic Code (FAC) by introducing the Second-Order Many-Body Perturbation Theory (MBPT) method. Additionally, the Breit interaction and leading quantum electrodynamic effects (QED) are included as perturbations in extensive relativistic configuration interaction (RCI) calculations. We signal that, our calculations for SnXLIX are made for the first time and they provide to date the most accurate and complete atomic data related to this ion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Daily grass reference evapotranspiration with Meteosat Second Generation shortwave radiation and reference ET products

Here we assess the accuracy of estimating daily grass reference evapotranspiration (PM-ET o ) using daily shortwave radiation (R s ) and reference evapotranspiration (ET REF ) products provided by the Meteosat Second Generation (MSG) geostationary satellite delivered by the Satellite Applications Facility on Land Surface Analysis (LSA-SAF) framework. The accuracy of using reanalysis ERA5 shortwave radiation data (R s ERA5 ) provided by the European Center for Medium-Range Weather Forecasts (ECMWF) is also evaluated. The assessments were performed using observed weather variables at 37 weather stations distributed across continental Portugal, where climate conditions range from semi-arid to humid, and 12 weather stations located in Azores islands, characterized by humid, windy and often cloudy conditions. This study&#x27;s use of data from a variety of climate conditions contributed to a unique and innovative assessment of the usability of LSA-SAF and ERA5 products for ET o estimation. The first assessment focused on comparing LSA-SAF estimates of R s (R s LSA-SAF ) against ground stations (R s ground ). The results showed a good matching between the two R s data sets for continental Portugal but a tendency for R s LSA-SAF to under-estimate R s ground in the cloudy islands of Azores. ET o values computed using R s LSA-SAF data and observed temperature, humidity and wind speed (ET o LSA-SAF ) were then compared with PM-ET o estimates with ground-based data, which were used as benchmark; input data of temperature and humidity needed for PM-ET o were quality checked for surface aridity effects. It was observed that ET o LSA-SAF is strongly correlated with PM-ET o (R 2 > 0.97) for most locations in continental Portugal, with regression coefficient of a linear regression forced to the origin ranging between 0.95 and 1.05, mean root mean square error (RMSE) of 0.13 mm d –1 , and Nash and Sutcliff efficiency of modeling (EF) above 0.95. For most Azores locations, ET o LSA-SAF over-estimated PM-ET o . This is likely a consequence of the high spatio-temporal heterogeneity of weather conditions that occur in these oceanic islands together with the different footprints of satellite (averaged over the pixel) and station observations. Reanalysis ERA5 shortwave radiation data presented similar behavior to the LSA-SAF products, however with slightly lower accuracy. The daily LSA-SAF ET REF product (ET REF LSA-SAF ) was assessed and results have shown a good accuracy of this product, with acceptable RMSE and high EF values, for continental Portugal but a low accuracy for the Azores islands. A simplified bias correction approach was shown to improve both ET o derived from the LSA-SAF products, namely for Azores stations, which seem to be representative of smaller areas. The use of the FAO-PM temperature approach (PMT) was also assessed using the R s LSA-SAF and R s ERA5 data, which showed a superiority of the LSA-SAF product for ET o estimations (ET o PMT LSA-SAF ). No significant differences (p < 0.05) were observed in terms of the median value of the RMSE when adopting ET o PMT and ET REF LSA-SAF . Differently, results showed that using the R s LSA-SAF in the PMT approach (ET o PMT LSA-SAF ) produces significantly better RMSE results than ET o PMT and ET REF LSA-SAF . Overall, the performed assessment allows concluding that the use of R s LSA-SAF, and to a lesser extent the use of the R s ERA5 , highly improves the accuracy of computation of ET o when R s observations are not available, including when only temperature data are accessible. The use of the ET REF LSA-SAF product is a good alternative when observed weather data are not available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

The second law of thermodynamics under unitary evolution and external operations

The von Neumann entropy cannot represent the thermodynamic entropy of equilibrium pure states in isolated quantum systems. The diagonal entropy, which is the Shannon entropy in the energy eigenbasis at each instant of time, is a natural generalization of the von Neumann entropy and applicable to equilibrium pure states. We show that the diagonal entropy is consistent with the second law of thermodynamics upon arbitrary external unitary operations. In terms of the diagonal entropy, thermodynamic irreversibility follows from the facts that quantum trajectories under unitary evolution are restricted by the Hamiltonian dynamics and that the external operation is performed without reference to the microscopic state of the system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

New implicitly solvable potential produced by second order shape invariance

The procedure proposed recently by Bougie et al. (2010) to study the general form of shape invariant potentials in one-dimensional Supersymmetric Quantum Mechanics (SUSY QM) is generalized to the case of Higher Order SUSY QM with supercharges of second order in momentum. A new shape invariant potential is constructed by this method. It is singular at the origin, it grows at infinity, and its spectrum depends on the choice of connection conditions in the singular point. The corresponding Schrödinger equation is solved explicitly: the wave functions are constructed analytically, and the energy spectrum is defined implicitly via the transcendental equation which involves Confluent Hypergeometric functions. - Highlights: • New potential with 2nd order irreducible shape invariance was constructed. • The connection conditions at the singularity of potential were obtained. • The explicit expressions for all wave functions were derived. • The implicit equation for the energy spectrum was obtained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Dynamical centrosymmetry breaking — A novel mechanism for second harmonic generation in graphene

We discover an unusual phenomenon that occurs when a graphene monolayer is illuminated by a short and intense pulse at normal incidence. Due to the pulse-induced oscillations of the Dirac cones, a dynamical breaking of the layer’s centrosymmetry takes place, leading to the generation of second harmonic waves. We prove that this result can only be found by using the full Dirac equation and show that the widely used semiconductor Bloch equations fail to reproduce this and some other important physics of graphene. Our results open new windows in the understanding of nonlinear light-matter interactions in a wide variety of new 2D materials with a gapped or ungapped Dirac-like dispersion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Efficiency enhancement of slow-wave electron-cyclotron maser by a second-order shaping of the magnetic field in the low-gain limit

Based on the anomalous Doppler effect, we put forward a proposal to enhance the conversion efficiency of the slow-wave electron cyclotron masers (ECM) under the resonance condition. Compared with previous studies, we add a second-order shaping term in the guild magnetic field. Theoretical analyses and numerical calculations show that it can enhance the conversion efficiency in the low-gain limit. The case of the initial velocity spread of electrons satisfying the Gaussian distribution is also analysed numerically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Development and validation of a second-order thermal network model for residential buildings

Heating, Ventilation, and Air Conditioning (HVAC) systems can maintain the space air temperature of residential buildings, either directly by heating/cooling the air, or indirectly via heat transfer to and from the building structure that acts as a thermal mass. Hence, HVAC systems can help achieve load shifting, peak load reduction, and/or energy cost saving, thus enabling grid-interactive HVAC operation. A home thermal model that can accurately reflect the dynamics of the space air and interior wall surface temperatures, is therefore valuable. This paper develops such a model using the standard RC (resistance-capacitance) approach. The model contains a virtual envelope node and an internal space node and is thus second-order. A hybrid parameter identification scheme, made up of the least-squares and optimal search methods, is also developed. The proposed model and scheme were validated using data collected from a test home. It was found that a modest amount of training data was sufficient to yield reliable parameter estimates and accurate prediction. It was also found that when making 24-hour-ahead prediction of the space air temperature, both methods had comparable performances when the training data began in a transition season. However, when they began in an HVAC season, the optimal search method performed better. Furthermore, the least-squares method is recommended during a transition season due to its lower computational burden, while the optimal search method is recommended during an HVAC season due to its better estimation performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2012

Structural evolution and tissue-specific expression of tetrapod-specific second isoform of secretory pathway Ca{sup 2+}-ATPase

Highlights: Black-Right-Pointing-Pointer Full-length secretory pathway Ca-ATPase (SPCA2) cloned from rat duodenum. Black-Right-Pointing-Pointer ATP2C2 gene (encoding SPCA2) exists only in genomes of Tetrapoda. Black-Right-Pointing-Pointer Rat and pig SPCA2 are expressed in intestines, lung and some secretory glands. Black-Right-Pointing-Pointer Subcellular localization of SPCA2 may depend on tissue type. Black-Right-Pointing-Pointer In rat duodenum, SPCA2 is localized in plasma membrane-associated compartments. -- Abstract: Secretory pathway Ca-ATPases are less characterized mammalian calcium pumps than plasma membrane Ca-ATPases and sarco-endoplasmic reticulum Ca-ATPases. Here we report analysis of molecular evolution, alternative splicing, tissue-specific expression and subcellular localization of the second isoform of the secretory pathway Ca-ATPase (SPCA2), the product of the ATP2C2 gene. The primary structure of SPCA2 from rat duodenum deduced from full-length transcript contains 944 amino acid residues, and exhibits 65% sequence identity with known SPCA1. The rat SPCA2 sequence is also highly homologous to putative human protein KIAA0703, however, the latter seems to have an aberrant N-terminus originating from intron 2. The tissue-specificity of SPCA2 expression is different from ubiquitous SPCA1. Rat SPCA2 transcripts were detected predominantly in gastrointestinal tract, lung, trachea, lactating mammary gland, skin and preputial gland. In the newborn pig, the expression profile is very similar with one remarkable exception: porcine bulbourethral gland gave the strongest signal. Upon overexpression in cultured cells, SPCA2 shows an intracellular distribution with remarkable enrichment in Golgi. However, in vivo SPCA2 may be localized in compartments that differ among various tissues: it is intracellular in epidermis, but enriched in plasma membranes of the intestinal epithelium. Analysis of SPCA2 sequences from various vertebrate species argue that ATP2C2 gene radiated from ATP2C1 (encoding SPCA1) during adaptation of tetrapod ancestors to terrestrial habitats.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Flow-driven spectral chaos (FSC) method for long-time integration of second-order stochastic dynamical systems

For decades, uncertainty quantification techniques based on the spectral approach have been demonstrated to be computationally more efficient than the Monte Carlo method for a wide variety of problems, particularly when the dimensionality of the probability space is relatively low. The time-dependent generalized polynomial chaos (TD-gPC) is one such technique that uses an evolving orthogonal basis to better represent the stochastic part of the solution space in time. Here in this paper, we present a new numerical method that uses the concept of enriched stochastic flow maps to track the evolution of the stochastic part of the solution space in time. The computational cost of this proposed flow-driven stochastic chaos (FSC) method is an order of magnitude lower than TD-gPC for comparable solution accuracy. This gain in computational cost is realized because, unlike most existing methods, the number of basis vectors required to track the stochastic part of the solution space does not depend upon the dimensionality of the probability space. Four representative numerical examples are presented to demonstrate the performance of the FSC method for long-time integration of second-order stochastic dynamical systems in the context of stochastic dynamics of structures.

FSC↗
DOE OSTI2018

Unconditionally stable, second-order schemes for gradient-regularized, non-convex, finite-strain elasticity modeling martensitic phase transformations

In the setting of continuum elasticity martensitic phase transformations are characterized by a non-convex free energy density function that possesses multiple wells in strain space and includes higher-order gradient terms for regularization. Metastable martensitic microstructures, defined as solutions that are local minimizers of the total free energy, are of interest and are obtained as steady state solutions to the resulting transient formulation of Toupin’s gradient elasticity at finite strain. This type of problem poses several numerical challenges including stiffness, the need for fine discretization to resolve microstructures, and following solution branches. Accurate time-integration schemes are essential to obtain meaningful solutions at reasonable computational cost. Here in this work we introduce two classes of unconditionally stable second-order time-integration schemes for gradient elasticity, each having relative advantages over the other. Numerical examples are shown highlighting these features.

42 ENGINEERING↗
DOE OSTI2020

The second-generation Shifted Boundary Method and its numerical analysis

Recently, the Shifted Boundary Method (SBM) was proposed within the class of unfitted (or immersed, or embedded) finite element methods. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we present enhanced variational SBM formulations for the Poisson and Stokes problems with improved flexibility and robustness. These simplified variational forms allow to relax some of the assumptions required by the mathematical proofs of stability and convergence of earlier implementations. First, we show that these new SBM implementations can be proved asymptotically stable and convergent even without the rather restrictive assumption that the inner product between the normals to the true and surrogate boundaries is positive. Second, we show that it is not necessary to introduce a stabilization term involving the tangential derivatives of the solution at Dirichlet boundaries, therefore avoiding the calibration of an additional stabilization parameter. Finally, we prove enhanced L 2 -estimates without the cumbersome assumption – of earlier proofs – that the surrogate domain is convex. Instead we rely on a conventional assumption that the boundary of the true domain is smooth, which can also be replaced by requiring convexity of the true domain. The aforementioned improvements open the way to a more general and efficient implementation of the Shifted Boundary Method, particularly in complex three-dimensional geometries. We complement these theoretical developments with numerical experiments in two and three dimensions.

42 ENGINEERING↗
DOE OSTI2021

Second-order invariant domain preserving approximation of the compressible Navier–Stokes equations

Here, we present a fully discrete approximation technique for the compressible Navier–Stokes equations that is second-order accurate in time and space, semi-implicit, and guaranteed to be invariant domain preserving. The restriction on the time step is the standard hyperbolic CFL condition, i.e. τ ≲ O ( h ) / V where V is some reference velocity scale and h the typical meshsize.

42 ENGINEERING↗
DOE OSTI2016

Massively parallel and linear-scaling algorithm for second-order Moller–Plesset perturbation theory applied to the study of supramolecular wires

Here, we present a scalable cross-platform hybrid MPI/OpenMP/OpenACC implementation of the Divide–Expand–Consolidate (DEC) formalism with portable performance on heterogeneous HPC architectures. The Divide–Expand–Consolidate formalism is designed to reduce the steep computational scaling of conventional many-body methods employed in electronic structure theory to linear scaling, while providing a simple mechanism for controlling the error introduced by this approximation. Our massively parallel implementation of this general scheme has three levels of parallelism, being a hybrid of the loosely coupled task-based parallelization approach and the conventional MPI +X programming model, where X is either OpenMP or OpenACC. We demonstrate strong and weak scalability of this implementation on heterogeneous HPC systems, namely on the GPU-based Cray XK7 Titan supercomputer at the Oak Ridge National Laboratory. Using the “resolution of the identity second-order Moller–Plesset perturbation theory” (RI-MP2) as the physical model for simulating correlated electron motion, the linear-scaling DEC implementation is applied to 1-aza-adamantane-trione (AAT) supramolecular wires containing up to 40 monomers (2440 atoms, 6800 correlated electrons, 24 440 basis functions and 91 280 auxiliary functions). This represents the largest molecular system treated at the MP2 level of theory, demonstrating an efficient removal of the scaling wall pertinent to conventional quantum many-body methods.

97 MATHEMATICS AND COMPUTING↗