Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical stability”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

NASA Langley FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

This work summarizes the results of FUN3D analyses conducted for the 1st AIAA Stability and Control Workshop on behalf of participants from the NASA Langley Research Center. The workshop was created to establish best practices for the prediction of stability and control derivatives using computational fluid dynamics and assess the limitations of these methods when those best practices are applied. The inaugural workshop considered the ONERA version of the NASA/Boeing Common Research Model, which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel data at small sideslip angles remain unpublished and served as ‘blind’ data for computational comparisons. The present research generated workshop test case data using the NASA FUN3D code, which is a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Steady- state numerical simulations were conducted for workshop test cases investigating the following: grid convergence, Mach number effect on static stability, wind tunnel sting increments, static stability-derivative calculations, and a sideslip angle sweep. Results were generated for two series of unstructured, mixed-element grids, one set provided by the workshop and another set created using the HeldenMesh grid generation software. The results provided include total- and component-level breakdowns of the force and moment coefficients, in addition to sectional pressure distributions for the wing and tail components for comparisons to wind tunnel data.

CFD↗
NASA NTRS1979

Analysis of solar radiation pressure induced coupled librations of a gravity stabilized axisymmetric satellite

This paper presents an analysis of solar radiation pressure induced coupled librations of gravity stabilized cylindrical spacecraft with a special reference to geostationary communication satellites. The Lagrangian approach is used to obtain the corresponding equations of motion. The solar induced torques are assumed to be free of librational angles and are represented by their Fourier expansion. The response and periodic solutions are obtained through linear and nonlinear analyses, using the method of harmonic balance in the latter case. The stability conditions are obtained using Routh-Hurwitz criteria. To establish the ranges of validity the analytic response is compared with the numerical solution. Finally, values of the system parameters are suggested to make the satellite behave as desired. Among these is a possible approach to subdue the solar induced roll resonance. It is felt that the approximate analysis presented here should significantly reduce the computational efforts involved in the design and stability analysis of the systems.

Shrivastava, S. K.↗
NASA NTRS1996

Development of Three-Dimensional Flow Code Package to Predict Performance and Stability of Aircraft with Leading Edge Ice Contamination

In the work reported herein, a simplified, uncoupled, zonal procedure is utilized to assess the capability of numerically simulating icing effects on a Boeing 727-200 aircraft. The computational approach combines potential flow plus boundary layer simulations by VSAERO for the un-iced aircraft forces and moments with Navier-Stokes simulations by NPARC for the incremental forces and moments due to iced components. These are compared with wind tunnel force and moment data, supplied by the Boeing Company, examining longitudinal flight characteristics. Grid refinement improved the local flow features over previously reported work with no appreciable difference in the incremental ice effect. The computed lift curve slope with and without empennage ice matches the experimental value to within 1%, and the zero lift angle agrees to within 0.2 of a degree. The computed slope of the un-iced and iced aircraft longitudinal stability curve is within about 2% of the test data. This work demonstrates the feasibility of a zonal method for the icing analysis of complete aircraft or isolated components within the linear angle of attack range. In fact, this zonal technique has allowed for the viscous analysis of a complete aircraft with ice which is currently not otherwise considered tractable.

Aircraft icing↗
NASA NTRS2018

Transition Delay via Vortex Generators in a Hypersonic Boundary Layer at Flight Conditions

The potential of realizable, stationary streaks undergoing non-modal growth to stabilize a hypersonic boundary-layer flow and, subsequently, delay the laminar-turbulent transition onset, is studied via numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment, namely, a 7-degree half-angle cone with 2.5 mm nose radius, freestream Mach number of 5.30, freestream unit Reynolds number equal to 13.42 x 10(exp 6)/m, and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the test article. This paper investigates flow modifications induced by wall-mounted vortex generators (VGs), followed by an analysis of the modal instability of the perturbed, streaky boundary-layer flow. Results are presented both for a single array of VGs that was designed on the basis of optimal growth theory and for a VG configuration involving two separate arrays with opposite orientations that ware designed to provide staged control of flow instabilities while simultaneously reducing the amplification of streak instabilities resulting from the control devices. Earlier research had shown that the onset of transition during the HIFiRE-1 flight experiment, which did not include any control devices, correlated with an amplification factor of N = 14.7 for the planar Mack modes. If one assumes that the transition N -factor is not affected by the introduction of the VGs, then the control configurations based on a single array of VGs and two separate arrays would result in a transition delay of 17% and 40%, respectively. These findings suggest a passive flow control s to induce streaks that would delay transition in hypersonic boundary dominated by Mack-mode instabilities.

Incompressible boundary layer↗
NASA NTRS1986

Stability and control of compressible flows over a surface with concave-conves curvature

The active control of spatially unstable disturbances in a laminar, two-dimensional, compressible boundary layer over a curved surface is numerically simulated. The control is effected by localized time-periodic surface heating. We consider two similar surfaces of different heights with concave-convex curvature. In one, the height is sufficiently large so that the favorable pressure gradient is sufficient to stabilize a particular disturbance. In the other case the pressure gradient induced by the curvature is destabilizing. It is shown that by using active control that the disturbance can be stabilized. The results demonstrate that the curvature induced mean pressure gradient significantly enhances the receptivity of the flow localized time-periodic surface heating and that this is a potentially viable mechanism in air.

Maestrello, L.↗
NASA NTRS2007

GRMHD/RMHD Simulations and Stability of Magnetized Spine-Sheath Relativistic Jets

A new general relativistic magnetohydrodynamics (GRMHD ) code "RAISHIN" used to simulate jet generation by rotating and non-rotating black holes with a geometrically thin Keplarian accretion disk finds that the jet develops a spine-sheath structure in the rotating black hole case. Spine-sheath structure and strong magnetic fields significantly modify the Kelvin-Helmholtz (KH) velocity shear driven instability. The RAISHIN code has been used in its relativistic magnetohydrodynamic (RMHD) configuration to study the effects of strong magnetic fields and weakly relativistic sheath motion, cl2, on the KH instability associated with a relativistic, Y = 2.5, jet spine-sheath interaction. In the simulations sound speeds up to ? c/3 and Alfven wave speeds up to ? 0.56 c are considered. Numerical simulation results are compared to theoretical predictions from a new normal mode analysis of the RMHD equations. Increased stability of a weakly magnetized system resulting from c/2 sheath speeds and stabilization of a strongly magnetized system resulting from d 2 sheath speeds is found.

Hardee, Philip↗
NASA NTRS1991

Nonlinear evolution of subsonic and supersonic disturbances on a compressible free shear layer

The effects of a nonlinear-nonequilibrium-viscous critical layer on the spatial evolution of subsonic and supersonic instability modes on a compressible free shear layer is considered. It is shown that the instability wave amplitude is governed by an integrodifferential equation with cubic-type nonlinearity. Numerical and asymptotic solutions to this equation show that the amplitude either ends in a singularity at a finite downstream distance or reaches an equilibrium value, depending on the Prandtl number, viscosity law, viscous parameter and a real parameter which is determined by the linear inviscid stability theory. A necessary condition for the existence of the equilibrium solution is derived, and whether or not this condition is met is determined numerically for a wide range of physical parameters including both subsonic and supersonic disturbances. it is found that no equilibrium solution exists for the subsonic modes unless the temperature ratio of the low-to-high-speed streams exceeds a critical value, while equilibrium solutions for the most rapidly growing supersonic mode exist over most of the parameter range examined.

Leib, S. J.↗
NASA NTRS2009

Re-evaluation of an Optimized Second Order Backward Difference (BDF2OPT) Scheme for Unsteady Flow Applications

Recent experience in the application of an optimized, second-order, backward-difference (BDF2OPT) temporal scheme is reported. The primary focus of the work is on obtaining accurate solutions of the unsteady Reynolds-averaged Navier-Stokes equations over long periods of time for aerodynamic problems of interest. The baseline flow solver under consideration uses a particular BDF2OPT temporal scheme with a dual-time-stepping algorithm for advancing the flow solutions in time. Numerical difficulties are encountered with this scheme when the flow code is run for a large number of time steps, a behavior not seen with the standard second-order, backward-difference, temporal scheme. Based on a stability analysis, slight modifications to the BDF2OPT scheme are suggested. The performance and accuracy of this modified scheme is assessed by comparing the computational results with other numerical schemes and experimental data.

Vatsa, Veer N.↗
NASA NTRS1996

Simulation of Nonlinear Instabilities in an Attachment-Line Boundary Layer

The linear and the nonlinear stability of disturbances that propagate along the attachment line of a three-dimensional boundary layer is considered. The spatially evolving disturbances in the boundary layer are computed by direct numerical simulation (DNS) of the unsteady, incompressible Navier-Stokes equations. Disturbances are introduced either by forcing at the in ow or by applying suction and blowing at the wall. Quasi-parallel linear stability theory and a nonparallel theory yield notably different stability characteristics for disturbances near the critical Reynolds number; the DNS results con rm the latter theory. Previously, a weakly nonlinear theory and computations revealed a high wave-number region of subcritical disturbance growth. More recent computations have failed to achieve this subcritical growth. The present computational results indicate the presence of subcritically growing disturbances; the results support the weakly nonlinear theory. Furthermore, an explanation is provided for the previous theoretical and computational discrepancy. In addition, the present results demonstrate that steady suction can be used to stabilize disturbances that otherwise grow subcritically along the attachment line.

Joslin, Ronald D.↗
NASA NTRS1981

Current research activities at the Texas Center for Orbital Mechanics

The development of new techniques and their applications to problems in satellite geodesy, oceanography, orbit determination, long-time prediction, regularization, potential determination, stability of natural and artificial satellites and planets, etc. are described. Applications to polar motion, earth rotation and to tectonic plate motion are discussed as related to the Lageos, Starlette and Seasat missions. Numerical averaging techniques are described for the determination of the long term motion of near-resonant gravitational systems with application to Jupiter's Galilean satellites. Regions of long-time stability and ergodicity are established in the restricted problem.

Szebehely, V. G.↗
NASA NTRS1976

An analytical method for ride quality of flexible airplanes

A new and easily used state variable method of aircraft ride quality analysis is developed and its use is illustrated by a numerical example using DC-8 airplane equations of motion. The method readily allows any type of stability augmentation system to be included for good handling qualities. As part of the ride quality analysis, a particularly powerful method of specifying and achieving all roots of a desired closed-loop characteristic equation by use of full state variable feedback is described and applied to the DC-8 example.

Swaim, R. L.↗
NASA NTRS1977

An analytical method for ride quality of flexible airplanes

A new and easily used method of aircraft ride quality analysis is developed, and its use is illustrated by a numerical example using DC-8 airplane equations of motion. The method readily allows any type of stability augmentation system to be included for good handling qualities. A particularly powerful method of specifying and achieving all roots of a desired closed-loop characteristic equation by use of full state variable feedback is described and applied to the DC-8 example.

Swaim, R. L.↗
NASA NTRS1984

Experimental study of bubble cavities attached to a rotating shaft in a reservoir

Bubble cavities formed by air entrainment and attached to a rotating shaft in an oil reservoir were studied. The cavities appear to the unaided eye as toroidal. High speed photography, however, reveals the individuality of the bubble cavities and their near solid body rotational characteristics. The cavities are distorted by the rotation effects but remain attached and tend to merge because of edge effects in the axial direction. The flow field within the reservoir is influenced by the unusual character of the two phase fluid found there; the vorticity is readily visualized. Other examples of vapor entrapment at the inlet of an eccentric rotor are also discussed. A simplified analytical method is provided, and a numerical analysis is being investigated. Vapor (void) entrainment and generation can significantly alter leakage rates and stability of seals, bearings, and dampers. Recognition of these effects in the component design systems will result only after detailed studies of the above phenomena.

Hendricks, R. C.↗
NASA NTRS1984

Relative motion characteristics of 2 near-Earth Satellites

The stability of the nonlinear dynamical system of two GRAVSAT - type satellites was investigated by performing several numerical experiments which provide the simulations of the relative motion characteristics between the two satellites for various specified time intervals. The simulations included the relative range, range-rate, and relative acceleration magnitude. These simulations were generated with respect to appropriate initial orbital elements which were obtained such that the instantaneous separation distance between the two satellites has small fluctuations from a specified constant separation distance. The simulation results indicate that the behavior of the relative motions is very sensitive to the initial orbital elements of the satellites and that for a specified time interval of interest. A stable behavior is possible only with the use of an appropriate set of initial orbital elements compatible with the gravity field used to derive them.

Schutz, B. E.↗
NASA NTRS1988

Cooling flows and the stability of radio jets

This paper considers what happens to radio jets which, propagating out from the nucleus of an elliptical galaxy, encounters a point at which the ambient gas pressure increases very rapidly. The pressure variation within the sonic radius of a cooling flow is calculated, and a semianalytical calculation of the stability of a jet which encounters a rapidly increasing pressure gradient is given. The results of numerical hydrodynamic simulations which test the results of the analytic calculations are shown. The case where the jet encounters the sonic region at an oblique angle is considered, and it is suggested that this may provide an explanation for jets which are bent on small scales. These results are applied to radio sources associated with central dominant cluster galaxies and with more normal elliptical galaxies.

Soker, Noam↗
NASA NTRS1988

Effects of subcritical crack growth on fracture toughness of ceramics assessed in chevron-notched three-point bend tests

A numerical computational study was carried out to assess the effects of subcritical crack growth on crack stability in the chevron-notched three-point bend specimens. A power-law relationship between the subcritical crack velocity and the applied stress intensity were used along with compliance and stress-intensity relationships for the chevron-notched bend specimen to calculate the load response under fixed deflection rate and a machine compliance. The results indicate that the maximum load during the test occurs at the same crack length for all the deflection rates; the maximum load, however, is dependent on the deflection rate for rates below the critical rate. The resulting dependence of the apparent fracture toughness on the deflection rate is compared to experimental results on soda-lime glass and polycrystalline alumina.

Chao, L. Y.↗
NASA NTRS1988

A stable, high power optically pumped far infrared laser system

The generation of 1.25 watts of CW laser power at the 119-micron (2522.8 GHz) methanol line is reported. The maximum frequency fluctuation of the free running laser is less than + or - 100 kHz per hour. This laser has also been tested on numerous other lines ranging from 403.7 GHz (HCOOH) to 5260 GHz (CH3OD) with improved power and stability.

Farhoomand, Jam↗