Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical stability”

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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Updates to Implicit Edge-Based Gradient Methods

In this paper, we report updates to the implicit edge-based gradient methods originally introduced in [H. Nishikawa, AIAA Paper 2020-3048, 2020]. First, we clarify the relationship between gradient accuracy and truncation error and show that the quadratic method involves a free parameter. Then, we provide a complete description and a simplified matrix form of the implicit gradient systems including a consistent boundary treatment, and derive a set of parameters for achieving fourth-order gradient accuracy on regular tetrahedral grids. A stability analysis is performed for a relaxation scheme used to solve the implicit gradient systems, and the result serves as a guide for choosing parameters. Numerical results are shown for accuracy verification and also for realistic inviscid flow problems in three dimensions, including flows with shock waves.

Weighted Least-Squares↗
DOE OSTI2025

Laser-Induced Breakdown Spectroscopy (LIBS) Sensing for Environmental and Subsurface Monitoring

Groundwater monitoring is essential to timely and accurately reflect the current situations and the trends of water quality. However, the long-term stability and survival of any potential in-situ monitoring method is threatened by the harsh conditions and the groundwater monitoring in downhole environments poses numerous challenges to the sensor community. Laser induced breakdown spectroscopy (LIBS) has been demonstrated as a promising technology for chemical monitoring in high-temperature high-pressure (HTHP) environments and hard to reach places. The technique demonstrates several advantages, such as rapid, real-time, in-situ, and simultaneous detection of multiple elements with simple or no sample preparation. This chapter includes a brief review of the field-portable LIBS systems, and development of a compact, robust, and simple LIBS instrument for downhole HTHP water quality monitoring. The fieldable prototype sensor is tested in an onsite monitoring well where trace elements’ concentrations are tracked over an extended period. The testing has verified that the fiber coupled design performs as desired. The system shows good calibration linearity for tested elements and collection times, and Limits of Detection (LODs) that are comparable to those of tabletop LIBS instruments. In addition to groundwater quality monitoring, the fabricated LIBS-based sensor could have widespread sub-surface detection applications.

Jain, Jinesh [NETL Site Support Contractor, Nation↗
NASA NTRS1999

The Effect of Viscosity on the Spherical Stability of Oscillating Gas Bubbles

Gas bubbles driven in radial oscillations are subject to an instability of the spherical shape that is opposed by surface tension and viscosity. An exact linear formulation for the study of the phenomenon has been available for many years, but its complexity has discouraged a detailed investigation. With the recent theory of sonoluminescence of Lohse and co-workers, there has been a renewed interest in the problem and new data have become available. This paper presents a numerical method for the solution of the pertinent equations and compares the theory with these new data. The coupling of the strong nonlinearity of the bubble radial oscillations with the parametric mechanism of the surface instability results in a very complex structure for the stability boundary. Nevertheless, a good agreement between theory and data is found. A comparison with earlier approximate models is also made.

Hao, Y.↗
NASA NTRS1985

Effects of stratospheric lapse rate on numerically simulated thunderstorm cloud top structure

The study of overshooting midlatitude severe thunderstorm tops has led to the detection of some intriguing patterns. This paper provides a summary of an investigation of the effects of stratospheric stability on thunderstorm cloud top structure, taking into account the utilization of a fully three-dimensional (3D) anelastic model discussed by Schlesinger (1975, 1978, 1984). The layout of model experiments is discussed, giving attention to initial environmental soundings, and initial forcing. The time variation of the vertical velocity extreme is examined, and a description of the mature storm structure is presented. Constant-altitude plane views throughout storm depth are considered along with cloud top height fields. With the stratospheric lapse, the warm region is small and breaks down into separate spots with secondary cold spots downshear of them.

Schlesinger, R. E.↗
DOE OSTI2024

Performance stability of plastics for neutron-gamma pulse shape discrimination

The introduction of the first commercially available plastic scintillators with pulse shape discrimination (PSD) offered by Eljen Technology marked progress towards potentially replacing liquid scintillators in neutron detection. However, use of these plastics over several recent years has revealed an important flaw in these materials: the eventual degradation of scintillation light output and PSD performance. Here, studies described in this paper considered possible reasons for this degradation. Experiments conducted with numerous lab-prepared and commercial EJ-276 samples showed that the main factor affecting the instability is oxidation that involves highly reactive radicals generated during the polymerization process or from the further breakdown of polymer chains under oxygen/air exposure. Based on the obtained results, the stability of scintillation performance for PPO (2,5-Diphenyloxazole)-based PSD plastics has been improved through modifications of the composition via the utilization of scintillation dyes and compounds with antioxidant properties that diminish the effects of oxidation. Elements of these studies were used in the commercial production of the most recent EJ-276D PSD plastic version for fast neutron detection and 6Li-loaded plastics for thermal neutron and antineutrino detection applications. Performance projections of the new PSD plastics indicate a likely degradation of less than 10 % over 5–10 years in comparison to previous EJ-276 that might lose up to 30–40 % of the scintillation light during 1–2 years of storage or deployment under ambient conditions.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Numerical simulation of conservation laws

A new numerical framework for solving conservation laws is being developed. This new approach differs substantially from the well established methods, i.e., finite difference, finite volume, finite element and spectral methods, in both concept and methodology. The key features of the current scheme include: (1) direct discretization of the integral forms of conservation laws, (2) treating space and time on the same footing, (3) flux conservation in space and time, and (4) unified treatment of the convection and diffusion fluxes. The model equation considered in the initial study is the standard one dimensional unsteady constant-coefficient convection-diffusion equation. In a stability study, it is shown that the principal and spurious amplification factors of the current scheme, respectively, are structurally similar to those of the leapfrog/DuFort-Frankel scheme. As a result, the current scheme has no numerical diffusion in the special case of pure convection and is unconditionally stable in the special case of pure diffusion. Assuming smooth initial data, it will be shown theoretically and numerically that, by using an easily determined optimal time step, the accuracy of the current scheme may reach a level which is several orders of magnitude higher than that of the MacCormack scheme, with virtually identical operation count.

Chang, Sin-Chung↗
DOE OSTI2022

Tri-Sectional Approximation of the Shortest Path to Long-Term Voltage Stability Boundary with Distributed Energy Resources

Ensuring long-term voltage stability is critical for reliable operations of power grids. High share of distributed energy resources (DERs) can create complicated system operation modes that may invalidate the traditional long-term voltage stability analysis based on typical operation modes. To address this challenge, this paper investigates how to compute the shortest path to the voltage stability boundary in the DER aggregated load space with large dispersion. Instead of working in the Euclidean space, we establish the analysis and computations on the algebraic power flow manifold to better capture the curvature change of the shortest path along the direction of losing stability. A modified optimal control framework is presented for obtaining the ground-truth of the smooth shortest path on the manifold. To efficiently and accurately solve for the shortest path, we further leverage the geometric features of the power flow manifold and propose a tri-sectional approximation model that is scalable for large-scale systems. Several numerical examples, up to the 1354-bus system, with different DER penetration levels and high dimensional renewable power injection variations are evaluated. The simulation results demonstrate that the tri-sectional approximation achieves high accuracy and efficiency to approximate the shortest path to the voltage stability boundary.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Physics constrained learning for data-driven inverse modeling from sparse observations

Deep neural networks (DNN) have been used to model nonlinear relations between physical quantities. Those DNNs are embedded in physical systems described by partial differential equations (PDE) and trained by minimizing a loss function that measures the discrepancy between predictions and observations in some chosen norm. This loss function often includes the PDE constraints as a penalty term when only sparse observations are available. As a result, the PDE is only satisfied approximately by the solution. However, the penalty term typically slows down the convergence of the optimizer for stiff problems. We present a new approach that trains the embedded DNNs while numerically satisfying the PDE constraints. We develop an algorithm that enables differentiating both explicit and implicit numerical solvers in reverse-mode automatic differentiation. This allows the gradients of the DNNs and the PDE solvers to be computed in a unified framework. We demonstrate that our approach enjoys faster convergence and better stability in relatively stiff problems compared to the penalty method. Furthermore, our approach allows for the potential to solve and accelerate a wide range of data-driven inverse modeling, where the physical constraints are described by PDEs and need to be satisfied accurately.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Long-term Dynamical Stability in the Outer Solar System. I. The Regular and Chaotic Evolution of the 34 Largest Trans-Neptunian Objects

We carried out an extensive analysis of the stability of the outer solar system, making use of the frequency analysis technique over short-term integrations of nearly 100,000 test particles, as well as a statistical analysis of 200 1 Gyr long numerical simulations, which consider the mutual perturbations of the giant planets and the 34 largest trans-Neptunian objects (we have called all 34 objects “dwarf planets,” DPs, even if probably only the largest of them are true DPs). From the frequency analysis, we produced statistical diffusion maps for a wide region of the a–e phase-space plane; we also present the average diffusion time for orbits as a function of perihelion. We later turned our attention to the 34 DPs, making an individualized analysis for each of them and producing a first approximation of their future stability. From the 200 distinct realizations of the orbital evolution of the 34 DPs, we classified the sample into three categories, including 17 stable, 11 unstable, and 6 resonant objects; we also found that, statistically, two objects from the sample will leave the trans-Neptunian region within the next gigayear, most likely being ejected from the solar system, but with a nonnegligible probability of going inside the orbit of Neptune, either to collide with a giant planet or even falling to the inner solar system, where our simulations are no longer able to resolve their continuous evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1975

A numerical method for the prediction of high-speed boundary-layer transition using linear theory

A method is described of estimating the location of transition in an arbitrary laminar boundary layer on the basis of linear stability theory. After an examination of experimental evidence for the relation between linear stability theory and transition, a discussion is given of the three essential elements of a transition calculation: (1) the interaction of the external disturbances with the boundary layer; (2) the growth of the disturbances in the boundary layer; and (3) a transition criterion. The computer program which carried out these three calculations is described. The program is first tested by calculating the effect of free-stream turbulence on the transition of the Blasius boundary layer, and is then applied to the problem of transition in a supersonic wind tunnel. The effects of unit Reynolds number and Mach number on the transition of an insulated flat-plate boundary layer are calculated on the basis of experimental data on the intensity and spectrum of free-stream disturbances. Reasonable agreement with experiment is obtained in the Mach number range from 2 to 4.5.

Mack, L. M.↗
DOE OSTI2025

Modeling gas, hydrates, and slope stability on the U.S. Atlantic margin during Pleistocene glacial cycles

Changes in temperature and sea level can cause dissociation of methane hydrates in shallow marine sediments, leading to seafloor destabilization. Along the U.S. Atlantic margin, there exists a well-documented history of slope failure and numerous recorded occurrences of gas seeps. Several studies have linked slope failure in the region to gas seepage and hydrate dissociation driven by glacial-interglacial transitions, but this linkage has not been quantitatively demonstrated. Along the shelf edge, in an area where shallow methane gas seeps have been identified, we modeled methane gas and hydrate formation using a one-dimensional fluid flow model. Methane gas formation was modeled over the last 120,000 years to simulate a glacial-interglacial cycle. Here, we ran this model at 16,044 individual locations in the region between 29° N – 45° N and 82° W – 66° W at a resolution of 1 x 1 arcminutes, focusing specifically on water depths between 200 and 1000 m that bracket the seafloor outcrop of the base of the hydrate stability zone. Using historic temperature and pressure records from the last 120,000 years, sediment properties in the area, and factor of safety calculations, we found that hydrate dissociation alone is unlikely to cause sediment failure in the region, implying that an additional driving force would be necessary for failure to occur.

Hydrate dissociation↗
DOE OSTI2023

Extended isogeometric analysis of multi-material and multi-physics problems using hierarchical B-splines

Here, this paper presents an immersed, isogeometric finite element framework to predict the response of multi-material, multi-physics problems with complex geometries using locally refined discretizations. To circumvent the need to generate conformal meshes, this work uses an extended finite element method (XFEM) to discretize the governing equations on non-conforming, embedding meshes. A flexible approach to create truncated hierarchical B-splines discretizations is presented. This approach enables the refinement of each state variable field individually to meet field-specific accuracy requirements. To obtain an immersed geometry representation that is consistent across all hierarchically refined B-spline discretizations, the geometry is immersed into a single mesh, the XFEM background mesh, which is constructed from the union of all hierarchical B-spline meshes. An extraction operator is introduced to represent the truncated hierarchical B-spline bases in terms of Lagrange shape functions on the XFEM background mesh without loss of accuracy. The truncated hierarchical B-spline bases are enriched using a generalized Heaviside enrichment strategy to accommodate small geometric features and multi-material problems. The governing equations are augmented by a formulation of the face-oriented ghost stabilization enhanced for locally refined B-spline bases. We present examples for two- and three-dimensional linear elastic and thermo-elastic problems. The numerical results validate the accuracy of our framework. The results also demonstrate the applicability of the proposed framework to large, geometrically complex problems.

42 ENGINEERING↗
DOE OSTI2022

Do ambient shear and thermal stratification impact wind turbine tip-vortex breakdown?

Modern wind turbines experience uneven inflow conditions across the rotor, due to the ambient flow’s shear and thermal stratification. Such conditions alter the shape and length of turbine wakes and thus impact the loads and power generation of downstream turbines. To this end, understanding the spatial evolution of the individual wakes under different atmospheric conditions is key to controlling and optimising turbine arrays. With this numerical study we aim to obtain a better understanding of the fundamental physics governing the near-wake dynamics of wind turbines under shear and thermal stability, by examining their tip-vortex breakup mechanisms. Our approach considers scale-resolving simulations of a single turbine wake under a linear shear profile as well as the application of harmonic tip perturbations to trigger flow instabilities. For the subsequent analysis we use the proper orthogonal decomposition (POD) method to extract coherent structures from the flow, and we also calculate mean kinetic energy fluxes to quantify each coherent structure’s contribution to wake recovery. The wake’s helical spiral is found to hinder wake recovery for all studied ambient flow conditions, whereas the mutual inductance instability has positive MKE flux leading to an enhanced wake recovery. Finally, the ambient shear has the largest impact on the local MKE flux with respect to downstream location by changing the shape of the curve and location of extrema, whereas thermal stratification has only a minimal impact on the magnitude of the near-wake local MKE flux distribution.

17 WIND ENERGY↗
DOE OSTI2024

Impact of hydrogenation on the stability and mechanical properties of amorphous boron nitride

Abstract Interconnect materials with ultralow dielectric constant, and good thermal and mechanical properties are crucial for the further miniaturization of electronic devices. Recently, it has been demonstrated that ultrathin amorphous boron nitride (aBN) films have a very low dielectric constant, high density (above 2.1 g cm −3 ), high thermal stability, and mechanical properties. The excellent properties of aBN derive from the nature and degree of disorder, which can be controlled at fabrication, allowing tuning of the physical properties for desired applications. Here, we report an improvement in the stability and mechanical properties of aBN upon hydrogen doping. With the introduction of a Gaussian approximation potential for atomistic simulations, we investigate the changing morphology of aBN with varying H doping concentrations. We found that for 8 at% of H doping, the concentration ofsp 3 -hybridized atoms reaches to a maximum which leads to an improvement of thermal stability and mechanical properties by 20%. These results will be a guideline for experimentalists and process engineers to tune the growth conditions of aBN films for numerous applications.

Materials Science↗
NASA NTRS1974

Influence of rotation on the maximum mass of pulsationally stable stars

The critical mass for stability against radial pulsations in rotating, homogeneous main-sequence stars is found to be greater than in the case of no rotation. Analytic and detailed numerical models show that the critical mass rises steeply with increasing concentration of angular momentum to the center of the star. For uniform rotation near breakup velocity at the star's equator the critical mass is about 850 solar masses if an electron-scattering opacity is used, or about 5000 solar masses if the opacities of Cox and Stewart are used. For nonuniform rotation with a constant ratio of centrifugal force to gravity in the star, the critical mass becomes 'infinite' long before breakup velocity is attained. The relevance of the present results to several observational problems is noted.

Stothers, R.↗
NASA NTRS1979

A note on a corrector formula for the numerical solution of ordinary differential equations

A new corrector formula for predictor-corrector methods for numerical solutions of ordinary differential equations is presented. Two considerations for choosing corrector formulas are given: (1) the coefficient in the error term and (2) its stability properties. The graph of the roots of an equation plotted against its stability region, of different values, is presented along with the tables that correspond to various corrector equations, including Hamming's and Milne and Reynolds'.

Chien, Y.-C.↗