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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Olsson & Oliger, Gerritsen & Olsson, and Yee et al. (15, 7, 37) is revisited. The objective of this paper is to prove for the first time that the entropy split scheme is an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for nonlinear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. The proof is to replace the spatial derivatives by summation-by-parts (SBP) difference operators in the entropy split form of the equations using the physical entropy of the Euler equations. The numerical boundary closure follows directly from the SBP operator. No additional numerical boundary procedure is required. In contrast, Tadmor-type entropy conserving schemes (31) using mathematical entropies and more recently in (35], do not naturally come with a numerical boundary closure and a generalized SBP operator has to be developed (18). Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grubber splitting.

Sjogreen, Bjorn↗
NASA NTRS1986

Nonlinear flexure and torsion of rotating beams, with application to helicopter rotor blades. I - Formulation. II - Response and stability results

The dynamic response and aeroelastic stability of rotating beams such as helicopter blades is investigated analytically. The Hamilton principle is used to formulate the equations of motion for extensional and inextensional beams with precone angles and variable pitch angles, taking higher-order nonlinearities into account. The derivation of the equations and their approximate solution by a Galerkin procedure are explained in detail, and numerical results of equilibrium solutions and stability analyses are presented graphically.

Crespo Da Silva, M. R. M.↗
NASA NTRS1981

Cylindrical annulus models

The stability diagram for the atmospheric general circulation experiment (AGCE) was numerically determined. The cylindrical code of Warn-Varnas modified for linear stability and adapted to spherical geometry. The code computes axisymmetric basic states and then test their stability to zonal disturbances. The Warn-Varnas code and some of the basic state results are reported.

Leslie, F.↗
NASA NTRS1990

The inviscid stability of supersonic flow past a sharp cone

The laminar boundary layer which forms on a sharp cone in a supersonic freestream, where lateral curvature plays a key role in the physics of the problem is considered. This flow is then analyzed from the point of view of linear, temporal, inviscid stability. The basic, non-axisymmetric disturbance equations are derived for general flows of this class, and a so called triply generalized inflexion condition is found for the existence of subsonic neutral modes of instability. This condition is analogous to the well-known generalized inflexion condition found in planar flows, although in the present case the condition depends on both axial and aximuthal wavenumbers. Extensive numerical results are presented for the stability problem at a freestream Mach number of 3.8, for a range of streamwise locations. These results reveal that a new mode of instability may occur, peculiar to flows of this type involving curvature. Additionally, asymptotic analyses valid close to the tip of the cone, far downstream of the cone are presented, and these give a partial (asymptotic) description of this additional mode of instability.

Duck, Peter W.↗
NASA NTRS2003

Technical Evaluation Report, Part A - Vortex Flow and High Angle of Attack

A symposium entitled Vortex Flow and High Angle of Attack was held in Loen, Norway, from May 7 through May 11, 2001. The Applied Vehicle Technology (AVT) panel, under the auspices of the Research and Technology Organization (RTO), sponsored this symposium. Forty-eight papers, organized into nine sessions, addressed computational and experimental studies of vortex flows pertinent to both aircraft and maritime applications. The studies also ranged from fundamental fluids investigations to flight test results, and significant results were contributed from a broad range of countries. The principal emphasis of this symposium was on "the understanding and prediction of separation-induced vortex flows and their effects on military vehicle performance, stability, control, and structural design loads." It was further observed by the program committee that "separation- induced vortex flows are an important part of the design and off-design performance of conventional fighter aircraft and new conventional or unconventional manned or unmanned advanced vehicle designs (UAVs, manned aircraft, missiles, space planes, ground-based vehicles, and ships)." The nine sessions addressed the following topics: vortical flows on wings and bodies, experimental techniques for vortical flows, numerical simulations of vortical flows, vortex stability and breakdown, vortex flows in maritime applications, vortex interactions and control, vortex dynamics, flight testing, and vehicle design. The purpose of this paper is to provide brief reviews of these papers along with some synthesizing perspectives toward future vortex flow research opportunities. The paper includes the symposium program. (15 refs.)

VORTEX BREAKDOWN↗
NASA NTRS1993

Mapping the stability region of the 3:2 Neptune-Pluto resonance

Pluto and Charon are most likely the remnants of a large number of objects that existed in the Uranus-Neptune region at early epochs of the solar system. Numerical integrations have shown that, in general, such objects were ejected from the planetary region on timescales of approximately 10(exp 7) years after Neptune and Uranus reached their current masses. It is thought that the Pluto-Charon system survived to current times without being dynamically removed in this way because it is trapped in a set of secular and mean motion resonances with Neptune. The best-known Pluto-Neptune orbit coupling is the 3:2 mean motion resonance discovered almost 30 years ago by C. Cohen and E. Hubbard. These workers showed that the resonance angle, delta is equivalent to 3(lambda(sub P)) - 2(lambda(sub N)) - omega-bar(sub P) where omega-bar(sub P) is the longitude of perihelion of the Pluto-Charon system, and lambda(sub N) and lambda(sub P) are the mean longitude of Neptune and Pluto-Charon respectively, librates about 180 deg with an amplitude, A(sub delta), of 76 deg. A numerical simulation project to map out the stability region of the 3:2 resonance is reported. The results of these simulations are important to understanding whether Pluto's long-term heliocentric stability requires only the 3:2 resonance, or whether it instead requires one or more of the other Pluto-Neptune resonances. Our study also has another important application. By investigating stability timescales as a function of orbital elements, we gain insight into the fraction of orbital phase space which the stable 3:2 resonance occupies. This fraction is directly related to the probability that the Pluto-Charon system (and possibly other small bodies) could have been captured into this resonance.

Levison, H. F.↗
NASA NTRS2010

Receptivity of Hypersonic Boundary Layers over Straight and Flared Cones

The effects of adverse pressure gradients on the receptivity and stability of hypersonic boundary layers were numerically investigated. Simulations were performed for boundary layer flows over a straight cone and two flared cones. The steady and the unsteady flow fields were obtained by solving the two-dimensional Navier-Stokes equations in axi-symmetric coordinates using the 5th order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The mean boundary layer profiles were analyzed using local stability and non-local parabolized stability equations (PSE) methods. After the most amplified disturbances were identified, two-dimensional plane acoustic waves were introduced at the outer boundary of the computational domain and time accurate simulations were performed. The adverse pressure gradient was found to affect the boundary layer stability in two important ways. Firstly, the frequency of the most amplified second-mode disturbance was increased relative to the zero pressure gradient case. Secondly, the amplification of first- and second-mode disturbances was increased. Although an adverse pressure gradient enhances instability wave growth rates, small nose-tip bluntness was found to delay transition due to the low receptivity coefficient and the resulting weak initial amplitude of the instability waves. The computed and measured amplitude-frequency spectrums in all three cases agree very well in terms of frequency and the shape except for the amplitude.

Balakumar, Ponnampalam↗
NASA NTRS1982

On the stability of Schwarzschild's triaxial galaxy model

The numerical model for a triaxial galaxy constructed by Schwartzschild (1979) has been tested by means of fully self-consistent three-dimensional n-body numerical experiments that use a representation of this model to start the integration. The integration, study models, initial load, and correcting potential of the experiment are detailed. Preliminary checks included density values, potential values, and integration in potential of an ideal model. The results showed some particle loss and no significant shape change. The model is rugged and robust, and no growing disturbances with growth rate in excess of 0.5 per crossing time could be detected. Runs were done at 50 percent excess halo and with tube orbits reversed.

Smith, B. F.↗
NASA NTRS1999

Experimental and Numerical Optimization of a High-Lift System to Improve Low-Speed Performance, Stability, and Control of an Arrow-Wing Supersonic Transport

An investigation was performed to evaluate leading-and trailing-edge flap deflections for optimal aerodynamic performance of a High-Speed Civil Transport concept during takeoff and approach-to-landing conditions. The configuration used for this study was designed by the Douglas Aircraft Company during the 1970's. A 0.1-scale model of this configuration was tested in the Langley 30- by 60-Foot Tunnel with both the original leading-edge flap system and a new leading-edge flap system, which was designed with modem computational flow analysis and optimization tools. Leading-and trailing-edge flap deflections were generated for the original and modified leading-edge flap systems with the computational flow analysis and optimization tools. Although wind tunnel data indicated improvements in aerodynamic performance for the analytically derived flap deflections for both leading-edge flap systems, perturbations of the analytically derived leading-edge flap deflections yielded significant additional improvements in aerodynamic performance. In addition to the aerodynamic performance optimization testing, stability and control data were also obtained. An evaluation of the crosswind landing capability of the aircraft configuration revealed that insufficient lateral control existed as a result of high levels of lateral stability. Deflection of the leading-and trailing-edge flaps improved the crosswind landing capability of the vehicle considerably; however, additional improvements are required.

Hahne, David E.↗
NASA NTRS2009

Mean Element Propagations Using Numerical Averaging

The long-term evolution characteristics (and stability) of an orbit are best characterized using a mean element propagation of the perturbed two body variational equations of motion. The averaging process eliminates short period terms leaving only secular and long period effects. In this study, a non-traditional approach is taken that averages the variational equations using adaptive numerical techniques and then numerically integrating the resulting EOMs. Doing this avoids the Fourier series expansions and truncations required by the traditional analytic methods. The resultant numerical techniques can be easily adapted to propagations at most solar system bodies.

orbit propagations↗
NASA NTRS1999

Numerical Study of Buoyancy and Differential Diffusion Effects on the Structure and Dynamics of Triple Flames

Triple flames arise in a number of practical configurations where fuel and oxidizer are partially premixed, such as in the base of a lifted jet flame. Past experimental studies, theoretical analyses, and numerical modeling of triple flames suggested the potential role of triple flames in stabilizing turbulent flames and in promoting flame propagation. From recent numerical simulations of laminar triple flames, a strong influence of differential diffusion among species and heat on the triple flame structure has been gradually appreciated. This paper reports preliminary numerical results on the influence of gravity and differential diffusion effects on the structure and dynamics of triple flames with a one-step global irreversible chemistry model.

Chen, J. -Y.↗
NASA NTRS1985

Energy Stability of Thermocapillary Convection in Models of the Float Zone Process

The energy-stability of thermocapillary convection in models of the float-zone, crystal-growing process was studied. Stability limits, as functions of pertinent parameters, that will identify conditions which will not allow the existence of an undesirable oscillatory flow instability were determined. Such instabilities may occur in the space processing of semiconductor materials. The determination of the stability limits will involve two sets of numerical computations: (1) solution of the nonlinear governing equations together with the appropriate boundary conditions to determine the basic state (in general, velocity, pressure and temperature fields and the displacement of free surfaces and interfaces); and (2) solution of a nonlinear Euler-Lagrange systems for the energy-stability limit. Both computations, while difficult, should be within the scope of available computer capability and available concepts in numerical analysis. Finite-element methods are attractive candidates for the numerical work.

Neitzel, G. P.↗
NASA NTRS1993

Bifurcation analysis of a simple analytic model of self-propagating star formation

We investigate the structure and stability of rotationally symmetric nonhomogeneous time-independent solutions derived from a simple analytic model of self-propagating star formation. For this purpose we employ two methodologies: We use bifurcation theoretical methods to prove the existence of nonhomogeneous axisymmetric stationary solutions of an appropriate nonlinear evolution equation for the stellar density. We show that the nonhomogeneous solution branch bifurcates from the homogeneous one at a critical parameter value of the star formation rate. Further, the analytical theory allows us to show that the new solution set is stable in the weakly nonlinear regime near the bifurcation point. To follow the solution branch further, we use numerical methods. The numerical calculation shows the structure and stability of these solutions. We conclude that no periodic time-dependent solutions of this special model exist, and no further bifurcations can be found. The same results have been found in simulations of stochastic self-propagating star formation based on similar models. Therefore, our findings provide a natural explanation, why long-lived large-scale structure have not been found in those simulations.

Neukirch, Thomas↗
NASA NTRS1986

The effect of beam-driven return current instability on solar hard X-ray bursts

The problem of electrostatic wave generation by a return current driven by a small area electron beam during solar hard X-ray bursts is discussed. The marginal stability method is used to solve numerically the electron and ion heating equations for a prescribed beam current evolution. When ion-acoustic waves are considered, the method appears satisfactory and, following an initial phase of Coulomb resistivity in which T sub e/T sub i rise, predicts a rapid heating of substantial plasma volumes by anomalous ohmic dissipation. This hot plasma emits so much thermal bremsstrahlung that, contrary to previous expectations, the unstable beam-plasma system actually emits more hard X-rays than does the beam in the purely collisional thick target regime relevant to larger injection areas. Inclusion of ion-cyclotron waves results in ion-acoustic wave onset at lower T sub e/T sub i and a marginal stability treatment yields unphysical results.

Cromwell, D.↗
NASA NTRS1989

The inviscid axisymmetric stability of the supersonic flow along a circular cylinder

The supersonic flow past a thin straight circular cylinder is investigated. The associated boundary layer flow (i.e., the velocity and temperature field) is computed; the asymptotic, far downstream solution is obtained, and compared with the full numerical results. The inviscid, linear, axisymmetric (temporal) stability of this boundary layer is also studied. A so called doubly generalized inflexion condition is derived, which is a condition for the existence of so called subsonic neutral modes. The eigenvalue problem (for the complex wavespeed) is computed for two freestream Mach numbers (2.8 and 3.8), and this reveals that curvature has a profound effect on the stability of the flow. The first unstable inviscid mode is seen to rapidly disappear as curvature is introduced, while the second (and generally the most important) mode suffers a substantially reduced amplification rate.

Duck, Peter W.↗
NASA NTRS2019

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Ols- son & Oliger (1994), Gerritsen & Olsson (1996), and Yee et al. (2000) is revisited. The entropy splitting is a form of skew-symmetric splitting of the nonlinear Euler flux derivatives. Central differencing applied to the entropy splitting form of the Euler flux derivatives together with SBP difference operators will, hereafter, be referred to as entropy split schemes. This study is prompted by the recent growing interest in numerical methods for which a discrete entropy conservation law holds, a discrete global entropy conservation can be proved and/or the numerical method possesses a stable entropy in the framework of SBP difference operators and L2-energy norm estimate. The objective of this paper is to recast the entropy split scheme as the re- cent definition of an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for non- linear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grub- ber splitting.

Sjogreen, Bjorn↗
NASA NTRS1990

The inviscid axisymmetric stability of the supersonic flow along a circular cylinder

The supersonic flow past a thin straight circular cylinder is investigated. The associated boundary-layer flow (i.e. the velocity and temperature field) is computed; the asymptotic, far downstream solution is obtained, and compared with the full numerical results. The inviscid, linear, axisymmetric (temporal) stability of this boundary layer is also studied. A so-called 'doubly generalized' inflexion condition is derived, which is a condition for the existence of so-called 'subsonic' neutral modes. The eigenvalue problem (for the complex wavespeed) is computed for two free-stream Mach numbers (2.8 and 3.8), and this reveals that curvature has a profound effect on the stability of the flow. The first unstable inviscid mode is seen to disappear rapidly as curvature is introduced, while the second (and generally the most important) mode suffers a substantially reduced amplification rate.

Duck, Peter W.↗