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 1,315 records · Page 73

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Stability derivatives for bodies of revolution at subsonic speeds

The paper considers a rigid pointed body of revolution in a steady uniform subsonic flow. The body performs harmonic small-amplitude pitching oscillations around its zero angle of attack position. The body is assumed to be smooth and sufficiently slender so that the small perturbation concept can be applied. The basis of the method used, following Revell (1960), is the relation of a body-fixed perturbation potential to the general velocity potential. Normal force distributions as well as total force and moment coefficients are calculated for parabolic spindles and the numerical results show good agreement between Revell's second-order slender body theory and the present theory for the static stability derivatives of the parabolic spindles.

Liu, D. D.↗

NASA NTRS1993

A new approach to mixed H2/H infinity controller synthesis using gradient-based parameter optimization methods

In the past few years, the mixed H(sub 2)/H-infinity control problem has been the object of much research interest since it allows the incorporation of robust stability into the LQG framework. The general mixed H(sub 2)/H-infinity design problem has yet to be solved analytically. Numerous schemes have considered upper bounds for the H(sub 2)-performance criterion and/or imposed restrictive constraints on the class of systems under investigation. Furthermore, many modern control applications rely on dynamic models obtained from finite-element analysis and thus involve high-order plant models. Hence the capability to design low-order (fixed-order) controllers is of great importance. In this research a new design method was developed that optimizes the exact H(sub 2)-norm of a certain subsystem subject to robust stability in terms of H-infinity constraints and a minimal number of system assumptions. The derived algorithm is based on a differentiable scalar time-domain penalty function to represent the H-infinity constraints in the overall optimization. The scheme is capable of handling multiple plant conditions and hence multiple performance criteria and H-infinity constraints and incorporates additional constraints such as fixed-order and/or fixed structure controllers. The defined penalty function is applicable to any constraint that is expressible in form of a real symmetric matrix-inequity.

Ly, Uy-Loi↗

NASA NTRS1993

Stability of outer planetary orbits around binary stars - A comparison of Hill's and Laplace's stability criteria

A comparison is made between the stability criteria of Hill and that of Laplace to determine the stability of outer planetary orbits encircling binary stars. The restricted, analytically determined results of Hill's method by Szebehely and coworkers and the general, numerically integrated results of Laplace's method by Graziani and Black (1981) are compared for varying values of the mass parameter mu. For mu = 0 to 0.15, the closest orbit (lower limit of radius) an outer planet in a binary system can have and still remain stable is determined by Hill's stability criterion. For mu greater than 0.15, the critical radius is determined by Laplace's stability criterion. It appears that the Graziani-Black stability criterion describes the critical orbit within a few percent for all values of mu.

Kubala, A.↗

NASA NTRS1993

Stability of the Halley Cometo-Sheath with Resistivity and Plasma Motion

The MHD stability of the cometary inner shealth determined by the balance between the inward Lorentz body force and the outward ion-neutral drag force is investigated by numerically solving the wave equations which include resistivity, plasma motion and plasma pressure with the help of two-point boundary value method.

Comets Cometosheath MHD Stability↗

NASA NTRS1976

A two-dimensional finite difference solution for the transient thermal behavior of tubular solar collector

A numerical procedure was established using the finite-difference technique in the determination of the time-varying temperature distribution of a tubular solar collector under changing solar radiancy and ambient temperature. Three types of spatial discretization processes were considered and compared for their accuracy of computations and for selection of the shortest computer time and cost. The stability criteria of this technique was analyzed in detail to give the critical time increment to ensure stable computations. The results of the numerical analysis were in good agreement with the analytical solution previously reported. The numerical method proved to be a powerful tool in the investigation of the collector sensitivity to two different flow patterns and several flow control mechanisms.

Lansing, F. L.↗

NASA NTRS1980

A two-dimensional finite-difference solution for the transient thermal behavior of a tubular solar collector

A numerical procedure was established using the finite-difference technique in the determination of the time-varying temperature distribution of a tubular solar collector under changing solar radiancy and ambient temperature. Three types of spatial discretization processes were considered and compared for their accuracy of computations and for selection of the shortest computer time and cost. The stability criteria of this technique were analyzed in detail to give the critical time increment to ensure stable computations. The results of the numerical analysis were in good agreement with the analytical solution previously reported. The numerical method proved to be a powerful tool in the investigation of the collector sensitivity to two different flow patterns and several flow control mechanisms.

Lansing, F. L.↗

DOE OSTI2021

Numerical modeling of a proton spin-flipping system in the spin transparency mode at an integer spin resonance in JINR's Nuclotron

In this paper we propose a lattice insertion for the Nuclotron ring called a “spin navigator” that can adjust any direction of the proton polarization in the orbital plane using weak solenoids. The polarization control is realized in the spin transparency mode at the energy of 108 MeV, which corresponds to the integer spin resonance γ G = 2. The requirements on the navigator solenoid fields are specified considering the criteria for stability of the spin motion during any manipulation of the polarization direction in an experiment. Additionally, this paper presents the results of numerical modeling of the proton spin dynamics in the Nuclotron ring operated in the spin transparency mode. The verified spin navigator is aimed at an experimental study of a spin-flipping system using the Nuclotron ring. The results are relevant to the NICA (JINR), EIC (BNL) and COSY (FZJ) facilities where the spin transparency mode can be applied for polarization control.

47 OTHER INSTRUMENTATION↗

NASA NTRS1991

The asymptotic spectra of banded Toeplitz and quasi-Toeplitz matrices

Toeplitz matrices occur in many mathematical, as well as, scientific and engineering investigations. This paper considers the spectra of banded Toeplitz and quasi-Toeplitz matrices with emphasis on non-normal matrices of arbitrarily large order and relatively small bandwidth. These are the type of matrices that appear in the investigation of stability and convergence of difference approximations to partial differential equations. Quasi-Toeplitz matrices are the result of non-Dirichlet boundary conditions for the difference approximations. The eigenvalue problem for a banded Toeplitz or quasi-Toeplitz matrix of large order is, in general, analytically intractable and (for non-normal matrices) numerically unreliable. An asymptotic (matrix order approaches infinity) approach partitions the eigenvalue analysis of a quasi-Toeplitz matrix into two parts, namely the analysis for the boundary condition independent spectrum and the analysis for the boundary condition dependent spectrum. The boundary condition independent spectrum is the same as the pure Toeplitz matrix spectrum. Algorithms for computing both parts of the spectrum are presented. Examples are used to demonstrate the utility of the algorithms, to present some interesting spectra, and to point out some of the numerical difficulties encountered when conventional matrix eigenvalue routines are employed for non-normal matrices of large order. The analysis for the Toeplitz spectrum also leads to a diagonal similarity transformation that improves conventional numerical eigenvalue computations. Finally, the algorithm for the asymptotic spectrum is extended to the Toeplitz generalized eigenvalue problem which occurs, for example, in the stability of Pade type difference approximations to differential equations.

Beam, Richard M.↗

NASA NTRS1989

Numerical algorithms for computations of feedback laws arising in control of flexible systems

Several continuous models will be examined, which describe flexible structures with boundary or point control/observation. Issues related to the computation of feedback laws are examined (particularly stabilizing feedbacks) with sensors and actuators located either on the boundary or at specific point locations of the structure. One of the main difficulties is due to the great sensitivity of the system (hyperbolic systems with unbounded control actions), with respect to perturbations caused either by uncertainty of the model or by the errors introduced in implementing numerical algorithms. Thus, special care must be taken in the choice of the appropriate numerical schemes which eventually lead to implementable finite dimensional solutions. Finite dimensional algorithms are constructed on a basis of a priority analysis of the properties of the original, continuous (infinite diversional) systems with the following criteria in mind: (1) convergence and stability of the algorithms and (2) robustness (reasonable insensitivity with respect to the unknown parameters of the systems). Examples with mixed finite element methods and spectral methods are provided.

Lasiecka, Irena↗

NASA NTRS1987

Nonparallel stability of boundary layers

The asymptotic formulations of the nonparallel linear stability of incompressible growing boundary layers are critically reviewed. These formulations can be divided into two approaches. The first approach combines a numerical method with either the method of multiple scales, or the method of averaging, of the Wentzel-Kramers-Brillouin (WKB) approximation; all these methods yield the same result. The second approach combined a multi-structure theory with the method of multiple scales. The first approach yields results that are in excellent agreement with all available experimental data, including the growth rates as well as the neutral stability curve. The derivation of the linear stability of the incompressible growing boundary layers is explained.

Nayfeh, Ali H.↗

NASA NTRS2003

Functional Coordination of a Full-Body Gaze Control Mechanisms Elicited During Locomotion

Control of locomotion requires precise interaction between several sensorimotor subsystems. Exposure to the microgravity environment of spaceflight leads to postflight adaptive alterations in these multiple subsystems leading to postural and gait disturbances. Countermeasures designed to mitigate these postflight gait alterations will need to be assessed with a new generation of functional tests that evaluate the interaction of various elements central to locomotor control. The goal of this study is to determine how the multiple, interdependent, full- body sensorimotor subsystems aiding gaze stabilization during locomotion are functionally coordinated. To explore this question two experiments were performed. In the first study (Study 1) we investigated how alteration in gaze tasking changes full-body locomotor control strategies. Subjects (n=9) performed two discreet gaze stabilization tasks while walking at 6.4 km/hr on a motorized treadmill: 1) focusing on a central point target; 2) reading numeral characters; both presented at 2m in front at eye level. The second study (Study 2) investigated the potential of adaptive remodeling of the full-body gaze control systems following exposure to visual-vestibular conflict. Subjects (n=14) walked (6.4 km/h) on the treadmill before and after they were exposed to 0.5X minifying lenses worn for 30 minutes during self-generated sinusoidal vertical head rotations performed while seated. In both studies we measured: temporal parameters of gait, full body sagittal plane segmental kinematics of the head, trunk, thigh, shank and foot, accelerations along the vertical axis at the head and the shank, and the vertical forces acting on the support surface. Results from Study 1 showed that while reading numeral characters as compared to the central point target: 1) compensatory head pitch movements were on average 22% greater 2) the peak acceleration measured at the head was significantly reduced by an average of 13% in four of the six subjects 3) the knee joint total movement was on average 11% greater during the period from the heel strike event to the peak knee flexion event in stance phase of the gait cycle. Results from Study 2 indicate that following exposure to visual-vestibular conflict changes in full-body strategies were observed consistent with the requirement to aid gaze stabilization during locomotion.

Bloomberg, Jacob J.↗

DOE OSTI2026

Thermo-hydro-mechanical analysis of subsurface ice-based thermal energy storage

Ice-based thermal energy storage systems are widely utilized for cooling and managing peak electrical demand globally, offering daily or weekly storage capabilities for both individual homes and larger office buildings. However, scaling these systems for district-level cooling or integrating them with renewable energy sources presents challenges, especially in accommodating larger volumes and addressing seasonal storage requirements in densely populated urban areas. This paper proposes a novel solution by evaluating subsurface ice-based thermal energy storage, in which the underground is subjected to seasonal freeze/thaw cycles. However, these cycles may influence ground behavior, affecting pore pressure and inducing ground movement. To systematically investigate these challenges, we enhance the TOUGH-FLAC simulator by integrating water/ice phase change capabilities and updating the effective stress–strain constitutive relation. Both modifications are validated against analytical solutions or experimental data. Through numerical simulations spanning a decade with ten seasonal freeze/thaw cycles, we evaluate the performance and long-term stability of a generic subsurface ice-based thermal energy storage system, considering factors such as ground permeability, freezing pipe spacing, freeze/thaw damage, and glycol solution temperature. The simulations indicate that ice formation induces pore pressure variations that drive seasonal surface heave and settlement, controlled by ground permeability, pipe spacing, and glycol solution temperature, along with tensile and localized shear deformation around freeze pipes. This highlights the need for accurate ground property characterization and geomechanical analysis for subsurface ice-based thermal energy storage.

Ground freezing↗

DOE OSTI2022

Geometrically stabilized skyrmionic vortex in FeGe tetrahedral nanoparticles

The concept of topology has dramatically expanded the research landscape of magnetism, leading to the discovery of numerous magnetic textures with intriguing topological properties. A magnetic skyrmion is an emergent topological magnetic texture with a string-like structure in three dimensions and a disk-like structure in one and two dimensions. Skyrmions in zero dimensions have remained elusive due to challenges from many competing orders. Here, by combining electron holography and micromagnetic simulations, we uncover the real-space magnetic configurations of a skyrmionic vortex structure confined in a B20-type FeGe tetrahedral nanoparticle. An isolated skyrmionic vortex forms at the ground state and this texture shows excellent robustness against temperature without applying a magnetic field. Furthermore, our findings shed light on zero-dimensional geometrical confinement as a route to engineer and manipulate individual skyrmionic metastructures.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Improvement of Numerical Methods in Pronghorn

During the fiscal year of 2021, the finite-volume method (FVM) was deployed in Pronghorn to provide improved efficiency, stability, and accuracy for coarse-mesh, thermal-hydraulics problems. While the main goals of fiscal year 2021 were met, several issues emerged from the early deployment of the finite-volume method in Pronghorn. These issues were: The compressible and incompressible formulations that were implemented are inadequate for many nuclear reactor flow problems; Omission of terms accounting for the porosity and Darcy-Forchheimer body force discontinuities in the Rhie-Chow interpolation lead to oscillations in pressure and velocity at these discontinuities; The FVM lacks a correction for non-orthogonal grids for computing accurate pressure gradient leading to loss of accuracy in regions with skewed elements; The FVM currently uses a monolithic solver. Monolithic solvers have issues dealing with the saddle-point nature of the discretized fluid equations. The result are bad convergence if direct factorization is not used and large memory consumption when direct factorization is used. These four issues are addressed in this report. In particular, we report the completion of the following task: Implementation of a weakly compressible formulation in the MOOSE Navier-Stokes module; Implementation of Moukalled’s method for including body forces in the RCI. Additionally, we identified the need to smooth the porosity using Moukalled’s face-cell smoothing operator; Implementation of a non-orthogonal correction for the Green-Gauss gradient computation; Preliminary implementation of a SIMPLE segregated solver.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1992

The effect of initial imperfections on shell stability: An updated review

The development of the Delft Interactive Shell Design Code (DISDECO) is described. The purpose is to make the accumulated theoretical, numerical, and practical knowledge of the last years readily accessible to users interested in the analysis of buckling sensitive structures. With this open ended, hierarchical, interactive computer code, the user can access workstation programs of increasing complexity. As a first update of DISDECO, the computational modules of the initial level were modified so as to handle combined loads consisting of axial compression, internal and external pressure, and clockwise or counter clockwise torque. The applied load is assumed to have a uniform spatial distribution and is divided into a fixed and a variable part. The magnitude of the variable part is allowed to vary in proportion to a load parameter Lambda. This leads to an eigenvalue problem for the critical lambda sub C.

Arbocz, J.↗

NASA NTRS2001

Driving Parameters for Distributed and Centralized Air Transportation Architectures

This report considers the problem of intersecting aircraft flows under decentralized conflict avoidance rules. Using an Eulerian standpoint (aircraft flow through a fixed control volume), new air traffic control models and scenarios are defined that enable the study of long-term airspace stability problems. Considering a class of two intersecting aircraft flows, it is shown that airspace stability, defined both in terms of safety and performance, is preserved under decentralized conflict resolution algorithms. Performance bounds are derived for the aircraft flow problem under different maneuver models. Besides analytical approaches, numerical examples are presented to test the theoretical results, as well as to generate some insight about the structure of the traffic flow after resolution. Considering more than two intersecting aircraft flows, simulations indicate that flow stability may not be guaranteed under simple conflict avoidance rules. Finally, a comparison is made with centralized strategies to conflict resolution.

Feron, Eric↗

DOE OSTI2021

Impact of Increased Inverter-Based Resources on Power System Small-Signal Stability

The transformation of the power system to include more distributed energy resources (DER) implies an increase in the number of inverter-based resources deployed on the grid. Envisioning future scenarios, this paper presents a small-signal stability analysis for a power grid comprising synchronous generators and inverter-based resources. Three types of inverter control are considered: grid following, droop-controlled grid forming, and virtual oscillator control grid forming. Although small-signal stability of power systems is a widely studied topic, systematic analysis of mixed machine-inverter systems with detailed control models at various inverter levels are limited. This paper addresses the gap with numerical simulations tailored to the IEEE 39-bus system. Results show the system may become unstable at high inverter level of grid-following inverters, and grid-forming inverter control can potentially improve system stability, thereby enabling very high level of DERs.

grid-following inverter↗

DOE OSTI2021

A scalable exponential-DG approach for nonlinear conservation laws: With application to Burger and Euler equations

In this work, we propose an Exponential DG framework for partial differential equations. We decompose 7 governing equations into linear and nonlinear parts to which we apply the discontinuous Galerkin 8 (DG) spatial discretization. In particular, we construct the linear part using Jacobian that effectively 9 capture stiff characteristics in the system. The former is integrated analytically, whereas the latter 10 is approximated. This approach i) is stable with a large Courant number (Cr > 1); ii) supports 11 high-order solutions both in time and space; iii) is computationally favorable compared to IMEX 12 DG methods with no preconditioner; iv) becomes comparable to explicit RKDG methods on uniform 13 mesh and beneficial on non-uniform grid for Euler equations; v) is scalable in a modern massively 14 parallel computing architecture due to its explicit nature of exponential time integrators and com15 pact communication stencil of DG method. Numerical results demonstrate the performance of our 16 proposed methods through various examples. We also discuss the stability and convergence analysis 17 for our exponential DG scheme in the context of Burgers equation.

42 ENGINEERING↗