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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

NASA NTRS1986

On compactness of admissible parameter sets: Convergence and stability in inverse problems for distributed parameter systems

A series of numerical examples is reported and several algorithms compared for estimation of coefficients in differential equation models. Unconstrained, constrained and Tikhonov regularization methods are tested for their behavior with regard to both convergence (of approximation methods for the states and parameters) and stability (continuity of the estimates with respect to perturbations in the data or observed states).

Banks, H. T.↗

DOE OSTI2023

ML-Based Power System Stability Assessment Considering Network Topology Changes: WECC 20,000+ Bus System Case Study

Modern power grids are fast-changing and thus require real-time monitoring and online stability assessment. With the rapid development of machine learning (ML) techniques, using data-driven models to provide fast and accurate estimations of power system stability marginal information, such as frequency nadir for frequency stability and critical clearing time (CCT) for transient stability, have become possible. However, despite the numerous research on ML-based methods for frequency nadir and CCT prediction, there is limited work on the impact of different network topology changes. Furthermore, most previous studies only focused on small or synthetic systems, and there is a lack of research on actual large power system models. In this paper, the above issues are addressed by studying the actual U.S. Western Electricity Coordinating Council (WECC) system model with more than 20,000 buses. Massive simulations are conducted in PowerWorld Simulator to study the impact of various topology change scenarios on both frequency stability and transient stability. System operating information is extracted from the success dispatch cases of various network topologies to generate a comprehensive dataset for ML-based models. Two ML methods, random forest (RF) and multilayer perceptron (MLP) neural network, are trained and tested for both frequency nadir prediction and CCT prediction. Test results have proven the models are capable of online stability assessment for large power networks such as the WECC system with sufficient accuracy.

critical clearing time↗

NASA NTRS1974

Linear viscous stability analysis of the stratified Bickley jet

The stability of the stably stratified Bickley jet to small disturbances in a viscous, heat conducting fluid is explored by numerical integration of the governing differential equation. Curves of neutral stability are presented in the alpha-R plane and the parametric dependence of these curves upon Richardson number and Prandtl number is examined. The results show a curious bimodal behavior of the curve of neutral stability for small Richardson numbers. In addition, closed neutral curves are found for Richardson numbers in excess of the critical value for the inviscid flow. The unstable regions in the interior of these curves increase with decreasing Prandtl number. Energy balances are computed which show the relative magnitude of the Reynolds stress production, buoyancy dissipation, and viscous dissipation terms across the flow. As the thermal diffusivity becomes large compared with the kinematic viscosity, the buoyancy term is suppressed and the thermal stratification becomes irrelevant.

Gage, K. S.↗

DOE OSTI2022

Embedded pairs for optimal explicit strong stability preserving Runge–Kutta methods

We construct a family of embedded pairs for optimal explicit strong stability preserving Runge–Kutta methods of order 2 ≤ p ≤ 4 to be used to obtain numerical solution of spatially discretized hyperbolic PDEs. In this construction, the goals include non-defective property, large stability region, and small error values as defined in Dekker and Verwer (1984) and Kennedy et al. (2000). The new family of embedded pairs offer the ability for strong stability preserving (SSP) methods to adapt by varying the step-size. Through several numerical experiments, we assess the overall effectiveness in terms of work versus precision while also taking into consideration accuracy and stability.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Numerical analyses of exponential time-differencing schemes for the solution of atmospheric models

In high-resolution numerical weather prediction models, fast-moving acoustic waves must be treated in a stable manner. Among the implicit–explicit (IMEX) schemes used for the solution of these models, horizontally explicit, vertically implicit (HEVI) methods show good stability and scalability on massively parallel machines. In this work, we present two classes of exponential time-differencing (ETD) methods for atmospheric models that use a HEVI splitting strategy, one being a three-stage method with the addition of artificial diffusion, the second based on a Strang splitting approach. Overall, the stability properties of the methods are analyzed and numerical examples are provided, which compare time-step restrictions and cost-to-accuracy ratios of the new methods with those for existing approaches.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1991

Resonant Triad in Boundary-Layer Stability: Composite Solution and Comparison with Observations - Part 2

Here, numerical results are computed from an asymptotic near-resonance triad analysis. The analysis considers a resonant triad of instability waves consisting of a plane fundamental wave and a pair of symmetrical oblique subharmonic waves. The relevant scaling ensures that nonlinearity is confined to a distinct critical layer. The analysis is first used to form a composite solution that accounts for both the flow divergence and nonlinear effects. It is shown that the backreaction on the plane Tollmien Schlichting (TS) fundamental wave, although fully accounted for, is of little significance. The observed enhancement at the fundamental frequency disturbance is not in the plane TS wave, but is caused by nonlinearly generated waves at the fundamental frequency that result from nonlinear interactions in the critical layer. The saturation of the oblique waves is caused by their self-interaction. The nonlinear phase-locking phenomenon, the location of resonance with respect to the neutral stability curve, low frequency effects, detuning in the streamwise wave numbers, and nonlinear distortion of the mode shapes are discussed. Nonlinearity modifies the initially two dimensional Blasius profile into a fuller one with spanwise periodicity. The interactions at a wide range of unstable spanwise wave numbers are considered, and the existence of a preferred spanwise wave number is explained by means of the vorticity distribution in the critical layer. Besides presenting novel features of the phenomena and explaining the delicate mechanisms of the interactions, the results of the theory are in excellent agreement with experimental and numerical observations for all stages of the development and for various input parameters.

Mankbadi, Reda R.↗

NASA NTRS1991

Numerical simulation of spatially-evolving instability in plane channel flow

The spatial stability of plane channel flow is analyzed using a three-dimensional, time-dependent spectral/finite difference code (Danabasoglu et al., 1990) which integrates numerically the Navier-Stokes equations. The study centers on inflow disturbance amplitudes effects on the secondary instability. The resolution requirements along the spacewise direction, which become critical before the breakdown stage, are of particular interest. A direct comparison is made with the experiments of Nishioka et al. (1980).

Danabasoglu, G.↗

NASA NTRS1977

A charging model for three-axis stabilized spacecraft

A charging model was developed for geosynchronous, three-axis stabilized spacecraft when under the influence of a geomagnetic substorm. The differential charging potentials between the thermally coated or blanketed outer surfaces and metallic structure of a spacecraft were determined when the spacecraft was immersed in a dense plasma cloud of energetic particles. The spacecraft-to-environment interaction was determined by representing the charged particle environment by equivalent current source forcing functions and by representing the spacecraft by its electrically equivalent circuit with respect to the plasma charging phenomenon. The charging model included a sun/earth/spacecraft orbit model that simulated the sum illumination conditions of the spacecraft outer surfaces throughout the orbital flight on a diurnal as well as a seasonal basis. Transient and steady-state numerical results for a three-axis stabilized spacecraft are presented.

Massaro, M. J.↗

NASA NTRS1990

Use of satellite data and modeling to assess the influence of stratospheric processes on the troposphere

Current research focused on the examination of the effects of ozone heating and Newtonian cooling on the linear stability of transient planetary waves. A finite difference numerical model was developed in which the basic state wind, temperature, and ozone distributions were constructed using ground-based and satellite data. Preliminary results indicate that ozone dynamics interaction may play a more important role in the stability and maintenance of planetary waves than previously thought, particularly in summer when the mean solar zenith angle is smallest and thus ozone heating largest. Examination of the finite amplitude interactions among radiation, ozone, and dynamics was carried out. Self-consistent, coupled equations governing the weakly nonlinear interactions are currently being analyzed to provide a better understanding of wave-mean flow interactions and ozone transport in a continuously stratified model of the troposphere-stratosphere coupled system.

Nathan, Terrence↗

NASA NTRS1993

Direct numerical simulation of hot jets

The ultimate motivation of this work is to investigate the stability of two dimensional heated jets and its implications for aerodynamic sound generation from data obtained with direct numerical simulations (DNS). As pointed out in our last report, these flows undergo two types of instabilities, convective or absolute, depending on their temperature. We also described the limits of earlier experimental and theoretical studies and explained why a numerical investigation could give us new insight into the physics of these instabilities. The aeroacoustical interest of these flows was also underlined. In order to reach this goal, we first need to succeed in the DNS of heated jets. Our past efforts have been focused on this issue which encountered several difficulties. Our numerical difficulties are directly related to the physical problem we want to investigate since these absolutely or almost absolutely unstable flows are by definition very sensitive to the smallest disturbances and are very likely to reach nonlinear saturation through a numerical feedback mechanism. As a result, it is very difficult to compute a steady laminar solution using a spatial DNS. A steady state was reached only for strongly co-flowed jets, but these flows are almost equivalent to two independent mixing layers. Thus they are far from absolute instability and have much lower growth rates.

Jacob, Marc C.↗

NASA NTRS2011

Effects of Nose Bluntness on Hypersonic Boundary-Layer Receptivity and Stability Over Cones

The receptivity to freestream acoustic disturbances and the stability properties of hypersonic boundary layers are numerically investigated for boundary-layer flows over a 5 straight cone at a freestream Mach number of 6.0. To compute the shock and the interaction of the shock with the instability waves, the Navier-Stokes equations in axisymmetric coordinates were solved. In the governing equations, inviscid and viscous flux vectors are discretized using a fifth-order accurate weighted-essentially-non-oscillatory scheme. A third-order accurate total-variation-diminishing Runge-Kutta scheme is employed for time integration. After the mean flow field is computed, disturbances are introduced at the upstream end of the computational domain. The appearance of instability waves near the nose region and the receptivity of the boundary layer with respect to slow mode acoustic waves are investigated. Computations confirm the stabilizing effect of nose bluntness and the role of the entropy layer in the delay of boundary-layer transition. The current solutions, compared with experimental observations and other computational results, exhibit good agreement.

Kara, Kursat↗

NASA NTRS2014

Entropy Stable Wall Boundary Conditions for the Compressible Navier-Stokes Equations

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite volume, finite difference, discontinuous Galerkin, and flux reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Parsani, Matteo↗

NASA NTRS1934

Supplemental data and calculations of the lateral stability of airplanes

In connection with the DVL Report 272 on the theory of the lateral stability of airplanes, the formal results are here amplified in some respects and their technical significance again briefly explained. Three numerical examples show how model tests for checking the lateral stability are to be evaluated and supplemented, if necessary, and how the stability limits depend on the design of the airplane and on the conditions of flight.

Mathias, Gotthold↗

NASA NTRS2010

Harmonic Balance Computations of Fan Aeroelastic Stability

A harmonic balance (HB) aeroelastic analysis, which has been recently developed, was used to determine the aeroelastic stability (flutter) characteristics of an experimental fan. To assess the numerical accuracy of this HB aeroelastic analysis, a time-domain aeroelastic analysis was also used to determine the aeroelastic stability characteristics of the same fan. Both of these three-dimensional analysis codes model the unsteady flowfield due to blade vibrations using the Reynolds-averaged Navier-Stokes (RANS) equations. In the HB analysis, the unsteady flow equations are converted to a HB form and solved using a pseudo-time marching method. In the time-domain analysis, the unsteady flow equations are solved using an implicit time-marching approach. Steady and unsteady computations for two vibration modes were carried out at two rotational speeds: 100 percent (design) and 70 percent (part-speed). The steady and unsteady results obtained from the two analysis methods compare well, thus verifying the recently developed HB aeroelastic analysis. Based on the results, the experimental fan was found to have no aeroelastic instability (flutter) at the conditions examined in this study.

Bakhle, Milind A.↗

NASA NTRS1990

Efficient solution of parabolic equations by Krylov approximation methods

Numerical techniques for solving parabolic equations by the method of lines is addressed. The main motivation for the proposed approach is the possibility of exploiting a high degree of parallelism in a simple manner. The basic idea of the method is to approximate the action of the evolution operator on a given state vector by means of a projection process onto a Krylov subspace. Thus, the resulting approximation consists of applying an evolution operator of a very small dimension to a known vector which is, in turn, computed accurately by exploiting well-known rational approximations to the exponential. Because the rational approximation is only applied to a small matrix, the only operations required with the original large matrix are matrix-by-vector multiplications, and as a result the algorithm can easily be parallelized and vectorized. Some relevant approximation and stability issues are discussed. We present some numerical experiments with the method and compare its performance with a few explicit and implicit algorithms.

Gallopoulos, E.↗

NASA NTRS2015

Entropy Stable Wall Boundary Conditions for the Three-Dimensional Compressible Navier-Stokes Equations

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators on unstructured grids are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite difference, finite volume, discontinuous Galerkin, and flux reconstruction/correction procedure via reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Parsani, Matteo↗

NASA NTRS2000

The Generation of Smooth High Speed Solar Wind From Plume-Interplume Mixing

Ulysses has shown that fast solar wind is extremely smooth, with a variance of less than 5%, in contrast to slow wind with a variance of approximately 30%. Now UVCS has produced the surprising result that the flow speed within coronal holes, the source of fast wind, is not at all smooth. Specifically, Giordano et al. (ApJ, v531, L79-L82, 2000) report that at 1.7 R(sub SUN) the interplume flow speed is typically more than twice the plume flow speed. Other less direct evidence supports this same result, with speeds from less than 300 to over 1000 km/s reported at approximately 5 R(sub SUN). This presents the paradox of how strongly differing plume and interplume flow speeds can exist near the Sun and be absent far from the Sun. The only answer is that plume and interplume material or momentum must be strongly mixed and that the mixing must occur mainly inside 0.3 AU to be consistent with Helios observations of smooth fast wind. Pressure balance structures (PBSs) and He abundance anomalies (Reisenfeld, et al., GRL, v26(13), 1805-1808, 1999) have been identified as interplanetary remnants of plumes, implying momentum mixing is the dominant coronal process. One possible source for plume/interplume momentum mixing is MHD Kelvin-Helmholtz (KH) shear instabilities occurring on the velocity shear interfaces. The velocity shear is a source of free energy and KH fluctuations could, through nonlinear cascade, provide the forcing required for the plasma oscillations (Cranmer, ApJ, v532, 1197-1208, 2000) reported to exist in coronal holes. The physical properties in coronal holes are now sufficiently well known that we can show plume/interplume shear interfaces become unstable to the KH instability at 5 - 10 R(sub SUN). The KH dispersion relation can be used to analyze marginal stability, the most unstable wavelengths, and linear growth rates. Numerical simulations can be used to verify results from the linear analysis and study the nonlinear development of KH modes. Here we will describe the marginal stability criterion for the KH instability, how and where this condition is met in the corona, and the possible character of the resultant fluctuations.

Suess, Steve↗

NASA NTRS1986

Entropy conservative finite element schemes

The question of entropy stability for discrete approximations to hyperbolic systems of conservation laws is studied. The amount of numerical viscosity present in such schemes is quantified and related to their entropy stability by means of comparison. To this end, two main ingredients are used: entropy variables and the construction of certain entropy conservative schemes in terms of piecewise-linear finite element approximations. It is then shown that conservative schemes are entropy stable, if and (for three-point schemes) only if, they contain more numerical viscosity than the abovementioned entropy conservation ones.

Tadmor, E.↗