Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear 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.

At least 1,081 records · Page 60

The ultrahigh pressure stability of silver: An experimental and theoretical study

We measured the atomic volume of Ag in a toroidal diamond anvil cell to a maximum pressure of 416 GPa and calculated the atomic volume and elastic constants of Ag up to 750 and 460 GPa, respectively. Our density functional theory calculations at 0 K utilize an all-electron fully relativistic method and agree well with our volume measurements, particularly at pressures above ~75 GPa. In this work, we corrected our experimental results for non-hydrostaticity using a line shift analysis, and the resulting Vinet equation of state (EOS) parameters are reported. We find that the uniaxial stress sustained by Ag increases linearly up to 4.5 GPa at a pressure of 416 GPa. Our experimental results indicate that the fcc structure of Ag remains stable to at least 416 GPa at room temperature. Our theoretical results show that C44 increases as pressure increases and reaches a maximum at ~100 GPa above which it begins to decrease, a sign that the fcc structure of Ag is becoming unstable, and at V/V 0 = 0.30, the bcc structure is lower in energy than fcc.

36 MATERIALS SCIENCE↗

Multi-scale and multi-technic microstructure analysis of a linear friction weld of the metastable-β titanium alloy Ti-5Al-2Sn-2Zr-4Mo-4Cr (Ti17) towards a new Post-Weld Heat Treatment

Highlights: • Soft joint showing sharp macro-zones and marked by asymmetrical process effects • Gradient of thermomechanical loads identified through progressive stages of transformation • Highly textured recrystallized β grains at the joint core with distorted β lattice • Lattices, phase fractions, microstructures and hardness homogenized after heat treatment - Abstract: Linear Friction Welding (LFW) is a solid-state joining process producing narrow joints mainly developed for the aircraft industry. The Thermo-Mechanical Processing (TMP) loads imposed on the material cause complex microstructure transformations. This study aimed at characterizing the local microstructures and crystallographic configurations to identify the mechanisms impacting the stability and the mechanical properties of a Ti17 LFW joint with the objective of establishing an innovative homogenizing Post-Weld Heat Treatment (PWHT). The process caused a narrow zone of ±2 mm from the friction interface to reach the isothermal α → β transformation temperature under a heating rate of 500 K/s followed by quenching (−100 K/s). The TMP resulted in the formation a ±3 mm wide soft Process Affected Zone (PAZ) delimited by three characteristic zones: the Welding Line (WL width ≈ 400 μm), The Thermo-Mechanically Affected Zone (TMAZ width ≈ 1 mm) and the Heat Affected Zone (HAZ). β TMP resulted in the total dissolution of the α phase in the joint core (WL + TMAZ) and formed strongly textured {110}<111> β recrystallized grains in the WL by continuous dynamic recrystallization of squeezed prior-β grains in the TMAZ. α + β TMP led to a gradual loss of hardness in the HAZ whose microstructure is similar to the one of the base material (BM) but is subject to gradual α laths dissolution due to the local temperatures experienced during LFW. A PWHT consisting in a β annealing followed by an α + β ageing resulted in a uniformly hardened weld with a homogenized Widmanstätten microstructure in the whole assembly.

36 MATERIALS SCIENCE↗

Properties and Autoignition Reactivity of Diesel Boiling Range Ethers Produced from Guerbet Alcohols

We examine the properties of diesel boiling range ethers made from coupling of alcohols produced by oligomerization of ethanol (Guerbet alcohols) for their utility as low-carbon liquid fuel blendstocks. Basic properties of boiling point, flash point, freezing point, density and viscosity are well suited for blending into diesel fuels. For the mixture of ethers the lightest component, di-n-butyl ether, can be present at up to 20 vol% while still having adequately high flashpoint for safe handling. Soot formation tendency (as yield sooting index) is well below that of conventional diesel. The ethers have similar compatibility with elastomers as conventional diesel, based on Hansen solubility parameter analysis. Oxidation stability was assessed for 30 vol% blends of individual ethers in a conventional diesel fuel using a long-term storage test. Over 6 weeks we observed no formation of peroxides or degradation. n-alkyl ethers with carbon number of 8 or higher have cetane number over 100, which is outside the defined range of cetane number, while branched ethers are over 70. The ethers also blend antagonistically into conventional diesel for cetane number, meaning that the blend cetane value is lower than predicted based on a linear by volume, mass, or mole model. We show that aromatics and naphthenes likely act as radical scavengers to slow or shut down autoignition of the highly reactive ethers at low to medium blend levels. Overall, diesel boiling range ethers show significant promise as high quality low-net carbon diesel blendstocks.

09 BIOMASS FUELS↗

Structural and physical properties of Ni 1- x V x alloys around and away from quantum critical point

We investigate the room temperature structure (global and local), temperature dependent magnetic and transport behaviour of Ni 1-x V x (0$&LessSlantEqual;$x$ &LessSlantEqual;$ $0.13$) alloys. Our Energy Dispersive Analysis of x-rays results show that the prepared compositions are stoichiometric. With increase in V doping, the compounds exhibit a quantum phase transition around x c = 0.12, where the ferromagnetic phase is suppressed. Further, our results show that all the compounds stabilize in face centred cubic structure at RT and the lattice parameter shows unusual behaviour close to x c . The magnetic and heat capacity studies show signature of Griffiths phase on either side of x c . From 25 K to the lowest collected temperature, we observe a linear T dependence of resistivity at x = 0.1 and around x c , which is separated by a Fermi-liquid region around x = 0.106. This suggests that the origin of the transport behaviour is different around the quantum critical point and away from it. Our Ni K-edge x-ray Absorption Spectroscopy results show that there is a significant reduction in the first coordination number around Ni central atom on doping. Further, with doping, there is distortion in the first coordination shell around Ni. This suggests, with V doping, the local structure around Ni is different from the global structure as obtained from the x-ray Diffraction results. Interestingly, with doping, we observe a direct connection between the extent of distortion at RT and the magnetic disorder obtained at 2 K. We believe our results will motivate the scientific community to further study the interplay between the structural disorder and quantum fluctuations with temperature at the local level.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermal oxidative degradation reactions of linear perfluoroalky lethers

Thermal and thermal oxidative stability studies were performed on linear perfluoro alkyl ether fluids. The effect on degradation by metal catalysts and degradation inhibitors are reported. The liner perfluoro alkylethers are inherently unstable at 316 C in an oxidizing atmosphere. The metal catalysts greatly increased the rate of degradation in oxidizing atmospheres. In the presence of these metals in an oxidizing atmosphere, the degradation inhibitors were highly effective in arresting degradation at 288 C. However, the inhibitors had only limited effectiveness at 316 C. The metals promote degradation by chain scission. Based on elemental analysis and oxygen consumption data, the linear perfluoro alkylether fluids have a structural arrangement based on difluoroformyl and tetrafluoroethylene oxide units, with the former predominating.

Jones, W. R., Jr.↗

Cosmic ray secular variations in terrestrial records and aurorae

The rediscovery that the Sun and the solar wind can undergo important changes on historical time scales has brought into question the stability of the cyclic behavior of past time series of solar and solar-terrestrial origin. It was found by Vector Fourier analysis that the solar 11 year cycle is present in the series of 10Be, delta 180, in ice cores and of thermoluminescence (TL) in sea sediments during the last Millennia with a frequency modulation, related to the Sun behavior, as tested by comparison with the Sunspot number R sub z series. It was shown that the cyclogram of the series of yearly Aurorae from 1721 to 1979 linear-regression-corrected-for-R sub z is straight for the periodicity zeta=11,1y, which indicates that such periodicity is constant in time corresponding to the only line present in the 11y band. The maxima of this component appear at the same time together with the high speed solar wind streams taking place in coronal holes situated in high heliolatitudes. It is evidenced that the 11 year cycle has undergone frequency oscillations on a time scale of two centuries, although it is very difficult to determine the periodicities with high accuracy.

Attolini, M. R.↗

Stochastic stability and performance robustness of linear multivariable systems

Stochastic robustness, a simple technique used to estimate the robustness of linear, time-invariant systems, is applied to a single-link robot arm control system. Concepts behind stochastic stability robustness are extended to systems with estimators and to stochastic performance robustness. Stochastic performance robustness measures based on classical design specifications are introduced, and the relationship between stochastic robustness measures and control system design parameters is discussed. The application of stochastic performance robustness and the relationship between performance objectives and design parameters are demonstrated by means of the example. The results prove stochastic robustness to be a good overall robustness analysis method that can relate robustness characteristics to control system design parameters.

Stengel, Robert F.↗

Reduced conservatism in stability robustness bounds by state transformation

This note addresses the issue of 'conservatism' in the time domain stability robustness bounds obtained by the Liapunov approach. A state transformation is employed to improve the upper bounds on the linear time-varying perturbation of an asymptotically stable linear time-invariant system for robust stability. This improvement is due to the variance of the conservatism of the Liapunov approach with respect to the basis of the vector space in which the Liapunov function is constructed. Improved bounds are obtained, using a transformation, on elemental and vector norms of perturbations (i.e., structured perturbations) as well as on a matrix norm of perturbations (i.e., unstructured perturbations). For the case of a diagonal transformation, an algorithm is proposed to find the 'optimal' transformation. Several examples are presented to illustrate the proposed analysis.

Yedavalli, R. K.↗

On the stability of a time dependent boundary layer

The aim of this article is to determine the stability characteristics of a Rayleigh layer, which is known to occur when the fluid above a flat plate has a velocity imparted to it (parallel to the plate). This situation is intrinsically unsteady, however, as a first approximation we consider the instantaneous stability of the flow. The Orr-Sommerfeld equation is found to govern fixed downstream wavelength linear perturbations to the basic flow profile. By the solution of this equation, we can determine the Reynolds numbers at which the flow is neutrally stable; this quasisteady approach is only formally applicable for infinite Reynolds numbers. We shall consider the large Reynolds number limit of the original problem and use a three deck mentality to determine the form of the modes. The results of the two calculations are compared, and the linear large Reynolds number analysis is extended to consider the effect of weak nonlinearity in order to determine whether the system is subcritical or supercritical.

Otto, S. R.↗

Nonlinear rotordynamics analysis

Three analytic consequences of the nonlinear Jeffcott equations are examined. The primary application of these analyses is directed toward understanding the excessive vibrations recorded in the Liquid Oxygen (LOX) pump of the Space Shuttle Main Engine (SSME) during hot firing ground testing. The first task is to provide bounds on the coefficients of the equations which delimit the two cases of numerical solution as a circle or an annulus. The second task examines the mathematical generalization to multiple forcing functions, which includes the special problems of mass imbalance, side force, rubbing, and combination of these forces. Finally, stability and boundedness of the steady-state solutions is discussed and related to the corresponding linear problem.

Day, W. B.↗

A detailed study of mean-flow solutions for stability analysis of transitional flows

A finite-volume upwind-difference parabolized Navier-Stokes code is utilized to obtain laminar mean-flow solutions at Mach 3.5 on a half-angle cone of 5 deg at an angle-of-attack of 2 deg. A detailed study is conducted on this configuration; the main focus is the velocity profiles in the leeward and windward symmetry planes at various axial locations. Comparisons of the solution profiles are made with both a central-difference code that incorporates scalar and matrix dissipation models and another state-of-the-art upwind-difference finitevolume code. The results obtained emphasize the importance of using matrix dissipation models for schemes that require explicit artificial dissipation. These results also illustrate the accuracy and efficiency of the planeby-plane marching procedure for computing mean-flow solutions for predicting the onset of transition with linear instability.

Ramakrishnan, R.↗

Magnetic field-temperature phase diagram of the spin-$\frac{1}{2}$ triangular lattice antiferromagnet KYbSe 2

A quantum spin liquid (QSL) is a state of matter characterized by fractionalized quasiparticle excitations, quantum entanglement, and a lack of long-range magnetic order. However, QSLs have evaded definitive experimental observation. Several Yb 3+ -based triangular lattice antiferromagnets with effective 𝑆 = $\frac{1}{2}$ have been suggested to stabilize the QSL state as the ground state. Here, in this work, we build a comprehensive magnetic temperature phase diagram of a high-quality single crystalline KYbSe 2 via heat capacity and magnetocaloric effect down to 30 mK with magnetic field applied along the 𝑎 axis. At zero magnetic field, we observe the magnetic long-range order at 𝑇 N =0.29 K entering 120 degrees ordered state in heat capacity, consistent with neutron scattering studies. Analysis of the low-temperature (𝑇) specific heat (𝐶) at zero magnetic field indicates linear 𝑇 dependence of 𝐶/𝑇 and a broad hump of 𝐶/𝑇 in the proximate QSL region above 𝑇 N . By applying magnetic field, we observe the up-up-down phase with 1/3 magnetization plateau and oblique phases, in addition to two new phases. These observations strongly indicate that while KYbSe 2 closely exhibits characteristics resembling an ideal triangular lattice, deviations may exist, such as the effect of the next-nearest-neighbor exchange interaction, calling for careful consideration for spin Hamiltonian modeling. Further investigations into tuning parameters, such as chemical pressure, could potentially induce an intriguing QSL phase in the material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Eddy heat fluxes and stability of planetary waves. I, II

The stability of baroclinic Rossby waves in a zonal shear flow was analyzed by a linear, quasigeostrophic, two-level, adiabatic, and frictionless midlatitude beta-plane model. The ratio of the basic wave scale and the radius of deformation together with two nondimensional parameters which describe the amplitudes of the barotropic and baroclinic components of the basic wave constitute the three parameters of the stability problem. The parameter space is partitioned according to the dominant energy source for instability; the Lorenz and Kim conditions are characterized by significant horizontal and vertical shears of the basic wave, while the Phillips regime has a strong zonal flow. The stability analysis is then applied to the atmosphere, with the primary motivation being to examine the midlatitude planetary scale (zonal wavenumbers 1, 2, 3) transient waves that transport heat. It is found that the most unstable mode consists of a spectrum of waves, with a maximum amplitude at wavenumber 3; the response is thus maximum at a zonal scale intermediate between the basic wave scale and the radius of deformation.

Lin, C. A.↗

On the instability of Goertler vortices to nonlinear travelling waves

Recent theoretical work by Hall and Seddougui (1989) has shown that strongly nonlinear, high wavenumber Goertler vortices developing within a boundary layer flow are susceptible to a secondary instability which takes the form of travelling waves confined to a thin region centered at the outer edge of the vortex. The case is considered in which the secondary mode could be satisfactorily described by a linear stability theory and herein the objective is to extend this investigation of Hall and Seddougui (1989) into the nonlinear regime. It was found that at this stage not only does the secondary mode become nonlinear but it also interacts with itself so as to modify the governing equations for the primary Goertler vortex. In this case then, the vortex and the travelling wave drive each other and, indeed, the whole flow structure is described by an infinite set of coupled, nonlinear differential equations. A Stuart-Watson type of weakly nonlinear analysis of these equations is undertaken and concluded, in particular, that on this basis there exist stable flow configurations in which the travelling mode is of finite amplitude. Implications of the findings for practical situations are discussed and it is shown that the theoretical conclusions drawn here are in good qualitative agreement with available experimental observations.

Seddougui, Sharon O.↗

High throughput exploration of the oxidation landscape in high entropy alloys

We report high entropy alloys (HEAs) have gained interest for structural applications in extreme environments. With a potentially vast chemical and phase space, there are significant opportunities to discover superior performing alloys. Crucial for most high-temperature applications is understanding and mitigating the oxidation behavior of these chemically complex alloys. Most experimental and computational HEA studies have focused on a limited set of compositions and only a fraction these compositions have been characterized for oxidation. We present a high-throughput framework that utilizes density-functional theory (DFT) in concert with a combined machine-learning model and grand-canonical linear programming for assessing phase stability, phase-fraction, chemical activity and high-temperature survivability of arbitrary HEAs. This framework considers temperature dependent contributions to the Gibbs energy of the competing phases arising from short-range order and vibrational entropy. We demonstrate the effectiveness of the framework by assessing the thermodynamic stability, oxidation behavior, chemical activity, and phase decomposition of body-centered cubic Mo-W-Ta-Ti-Zr refractory HEAs. A total of 51 compositions were analyzed and ranked in order of their survivability based on the Pareto-front analysis. Oxidation was performed at 1373 K on fours samples in air showing the difference in oxidation behavior determined experimentally through scale thickness and their mass changes. The insights on oxidation behavior presented in this work will enable the fast assessment of technologically useful HEAs needed for future structural application in extreme conditions

36 MATERIALS SCIENCE↗

Toroidal modeling of energetic passing particle drift kinetic effects on tearing mode stability

Abstract Drift kinetic effects of the neutral beam injection induced passing energetic particles (EPs) on the linear stability of the n = 1 tearing mode (TM) (with the dominant poloidal harmonic of m = 2) are numerically investigated utilizing the MARS-K code (Liu et al 2008 Phys. Plasmas 15 112503), in a tokamak plasma with finite equilibrium pressure and anisotropic thermal transport. In the low plasma pressure regime, it is found that co- (counter-) passing EPs stabilize (destabilize) the TM, agreeing with previous studies. However, as the plasma pressure increases beyond a critical value, it is found that co-passing EPs also destabilize the mode. An in-depth analysis reveals that the net effect of co-passing EPs is a result of competition between the stabilizing contribution from the non-adiabatic drift kinetic terms and the destabilizing contribution associated with adiabatic terms, with the latter becoming more dominant at higher equilibrium pressure. Non-perturbative magnetohydrodynamic-kinetic hybrid modeling also finds that co- and counter-passing EPs modify the TM eigenfunction differently, with the counter-passing EPs enhancing the sideband harmonics. Furthermore, effects of the plasma resistivity and toroidal rotation, as well as that of the equilibrium distribution of EPs in the particle pitch angle space, are also investigated, showing asymmetric results on the TM stability between the co- and counter-passing EPs. The first order finite orbit width correction is found to be stabilizing with co-passing EPs and destabilizing with counter-passing particles. Finally, drift resonances between passing EPs and the TM induce finite frequency to the mode and generate finite net torques inside the plasma, due to the neoclassical toroidal viscosity and the Reynolds stress associated with 3D perturbations.

Physics↗

LIFT CONTROL DURING ATMOSPHERE ENTRY FROM SUPERCIRCULAR VELOCITY

This paper presents an analysis of a guidance method which uses a reference trajectory. The four state variables needed to prescribe the trajectory are used as follows : Velocity is made the independent variable, and the errors in the rate-of-climb, acceleration, and range variables away from the reference are used to govern the lift. A linearized form of the motion equations is used to show that this represents a third-order control system. First- and second- order control terms (rate of climb and acceleration inputs) are shown to determine the entry corridor depth by stabilizing the trajectory so that the vehicle does not skip back out of the atmosphere or does not exceed a specified acceleration limit. The destabilizing effect that range input (the third-order control term) can have is illustrated and the results indicate that a low value of range input gain must be used at the high supercircular velocities while larger values of range input gain can be used at lower velocities. The usable corridor depth and range capability with this guidance system are demonstrated for a lifting capsule (L/D = 0.5). The practical applications of this system are illustrated with a fixed trim configuration wherein roll angle is used to command the desired lift. The results show that the guidance system requires only one reference trajectory for abort entry conditions as well as for entry conditions near the design values.

Atmosphere↗

Transformation of two and three-dimensional regions by elliptic systems

Finite difference methods for composite grids were analyzed. It was observed that linear interpolation between grids would suffice only where low order accuracy was required. In the context of fluid flow, this would be in regions where the flow was essentially free stream. Higher order interpolation schemes were also investigated. The well known quadratic and cubic interpolating polynomials would increase the formal accuracy of the overall numerical algorithm. However, it can also be shown that the stability of the algorithm may be adversely affected. Further numerical results are needed in order to assess the nature of this instability induced by the interpolation procedure. Finally, error analysis and the order of difference expressions on general curvilinear coordinates are discussed.

Thompson, J. F.↗