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.

SourceResearch connectionFrom yearThrough year

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

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

NASA NTRS1990

The inviscid stability of supersonic flow past a sharp cone

The effects of lateral curvature on the development of supersonic laminar inviscid boundary-layer flow on a sharp cone with adiabatic wall conditions are investigated analytically, with a focus on the linear temporal inviscid stability properties. The derivation of the governing equations and of a 'triply generalized' inflexion condition is outlined, and numerical results for freestream Mach number 3.8 are presented in extensive graphs and characterized in detail. A third instability mode related to the viscous mode observed by Duck and Hall (1990) using triple-deck theory is detected and shown to be more unstable and to have larger growth rates than the second mode in some cases. It is found that the 'sonic' neutral mode is affected by the lateral curvature and becomes a supersonic neutral mode.

Duck, Peter W.↗

NASA NTRS1992

The three-dimensional evolution of a plane mixing layer. Part 2: Pairing and transition to turbulence

The evolution of three-dimensional temporally evolving plane mixing layers through as many as three pairings was simulated numerically. Initial conditions for all simulations consisted of a few low-wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity. Three-dimensional perturbations were used with amplitudes ranging from infinitesimal to large enough to trigger a rapid transition to turbulence. Pairing is found both to inhibit the growth of infinitesimal three-dimensional disturbances and to trigger the transition to turbulence in highly three dimensional flows. The mechanisms responsible for the growth of three-dimensionality as well as the initial phases of the transition to turbulence are described. The transition to turbulence is accompanied by the formation of thin sheets of span wise vorticity, which undergo a secondary roll up. Transition also produces an increase in the degree of scalar mixing, in agreement with experimental observations of mixing transition. Simulations were also conducted to investigate changes in span wise length scale that may occur in response to the change in stream wise length scale during a pairing. The linear mechanism for this process was found to be very slow, requiring roughly three pairings to complete a doubling of the span wise scale. Stronger three-dimensionality can produce more rapid scale changes but is also likely to trigger transition to turbulence. No evidence was found for a change from an organized array of rib vortices at one span wise scale to a similar array at a larger span wise scale.

Moser, Robert D.↗

NASA NTRS1992

Investigation of the flight mechanics simulation of a hovering helicopter

The flight mechanics simulation of a hovering helicopter is investigated by comparing the results of two different numerical models with flight test data for a hovering AH-64 Apache. The two models are the U.S. Army BEMAP and the Technion model. These nonlinear models are linearized by applying a numerical linearization procedure. The results of the linear models are compared with identification results in terms of eigenvalues, stability and control derivatives, and frequency responses. Detailed time histories of the responses of the complete nonlinear models, as a result of various pilots' inputs, are compared with flight test results. In addition the sensitivity of the models to various effects are also investigated. The results are discussed and problematic aspects of the simulation are identified.

Chaimovich, M.↗

NASA NTRS1993

On the generalized VIP time integral methodology for transient thermal problems

The paper describes the development and applicability of a generalized VIrtual-Pulse (VIP) time integral method of computation for thermal problems. Unlike past approaches for general heat transfer computations, and with the advent of high speed computing technology and the importance of parallel computations for efficient use of computing environments, a major motivation via the developments described in this paper is the need for developing explicit computational procedures with improved accuracy and stability characteristics. As a consequence, a new and effective VIP methodology is described which inherits these improved characteristics. Numerical illustrative examples are provided to demonstrate the developments and validate the results obtained for thermal problems.

Mei, Youping↗

NASA NTRS1992

Multi-dimensional high order essentially non-oscillatory finite difference methods in generalized coordinates

The nonlinear stability of compact schemes for shock calculations is investigated. In recent years compact schemes were used in various numerical simulations including direct numerical simulation of turbulence. However to apply them to problems containing shocks, one has to resolve the problem of spurious numerical oscillation and nonlinear instability. A framework to apply nonlinear limiting to a local mean is introduced. The resulting scheme can be proven total variation (1D) or maximum norm (multi D) stable and produces nice numerical results in the test cases. The result is summarized in the preprint entitled 'Nonlinearly Stable Compact Schemes for Shock Calculations', which was submitted to SIAM Journal on Numerical Analysis. Research was continued on issues related to two and three dimensional essentially non-oscillatory (ENO) schemes. The main research topics include: parallel implementation of ENO schemes on Connection Machines; boundary conditions; shock interaction with hydrogen bubbles, a preparation for the full combustion simulation; and direct numerical simulation of compressible sheared turbulence.

Shu, Chi-Wang↗

NASA NTRS1994

Numerical simulation of supersonic boundary layer transition

The present contribution reviews some of the recent progress obtained at our group in the direct numerical simulation (DNS) of compressible boundary layer transition. Elements of the different simulation approaches and numerical techniques employed are surveyed. Temporal and spatial simulations, as well as comparisons with results obtained from Parabolized Stability Equations, are discussed. DNS results are given for flat plate boundary layers in the Mach number range 1.6 to 4.5. A temporal DNS at Mach 4.5 has been continued through breakdown all the way to the turbulent stage. In addition results obtained with a recently developed extended temporal DNS approach are presented, which takes into account some nonparallel effects of a growing boundary layer. Results from this approach are quite close to those of spatial DNS, while preserving the efficiency of the temporal DNS.

Guo, Y.↗

NASA NTRS2013

Crystal Shape Evolution in Detached Bridgman Growth

Detached (or dewetted) Bridgman crystal growth defines that process in which a gap exists between a growing crystal and the crucible wall. Existence of the gap provides several advantages, including no sticking of the crystal to the crucible wall, reduced thermal and mechanical stresses, reduced dislocations, and no heterogeneous nucleation by the crucible. Numerical calculations are used to determine the conditions in which a gap can exist. According to crystal shape stability theory, only some of these gap widths will be dynamically stable. Beginning with a crystal diameter that differs from stable conditions, the transient crystal growth process is analyzed. In microgravity, dynamic stability depends only on capillary effects and is decoupled from heat transfer. Depending on the initial conditions and growth parameters, the crystal shape will evolve towards the crucible wall, towards a stable gap width, or towards the center of the crucible, collapsing the meniscus. The effect of a tapered crucible on dynamic stability is also described

Volz, M. P.↗

NASA NTRS2013

Crystal Shape Evolution in Detached Bridgman Growth

Volz, M. P.↗

NASA NTRS2014

Design and Analysis of Morpheus Lander Flight Control System

The Morpheus Lander is a vertical takeoff and landing test bed vehicle developed to demonstrate the system performance of the Guidance, Navigation and Control (GN&C) system capability for the integrated autonomous landing and hazard avoidance system hardware and software. The Morpheus flight control system design must be robust to various mission profiles. This paper presents a design methodology for employing numerical optimization to develop the Morpheus flight control system. The design objectives include attitude tracking accuracy and robust stability with respect to rigid body dynamics and propellant slosh. Under the assumption that the Morpheus time-varying dynamics and control system can be frozen over a short period of time, the flight controllers are designed to stabilize all selected frozen-time control systems in the presence of parametric uncertainty. Both control gains in the inner attitude control loop and guidance gains in the outer position control loop are designed to maximize the vehicle performance while ensuring robustness. The flight control system designs provided herein have been demonstrated to provide stable control systems in both Draper Ares Stability Analysis Tool (ASAT) and the NASA/JSC Trick-based Morpheus time domain simulation.

Jang, Jiann-Woei↗

NASA NTRS2013

Evolution of the Shape of Detached GeSi Crystals in Microgravity

A series of GeSi crystal growth experiments are planned to be conducted in the Low Gradient Furnace (LGF) onboard the International Space Station. An objective of these experiments is to understand the mechanisms of detached Bridgman growth, a process in which a gap exists between the growing semiconductor crystal and the crucible wall. Crystals grown without wall contact have superior quality to otherwise similar crystals grown in direct contact with a container, especially with respect to impurity incorporation, formation of dislocations, and residual stress in crystals. Numerical calculations are used to determine the conditions in which a gap can exist. According to crystal shape stability theory, only some of these gap widths will be dynamically stable. Beginning with a crystal diameter that differs from stable conditions, the transient crystal growth process is analyzed. In microgravity, dynamic stability depends only on capillary effects and is decoupled from heat transfer. Depending on the initial conditions and growth parameters, the crystal shape will evolve towards the crucible wall, towards a stable gap width, or towards the center of the crucible, collapsing the meniscus.

Volz, M. P.↗

NASA NTRS2013

Evolution of the Shape of Detached GeSi Crystals in Microgravity

Volz, M. P.↗

NASA NTRS1961

Estimation of Rotary Stability Derivatives at Subsonic and Transonic Speeds

The first part of this paper pertains to the estimation of subsonic rotary stability derivatives of wings. The unsteady potential flow problem is solved by a superposition of steady flow solutions. Numerical results for the damping coefficients of triangular wings are presented as functions of aspect ratio and Mach number, and are compared with experimental results over the Mach number range 0 to 1. In the second part, experimental results are used. to point out a close correlation between the nonlinear variations with angle of attack of the static pitching-moment curve slope and the damping-in-pitch coefficient. The underlying basis for the correlation is found as a result of an analysis in which the indicial function concept and. the principle of super-position are adapted to apply to the nonlinear problem. The form of the result suggests a method of estimating nonlinear damping coefficients from results of static wind-tunnel measurements.

Tobak, Murray↗

NASA NTRS2016

Nonlinear Transient Growth and Boundary Layer Transition

Parabolized stability equations (PSE) are used in a variational approach to study the optimal, non-modal disturbance growth in a Mach 3 at plate boundary layer and a Mach 6 circular cone boundary layer. As noted in previous works, the optimal initial disturbances correspond to steady counter-rotating streamwise vortices, which subsequently lead to the formation of streamwise-elongated structures, i.e., streaks, via a lift-up effect. The nonlinear evolution of the linearly optimal stationary perturbations is computed using the nonlinear plane-marching PSE for stationary perturbations. A fully implicit marching technique is used to facilitate the computation of nonlinear streaks with large amplitudes. To assess the effect of the finite-amplitude streaks on transition, the linear form of plane- marching PSE is used to investigate the instability of the boundary layer flow modified by spanwise periodic streaks. The onset of bypass transition is estimated by using an N- factor criterion based on the amplification of the streak instabilities. Results show that, for both flow configurations of interest, streaks of sufficiently large amplitude can lead to significantly earlier onset of transition than that in an unperturbed boundary layer without any streaks.

Laminar turbulent transition↗

NASA NTRS2016

Influence of Stationary Crossflow Modulation on Secondary Instability

A likely scenario for swept wing transition on subsonic aircraft with natural laminar flow involves the breakdown of stationary crossflow vortices via high frequency secondary instability. A majority of the prior research on this secondary instability has focused on crossflow vortices with a single dominant spanwise wavelength. This paper investigates the effects of the spanwise modulation of stationary crossflow vortices at a specified wavelength by a subharmonic stationary mode. Secondary instability of the modulated crossflow pattern is studied using planar, partial-differential-equation based eigenvalue analysis. Computations reveal that weak modulation by the first subharmonic of the input stationary mode leads to mode splitting that is particularly obvious for Y-type secondary modes that are driven by the wall-normal shear of the basic state. Thus, for each Y mode corresponding to the fundamental wavelength of results in unmodulated train of crossflow vortices, the modulated flow supports a pair of secondary modes with somewhat different amplification rates. The mode splitting phenomenon suggests that a more complex stationary modulation such as that induced by natural surface roughness would yield a considerably richer spectrum of secondary instability modes. Even modest levels of subharmonic modulation are shown to have a strong effect on the overall amplification of secondary disturbances, particularly the Z-modes driven by the spanwise shear of the basic state. Preliminary computations related to the nonlinear breakdown of these secondary disturbances provide interesting insights into the process of crossflow transition in the presence of the first subharmonic of the dominant stationary vortex.

Vortex structure↗

NASA NTRS2026

Improvement in the Thermal-to-Structural Model Mapping Process for Integrated Modeling for the Roman Space Telescope

Integrated Modeling has been a key component of verifying optical requirements for the Nancy Grace Roman Space Telescope (RST) that are either impossible or impractical to verify exclusively through ground testing. Two major areas for integrated Modeling are Jitter and Thermal Distortion that require the exchanges of model performance predictions across disciplines. In both cases, distortions are impressed on optical models to evaluate the impact on boresight alignment and wave front error. In the case of Jitter, the disturbances are driven by reactions to motions most often from actuators; however, in the case of thermal distortion, the motions are driven by thermal expansion or contraction as a result of changing temperatures. This then requires a link further upstream to the thermal model, which is used to predict the thermal performance and temperature gradients and stability. The process for mapping temperatures from a thermal model to a corresponding structural model has been performed numerous times through the RST project lifecycle, with improvements in the accuracy, verification, and effort sought throughout. This paper describes some of the recent improvements to the process, including: capture of the visualization parameters, automatic generation of the mapped images for both the thermal and structural model groupings, and reduction in the effort to assemble the full set of mapped temperatures. These upgrades have greatly reduced the manual effort associated with thermal mapping and allowed for faster turn-around of Integrated Modeling predictions.

Thermal Mapping↗

NASA NTRS1994

Stabilization Mechanisms and Burning Rates of Cylindrical Burner Flames

A study is conducted of the structure and response of curved (but unstretched), cylindrically-symmetric 1D premixed flames from a cylindrical porous burner. The study has employed (1) activation-energy asymptotics with one-step reaction constant and constant properties; (2) a numerical computation which encompassed detailed chemistry and transport behavior, and (3) drop-tower microgravity tests. Attention was given to the relative importance of heat loss vs. flow divergence as the dominant mechanism for flame stabilization; the results show that, with increasing flow discharge rate, the dominant flame stabilization mechanism changes from heat loss to flow divergence.

Eng, J. A.↗

NASA NTRS2002

Identifying head-trunk and lower limb contributions to gaze stabilization during locomotion

The goal of the present study was to determine how the multiple, interdependent full-body sensorimotor subsystems respond to a change in gaze stabilization task constraints during locomotion. Nine subjects performed two gaze stabilization tasks while walking at 6.4 km/hr on a motorized treadmill: 1) focusing on a central point target; 2) reading numeral characters; both presented at 2 m in front at the level of their eyes. While subjects performed the tasks we measured: temporal parameters of gait, full body sagittal plane segmental kinematics of the head, trunk, thigh, tibia and foot, accelerations along the vertical axis at the head and the tibia, and the vertical forces acting on the support surface. We tested the hypothesis that with the increased demands placed on visual acuity during the number recognition task, subjects would modify full-body segmental kinematics in order to reduce perturbations to the head in order to successfully perform the task. We found that while reading numeral characters as compared to the central point target: 1) compensatory head pitch movement was on average 22% greater despite the fact that the trunk pitch and trunk vertical translation movement control were not significantly changed; 2) coordination patterns between head and trunk as reflected by the peak cross correlation between the head pitch and trunk pitch motion as well as the peak cross correlation between the head pitch and vertical trunk translation motion were not significantly changed; 3) knee joint total movement was on average 11% greater during the period from the heel strike event to the peak knee flexion event in stance phase of the gait cycle; 4) peak acceleration measured at the head was significantly reduced by an average of 13% in four of the six subjects. This was so even when the peak acceleration at the tibia and the transmission of the shock wave at heel strike (measured by the peak acceleration ratio of the head/tibia and the time lag between the tibial and head peak accelerations) remained unchanged. Taken together these results provide further evidence that the full body contributes to gaze stabilization during locomotion, and that its different functional elements can be modified online to contribute to gaze stabilization for different visual task constraints.

NASA Center JSC↗

NASA NTRS1980

Nonlinear grid error effects on numerical solution of partial differential equations

Finite difference solutions of nonlinear partial differential equations require discretizations and consequently grid errors are generated. These errors strongly affect stability and convergence properties of difference models. Previously such errors were analyzed by linearizing the difference equations for solutions. Properties of mappings of decadence were used to analyze nonlinear instabilities. Such an analysis is directly affected by initial/boundary conditions. An algorithm was developed, applied to nonlinear Burgers equations, and verified computationally. A preliminary test shows that Navier-Stokes equations may be treated similarly.

Dey, S. K.↗