Finite time stochastic stability and the analysis of tracking systems.
Liapunov functions for finite time stochastic stability and analysis of tracking system
SYSTEM STABILITY↗
SEARCH · Engineering Papers
Results for “stability analysis”
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.
A describing function analysis in closed loop systems with pulse frequency modulators
Describing functions for stability analysis of integral pulse frequency modulated unity feedback closed loop system
Fitzgerald, J. F.↗
Computer analysis of imperfection sensitivity of ring-stiffened orthotropic shells of revolution
Stability analysis of ring-stiffened orthotropic, multilayered shells of revolution under axisymmetric torsionless loads by digital computer
Cohen, G. A.↗
Analysis of a control system containing two nonlinearities
Parameter plane method of stability analysis applied to spacecraft attitude control system with two nonlinearities related by nonlinear differential equation
Seltzer, S. M.↗
Impact of equilibrium radial electric field on energy loss process after pedestal collapse
An impact of the equilibrium radial electric field on energy loss processes after pedestal collapse is numerically investigated using the BOUT++ framework. Using linear stability analysis, the resistive ballooning mode is shown to be stabilized by the radial shear of the equilibrium radial electric field. On the other hand, the energy loss level after the pedestal collapse increases if the equilibrium radial electric field is taken into account. The spatio-temporal and phase diagram analyses show that the equilibrium radial electric field partially cancels the fluctuation-driven toroidally axisymmetric radial electric field and weakens the E × B shearing rate after pedestal collapse, weakening the turbulence suppression by vortex shearing. Finally, the equilibrium radial electric field therefore increases turbulence intensity in nonlinear cyclic oscillations among pressure gradient, E × B shearing rate, and turbulence intensity, which gives rise to subsequent bursts of turbulent transport and increases the energy loss level.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Computation of linear MHD instabilities with the multi-region relaxed MHD energy principle
Over the last decade, a variational principle based on a generalisation of Taylor's relaxation, referred to as multi-region relaxed magnetohydrodynamics (MRxMHDs) has been developed. The numerical solutions of the MRxMHD equilibria have been constructed using the Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502). In principle, SPEC could also be established to describe the MRxMHD stability of an equilibrium. In this work, a theoretical framework is developed to relate the second variation of the energy functional to the so-called Hessian matrix, enabling the prediction of MHD linear instabilities of cylindrical plasmas, and is implemented in SPEC. The negative and positive eigenvalues of the Hessian matrix predict the stability of an equilibrium. Verification studies of SPEC stability results with the M3D-C 1 code and the tearing mode $\Delta^{^{\prime}}$ criterion have been conducted for ideal and resistive MHD instabilities, respectively, in a pressureless cylindrical tokamak, and have shown good agreement. Overall, our stability analysis is a critical step towards understanding the MHD stability of three-dimensional MHDs where nested flux surfaces, magnetic islands and stochastic regions co-exist.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Two-Level Sketching Alternating Anderson Acceleration for Complex Physics Applications
We present a novel two-level sketching extension of the Alternating Anderson–Picard (AAP) method for accelerating fixed-point iterations in challenging single- and multiphysics simulations governed by discretized PDEs. Our approach combines a static, physics-based projection that reduces the least-squares (LS) problem to the most informative field (e.g., via Schur-complement insight) with a dynamic, algebraic sketching stage driven by a backward stability analysis under Lipschitz continuity. We introduce inexpensive estimators for stability thresholds and cache-aware randomized selection strategies to balance computational cost against memory access overhead. The resulting algorithm solves reduced LS systems in place, minimizes memory footprints, and seamlessly alternates between low-cost Picard updates and Anderson mixing. Implemented in Julia, our two-level sketching AAP achieves up to 50% time-to-solution reductions compared to standard Anderson acceleration—without degrading convergence rates—on benchmark problems including Stokes, 𝑝-Laplacian, bidomain, and Navier–Stokes formulations at varying problem sizes. These results demonstrate the method’s robustness, scalability, and potential for integration into high-performance scientific computing frameworks. Our implementation is available open source in the AAP.jl library.
Barnafi, Nicolas [University of Chile, Santiago]↗
Existence and Stability of Menisci in Detached Bridgman Growth
Detached growth, also referred to as dewetted growth, is a Bridgman crystal growth process in which the melt is in contact with the crucible wall but the crystal is not. A meniscus bridges the gap between the top of the crystal and the crucible wall. The Young-Laplace capillary equation was used to calculate the crystal radii of detached states as a function of the pressure differential across the meniscus. The detached states depend on the contact angle of the melt with the crucible wall, the growth angle of the melt with respect to the solidifying crystal, and the Bond number. A static stability analysis was performed on the calculated detached states. The stability criterion was the sign of the second variation of the potential energy upon admissible meniscus shape perturbations. The conditions considered corresponded to the growth of Ge and InSb, in both terrestrial and microgravity conditions. Stability was found to depend significantly on whether the interior surface was considered to be microscopically rough or smooth, corresponding to pinned or unpinned states. It was also found that all meniscus shapes which are single-valued functions of the radius are statically stable in a microgravity environment.
Volz, M. P.↗
Experimental Determination of the Dynamic Hydraulic Transfer Function for the J-2X Oxidizer Turbopump: Methodology - Part One
An advanced methodology for extracting the hydraulic dynamic pump transfer matrix (Yp) for a cavitating liquid rocket engine turbopump inducer+impeller has been developed. The transfer function is required for integrated vehicle pogo stability analysis as well as optimization of local inducer pumping stability. Laboratory pulsed subscale waterflow test of the J-2X oxygen turbo pump is introduced and our new extraction method applied to the data collected. From accurate measures of pump inlet and discharge perturbational mass flows and pressures, and one-dimensional flow models that represents complete waterflow loop physics, we are able to derive Yp and hence extract the characteristic pump parameters: compliance, pump gain, impedance, mass flow gain. Detailed modeling is necessary to accurately translate instrument plane measurements to the pump inlet and discharge and extract Yp. We present the MSFC Dynamic Lump Parameter Fluid Model Framework and describe critical dynamic component details. We report on fit minimization techniques, cost (fitness) function derivation, and resulting model fits to our experimental data are presented. Comparisons are made to alternate techniques for spatially translating measurement stations to actual pump inlet and discharge.
Zoladz, Tom↗
Experimental Waterflow Determination of the Dynamic Hydraulic Transfer Function for the J-2X Oxidizer Turbopump: Results and Interpretation - Part Two
Experimental results describing the hydraulic dynamic pump transfer matrix (Yp) for a cavitating J-2X oxidizer turbopump inducer+impeller tested in subscale waterflow are presented. The transfer function is required for integrated vehicle pogo stability analysis as well as optimization of local inducer pumping stability. Dynamic transfer functions across widely varying pump hydrodynamic inlet conditions are extracted from measured data in conjunction with 1D-model based corrections. Derived Dynamic transfer functions are initially interpreted relative to traditional Pogo pump equations. Water-to-liquid oxygen scaling of measured cavitation characteristics are discussed. Comparison of key dynamic transfer matrix terms derived from waterflow testing are made with those implemented in preliminary Ares Upper Stage Pogo stability modeling. Alternate cavitating pump hydraulic dynamic equations are suggested which better reflect frequency dependencies of measured transfer matrices.
Zoladz, Tom↗
Study of Anti-Vortex Baffle Effect in Suppressing Swirling Flow in LOX Tank
Experimental results describing the hydraulic dynamic pump transfer matrix (Yp) for a cavitating J-2X oxidizer turbopump inducer+impeller tested in subscale waterflow are presented. The transfer function is required for integrated vehicle pogo stability analysis as well as optimization of local inducer pumping stability. Dynamic transfer functions across widely varying pump hydrodynamic inlet conditions are extracted from measured data in conjunction with 1D-model based corrections. Derived Dynamic transfer functions are initially interpreted relative to traditional Pogo pump equations. Water-to-liquid oxygen scaling of measured cavitation characteristics are discussed. Comparison of key dynamic transfer matrix terms derived from waterflow testing are made with those implemented in preliminary Ares Upper Stage Pogo stability modeling. Alternate cavitating pump hydraulic dynamic equations are suggested which better reflect frequency dependencies of measured transfer matrices.
Yang, H. Q.↗
Transition Within Leeward Plane of Axisymmetric Bodies at Incidence in Supersonic Flow
Boundary layer transition along the leeward symmetry plane of axisymmetric bodies at nonzero angle of attack in supersonic flow was investigated experimentally and numerically as part of joint research between the Japan Aerospace Exploration Agency (JAXA) and National Aeronautics and Space Administration (NASA). Transition over four axisymmetric bodies (namely, Sears-Haack body, semi-Sears-Haack body, straight cone and flared cone) with different axial pressure gradients was measured in two different facilities with different unit Reynolds numbers. The semi-Sears-Haack body and flared cone were designed at JAXA to broaden the range of axial pressure distributions. For a body shape with an adverse pressure gradient (i.e., flared cone), the experimentally measured transition patterns show an earlier transition location along the leeward symmetry plane in comparison with the neighboring azimuthal locations. For nearly zero pressure gradient (i.e.,straight cone), this feature is only observed at the larger unit Reynolds number. Later transition along the leeward plane was observed for the remaining two body shapes with a favorable pressure gradient. The observed transition patterns are only partially consistent with the numerical predictions based on linear stability analysis. Additional measurements are used in conjunction with the stability computations to explore the phenomenon of leeward line transition and the underlying transition mechanism in further detail.
Tokugawa, Naoko↗
Design optimization of natural laminar flow fuselages in compressible flow
An optimization method has been developed to design axisymmetric body shapes such as fuselages, nacelles, and external fuel tanks with increased transition Reynolds numbers in subsonic compressible flow. The new design method involves a constraint minimization procedure coupled with analysis of the inviscid and viscous flow regions, and linear stability analysis of the compressible boundary layer. Boundary-layer transition is predicted by a 'hybrid' transition criterion based on Granville's transition criterion and a criterion using linear stability theory coupled with the e(n)-method. A tiptank of a business-jet is used as an example to illustrate that the method can be utilized to design an axisymmetric body shape with extensive natural laminar flow. Boundary layer transition is predicted to occur at a transition Reynolds number of 6.04 x 10 to the 6th on the original tiptank. On the designed body shape, a transition Reynolds number of 7.22 x 10 to the 6th is predicted using compressible linear stability theory coupled with e(n)-method.
Dodbele, Simha S.↗
Design of laminar flow bodies in compressible flow
An optimization method was developed to design axisymmetric body shapes such as fuselages, nacelles, and external fuel tanks with increased transition Reynolds numbers in subsonic compressible flow. The method involves a constraint minimization procedure coupled with analysis of the inviscid and viscous flow regions, and linear stability analysis of the compressible boundary-layer. Boundary-layer transition is predicted by a 'hybrid' transition criterion based on Granville's transition criterion and a criterion using linear stability theory coupled with the e(sup n)-method. A tiptank of a business jet is used as an example to illustrate that the method can be utilized to design an axisymmetric body shape with extensive natural laminar flow. On the original tiptank boundary layer transition is predicted to occur at a transition Reynolds number of 6.04 x 10 6 on the original tiptank. On the designed body shape a transition Reynolds number of 7.22 x 10 6 is predicted using compressible linear stability theory coupled with e n -method.
Simha S. Dodbele↗
Design optimization of natural laminar flow bodies in compressible flow
An optimization method has been developed to design axisymmetric body shapes such as fuselages, nacelles, and external fuel tanks with increased transition Reynolds numbers in subsonic compressible flow. The new design method involves a constraint minimization procedure coupled with analysis of the inviscid and viscous flow regions and linear stability analysis of the compressible boundary-layer. In order to reduce the computer time, Granville's transition criterion is used to predict boundary-layer transition and to calculate the gradients of the objective function, and linear stability theory coupled with the e(exp n)-method is used to calculate the objective function at the end of each design iteration. Use of a method to design an axisymmetric body with extensive natural laminar flow is illustrated through the design of a tiptank of a business jet. For the original tiptank, boundary layer transition is predicted to occur at a transition Reynolds number of 6.04 x 10(exp 6). For the designed body shape, a transition Reynolds number of 7.22 x 10(exp 6) is predicted using compressible linear stability theory coupled with the e(exp n)-method.
Dodbele, Simha S.↗
Numerical Stability and Control Analysis Towards Falling-Leaf Prediction Capabilities of Splitflow for Two Generic High-Performance Aircraft Models
Aerodynamic analysis are performed using the Lockheed-Martin Tactical Aircraft Systems (LMTAS) Splitflow computational fluid dynamics code to investigate the computational prediction capabilities for vortex-dominated flow fields of two different tailless aircraft models at large angles of attack and sideslip. These computations are performed with the goal of providing useful stability and control data to designers of high performance aircraft. Appropriate metrics for accuracy, time, and ease of use are determined in consultations with both the LMTAS Advanced Design and Stability and Control groups. Results are obtained and compared to wind-tunnel data for all six components of forces and moments. Moment data is combined to form a "falling leaf" stability analysis. Finally, a handful of viscous simulations were also performed to further investigate nonlinearities and possible viscous effects in the differences between the accumulated inviscid computational and experimental data.
Charlton, Eric F.↗
A study analysis of cable-body systems totally immersed in a fluid stream
A general stability analysis of a cable-body system immersed in a fluid stream is presented. The analytical portion of this analysis treats the system as being essentially a cable problem, with the body dynamics giving the end conditions. The mathematical form of the analysis consists of partial differential wave equations, with the end and auxiliary conditions being determined from the body equations of motion. The equations uncouple to give a lateral problem and a longitudinal problem as in first order airplane dynamics. A series of tests on a tethered wind tunnel model provide a comparison of the theory with experiment.
Delaurier, J. D.↗
Multidisciplinary aeroelastic analysis of a generic hypersonic vehicle
This paper presents details of a flutter and stability analysis of aerospace structures such as hypersonic vehicles. Both structural and aerodynamic domains are discretized by the common finite element technique. A vibration analysis is first performed by the STARS code employing a block Lanczos solution scheme. This is followed by the generation of a linear aerodynamic grid for subsequent linear flutter analysis within subsonic and supersonic regimes of the flight envelope; the doublet lattice and constant pressure techniques are employed to generate the unsteady aerodynamic forces. Flutter analysis is then performed for several representative flight points. The nonlinear flutter solution is effected by first implementing a CFD solution of the entire vehicle. Thus, a 3-D unstructured grid for the entire flow domain is generated by a moving front technique. A finite element Euler solution is then implemented employing a quasi-implicit as well as an explicit solution scheme. A novel multidisciplinary analysis is next effected that employs modal and aerodynamic data to yield aerodynamic damping characteristics. Such analyses are performed for a number of flight points to yield a large set of pertinent data that define flight flutter characteristics of the vehicle. This paper outlines the finite-element-based integrated analysis procedures in detail, which is followed by the results of numerical analyses of flight flutter simulation.
Gupta, K. K.↗