Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical stabilization”

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 613 records · Page 34

On the numerical computation of nonlinear force-free magnetic fields

An algorithm has been developed to extrapolate nonlinear force-free magnetic fields from the photosphere, given the proper boundary conditions. This paper presents the results of this work, describing the mathematical formalism that was developed, the numerical techniques employed, and comments on the stability criteria and accuracy developed for these numerical schemes. An analytical solution is used for a benchmark test; the results show that the computational accuracy for the case of a nonlinear force-free magnetic field was on the order of a few percent (less than 5 percent). This newly developed scheme was applied to analyze a solar vector magnetogram, and the results were compared with the results deduced from the classical potential field method. The comparison shows that additional physical features of the vector magnetogram were revealed in the nonlinear force-free case.

Wu, S. T.↗

The inviscid stability of supersonic flow past axisymmetric bodies

The supersonic flow past a sharp cone is studied. The associated boundary layer flow (i.e., the velocity and temperature field) is computed. The inviscid linear temporal stability of axisymmetric boundary layers in general is considered, and in particular, a so-called 'triply generalized' inflection condition for 'subsonic' nonaxisymmetric neutral modes is presented. Preliminary numerical results for the stability of the cone boundary layer are presented for a freestream Mach number of 3.8. In particular, a new inviscid mode of instability is seen to occur in certain regimes, and this is shown to be related to a viscous mode found by Duck and Hall (1988).

Duck, Peter W.↗

Calculation of the lateral-dynamic stability of aircraft

Graphs and formulas are given with the aid of which all the aerodynamic coefficients required for computing the lateral dynamic stability can be determined. A number of numerical examples are given for obtaining the stability derivatives and solving the characteristic-stability equation. Approximate formulas are derived with the aid of which rapid preliminary computations may be made and the stability coefficients corrected for certain modifications of the airplane. A derivation of the lateral-dynamic-stability equations is included.

STABILITY, LATERAL - DYNAMIC↗

Calculated Low-Speed Steady and Time-Dependent Aerodynamic Derivatives for Several Different Wings Using a Discrete Vortex Method

Calculated numerical values for some aerodynamic terms and stability Derivatives for several different wings in unseparated inviscid incompressible flow were made using a discrete vortex method involving a limited number of horseshoe vortices. Both longitudinal and lateral-directional derivatives were calculated for steady conditions as well as for sinusoidal oscillatory motions. Variables included the number of vortices used and the rotation axis/moment center chordwise location. Frequencies considered were limited to the range of interest to vehicle dynamic stability (kb <.24 ). Comparisons of some calculated numerical results with experimental wind-tunnel measurements were in reasonable agreement in the low angle-of-attack range considering the differences existing between the mathematical representation and experimental wind-tunnel models tested. Of particular interest was the presence of induced drag for the oscillatory condition.

Riley, Donald R.↗

On the Stability of Strange Dwarf Hybrid Stars

We investigate the stability of “strange dwarfs”: white-dwarf-sized stars with a density discontinuity between a small dense core of quark matter and a thick low-density mantle of degenerate electrons. Previous work on strange dwarfs suggested that such a discontinuity could stabilize stars that would have been classified as unstable by the conventional criteria based on extrema in the mass–radius relation. We investigate the stability of such stars by numerically solving the Sturm–Liouville equations for the lowest-energy modes of the star. We find that the conventional criteria are correct, and strange dwarfs are not stable.

79 ASTRONOMY AND ASTROPHYSICS↗

High-order finite-volume modeling of drift waves

The study discusses high-order finite-volume numerical modeling of drift waves, which is an ubiquitous phenomenon in magnetized plasmas. It is found that some standard discretization methods applied to the conservative form of the governing equations can lead to a numerical instability. A method to stabilize high-order discretization is proposed and demonstrated to work in numerical simulations performed with the fourth-order finite-volume code COGENT. As practical examples, a stable drift-wave solution with adiabatic electrons and the collisionless (universal) drift-wave instability driven by electron kinetic effects are considered. Finally, application of the present analysis to a broader range of computational fluid dynamics systems is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Future of Land-Based Wind Turbine Rotor Technology: The Perspective from NREL

The Big Adaptive Rotor (BAR) Project looks at potential innovative pathways for future land-based wind turbine technology. BAR is led by the National Renewable Energy Laboratory and is sponsored by the Wind Energy Technology Office of the US Department of Energy. This talk will describe innovations and challenges faced by the teams working on BAR. The innovations include numerical models predicting the stability and performance of highly flexible blades, downwind rotors to maximize wind farm power, and controlled bending of 100-meter-long blades during rail transport. The challenges of developing predictive numerical models and formulating successful value propositions supporting new technologies such as distributed aero control devices will also be discussed.

BAR↗

CFD validation of a supersonic laminar flow control concept

A three-dimensional Navier-Stokes code is used in conjunction with a linear compressible stability analysis code to develop a numerical procedure for prediction of laminar flow transition. The procedure is applied to a modified F-16XL fighter with a laminar flow control glove at supersonic speed. Details of boundary layer stability analysis indicate that, computationally, laminar flow could be realized on the highly swept wing in the absence of the leading edge attachment-line contamination. Effects of the three-dimensionality of the flow were shown to be important in the boundary layer stability analysis. The numerically predicted surface pressures compare favorably with the flight test data.

Woan, Chung-Jin↗

Jupiter's winds and Arnol'd's second stability theorem: Slowly moving waves and neutral stability

Since the Voyager encounters in 1979, it has been known that Jupiter's cloud-top zonal winds violate the barotropic stability criterion. A vortex-tube stretching analysis of the Voyager wind data indicates that the more general Charney-Stern stability criterion is also violated. On the other hand, the zonal winds determined by tracking cloud features in Hubble Space Telescope images taken in 1991 precisely match the zonal winds determined by tracking cloud features in Voyager images, and it is hard to understand how a complicated zonal wind profile like Jupiter's could be unstable and yet not change at all in 12 years. In fact, there are at least two unknown ways to violate the Charney-Stern stability criterion and still have a stable flow. The better known of these is called Fjortoft's theorem, or Arnol'd's 1st theorem for the case of large-amplitude perturbations. Although the Fjortoft-Arnol'd theorem has been extended from the quasi-geostrophic equations to the primitive equations, the basic requirement that the potential vorticity be an increasing function of streamfunction is opposite to the case found in Jupiter, where the Voyager data indicate that the potential vorticity is a decreasing function of streamfunction. But this second case is precisely that which is covered by Arnol'd's 2nd stability theorem. In fact, the Voyager data suggest that Jupiter's zonal winds are neutrally stable with respect to Arnol'd's 2nd stability theorem. Here, we analyze the linear stability problem of a one-parameter family of sinusoidal zonal wind profiles that are close to neutral stability with respect to Arnol'd's 2nd stability theorem. We find numerically that the most unstable mode is always stationary, which may help to explain the slowly moving mode 10 waves observed on Jupiter. We find that violation of Arnol'd's 2nd stability theorem is both necessary and sufficient for instability of sinusoidal profiles. However, there appears to be no simple extension of Arnol'd's 2nd stability theorem to the primitive equations. Nevertheless, the primitive growth rates are small, and the primitive system is still governed by the quasi-geostrophic neutral-stability configuration.

Stamp, Andrew P.↗

Stability Limits and Dynamics of Nonaxisymmetric Liquid Bridges

Theoretical and experimental investigation of the stability of nonaxisymmetric and nonaxisymmetric bridges contained between equal and unequal radii disks as a function of Bond and Weber number with emphasis on the transition from unstable axisymmetric to stable nonaxisymmetric shapes. Numerical analysis of the stability of nonaxisymmetric bridges for various orientations of the gravity vector for equal and unequal disks. Experimental and theoretical investigation of large (nonaxisymmetric) oscillations and breaking of liquid bridges. This project involves both experimental and theoretical work. Static and dynamic experiments are conducted in a Plateau tank which makes a range of static Bond numbers accessible.

Alexander, J. Iwan D.↗

Dynamics and Statics of Nonaxisymmetric Liquid Bridges

Theoretical and experimental investigation of the stability of nonaxisymmetric and nonaxisymmetric bridges contained between equal and unequal radii disks as a function of Bond and Weber number with emphasis on the transition from unstable axisymmetric to stable nonaxisymmetric shapes, are conducted. Numerical analysis of the stability of nonaxisymmetric bridges between unequal disks for various orientations of the gravity vector is performed. Experimental and theoretical investigation of large (nonaxisymmetric) oscillations and breaking of liquid bridges are also conducted.

Alexander, J. Iwan D.↗

Finite volume computation of unsteady inviscid rotational transonic flows past airfoils in rigid body motion

Unsteady inviscid transonic flow over airfoils in arbitrary rigid body motion is analyzed numerically by solving the two-dimensional unsteady Euler equations in integral form using a finite volume scheme. The solution procedure is based on an explicit Runge-Kutta time-stepping scheme wherein the spatial terms are central-differenced and a combination of second- and fourth-differences in the flow variables are used to form the numerical dissipation terms to stabilize the scheme. Unsteady calculations are started from converged steady-state solutions as initial conditions. Nonreflective boundary conditions are imposed on the far-field boundaries. Results are presented and, where possible, validated against available numerical and experimental data for airfoils subjected to a step change in angle of attack, airfoils oscillating and plunging in transonic flow, and airfoils immersed in a time-varying free stream.

Damodaran, Murali↗

The effect of non-uniform mass loading on the linear, temporal development of particle-laden shear layers

The effect of non-uniformity in bulk particle mass loading on the linear development of a particle-laden shear layer is analyzed by means of a stochastic Eulerian-Eulerian model. From the set of governing equations of the two-fluid model, a modified Rayleigh equation is derived that governs the linear growth of a spatially periodic disturbance. Eigenvalues for this Rayleigh equation are determined numerically using proper conditions at the co-flowing gas and particle interface locations. For the first time, it is shown that non-uniform loading of small-inertia particles (Stokes number (St) <0.2) may destabilize the inviscid mixing layer development as compared to the pure-gas flow. The destabilization is triggered by an energy transfer rate that globally flows from the particle phase to the gas phase. For intermediate St (1 < St < 10), a maximum stabilizing effect is computed, while at larger St, two unstable modes may coexist. The growth rate computations from linear stability analysis are verified numerically through simulations based on an Eulerian-Lagrangian (EL) model based on the inviscid Euler equations and a point particle model. The growth rates found in numerical experiments using the EL method are in very good agreement with growth rates from the linear stability analysis and validate the destabilizing effect induced by the presence of particles with low St.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A numerical and analytical study of nonlinear bifurcations associated with the morphological stability of two-dimensional single crystals

The nonlinear stability of a two-dimensional single crystal of pure material in an undercooled melt is studied both analytically and numerically. The quasi-steady state approximation is used for the thermal fields, and the effects of different solid and liquid thermal conductivities, isotropic interfacial growth kinetics, and isotropic surface tension are included. The bifurcation analysis is performed by calculating the instantaneous value of the fundamental component of the local normal growth speed for an interface perturbed by a single Fourier shape component. Numerically, the fundamental component of the interfacial growth speed is found by Fourier analysis of the solution to an integrodifferential equation obeyed at the interface. Analytically, an expansion technique is used to derive a solvability condition defining each of these bifurcation points. The analytical and numerical results are in very close agreement. Almost all of the bifurcations are subcritical, and the results are presented by giving values of the Landau coefficient as a function of the different dimensionless parameters used in the model.

Brush, L. N.↗

Libration point staging options for SEI lunar missions - Station keeping implications

Libration point staging is one alternative to lunar orbit rendezvous which is being examined for the Space Exploration Initiative lunar missions. The libration point staging strategy being considered would eliminate the lunar mid-latitude accessibility constraint for long duration missions. Staging options for such a strategy are discussed with particular emphasis on the impact of the libration points' stability. The equations of motion for a satellite at the libration points are numerically integrated to investigate the stability and to determine the station-keeping requirements. These requirements are presented, and overall requirements for libration point staging briefly summarized. The results indicate that the slight eccentricity of the Moon's orbit would result in excessive station-keeping requirements for long duration missions.

Haynes, Davy A.↗

Numerical Simulation of an Enclosed Laminar Jet Diffusion Flame in Microgravity Environment: Comparison with ELF Data

Enclosed diffusion flames are commonly found in practical combustion systems, such as the power-plant combustor, gas turbine combustor, and jet engine after-burner. In these systems, fuel is injected into a duct with a co-flowing or cross-flowing air stream. The diffusion flame is found at the surface where the fuel jet and oxygen meet, react, and consume each other. In combustors, this flame is anchored at the burner (i.e., fuel jet inlet) unless adverse conditions cause the flame to lift off or blow out. Investigations of burner stability study the lift off, reattachment, and blow out of the flame. Flame stability is strongly dependent on the fuel jet velocity. When the fuel jet velocity is sufficiently low, the diffusion flame anchors at the burner rim. When the fuel jet velocity is increased, the flame base gradually moves downstream. However, when the fuel jet velocity increases beyond a critical value, the flame base abruptly jumps downstream. When this "jump" occurs, the flame is said to have reached its lift-off condition and the critical fuel jet velocity is called the lift-off velocity. While lifted, the flame is not attached to the burner and it appears to float in mid-air. Flow conditions are such that the flame cannot be maintained at the burner rim despite the presence of both fuel and oxygen. When the fuel jet velocity is further increased, the flame will eventually extinguish at its blowout condition. In contrast, if the fuel jet velocity of a lifted flame is reduced, the flame base moves upstream and abruptly returns to anchor at the burner rim. The fuel jet velocity at reattachment can be much lower than that at lift off, illustrating the hysteresis effect present in flame stability. Although there have been numerous studies of flame stability, the controlling mechanisms are not well understood. This uncertainty is described by Pitts in his review of various competing theories of lift off and blow out in turbulent jet diffusion flames. There has been some research on the stability of laminar flames, but most studies have focused on turbulent flames. It is also well known that the airflow around the fuel jet can significantly alter the lift off, reattachment and blow out of the jet diffusion flame. Buoyant convection is sufficiently strong in 1-g flames that it can dominate the flow-field, even at the burner rim. In normal-gravity testing, it is very difficult to delineate the effects of the forced airflow from those of the buoyancy-induced flow. Comparison of normal-gravity and microgravity flames provides clear indication of the influence of forced and buoyant flows on the flame stability. The overall goal of the Enclosed Laminar Flames (ELF) investigation (STS-87/USMP-4 Space Shuttle mission, November to December 1997) is to improve our understanding of the effects of buoyant convection on the structure and stability of co-flow diffusion flame, e.g., see http://zeta.lerc.nasa.gov/expr/elf.htm. The ELF hardware meets the experiment hardware limit of the 35-liter interior volume of the glovebox working area, and the 180x220-mm dimensions of the main door. The ELF experiment module is a miniature, fan-driven wind tunnel, equipped with a gas supply system. A 1.5-mm diameter nozzle is located on the duct's flow axis. The cross section of the duct is nominally a 76-mm square with rounded corners. The forced air velocity can be varied from about 0.2 to 0.9 m/s. The fuel flow can be set as high as 3 std. cubic centimeter (cc) per second, which corresponds to a nozzle exit velocity of up to 1.70 m/s. The ELF hardware and experimental procedure are discussed in detail in Brooker et al. The 1-g test results are repeated in several experiments following the STS-87 Mission. The ELF study is also relevant to practical systems because the momentum-dominated behavior of turbulent flames can be achieved in laminar flames in microgravity. The specific objectives of this paper are to evaluate the use reduced model for simulation of flame lift-off and blowout.

Jia, Kezhong↗