Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical stabilization”

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,099 records · Page 61

Thermal and capillary effects on the caprock mechanical stability at In Salah, Algeria

Thermo-mechanical effects are important in geologic carbon storage because CO 2 will generally reach the storage formation colder than the rock, inducing thermal stresses. Capillary functions, i.e., retention and relative permeability curves, control the CO 2 plume shape, which may affect overpressure and thus, caprock stability. To analyze these thermal and capillary effects, we numerically solve non-isothermal injection of CO 2 in deformable porous media considering the In Salah, Algeria, CO 2 storage site. Here, we find that changes in the capillary functions have a negligible effect on overpressure and thus, caprock stability is not affected by capillary effects. But, we show that for the strike slip stress regime prevalent at In Salah, stability decreases in the lowest parts of the caprock during injection due to cooling-induced thermal stresses. Simulations show that shear slip along pre-existing fractures may take place in the cooled region, whereas tensile failure is less likely to occur. Indeed, only the injection zone and the lowest tens of meters of the 900-m-thick caprock at In Salah might be affected by cooling effects, which would thus not jeopardize the overall sealing capacity of the caprock. Furthermore, faults are likely to remain stable far away from the injection well because outside the cooled region the injection-induced stress changes are not sufficient to exceed the anticipated shear strength of minor faults. Nonetheless, we recommend that thermal effects should be considered in the site characterization and injection design of future CO 2 injection sites to assess caprock stability and guarantee a permanent CO 2 storage.

54 ENVIRONMENTAL SCIENCES↗

Finite element or Galerkin type semidiscrete schemes

A finite element of Galerkin type semidiscrete method is proposed for numerical solution of a linear hyperbolic partial differential equation. The question of stability is reduced to the stability of a system of ordinary differential equations for which Dahlquist theory applied. Results of separating the part of numerical solution which causes the spurious oscillation near shock-like response of semidiscrete scheme to a step function initial condition are presented. In general all methods produce such oscillatory overshoots on either side of shocks. This overshoot pathology, which displays a behavior similar to Gibb's phenomena of Fourier series, is explained on the basis of dispersion of separated Fourier components which relies on linearized theory to be satisfactory. Expository results represented.

Durgun, K.↗

Numerical advection by conservation of second-order moments

A new, accurate, and nondiffusive method for three-dimensional advection of trace species is presented. The method preserves tracer structures by conserving the second-order moments of the spatial distribution of tracer during advection. Upstream transport and second-order tracer distribution are described, and the moments of the tracer distribution about the center of a grid box are formally defined and related to the polynomial distribution. Formulas are presented which describe how the moments of a grid box are decomposed into a unique set of moments centered about each subbox and how they are reassembled into new grid boxes. A one-dimensional example if tracer transport is given, and limits necessary to maintain positive tracer concentrations are derived. The accuracy and stability of this method are analytically examined, and numerical experiments testing the effective resolution are presented. The method is compared with other methods for numerical advection of tracers.

Prather, M. J.↗

Response studies of rotors and rotor blades with application to aeroelastic tailoring

Various tools for the aeroelastic stability and response analysis of rotor blades in hover and forward flight were developed and incorporated in a comprehensive package capable of performing aeroelastic tailoring of rotor blades in forward flight. The results indicate that substantial vibration reductions, of order 15-40%, in the vibratory hub shears can be achieved by relatively small modifications of the initial design. Furthermore the optimized blade can be up to 20% lighter than the original design. Accomplishments are reported for the following tasks: (1) finite element modeling of rotary-wing aeroelastic problems in hover and forward flight; (2) development of numerical methods for calculating the aeroelastic response and stability of rotor blades in forward fight; (3) formulation of the helicopter air resonance problem in hover with active controls; and (4) optimum design of rotor blades for vibration reduction in forward flight.

Friedmann, P. P.↗

Optimal subinterval selection approach for power system transient stability simulation

Power system transient stability analysis requires an appropriate integration time step to avoid numerical instability as well as to reduce computational demands. For fast system dynamics, which vary more rapidly than what the time step covers, a fraction of the time step, called a subinterval, is used. However, the optimal value of this subinterval is not easily determined because the analysis of the system dynamics might be required. This selection is usually made from engineering experiences, and perhaps trial and error. This paper proposes an optimal subinterval selection approach for power system transient stability analysis, which is based on modal analysis using a single machine infinite bus (SMIB) system. Fast system dynamics are identified with the modal analysis and the SMIB system is used focusing on fast local modes. An appropriate subinterval time step from the proposed approach can reduce computational burden and achieve accurate simulation responses as well. As a result, the performance of the proposed method is demonstrated with the GSO 37-bus system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Stability, accuracy, and efficiency of some underintegrated methods in finite element computations

In an attempt to increase computational efficiency in the numerical solution of highly nonlinear problems in solid and fluid mechanics, underintegrated finite element methods have been employed by many analysts. Underintegration refers to the use of a rule of an order lower than that required to integrate polynomial integrands exactly. The main drawback of this technique is related to the production of rank-deficient stiffness matrices, or equivalently an expanded kernel of the governing linear momentum operators. Such a development can introduce numerical instabilities. In order to overcome this difficulty, artificial stiffness or viscosity methods, or other stabilization methods have been proposed. One approach involves the elimination of spurious modes in a postprocessing operation. The present study is concerned with this a posteriori elimination method, taking into account the results which can be expected from it, and some of its possible extensions.

Jacquotte, O.-P.↗

Numerical Continuation Methods for Intrusive Uncertainty Quantification Studies

Rigorous modeling of engineering systems relies on efficient propagation of uncertainty from input parameters to model outputs. In recent years, there has been substantial development of probabilistic polynomial chaos (PC) Uncertainty Quantification (UQ) methods, enabling studies in expensive computational models. One approach, termed ”intrusive”, involving reformulation of the governing equations, has been found to have superior computational performance compared to non-intrusive sampling-based methods in relevant large-scale problems, particularly in the context of emerging architectures. However, the utility of intrusive methods has been severely limited due to detrimental numerical instabilities associated with strong nonlinear physics. Previous methods for stabilizing these constructions tend to add unacceptably high computational costs, particularly in problems with many uncertain parameters. In order to address these challenges, we propose to adapt and improve numerical continuation methods for the robust time integration of intrusive PC system dynamics. We propose adaptive methods, starting with a small uncertainty for which the model has stable behavior and gradually moving to larger uncertainty where the instabilities are rampant, in a manner that provides a suitable solution.

97 MATHEMATICS AND COMPUTING↗

Stability and instability of thermocapillary convection in models of the float-zone crystal-growth process

This project was concerned with the determination of conditions of guaranteed stability and instability for thermocapillary convection in a model of the float-zone crystal-growth process. This model, referred to as the half-zone, was studied extensively, both experimentally and theoretically. Our own earlier research determined, using energy-stability theory, sufficient conditions for stability to axisymmetric disturbances. Nearly all results computed were for the case of a liquid with Prandtl Number Pr = 1. Attempts to compute cases for higher Prandtl numbers to allow comparison with the experimental results of other researchers were unsuccessful, but indicated that the condition guaranteeing stability against axisymmetric disturbances would be a value of the Marangoni number (Ma), significantly higher than that at which oscillatory convection was observed experimentally. Thus, additional results were needed to round out the stability picture for this model problem. The research performed under this grant consisted of the following: (1) computation of energy-stability limits for non-axisymmetric disturbances; (2) computation of linear-stability limits for axisymmetric and non-axisymmetric disturbances; (3) numerical simulation of the basic state for half- and full-zones with a deformable free surface; and (4) incorporation of radiation heat transfer into a model energy-stability problem. Each of these is summarized briefly below.

Neitzel, G. P.↗

Modeling of Macroscopic Dynamics in Three-Dimensional Plasma Configurations (Final Technical Report)

A numerical simulation-based study to investigate the stability properties of inherently three-dimensional plasma magnetic configurations is presented. These configurations, called stellarators and heliotrons, use external helical coils to provide the twisting magnetic field that is required to confine ionized plasma particles. In contrast to the tokamak configuration, which nominally has the property of axisymmetry, electrical current does not need to be induced in stellarator and heliotron plasma. Experiments have shown that these three-dimensional configurations are remarkably robust to macroscopic instability, even when driven to relatively high levels of plasma pressure. Our understanding of how the plasma is able to withstand macroscopic instability will be enhanced by the work proposed here. The effort will generalize existing numerical modeling capability, and it will apply the new capability to model macroscopic dynamics in these three-dimensional configurations. The development work has been applied to the NIMROD simulation code (https://nimrodteam.org) that was originally based on axisymmetric problem domains. It has been generalized to accommodate the three-dimensional shaping of the plasma volume in stellarators and heliotrons. Its use of equilibrium data has also been generalized, and the implementation allows different numerical representation of the magnetic field. Analysis has been performed to inform the selection of numerical methods, and a method to accelerate the iterative solution of the necessary algebraic systems has been developed and implemented. Verification of the new NIMSTELL code on benchmark problems has also been an important part of the effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations: Global Asymptotic Behavior of Time Discretizations - 2

The global asymptotic nonlinear behavior of 1 1 explicit and implicit time discretizations for four 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODES) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how 'numerical' basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DEs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless of the stability of the spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDES.

Yee, H. C.↗

Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations: Global Asymptotic Behavior of Time Discretizations - Part 2

The global asymptotic nonlinear behavior of 11 explicit and implicit time discretizations for four 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how 'numerical' basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DEs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless of the stability of the spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDEs.

Yee, H. C.↗

PMU-Aided Voltage Security Assessment for a Wind Power Plant

Because wind power penetration levels in electric power systems are continuously increasing, voltage stability is a critical issue for maintaining power system security and operation. The traditional methods to analyze voltage stability can be classified into two categories: dynamic and steady-state. Dynamic analysis relies on time-domain simulations of faults at different locations; however, this method needs to exhaust faults at all locations to find the security region for voltage at a single bus. With the widely located phasor measurement units (PMUs), the Thevenin equivalent matrix can be calculated by the voltage and current information collected by the PMUs. This paper proposes a method based on a Thevenin equivalent matrix to identify system locations that will have the greatest impact on the voltage at the wind power plant's point of interconnection. The number of dynamic voltage stability analysis runs is greatly reduced by using the proposed method. The numerical results demonstrate the feasibility, effectiveness, and robustness of the proposed approach for voltage security assessment for a wind power plant.

17 WIND ENERGY↗

Analysis of ELM stability with extended MHD models in JET, JT-60U and future JT-60SA tokamak plasmas

The stability with respect to a peeling–ballooning mode (PBM) was investigated numerically with extended MHD simulation codes in JET, JT-60U and future JT-60SA plasmas. The MINERVA-DI code was used to analyze the linear stability, including the effects of rotation and ion diamagnetic drift (${\omega }_{* {\rm{i}}}$), in JET-ILW and JT-60SA plasmas, and the JOREK code was used to simulate nonlinear dynamics with rotation, viscosity and resistivity in JT-60U plasmas. It was validated quantitatively that the ELM trigger condition in JET-ILW plasmas can be reasonably explained by taking into account both the rotation and ${\omega }_{* {\rm{i}}}$ effects in the numerical analysis. When deuterium poloidal rotation is evaluated based on neoclassical theory, an increase in the effective charge of plasma destabilizes the PBM because of an acceleration of rotation and a decrease in ${\omega }_{* {\rm{i}}}$. The difference in the amount of ELM energy loss in JT-60U plasmas rotating in opposite directions was reproduced qualitatively with JOREK. By comparing the ELM affected areas with linear eigenfunctions, it was confirmed that the difference in the linear stability property, due not to the rotation direction but to the plasma density profile, is thought to be responsible for changing the ELM energy loss just after the ELM crash. A predictive study to determine the pedestal profiles in JT-60SA was performed by updating the EPED1 model to include the rotation and ${\omega }_{* {\rm{i}}}$ effects in the PBM stability analysis. It was shown that the plasma rotation predicted with the neoclassical toroidal viscosity degrades the pedestal performance by about 10% by destabilizing the PBM, but the pressure pedestal height will be high enough to achieve the target parameters required for the ITER-like shape inductive scenario in JT-60SA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Current driven instabilities of an electromagnetically accelerated plasma

A plasma instability that strongly influences the efficiency and lifetime of electromagnetic plasma accelerators was quantitatively measured. Experimental measurements of dispersion relations (wave phase velocities), spatial growth rates, and stability boundaries are reported. The measured critical wave parameters are in excellent agreement with theoretical instability boundary predictions. The instability is current driven and affects a wide spectrum of longitudinal (electrostatic) oscillations. Current driven instabilities, which are intrinsic to the high-current-carrying magnetized plasma of the magnetoplasmadynmic (MPD) accelerator, were investigated with a kinetic theoretical model based on first principles. Analytical limits of the appropriate dispersion relation yield unstable ion acoustic waves for T(i)/T(e) much less than 1 and electron acoustic waves for T(i)/T(e) much greater than 1. The resulting set of nonlinear equations for the case of T(i)/T(e) = 1, of most interest to the MPD thruster Plasma Wave Experiment, was numerically solved to yield a multiparameter set of stability boundaries. Under certain conditions, marginally stable waves traveling almost perpendicular to the magnetic field would travel at a velocity equal to that of the electron current. Such waves were termed current waves. Unstable current waves near the upper stability boundary were observed experimentally and are in accordance with theoretical predictions. This provides unambiguous proof of the existence of such instabilites in electromagnetic plasma accelerators.

Chouetri, E. Y.↗

A dynamic analysis of the motion of a low-wing general aviation aircraft about its calculated equilibrium flat spin mode

The properties of the flat spin mode of a general aviation configuration have been studied through analysis of rotary balance data, numerical simulation, and analytical study of the equilibrium state. The equilibrium state is predicted well from rotary balance data. The variations of yawing moment and pitching moment as functions of sideslip have been shown to be of great importance in obtaining accurate modeling. These dependencies are not presently available with sufficient accuracy from previous tests or theories. The stability of the flat spin mode has been examined extensively using numerical linearization, classical perturbation methods, and reduced order modeling. The stability exhibited by the time histories and the eigenvalue analyses is shown to be strongly dependent on certain static cross derivatives and more so on the dynamic derivatives. Explicit stability criteria are obtained from the reduced order models.

Tischler, M. B.↗

A method for the early detection of instabilities in rotordynamics

For a simple Jeffcott model with deadband, a method is developed to determine stability margins by an analysis of the Discrete Fourier transform of the system response. The viability of this method is demonstrated by means of numerical simulations. The circular behavior of the system response on the stability boundary is also explained.

Zalik, Richard A.↗

Creep Resistance of ZrO2 Ceramic Improved by the Addition of a Small Amount of Er2O3

Zirconia (ZrO2) has great technological importance in structural, electrical, and chemical applications. It is the crucial component for state-of-the art thermal barrier coatings and an enabling component as a solid electrolyte for solid-oxide fuel cell systems. Pure ZrO2 is of limited use for industrial applications because of the phase transformations that occur. Upon the addition of stabilizers, cubic (c-ZrO2) and tetragonal (t-ZrO2) forms can be preserved. It is the stabilized and partially stabilized forms of zirconia that function as thermal barrier coatings, solid electrolytes, and oxygen sensors and that have numerous applications in the electrochemical industry. The cubic form of ZrO2 is typically stabilized through Y2O3 additions. However, Y2O3-stabilized zirconia is susceptible to deformation at high temperatures (greater than 900 C) because of the large number of slip systems and the high oxygen diffusion rates, which result in high creep rates at high temperatures. Successful use of ZrO2 at high temperatures requires that new dopant additives be found that will retain or enhance the desirable properties of cubic ZrO2 and yet produce a material with lower creep rates. At the NASA Glenn Research Center, erbium oxide (Er2O3) was identified as a promising dopant for improving the creep resistance of. ZrO2. The selection of Er2O3 was based on the strong interactions of point defects and dislocations. Single crystals of 5 mol% Er2O3- doped ZrO2 rods (4 mm in diameter) and monofilaments (200 to 300 mm in diameter and 30 cm long) were grown using the laser-heated float zone technique, and their creep behavior was measured as a function of temperature. The addition of 5 mol% Er2O3 to single-crystal ZrO2 improved its creep resistance at high temperatures by 2 to 3 orders of magnitude over state-of-the-art Y2O3-doped crystals. Detailed microstructural characterization of ZrO2-Er2O3 single crystals has identified new mechanisms for improving the creep resistance of this class of materials. Adding Er2O3 to ZrO2 results in microstructure of stable and metastable tetragonal precipitates that with thermal treatment evolve to a tweed structure of nanosize tetragonal lamellae. The superior high-temperature creep resistance of Er2O3-doped ZrO2 is attributed to nanoscale precipitation hardening. Doping with Er2O3 will significantly increase the upper-use temperature limit of ZrO2. Potential applications include using Er2O3-doped ZrO2 as a high-temperature fiber for structural applications and adding Er2O3 to reduce the sintering rates of ZrO2 thermal barrier coatings. This work was conducted at Dpto. de F sica de la Materia Condensada, Universidad de Sevilla, Spain, and at NASA Glenn.

Martinez-Fernandez, Julian↗

Dynamical Instabilities in Systems of Multiple Short-period Planets Are Likely Driven by Secular Chaos: A Case Study of Kepler-102

We investigated the dynamical stability of high-multiplicity Kepler and K2 planetary systems. Our numerical simulations find instabilities in ∼20% of the cases on a wide range of timescales (up to 5 × 10{sup 9} orbits) and over an unexpectedly wide range of initial dynamical spacings. To identify the triggers of long-term instability in multiplanet systems, we investigated in detail the five-planet Kepler-102 system. Despite having several near-resonant period ratios, we find that mean-motion resonances are unlikely to directly cause instability for plausible planet masses in this system. Instead, we find strong evidence that slow inward transfer of angular momentum deficit (AMD) via secular chaos excites the eccentricity of the innermost planet, Kepler-102 b, eventually leading to planet–planet collisions in ∼80% of Kepler-102 simulations. Kepler-102 b likely needs a mass ≳0.1 M {sub ⊕}, hence a bulk density exceeding about half Earth’s, in order to avoid dynamical instability. To investigate the role of secular chaos in our wider set of simulations, we characterize each planetary system’s AMD evolution with a “spectral fraction” calculated from the power spectrum of short integrations (∼5 × 10{sup 6} orbits). We find that small spectral fractions (≲0.01) are strongly associated with dynamical stability on long timescales (5 × 10{sup 9} orbits) and that the median time to instability decreases with increasing spectral fraction. Our results support the hypothesis that secular chaos is the driver of instabilities in many nonresonant multiplanet systems and also demonstrate that the spectral analysis method is an efficient numerical tool to diagnose long-term (in)stability of multiplanet systems from short simulations.

79 ASTRONOMY AND ASTROPHYSICS↗