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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

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.↗
DOE OSTI2015

A Pareto analysis approach to assess relevant marginal CO{sub 2} footprint for petroleum products

Recently, linear programing (LP) models have been extended to track the marginal CO{sub 2} intensity of automotive fuels at the refinery gate. The obtained CO{sub 2} data are recommended for policy making because they capture the economic and environmental tensions as well as the processing effects related to oil products. However, they are proven to be extremely sensitive to small perturbations and therefore useless in practice. In this paper, we first investigate the theoretical reasons of this drawback. Then, we develop a multiple objective LP framework to assess relevant marginal CO{sub 2} footprints that preserve both defensibility and stability at a satisfactory level of acceptance. A case study illustrates this new methodology. - Highlights: • Refining LP models have limitations to provide useful marginal CO{sub 2} footprints. • A multi objective optimization framework is developed to assess relevant CO{sub 2} data. • Within a European Refinig industry, diesel is more CO{sub 2} intensive than gasoline.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1986

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.↗
NASA NTRS1993

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.↗
DOE OSTI2024

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↗
NASA NTRS1980

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.↗
NASA NTRS1990

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.↗
DOE OSTI2016

Property Analysis of Ethanol--Natural Gasoline--BOB Blends to Make Flex Fuel

Ten natural gasolines were analyzed for a wide range of properties, including Reid vapor pressure (RVP), benzene, sulfur, distillation, stability, metals, and aromatic content, to determine their quality. Benzene and sulfur content were sufficiently low in all but one of the samples that they could be blended without further upgrading. Four of these samples were selected to blend with blendstock for oxygenate blending (BOB) and ethanol to produce E51, E70, and E83 blends, targeting 7.8 and 9.0-psi finished fuels. The volume of each component in the blend was estimated using the Reddy model, with the assumption that the BOB and natural gasoline blend linearly and behave as a single component in the model calculations. Results show that the Reddy model adequately predicts the RVP of the finished blend for E51 and E70, but significantly underpredicts the RVP of E83 blends by nearly 2 psi. It is hypothesized that the underprediction is a function of the very low aromatic content of the E83 blends, even compared to the E51 and E70 blends.

03 NATURAL GAS↗
DOE OSTI2016

Effects of equivalence ratio variations on turbulent flame speed in lean methane/air mixtures under lean-burn natural gas engine operating conditions

Direct numerical simulations of turbulent premixed air flames were carried out using an inflow-outflow configuration to study the effects of equivalence ratio on the turbulent flame speed in lean mixtures. The inflow velocity was dynamically adjusted at a run-time to stabilize the flame brush location within the computational domain. Linear forcing was applied to maintain the turbulent intensities at desired levels. Numerous equivalence ratios near the lean limit were selected for the same turbulence properties and the normalized turbulent flame speed was shown to be a function of the equivalence ratio. Simulations were performed for over 80 eddy turnover times with the turbulent flame speed obtained by averaging the inflow velocity. Results revealed that the equivalence ratio does not have an explicit effect on the normalized turbulent flame speed over the lean limit. Analysis of flame surface area showed that the surface wrinkling produced by eddies of varying scales was not influenced by the change in equivalence ratios when the Karlovitz and Damkohler numbers are fixed. Lastly based on the flame surface statistics flame surface normal is preferentially parallel to the most compressive strain rate direction for all equivalence ratios.

03 NATURAL GAS↗
DOE OSTI2022

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↗
DOE OSTI2022

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↗
NASA NTRS1962

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↗
NASA NTRS1983

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.↗
DOE OSTI2023

Study of Microphonic Effects on the C100 Cryomodule for High Energy Electron Beam Accelerators

The Continuous Electron Beam Accelerator Facility (CEBAF) at Thomas Jefferson National Laboratory (JLab) is a particle accelerator which can accelerate an electron beam to relativistic speeds and apply the beam onto target samples. The C100 superconducting radio frequency (SRF) cavity is the primary accelerating structure of the C100 cryomodule, one of the many cryomodules which compose the CEBAF linear accelerator. SRF cavities are particularly sensitive to internal and external vibrations that can result in a phenomenon called microphonics which degrade the operational stability of a cryomodule. The purpose of this thesis is to investigate the significance of mechanical disturbances on the electromagnetic resonant frequency of a C100 SRF cavity. Knowledge of the mechanical resonance of the cavities and cryomodule sheds light into how these disturbances are most easily realized as deformation which causes radio frequency (RF) detuning. Three studies were conducted: the development and hammer test calibration of a Finite Element Analysis (FEA) model of a C100 cavity, the development and hammer test calibration of an FEA model of a C100 cavity string, and the hammer test of the C100-10R cryomodule at the Cryomodule Test Facility (CMTF). The cavity FEA model was found to accurately predict two modes found in two real cavities in a simply supported configuration. The cavity string FEA model leveraged the calibrated cavity FEA model but was not found to accurately predict the modal behavior of a real cavity string. Even so, the modal behavior of the cavity string inside the C100-10R cryomodule was captured during a hammer test while it was partially assembled. Finally, the C100-10R cryomodule was placed in the CMTF to study RF detuning. The RF detuning spectra during hammer hits and background noise was captured. The results of the hammer testing indicate two strong peaks at low frequencies (9-10 Hz and 22-23 Hz). These two frequencies were found to be nearly coincident to four instances of mechanical resonance found during the hammer testing done on the partially-assembled C100-10R. Because of this, these two modes are believed to contribute to RF detuning of the cryomodule. This test event also included the testing of the effectiveness of a configuration of BNNT canisters designed to act as dampers. While these tests show promising results, the lurking variables render these tests somewhat inconclusive.

Hull, Caleb James↗
DOE OSTI2021

Theoretical analysis of the saturation phase of the 1/1 energetic-ion-driven resistive interchange mode

The aim of the present study is to analyze the saturation regime of the energetic-ion-driven resistive interchange mode (EIC) in the LHD plasma. In this work, a set of nonlinear simulations are performed by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particle (EP) species. The hellically trapped EP component is introduced through a modification of the averaged drift velocity operator to include their precessional drift. The nonlinear simulation results show similar 1/1 EIC saturation phases with respect to the experimental observations, reproducing the enhancement of the n/m = 1/1 resistive interchange modes (RIC) amplitude and width as the EP β increases, the EP β threshold for the 1/1 EIC excitation, the further destabilization of the 1/1 EIC as the population of the helically trapped EP increases and the triggering of burst events. The frequency of the 1/1 EIC calculated during the burst event is 9.4 kHz and the 2/2 and 3/3 overtones are destabilized, consistent with the frequency range and the complex mode structure measured in the experiment. In addition, the simulation shows the inward propagation of the 1/1 EIC due to the nonlinear destabilization of the 3/4 and 2/3 energetic particle modes, leading to the partial overlapping between resonances during the burst event. Finally, the analysis of the 1/1 EIC stabilization phase shows the excitation of the 1/1 RIC as soon as the flattening induced by the 1/1 EIC in the pressure profile vanishes, leading to the retrieval of the pressure gradient at the plasma periphery and the overcoming of the RIC stability limit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Interim Evaluation of State-of-the-Art Controllers for Low-Power Stirling Generators for Space Applications

Stirling convertors use a linear alternator to extract power from the piston motion; this requires a coil and a magnet (Faraday’s law). The coil introduces unwanted phase shift to degrade power output as well as introduce instability; i.e., low power factor. Tuning capacitor is added to increase power output and stability; widely accepted that tuning capacitor value will affect both the power output and stability.

Free-Piston Stirling Convertor Controller↗
NASA NTRS1991

Equilibria, stability and bifurcations of rotating columns of fluid subjected to planar disturbances

Long gyrostatically rotating drops bonded by surface tension are amenable to conventional bifurcation analysis and newer, computer-aided analytical methods, and therefore are useful prototypes of three-dimensional drops. A study is conducted by setting aside instability to Rayleigh's axisymmetric mode and investigating the effects of translationally symmetric (planar) disturbances. The disadvantage of employing single-coordinate representation of drop shapes close to break-up is brought out. It is shown that a family of symmetric two-lobed shapes bifurcates from the main family of perfectly cylindrical shapes when the rotation rate reaches a critical value, in accord with the linearized hydrodynamic analysis of Hocking.

Benner, R. E., Jr.↗
DOE OSTI2019

Leapfrog Time-Stepping for Hermite Methods

We introduce here Hermite-leapfrog methods for first order linear wave systems. The new Hermite-leapfrog methods pair leapfrog time-stepping with the Hermite methods of Goodrich and co-authors et al. (Math Comput 75(254):595–630, 2006). The new schemes stagger field variables in both time and space and are high-order accurate for equations with smooth solutions and coefficients. We provide a detailed description of the method and demonstrate that the method conserves variable quantities. Higher dimensional versions of the method are constructed via tensor products. Numerical evidence and rigorous analysis in one space dimension establish stability and high-order convergence. Experiments demonstrating efficient implementations on a graphics processing unit are also presented.

97 MATHEMATICS AND COMPUTING↗