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.

SourceResearch connectionFrom yearThrough year

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

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Suppression of pressure oscillations in an open cavity by passive pneumatic control

A computational investigation has been conducted to determine the effectiveness of a passive control technique in suppressing the oscillations in an open cavity exposed to a supersonic flow. Time-accurate solutions of the unsteady, Reynolds-averaged, Navier-Stokes equations were obtained with an explicit predictor-corrector algorithm; the passive control was implemented through the use of a simple linear pressure-velocity law along the porous cavity floor. The computational code was validated by comparisons with experimental data for the cavity flow without control. The computational results with control demonstrate that the oscillations are suppressed; the fluid dynamic mechanism of the control is seen to be a stabilization of the free-shear layer. Spectral analysis of the unsteady data shows that the resonant frequencies are essentially unchanged by the control. Beneficial reductions in the cavity resonant drag are also observed.

Chokani, N.↗

NASA NTRS1950

Aerodynamic Characteristics of NACA RM-10 Missile in 8- by 6-foot Supersonic Wind Tunnel at Mach Numbers from 1.49 to 1.98 II : Presentation and Analysis of Force Measurements

Experimental investigation of aerodynamic forces acting on body of revolution (NACA RM-10 missile) with and without stabilizing fins was conducted at Mach numbers from 1.49 to 1.98 at angles of attack from 0 to 9 degrees and at Reynolds number of approximately 30,000,000. Comparison of experimental lift, drag, and pitching-moment coefficients and center of pressure location for body alone is made with linearized potential theory and a semiempirical method. Results indicate that aerodynamic characteristics were predicted more accurately by semiempirical method than by potential theory. Breakdown of measured drag coefficients into components of friction, pressure, and base-pressure drag is presented for body alone at zero angle of attack.

TAIL-BODY COMBINATIONS - MISSILES↗

NASA NTRS1993

The onset of electrohydrodynamic instability in isoelectric focusing

The onset of electrohydrodynamic motion associated with the imposition of an electric field across a thin layer of liquid has been investigated for the case in which the electrical conductivity varies linearly over the depth of the layer. The variation of the conductivity is due to concentration gradients in the charge-carrying solutes and its spatio-temporal evolution is represented by a convective-diffusion equation. When the viscous relaxation time is long compared to the time for charge relaxation, the analysis reveals that the neutral stability curves for the layer can be characterized by three dimensionless parameters: Ra(sub e) is equivalent to d(epsilon)E(sup 2, sub 0) Delta sigma/mu K(sub eff) sigma(sub 0), an electrical Rayleigh number; delta sigma/sigma(sub 0), the relative conductivity increment; and alpha, the transverse wave number of the disturbance. Here d is the thickness, epsilon is the dielectric constant, and mu is the viscosity of layer, E(sub 0) is the applied field strength at the lower conductivity boundary, and K(sub eff) is an effective diffusivity associated with the Brownian motion of the charge-carrying solutes. With viscous-stress-free boundaries, at which the electrical conductivity and the normal component of the electric field are prescribed, the critical Ra(sub e) is 1.504 x 10(exp 4) at a critical transverse wave number of 1.97 when Delta sigma/sigma(sub 0) is 10. As Delta sigma/sigma(sub 0) increases, the critical Ra(sub e) increases and shifts to shorter wavelength disturbances; the critical imposed field strength, however, passes through a minimum because the lower-conductivity boundary exerts a considerable stabilizing influence in the presence of steep conductivity gradients. Similar trends were obtained for liquid layers with rigid boundaries.

Baygents, James C.↗

NASA NTRS2010

Use of ILTV Control Laws for LaNCETS Flight Research

A report discusses the Lift and Nozzle Change Effects on Tail Shock (LaNCETS) test to investigate the effects of lift distribution and nozzle-area ratio changes on tail shock strength of an F-15 aircraft. Specific research objectives are to obtain inflight shock strength for multiple combinations of nozzle-area ratio and lift distribution; compare results with preflight prediction tools; and update predictive tools with flight results. The objectives from a stability and control perspective are to ensure adequate aircraft stability for the changes in lift distribution and plume shape, and ensure manageable transient from engaging and disengaging the ILTV research control laws. In order to change the lift distribution and plume shape of the F-15 aircraft, a decade-old Inner Loop Thrust Vectoring (ILTV) research control law was used. Flight envelope expansion was performed for the test configuration and flight conditions prior to the probing test points. The approach for achieving the research objectives was to utilize the unique capabilities of NASA's NF-15B-837 aircraft to allow the adjustment of the nozzle-area ratio and/or canard positions by engaging the ILTV research control laws. The ILTV control laws provide the ability to add trim command biases to canard positions, nozzle area ratios, and thrust vectoring through the use of datasets. Datasets consist of programmed test inputs (PTIs) that define trims to change the nozzle-area ratio and/or canard positions. The trims are applied as increments to the normally commanded positions. A LaNCETS non-linear, six-degrees-of-freedom simulation capable of realtime pilot-in-the-loop, hardware-in-the-loop, and non-real-time batch support was developed and validated. Prior to first flight, extensive simulation analyses were performed to show adequate stability margins with the changes in lift distribution and plume shape. Additionally, engagement/disengagement transient analysis was also performed to show manageable transients.

Moua, Cheng↗

DOE OSTI2023

Rapid Polyolefin Plastic Hydrogenolysis Mediated by Single-Site Heterogeneous Electrophilic/Cationic Organo-group IV Catalysts

A homologous series of cationic electrophilic group IV metal hydrides (M = Ti, Zr, Hf) created by chemisorption of the corresponding MNp 4 precursors on highly Brønsted acidic sulfated alumina (AlS) to yield loosely coordinated surface AlS/MNp 2 (Np = neopentyl) species are systematically characterized by ICP, EXAFS/XANES, DRIFTS, and solid-state high-resolution multidimensional NMR spectroscopy (SSNMR), as well as by energy span DFT computation. With effective stirring, these complexes readily undergo reaction with H 2 to yield AlS/M(alkyl)H species, which are highly active for the hydrogenolysis of diverse commercial polyethylenes, α-olefin-ethylene copolymers, isotactic polypropylene, and postconsumer polyolefins including high-density polyethylenes, yielding medium and small linear and branched hydrocarbons at turnover frequencies as high as 36,300 h –1 at 200 °C/17 atm H 2 for M = Zr. For given polyolefin and reaction conditions, turnover frequencies scale approximately as M = Zr > Hf > Ti, whereas catalyst thermal stability scales approximately as M = Hf ≈ Zr > Ti, and these trends are qualitatively understandable from the DFT analysis. These catalytic results reveal that the AlS/Hf(R)H-mediated hydrogenolysis favors wax-like and liquid products, whereas the AlS/Zr(R)H-mediated hydrogenolysis can be tuned between gases and liquids. Finally, DFT analysis identifies β-alkyl elimination as the turnover-limiting C–C scission process, which is particularly facile in these cationic d 0 complexes but not so in the neutrally charged analogues.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Predictions of core plasma performance for the SPARC tokamak

SPARC is designed to be a high-field, medium-size tokamak aimed at achieving net energy gain with Ion Cyclotron Range-of-Frequencies (ICRF) as its primary auxiliary heating mechanism. Empirical predictions with conservative physics indicate that SPARC baseline plasmas would reach Q~11, well above its mission objective of Q>2. To build confidence that SPARC will be successful, physics-based integrated modeling has also been performed. The TRANSP code coupled with the theory-based TGLF turbulence model and EPED predictions for pedestal stability find that Q~9 is attainable in standard H-mode operation and confirms Q>2 operation is feasible even with adverse assumptions. In this analysis, ion cyclotron waves are simulated with the full wave TORIC code and alpha heating is modeled with the Monte-Carlo fast ion NUBEAM module. Detailed analysis of expected turbulence regimes with linear and nonlinear CGYRO simulations is also presented, demonstrating that profile predictions with the TGLF reduced model are in reasonable agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1984

Effects of nozzle design parameters on the extent of quiet test flow at Mach 3.5

Tests results at the NASA Langley Research Center, involving a Mach 3.5 pilot quiet tunnel, have shown that laminar-layered nozzle walls improve boundary layer stability and reduce stream disturbance levels caused by eddy Mach wave radiation. This type of wall design is required to obtain transition Reynolds numbers on tests models as high as those previously observed in supersonic flight vehicles. The Mach 3.5 pilot nozzle wall boundary layers were tested for Tollmein-Schlichting and Goertler linear amplification, and, in an analysis of Goertler vortices in two axisymmetric Mach 5 nozzles, transition values were found to vary. These values were applied to several nozzles with similar throat heights but different expansion rates. Among the nozzles included in the study, a flat-wall radial flow nozzle and a proposed rod-wall nozzle were tested. For the highest test unit Reynolds number, it was determined that the nozzle wall surface finish should not exceed 0.3 micron. Oil flow studies have indicated that Goertler vortex disturbances were the dominant mechanism causing transition on the walls of the pilot nozzle.

Beckwith, I. E.↗

NASA NTRS2020

Streak Instabilities on HIFiRE-5 Elliptic Cone

The Hypersonic International Flight Research and Experimentation (HIFiRE) 5 flight experiment by the Air Force Research Laboratories and Australian Defense Science and Technology Organization yielded inflight, hypersonic transition data for a canonical, 3D configuration in the form of an elliptic cone at various combinations of angle-of-attack and yaw. Previous analysis has shown that boundary-layer stability calculations for the as-flown, zero-yaw flight conditions are able to capture the second mode dominated transition along the attachment line as well as the crossflow induced transition in between the attachment line and the minor axis symmetry plane. However, there is still a need to understand the transition in the vicinity of the minor axis plane of symmetry, where the boundary layer thickness exhibits an azimuthally compact bulge due to the collision of the secondary flow approaching from the attachment lines on both sides. This paper reports a computational analysis of the stability of the boundary layer flow along the bulge regions on both the leeward and windward sides of the vehicle surface. At a representative flight condition corresponding to a time instant of t = 514.813 seconds, computations are performed to identify several different unstable modes within the bulge regions, with the eventual goal of developing a linear stability correlation for the transition onset near the minor axis symmetry plane. The most amplified flow instabilities in the pretransitional region correspond to a “stem” mode that is concentrated within a pair of nearly vertical shear layers below the cap of the mushroom shaped bulge in the mean velocity contours. An excellent agreement is found between the predictions of instability amplification based on the plane marching parabolized stability equations (PSE) and the direct numerical solutions (DNS). Given the large number of unstable modes within the bulge region, the observed agreement demonstrates, for first time, that the PSE can provide accurate predictions for the development of flow instabilities in streak-like flow features over three-dimensional bodies. The peak N-factor of approximately 15 at the transition location estimated from the thermocouple measurements along the leeward centerline is close to the transition N-factor for Mack mode transition during the HIFiRE-1 circular cone flight experiment with nearly zero degrees angle of attack.

Meelan Choudhari↗

NASA NTRSDate not supplied

Electromagnetic Modeling of Wind tunnel Magnetic Suspension and Balance Systems

An analytical framework for the open-loop behavior of a permanent magnet element levitated within an applied magnetic field are shown, in the context of application to wind tunnel Magnetic Suspension and Balance Systems (MSBS). Various modes of motion are identified, which correlate with observed behavior of the NASA/ODU 6-inch MSBS. This system is being developed as a technology demonstrator for dynamic stability testing of atmospheric entry capsules, with possible application to a supersonic wind tunnel. The analytic equations are populated with field information from experimental measurements and a finite element model. Finally, a system simulation can provide alternative estimates of system dynamic behavior via model linearization, with preliminary comparisons presented. Taken together, the analytical framework, experimental measurement, finite element analysis, and dynamic simulation provide a complete understanding of the characteristics and behavior of the MSBS.

magnetic suspension↗

DOE OSTI2020

Predictions of core plasma performance for the SPARC tokamak

SPARC is designed to be a high-field, medium-size tokamak aimed at achieving net energy gain with ion cyclotron range-of-frequencies (ICRF) as its primary auxiliary heating mechanism. Empirical predictions with conservative physics indicate that SPARC baseline plasmas would reach $\textit{Q}$ ≈ 11, which is well above its mission objective of $\textit{Q}$ > 2. To build confidence that SPARC will be successful, physics-based integrated modelling has also been performed. The TRANSP code coupled with the theory-based trapped gyro-Landau fluid (TGLF) turbulence model and EPED predictions for pedestal stability find that $\textit{Q}$ ≈ 9 is attainable in standard H-mode operation and confirms $\textit{Q}$ > 2 operation is feasible even with adverse assumptions. In this analysis, ion cyclotron waves are simulated with the full wave TORIC code and alpha heating is modelled with the Monte–Carlo fast ion NUBEAM module. Detailed analysis of expected turbulence regimes with linear and nonlinear CGYRO simulations is also presented, demonstrating that profile predictions with the TGLF reduced model are in reasonable agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Effect of polycyclosilane microstructure on thermal properties

Thermal characterization of polysilanes has focused on the influence of organic side chains, whereas little is understood about the influence of silane backbone microstructure on thermal stability, phase properties, and pyrolysis. To address this knowledge gap, we prepared three distinct polycyclosilanes: linear polymers synthesized from the cyclosilane building blocks 1,4Si 6 and 1,3Si 6 , as well as a cyclic polymer of 1,3Si 6 . Thermal properties across the temperature range 25 to 600 °C were investigated using differential scanning calromietry (DSC) and thermogravimetric analysis (TGA). We found differences between linear and cyclic materials, including a phase transition unique to the cyclic polymer and lower rates of mass loss during pyrolysis. Density functional theory (DFT) calculations provided insight into microstructure-dependent pyrolysis.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Undistorted linear Bi chains with hypervalent bonding in La 3 TiBi 5 from single-crystal X-ray diffraction

The crystal structure of the lanthanum titanium bismuthide La 3 TiBi 5 (Pearson code hP18, Wyckoff sequence b d g2) has been established from single-crystal X-ray diffraction data and analyzed in detail using first-principles calculations. There are no anomalies pertaining to the atomic displacement parameter of the Ti site, previously reported based on a powder X-ray diffraction analysis of this compound. The anionic substructure contains columns of face-sharing TiBi 6 octahedra and linear Bi chains. Due to a significant La(5d) and Bi(6p) orbital mixing, a perfectly one-dimensional character of the Bi chains is not realised, while a three-dimensional electronic structure is established instead. The latter fact explains the stability of the polyanionic pnictide units against Peierls distortions. The hypervalent bonding in the Bi chains is reflected in a rather long Bi—Bi distance of 3.2264 (4) Å and a typical pattern of bonding and antibonding interactions, as revealed by electronic structure calculations.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Estimating and Calibrating DER Model Parameters Using Levenberg–Marquardt Algorithm in Renewable Rich Power Grid

The proliferation of inverter-based distributed energy resources (IBDERs) has increased the number of control variables and dynamic interactions, leading to new grid control challenges. For stability analysis and designing appropriate protection controls, it is important that IBDER models are accurate. This paper focuses on the accurate estimation and parameter calibration of DER_A, a recently proposed aggregated IBDER model. In particular, we focus on the parameters of the reactive power–voltage regulation module. We formulate the problem of parameter tuning as a non-linear least square minimization problem and solve it using the Levenberg–Marquardt (LM) method. The LM method is primarily chosen due to its flexibility in adaptively selecting between the steepest descent and Gauss–Newton methods through a damping parameter. The LM approach is used to minimize the error between the actual measurements and the estimated response of the model. Further, the computational challenges posed by the numerical calculation of the Jacobian are tackled using a quasi-Newton root-finding approach. The proposed method is validated on a real feeder model in the northeastern part of the United States. The feeder is modeled in OpenDSS and the measurements thus obtained are fed to the DER_A model for calibration. The simulation results indicate that our approach is able to successfully calibrate the relevant model parameters quickly and with high accuracy, with a total sum of square error of 3.57 × 10 –7 .

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

NASA NTRS1994

The buckling response of symmetrically laminated composite plates having a trapezoidal planform area

The focus of this work is the buckling response of symmetrically laminated composite plates having a planform area in the shape of an isosceles trapezoid. The loading is assumed to be inplane and applied perpendicular to the parallel ends of the plate. The tapered edges of the plate are assumed to have simply supported boundary conditions, while the parallel ends are assumed to have either simply supported or clamped boundary conditions. A semi-analytic closed-form solution based on energy principles and the Trefftz stability criterion is derived and solutions are obtained using the Rayleigh-Ritz method. Intrinsic in this solution is a simplified prebuckling analysis which approximates the inplane force resultant distributions by the forms Nx=P/W(x) and Ny=Nxy=0, where P is the applied load and W(x) is the plate width which, for the trapezoidal planform, varies linearly with the lengthwise coordinate x. The out-of-plane displacement is approximated by a double trigonometric series. This analysis is posed in terms of four nondimensional parameters representing orthotropic and anisotropic material properties, and two nondimensional parameters representing geometric properties. For comparison purposes, a number of specific plate geometry, ply orientation, and stacking sequence combinations are investigated using the general purpose finite element code ABAQUS. Comparison of buckling coefficients calculated using the semi-analytical model and the finite element model show agreement within 5 percent, in general, and within 15 percent for the worst cases. In order to verify both the finite element and semi-analytical analyses, buckling loads are measured for graphite/epoxy plates having a wide range of plate geometries and stacking sequences. Test fixtures, instrumentation system, and experimental technique are described. Experimental results for the buckling load, the buckled mode shape, and the prebuckling plate stiffness are presented and show good agreement with the analytical results regarding the buckling load and the prebuckling plate stiffness. However, the experimental results show that for some cases the analysis underpredicts the number of halfwaves in the buckled mode shape. In the context of the definitions of taper ratio and aspect ratio used in this study, it is concluded that the buckling load always increases as taper ratio increases for a given aspect ratio for plates having simply supported boundary conditions on the parallel ends. There are combinations of plate geometry and ply stackling sequences, however, that reverse this trend for plates having clamped boundary conditions on the parallel ends such that an increase in the taper ratio causes a decrease in the buckling load. The clamped boundary conditions on the parallel ends of the plate are shown to increase the buckling load compared to simply supported boundary conditions. Also, anisotropy (the D16 and D26 terms) is shown to decrease the buckling load and skew the buckled mode shape for both the simply supported and clamped boundary conditions.

Radloff, H. D., II↗

NASA NTRSDate not supplied

Boundary-Layer Transition Prediction Through Loose Coupling of OVERFLOW and LASTRAC

Transition prediction based on linear stability theory is expected to more accurately reflect the causality of transition onset than phenomenological transition models based on RANS-like transport equations. To help achieve the CFD vision 2030 aim of building a CFD tool chain with automated prediction of boundary layer transition, a technique to loosely tie the NASA OVERFLOW CFD solver with the LASTRAC stability analysis tool is described. The coupled solver is then used to compute transition over selected over a flat plate in a freestream with sufficiently low levels of turbulence, NLF(1)-0416 airfoil, the 6:1 prolate spheroid at an angle of attack, and a NASA juncture flow model with symmetric wing configuration. The findings show that the loosely coupled approach can reliably predict the transition location accurately in scenarios that are dominated by a single transition mechanism involving Tollmien-Schlichting instabilities, crossflow instabilities, or separation bubble-induced transition, or include a mixture of selected mechanisms. The toolset presents here appears to be robust to the prescription of the initial transition location, and it can lead to a converged solution in four or five rounds of the mean flow calculation and stability analysis, with minimal input from the user.

boundary layer transition↗

NASA NTRS1988

Rotordynamic analysis of a bearing tester

The properties of the solutions of a system of four coupled nonlinear differential equations that model the behavior of the rotating shaft of a bearing tester are studied. In particular, it is shown how the bounds for the rotations of these equations can be obtained from bounds for the solutions of the linearized equations. By studying the behavior of the Fourier transforms of the solution, the approach to the stability boundary can also be predicted. These conclusions are verified by means of numerical solutions of the equations, and of power spectrum density (PSD) plots.

Zalik, Richard A.↗

NASA NTRS1990

The effect of an electric field on the morphological stability of the crystal-melt interface of a binary alloy. III - Weakly nonlinear theory

The effect of a constant electric current on the crystal-melt interface morphology during directional solidification at constant velocity of a binary alloy is considered. A linear temperature field is assumed, and thermoelectric effects and Joule heating are neglected; electromigration and differing electrical conductivities of crystal and melt are taken into account. A two-dimensional weakly nonlinear analysis is carried out to third order in the interface amplitude, resulting in a cubic amplitude equation that describes whether the bifurcation from the planar state is supercritical or subcritical. For wavelengths corresponding to the most dangerous mode of linear theory, the demarcation between supercritical and subcritical behavior is calculated as a function of processing conditions and material parameters. The bifurcation behavior is a sensitive function of the magnitude and direction of the electric current and of the electrical conductivity ratio.

Wheeler, A. A.↗

DOE OSTI2018

Semi-numerical solution to a fractal telegraphic dual-porosity fluid flow model

We present a semi-numerical solution to a fractal telegraphic dual-porosity fluid flow model. It combines Laplace transform and finite difference schemes. The Laplace transform handles the time variable whereas the finite difference method deals with the spatial coordinate. This semi-numerical scheme is not restricted by space discretization and allows the computation of a solution at any time without compromising numerical stability or the mass conservation principle. Our formulation results in a non-analytically-solvable second-order differential equation whose spatial discretization yields a tridiagonal linear algebraic system. Moreover, we describe comparisons between semi-numerical and semi-analytical solutions for particular cases. Results agree well with those from semi-analytic solutions. Furthermore, we expose a parametric analysis from the coupled model in order to show the effects of relevant parameters on pressure profiles and flow rates for the case where neither analytic nor semi-analytic solutions are available.

97 MATHEMATICS AND COMPUTING↗