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.

SourceResearch connectionFrom yearThrough year

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

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS1976

The effects of spherical geometry on baroclinic instability

A baroclinic stability analysis is performed for a simple family of zonal shear profiles over a sphere, using a two-layer, quasi-geostrophic model. The stability properties and the structure of the most unstable waves are qualitatively similar to those on a beta-plane. However, the spherical geometry plays a major role in locating some of the important features of the most unstable waves. In particular, the locations of the maximum wave amplitude, maximum eddy heat fluxes, and maximum convergence of the eddy angular momentum flux are all well correlated with the location of the maximum excess of the vertical shear over the minimum value necessary for local instability on a sphere. Consequently the eddy momentum flux tends to generate a mid-latitude jet even if there is no preexisting mid-latitude jet in the basic state zonal flow. These findings suggest some of the elements needed for parameterizing the meridional variations of baroclinic eddy fluxes accurately.

Moura, A. D.↗

NASA NTRS1990

Morphological instability in rapid directional solidification

Mullins and Sekerka (1964) showed for fixed temperature gradient that the planar interface is linearly stable for all pulling speeds V above some critical value, the absolute stability limit. Near this limit, where solidification rates are rapid, the assumption of local equilibrium at the interface may be violated. Here, nonequilibrium effects are incorporated into a linear stability analysis of the planar front by allowing the segregation coefficient and interface temperature to depend on V in a thermodynamically consistent way. The absolute stability limit of the cellular mode is modified. A new oscillatory state is formed which, in the absence of latent heat, has a critical wavenumber of zero; by itself this instability would lead to the formation of solute bands in the solid. This mode has its own absolute-stability limit determined by solute trapping and kinetics. Under certain conditions, there exists a window of stability above the steady absolute-stability boundary and below the oscillatory-stability boundary; here the planar segregation-free state is restabilized.

Merchant, G. J.↗

DOE OSTI2023

Multimode theory of electron hole transverse instability

We present Vlasov–Poisson three-dimensional linear stability analysis of an initially planar electron hole structure, solving for the distribution function by integration along unperturbed orbits. The non-sinusoidal potential perturbation shape (parallel to $B$ ) is expanded in eigenfunctions of the adiabatic Poisson operator, generalizing the prior assumption of a rigid shift of the equilibrium. We show that the shiftmode is then modified by a second discrete mode plus an integral over a continuum of wave-like modes. A rigorous treatment shows that the continuum can be approximated effectively by a single mode that satisfies the external wave dispersion relation, thus making the perturbation a weighted sum of three modes. We find numerically the solution for the complex instability frequency, and the corresponding three mode amplitudes determining the perturbation eigenmode. This multimode analysis refines the accuracy of the prior single-mode results, giving slightly higher growth rates at most parameters, as expected from the extra mode shape freedom. Oscillating modes near stability boundaries have larger mode distortions which help explain particle-in-cell simulations that observe instability up to ${\sim }20$ % beyond the prior shiftmode thresholds, and narrowing of the perturbation. At high magnetic field, the multimode analysis predicts a reduction of the already small growth rate.

Physics↗

NASA NTRS1994

Use of rotation to suppress thermosolutal convection in directionally solidified binary alloys

Effects of rotation on onset of convection during plane-front directional solidification of Pb-Sn and the pseudobinary system mercury cadmium telluride (Hg(1-x)Cd(x)Te), and on dendritic solidification of Pb-Sn have been studied by means of linear stability analysis. Incorporating Coriolis and centrifugal accelerations into the momentum equation of Coriell et al., we find that under realistic processing conditions, a large degree of stabilization can be achieved using modest rotation rates for both Pb-Sn and mercury cadmium telluride. At a growth velocity of 5 micron/sec and nominal liquid-side temperature gradient of 200 K/cm in Pb-Sn, rotation at 500 rpm results in a hundredfold increase in the critical Sn concentration. Large increases in the maximum allowable growth velocity at fixed melt composition are also attainable with modest rotation rates. The effect is amplified under conditions of reduced gravitational acceleration. For Hg(1-x)Cd(x)Te, we have also studied the nonrotating case. The key differences are due to the existence of a composition range for Hg(1-x)Cd(x)Te in which the melt density has a local maximum as a function of temperature. When the melt solidifies by cooling from below, the liquid density may initially increase with distance above the interface, before ultimately decreasing as the melt temperature increases above the value at which the local density maximum occurs. In contrast to the Pb-Sn case where density depends monotonically on temperature and composition, for Hg(1-x)Cd(x)Te there exists a critical value of the growth velocity above which plane-front solidification is unstable for all bulk CdTe mole fractions. Again, rotation leads to significant inhibition of onset. We identify the predicted stabilization with the Taylor-Proudman mechanism by which rotation inhibits thermal convection in a single-component fluid heated from below. In a binary liquid undergoing solidification, rotation inhibits the onset of buoyancy-driven convection, and has no effect on the short-wavelength morphological instability. At large growth velocities, the plane-front interface between liquid and solid becomes unstable with respect to a morphological instability and solidification occurs dendritically, with a mushy zone of dendrites and interdendritic fluid separating the solid from the melt. For the Pb-Sn system, rotation substantially suppresses the onset of convection in the mushy zone and in the overlying liquid, holding open the promise that rotation can suppress freckling and other macrosegregation defects.

Pearlstein, Arne J.↗

NASA NTRSDate not supplied

Aerodynamic Design Optimization for Natural Laminar Flow Airfoils

Natural laminar flow technology is a passive laminar flow control (LFC) strategy that seeks to delay the onset of boundary-layer transition (BLT) through shape optimization to reduce the drag of the aerodynamic vehicle. Adjoint-based design optimization for LFC is proposed in an integrated multidisciplinary framework, which includes the computational fluid dynamics (CFD), geometry and grid deformation, and linear stability analysis (LSA) for transition prediction. In particular, the BLT location is predicted using the dual N-factor method that is based on a linear stability theory (LST) eigenvalue problem. The dual N-factor criterion accounts for the amplification of planar Tollmien-Schlichting (TS) and stationary crossflow (CF) boundary-layer instabilities to predict the transition location in three-dimensional boundary-layer flows. The adjoint-based shape optimization procedure is based on an iteratively coupled CFD and LSA methodology to converge the transition location and flow solutions, as well as to calculate the sensitivities of the aerodynamic metrics of interest with respect to the flow and shape design parameters. The RAE 2822 airfoil at 0 and 30 degrees yaw angles, an angle of attack of 0.72 degrees, and subsonic conditions (M∞ = 0.19, Rec = 5.6 × 106 ) are used as baseline configurations for design optimization. The angle of attack and the vertical displacement of free-form-deformation control points are used as design variables to reduce the drag coefficient while reaching a specified lift coefficient. The optimized unswept airfoil designs achieve a 30% drag reduction accompanied by a downstream shift of the transition locations over both suction and pressure sides of the airfoil. The initial design iterations for the swept case also show a favorable trend in the drag reduction with transition delay over both sides.

Transition↗

DOE OSTI2024

Architecture of TOI-561 planetary system

ABSTRACT We present new observations from CHEOPS (CHaracterising ExOPlanet Satellite) and TESS (Transiting Exoplanet Survey Satellite) to clarify the architecture of the planetary system hosted by the old Galactic thick disc star TOI-561. Our global analysis, which also includes previously published photometric and radial velocity data, incontrovertibly proves that TOI-561 is hosting at least four transiting planets with periods of 0.44 d (TOI-561 b), 10.8 d (TOI-561 c), 25.7 d (TOI-561 d), and 77.1 d (TOI-561 e) and a fifth non-transiting candidate, TOI-561f with a period of 433 d. The precise characterization of TOI-561’s orbital architecture is interesting since old and metal-poor thick disc stars are less likely to host ultrashort-period super-Earths like TOI-561 b. The new period of planet -e is consistent with the value obtained using radial velocity alone and is now known to be $77.14399\pm 0.00025$ d, thanks to the new CHEOPS and TESS transits. The new data allowed us to improve its radius ($R_p = 2.517 \pm 0.045\,\mathrm{ R}_{\rm{\oplus }}$ from 5 per cent to 2 per cent precision) and mass ($M_p = 12.4 \pm 1.4\, \mathrm{ M}_{\rm{\oplus }}$) estimates, implying a density of $\rho _p = 0.778 \pm 0.097\, \rho _{\rm{\oplus }}$. Thanks to recent TESS observations and the focused CHEOPS visit of the transit of TOI-561 e, a good candidate for exomoon searches, the planet’s period is finally constrained, allowing us to predict transit times through 2030 with 20-min accuracy. We present an updated version of the internal structure of the four transiting planets. We finally performed a detailed stability analysis, which confirmed the long-term stability of the outer planet TOI-561 f.

Piotto, G. (ORCID:0000000299376387)↗

NASA NTRS1974

The modified equation approach to the stability and accuracy analysis of finite-difference methods

The stability and accuracy of finite-difference approximations to simple linear partial differential equations are analyzed by studying the modified partial differential equation. Aside from round-off error, the modified equation represents the actual partial differential equation solved when a numerical solution is computed using a finite-difference equation. The modified equation is derived by first expanding each term of a difference scheme in a Taylor series and then eliminating time derivatives higher than first order by certain algebraic manipulations. The connection between 'heuristic' stability theory based on the modified equation approach and the von Neumann (Fourier) method is established. In addition to the determination of necessary and sufficient conditions for computational stability, a truncated version of the modified equation can be used to gain insight into the nature of both dissipative and dispersive errors.

Warming, R. F.↗

DOE OSTI2020

Energy-Casimir, dynamically accessible, and Lagrangian stability of extended magnetohydrodynamic equilibria

The formal stability analysis of Eulerian extended magnetohydrodynamics (XMHD) equilibria is considered within the noncanonical Hamiltonian framework by means of the energy-Casimir variational principle and the dynamically accessible stability method. Specifically, we find explicit sufficient stability conditions for axisymmetric XMHD and Hall MHD (HMHD) equilibria with toroidal flow and for equilibria with arbitrary flow under constrained perturbations. The dynamically accessible, second-order variation of the Hamiltonian, which can potentially provide explicit stability criteria for generic equilibria, is also obtained. Moreover, we examine the Lagrangian stability of the general quasineutral two-fluid model written in terms of MHD-like variables, by finding the action and the Hamiltonian functionals of the linearized dynamics, working within a mixed Lagrangian-Eulerian framework. Upon neglecting electron mass, we derive a HMHD energy principle, and in addition, the perturbed induction equation arises from Hamilton’s equations of motion in view of a consistency condition for the relation between the perturbed magnetic potential and the canonical variables.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1978

An extension of A-stability to alternating direction implicit methods

An alternating direction implicit (ADI) scheme was constructed by the method of approximate factorization. An A-stable linear multistep method (LMM) was used to integrate a model two-dimensional hyperbolic-parabolic partial differential equation. Sufficient conditions for the A-stability of the LMM were determined by applying the theory of positive real functions to reduce the stability analysis of the partial differential equations to a simple algebraic test. A linear test equation for partial differential equations is defined and then used to analyze the stability of approximate factorization schemes. An ADI method for the three-dimensional heat equation is also presented.

Warming, R. F.↗

NASA NTRS1987

Convenient stability criteria for difference approximations of hyperbolic initial-boundary value problems. II

The results of Goldberg and Tadmor (1985) are extended to achieve improved stability criteria for a large class of approximations to the initial boundary value problem associated with a particular hyperbolic system in a quarter plane. In a stability analysis, it is shown that the entire approximation is stable if and only if the scalar outflow components of its principal part are stable. Thus, the global stability question is reduced to that of a scalar, homogeneous outflow problem. The stability criteria for the reduced problem, which depend both on the basic scheme and the boundary conditions, but very little on the interaction between the two, are stated and used to establish previous examples and new ones, including a host of dissipative and nondissipative examples. There is no difficulty in extending the stability criteria to two-boundary problems and initial-boundary value problems with variable coefficients.

Goldberg, Moshe↗

NASA NTRS1992

Cellular instability in rapid directional solidification - Bifurcation theory

Merchant and Davis performed a linear stability analysis on a model for the directional solidification of a dilute binary alloy valid for all speeds. The analysis revealed that nonequilibrium segregation effects modify the Mullins and Sekerka cellular mode, whereas attachment kinetics has no effect on these cells. In this paper, the nonlinear stability of the steady cellular mode is analyzed. A Landau equation is obtained that determines the amplitude of the cells. The Landau coefficient here depends on both nonequilibrium segregation effects and attachment kinetics. This equation gives the ranges of parameters for subcritical bifurcation (jump transition) or supercritical bifurcation (smooth transition) to cells.

Braun, R. J.↗

NASA NTRS1993

Thermohydrodynamic analysis of cryogenic liquid turbulent flow fluid film bearings

A thermohydrodynamic analysis is presented and a computer code developed for prediction of the static and dynamic force response of hydrostatic journal bearings (HJB's), annular seals or damper bearing seals, and fixed arc pad bearings for cryogenic liquid applications. The study includes the most important flow characteristics found in cryogenic fluid film bearings such as flow turbulence, fluid inertia, liquid compressibility and thermal effects. The analysis and computational model devised allow the determination of the flow field in cryogenic fluid film bearings along with the dynamic force coefficients for rotor-bearing stability analysis.

Andres, Luis San↗

NASA NTRS1990

Ballooning mode instability at the plasmapause

The ballooning mode instability, which can excite hydromagnetic waves at the plasmapause, is studied in the presence of azimuthal plasma flows induced during geomagnetically disturbed periods. A general sufficient criterion for the ballooning mode stability is derived, for a change in the potential energy greater than or equal to 0, which involves the integration over an entire field line. A local stability analysis at the equatorial plasmapause region shows that the ballooning modes could be spontaneously generated via instability under at least two conditions: one is similar to the usual interchange condition, and the second to the quasi-interchange modes. Both of these local instability conditions can be derived from the general stability criterion. Finally an exact solution for the equilibrium state with flow is derived analytically, and the change in the potential energy is computed numerically. It is found that, in the cases studied, the flow does not spontaneously excite the ballooning modes; it only further stabilizes (or destabilizes) the ballooning spectrum if originally the system is stable (or unstable). The analysis would be useful for the interpretation of some of the low-frequency modes observed at the ground and near the equatorial plasmapause.

Lakhina, G. S.↗

DOE OSTI2025

Small-Signal Stability of Grid-Forming Inverters Using Current-Limiting and Frequency Stabilization

This paper presents a small-signal stability analysis of grid-forming (GFM) inverters under current-limiting conditions. It examines how adjustments in virtual impedance angles, implemented through advanced current-limiting and frequency stabilization techniques, influence small-signal stability. This paper studies a GFM inverter control integrating a fictitious power technique stabilizing primary control by adding a virtual power term and a hybrid current limiter integrating virtual impedance in the anti-wind-up feedback with current reference saturation limiting. A small-signal model is developed to assess the impact of virtual impedance angles on GFM inverter dynamics during grid disturbances, such as voltage drops. The findings indicate that although increasing the virtual impedance angle (to make it more inductive) enhances large-signal stability and voltage support during faults, it can induce oscillations and lead to instability if the angle exceeds certain thresholds. Based on the small-signal models, this paper provides design considerations for the current-limiter impedances to ensure reliable GFM inverter behavior under grid disturbances while maintaining small-signal stability.

current limiting↗

NASA NTRS1976

Excitation of the ordinary electromagnetic mode in low-beta plasmas

The excitation of the ordinary mode propagating perpendicular to an external magnetic field in a plasma is discussed for the case where the ratio of the kinetic energy associated with electron motion parallel to the external magnetic field to the energy density of the magnetic field is low. The analysis is performed for plasma in which in addition to the thermal electrons there is a hot anisotropic electron component. Relativistic effects are taken into account in the stability analysis. Both a ring distribution and a loss cone distribution are used to model the hot electron species. For the ring distribution typical growth rates are larger than that for the case of the loss cone distribution. Excitation of the ordinary mode is found at numerous harmonics of the cyclotron frequency for the ring distribution, however, the loss cone distribution is found to give instability primarily at the fundamental harmonic.

Freund, H. P.↗

NASA NTRS1989

Solution and sensitivity analysis of a complex transcendental eigenproblem with pairs of real eigenvalues

This paper considers complex transcendental eigenvalue problems where one is interested in pairs of eigenvalues that are restricted to take real values only. Such eigenvalue problems arise in dynamic stability analysis of nonconservative physical systems, i.e., flutter analysis of aeroelastic systems. Some available solution methods are discussed and a new method is presented. Two computational approaches are described for analytical evaluation of the sensitivities of these eigenvalues when they are dependent on other parameters. The algorithms presented are illustrated through examples.

Murthy, Durbha V.↗

NASA NTRS1989

Solution and sensitivity analysis of a complex transcendental eigenproblem with pairs of real eigenvalues

Murthy, D. V.↗

NASA NTRS1992

Effect of surface tension on the onset of convection in a double-diffusive layer

The effect of surface tension on the stability of a double-diffusive layer is considered using linear stability analysis. The surface tension is assumed to vary linearly with temperature and solute concentration. The eigenvalue problem is solved by the Galerkin method. Results show that the predicted stability boundary based on Marangoni effects alone is completely altered in the presence of buoyancy effects induced by low gravity levels (about 10 exp -5 g). At reduced gravity levels, salt-finger instability may onset in the overstable mode due to the stabilizing effect of surface tension. Fluid properties in terms of the Prandtl and the Lewis numbers have a profound effect on the stability conditions; opposite stability characteristics are found in salt solutions and in molten metals.

Chen, C. F.↗