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 577 records · Page 32

All sources NASA NTRS DOE OSTI

NASA NTRS1977

A charging model for three-axis stabilized spacecraft

A charging model was developed for geosynchronous, three-axis stabilized spacecraft when under the influence of a geomagnetic substorm. The differential charging potentials between the thermally coated or blanketed outer surfaces and metallic structure of a spacecraft were determined when the spacecraft was immersed in a dense plasma cloud of energetic particles. The spacecraft-to-environment interaction was determined by representing the charged particle environment by equivalent current source forcing functions and by representing the spacecraft by its electrically equivalent circuit with respect to the plasma charging phenomenon. The charging model included a sun/earth/spacecraft orbit model that simulated the sum illumination conditions of the spacecraft outer surfaces throughout the orbital flight on a diurnal as well as a seasonal basis. Transient and steady-state numerical results for a three-axis stabilized spacecraft are presented.

Massaro, M. J.↗

NASA NTRS1990

Use of satellite data and modeling to assess the influence of stratospheric processes on the troposphere

Current research focused on the examination of the effects of ozone heating and Newtonian cooling on the linear stability of transient planetary waves. A finite difference numerical model was developed in which the basic state wind, temperature, and ozone distributions were constructed using ground-based and satellite data. Preliminary results indicate that ozone dynamics interaction may play a more important role in the stability and maintenance of planetary waves than previously thought, particularly in summer when the mean solar zenith angle is smallest and thus ozone heating largest. Examination of the finite amplitude interactions among radiation, ozone, and dynamics was carried out. Self-consistent, coupled equations governing the weakly nonlinear interactions are currently being analyzed to provide a better understanding of wave-mean flow interactions and ozone transport in a continuously stratified model of the troposphere-stratosphere coupled system.

Nathan, Terrence↗

NASA NTRS1993

Direct numerical simulation of hot jets

The ultimate motivation of this work is to investigate the stability of two dimensional heated jets and its implications for aerodynamic sound generation from data obtained with direct numerical simulations (DNS). As pointed out in our last report, these flows undergo two types of instabilities, convective or absolute, depending on their temperature. We also described the limits of earlier experimental and theoretical studies and explained why a numerical investigation could give us new insight into the physics of these instabilities. The aeroacoustical interest of these flows was also underlined. In order to reach this goal, we first need to succeed in the DNS of heated jets. Our past efforts have been focused on this issue which encountered several difficulties. Our numerical difficulties are directly related to the physical problem we want to investigate since these absolutely or almost absolutely unstable flows are by definition very sensitive to the smallest disturbances and are very likely to reach nonlinear saturation through a numerical feedback mechanism. As a result, it is very difficult to compute a steady laminar solution using a spatial DNS. A steady state was reached only for strongly co-flowed jets, but these flows are almost equivalent to two independent mixing layers. Thus they are far from absolute instability and have much lower growth rates.

Jacob, Marc C.↗

NASA NTRS2011

Effects of Nose Bluntness on Hypersonic Boundary-Layer Receptivity and Stability Over Cones

The receptivity to freestream acoustic disturbances and the stability properties of hypersonic boundary layers are numerically investigated for boundary-layer flows over a 5 straight cone at a freestream Mach number of 6.0. To compute the shock and the interaction of the shock with the instability waves, the Navier-Stokes equations in axisymmetric coordinates were solved. In the governing equations, inviscid and viscous flux vectors are discretized using a fifth-order accurate weighted-essentially-non-oscillatory scheme. A third-order accurate total-variation-diminishing Runge-Kutta scheme is employed for time integration. After the mean flow field is computed, disturbances are introduced at the upstream end of the computational domain. The appearance of instability waves near the nose region and the receptivity of the boundary layer with respect to slow mode acoustic waves are investigated. Computations confirm the stabilizing effect of nose bluntness and the role of the entropy layer in the delay of boundary-layer transition. The current solutions, compared with experimental observations and other computational results, exhibit good agreement.

Kara, Kursat↗

NASA NTRS2014

Entropy Stable Wall Boundary Conditions for the Compressible Navier-Stokes Equations

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite volume, finite difference, discontinuous Galerkin, and flux reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Parsani, Matteo↗

NASA NTRS1934

Supplemental data and calculations of the lateral stability of airplanes

In connection with the DVL Report 272 on the theory of the lateral stability of airplanes, the formal results are here amplified in some respects and their technical significance again briefly explained. Three numerical examples show how model tests for checking the lateral stability are to be evaluated and supplemented, if necessary, and how the stability limits depend on the design of the airplane and on the conditions of flight.

Mathias, Gotthold↗

NASA NTRS2010

Harmonic Balance Computations of Fan Aeroelastic Stability

A harmonic balance (HB) aeroelastic analysis, which has been recently developed, was used to determine the aeroelastic stability (flutter) characteristics of an experimental fan. To assess the numerical accuracy of this HB aeroelastic analysis, a time-domain aeroelastic analysis was also used to determine the aeroelastic stability characteristics of the same fan. Both of these three-dimensional analysis codes model the unsteady flowfield due to blade vibrations using the Reynolds-averaged Navier-Stokes (RANS) equations. In the HB analysis, the unsteady flow equations are converted to a HB form and solved using a pseudo-time marching method. In the time-domain analysis, the unsteady flow equations are solved using an implicit time-marching approach. Steady and unsteady computations for two vibration modes were carried out at two rotational speeds: 100 percent (design) and 70 percent (part-speed). The steady and unsteady results obtained from the two analysis methods compare well, thus verifying the recently developed HB aeroelastic analysis. Based on the results, the experimental fan was found to have no aeroelastic instability (flutter) at the conditions examined in this study.

Bakhle, Milind A.↗

NASA NTRS1990

Efficient solution of parabolic equations by Krylov approximation methods

Numerical techniques for solving parabolic equations by the method of lines is addressed. The main motivation for the proposed approach is the possibility of exploiting a high degree of parallelism in a simple manner. The basic idea of the method is to approximate the action of the evolution operator on a given state vector by means of a projection process onto a Krylov subspace. Thus, the resulting approximation consists of applying an evolution operator of a very small dimension to a known vector which is, in turn, computed accurately by exploiting well-known rational approximations to the exponential. Because the rational approximation is only applied to a small matrix, the only operations required with the original large matrix are matrix-by-vector multiplications, and as a result the algorithm can easily be parallelized and vectorized. Some relevant approximation and stability issues are discussed. We present some numerical experiments with the method and compare its performance with a few explicit and implicit algorithms.

Gallopoulos, E.↗

NASA NTRS2015

Entropy Stable Wall Boundary Conditions for the Three-Dimensional Compressible Navier-Stokes Equations

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators on unstructured grids are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite difference, finite volume, discontinuous Galerkin, and flux reconstruction/correction procedure via reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Parsani, Matteo↗

NASA NTRS2000

The Generation of Smooth High Speed Solar Wind From Plume-Interplume Mixing

Ulysses has shown that fast solar wind is extremely smooth, with a variance of less than 5%, in contrast to slow wind with a variance of approximately 30%. Now UVCS has produced the surprising result that the flow speed within coronal holes, the source of fast wind, is not at all smooth. Specifically, Giordano et al. (ApJ, v531, L79-L82, 2000) report that at 1.7 R(sub SUN) the interplume flow speed is typically more than twice the plume flow speed. Other less direct evidence supports this same result, with speeds from less than 300 to over 1000 km/s reported at approximately 5 R(sub SUN). This presents the paradox of how strongly differing plume and interplume flow speeds can exist near the Sun and be absent far from the Sun. The only answer is that plume and interplume material or momentum must be strongly mixed and that the mixing must occur mainly inside 0.3 AU to be consistent with Helios observations of smooth fast wind. Pressure balance structures (PBSs) and He abundance anomalies (Reisenfeld, et al., GRL, v26(13), 1805-1808, 1999) have been identified as interplanetary remnants of plumes, implying momentum mixing is the dominant coronal process. One possible source for plume/interplume momentum mixing is MHD Kelvin-Helmholtz (KH) shear instabilities occurring on the velocity shear interfaces. The velocity shear is a source of free energy and KH fluctuations could, through nonlinear cascade, provide the forcing required for the plasma oscillations (Cranmer, ApJ, v532, 1197-1208, 2000) reported to exist in coronal holes. The physical properties in coronal holes are now sufficiently well known that we can show plume/interplume shear interfaces become unstable to the KH instability at 5 - 10 R(sub SUN). The KH dispersion relation can be used to analyze marginal stability, the most unstable wavelengths, and linear growth rates. Numerical simulations can be used to verify results from the linear analysis and study the nonlinear development of KH modes. Here we will describe the marginal stability criterion for the KH instability, how and where this condition is met in the corona, and the possible character of the resultant fluctuations.

Suess, Steve↗

NASA NTRS1986

Entropy conservative finite element schemes

The question of entropy stability for discrete approximations to hyperbolic systems of conservation laws is studied. The amount of numerical viscosity present in such schemes is quantified and related to their entropy stability by means of comparison. To this end, two main ingredients are used: entropy variables and the construction of certain entropy conservative schemes in terms of piecewise-linear finite element approximations. It is then shown that conservative schemes are entropy stable, if and (for three-point schemes) only if, they contain more numerical viscosity than the abovementioned entropy conservation ones.

Tadmor, E.↗

NASA NTRS2022

Interaction of a Tunnel-like Acoustic Disturbance Field with a Blunt Cone Boundary Layer at Mach 8

The existing measurements of laminar-to-turbulent transition over circular cones in conventional (i.e., “noisy”) hypersonic wind tunnels have established that the transition location moves downstream when the nose radius is increased from zero. However, this initially downstream movement slows down and ultimately reverses beyond a critical value of the nose radius, and may be related to external forcing in the form of freestream disturbances and/or surface roughness. To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, hypersonic boundary-layer receptivity to broadband freestream acoustic disturbances from the nozzle wall of a digital conventional wind tunnel is investigated by both direct numerical simulations (DNS) and modal and nonmodal stability analysis. A Mach 8 flow over a 7 deg half-angle cone with a nose radius of 𝑅 𝑛 = 5.2mm and freestream Reynolds number of 12.2 × 10 6 m -1 is considered. The results show that the broadband tunnel noise in the free stream of a convectional hypersonic wind tunnel (i.e., outside of the nozzle-wall turbulent boundary layer) can be well represented by an acoustic model with an ansatz of slow acoustic waves. With successful calibration of the model parameters against the precursor tunnel DNS, such an acoustic ansatz can successfully reproduce both the frequency-wavenumber spectra and the temporal evolution of the broadband tunnel noise radiated from the nozzle wall. Additionally, the DNS of the Mach 8 blunt cone with tunnel-like acoustic input above the bow shock showed that the spectra of wall-pressure and heat-transfer fluctuations recovers the signature of the axisymmetric waves predicted by the nonmodal analysis. Furthermore, the azimuthal wavenumber and frequency spectrum of the temperature fluctuations as a function of the wall-normal distance show higher amplitudes for three-dimensional waves above the boundary-layer edge. The numerical schlieren contours show the inclined structures commonly observed in blunt cone experiments, demonstrating that they correspond to three-dimensional structures due to freestream disturbances in the presence of an entropy layer.

Boundary-layer transition↗

NASA NTRS2022

Interaction of a Tunnel-like Acoustic Disturbance Field with a Blunt Cone Boundary Layer at Mach 8

Boundary layer transition↗

NASA NTRS1987

Numerical tests of a modified full implicit continuous Eulerian (FICE) scheme with projected normal characteristic boundary conditions for MHD flows

A numerical method has been developed based on a modified full implicit continuous Eulerian (FICE) scheme and projected normal characteristic boundary conditions for simulating MHD flows which undergo a long process of evolution. An astrophysical flow is chosen for illustration of this procedure, and numerical tests are made to verify the computational stability and physically realistic solution. Three computational tests have been accomplished; they are tests of solving methods, characteristic boundary condition, and time steps. The tests show that the program from the modified FICE scheme with proper boundary conditions and time steps can be made numerically stable for a time long enough to obtain physically plausible solutions.

Wu, S. T.↗

NASA NTRS1992

On the instability of boundary layers on heated flat plates

The stability of a boundary layer on a heated flat plate is investigated in the linear regime. The flow is shown to be unstable to longitudinal vortex structures which in general develop in a nonparallel manner in the streamwise direction. Solutions of the nonparallel equations are obtained numerically at O(1) values of the appropriate stability parameter, i.e., the Grashof number. The particular cases investigated relate to the situations when the instability is induced by localized or distributed wall roughness or nonuniform wall heating. The case when the vortices are induced by freestream disturbances is also considered. The fastest growing mode is found to be governed by a quasi-parallel theory at high wavenumbers. The wavenumber and growth rate of the fastest growing mode are found in closed form. At low wavenumbers the vortex instability is shown to be closely related to Tollmein-Schlichting waves. The effect of wall heating or cooling on the latter type of instability is discussed.

Hall, Philip↗

NASA NTRS1934

A study of the pitching moments and the stability characteristics of monoplanes

This note presents a study of the pitching moments and the stability characteristics of monoplanes. Expressions for the pitching-moment coefficient and the Diehl stability coefficient for the monoplane are developed, suitable for the use of airplane designers. The effective difference between the high-wing and low-wing types is portrayed and discussed. Comparisons between experimental and computed values are made. Charts for use in the solution of numerical values of the pitching-moment and stability coefficients are presented.

Higgins, George J↗

NASA NTRS2012

Large-Eddy / Reynolds-Averaged Navier-Stokes Simulations of a Dual-Mode Scramjet Combustor

Numerical simulations of reacting and non-reacting flows within a scramjet combustor configuration experimentally mapped at the University of Virginia s Scramjet Combustion Facility (operating with Configuration A ) are described in this paper. Reynolds-Averaged Navier-Stokes (RANS) and hybrid Large Eddy Simulation / Reynolds-Averaged Navier-Stokes (LES / RANS) methods are utilized, with the intent of comparing essentially blind predictions with results from non-intrusive flow-field measurement methods including coherent anti-Stokes Raman spectroscopy (CARS), hydroxyl radical planar laser-induced fluorescence (OH-PLIF), stereoscopic particle image velocimetry (SPIV), wavelength modulation spectroscopy (WMS), and focusing Schlieren. NC State's REACTMB solver was used both for RANS and LES / RANS, along with a 9-species, 19- reaction H2-air kinetics mechanism by Jachimowski. Inviscid fluxes were evaluated using Edwards LDFSS flux-splitting scheme, and the Menter BSL turbulence model was utilized in both full-domain RANS simulations and as the unsteady RANS portion of the LES / RANS closure. Simulations were executed and compared with experiment at two equivalence ratios, PHI = 0.17 and PHI = 0.34. Results show that the PHI = 0.17 flame is hotter near the injector while the PHI = 0.34 flame is displaced further downstream in the combustor, though it is still anchored to the injector. Reactant mixing was predicted to be much better at the lower equivalence ratio. The LES / RANS model appears to predict lower overall heat release compared to RANS (at least for PHI = 0.17), and its capability to capture the direct effects of larger turbulent eddies leads to much better predictions of reactant mixing and combustion in the flame stabilization region downstream of the fuel injector. Numerical results from the LES/RANS model also show very good agreement with OH-PLIF and SPIV measurements. An un-damped long-wave oscillation of the pre-combustion shock train, which caused convergence problems in some RANS simulations, was also captured in LES / RANS simulations, which were able to accommodate its effects accurately.

Fulton, Jesse A.↗

NASA NTRS2003

Finite Element Method for Capturing Ultra-relativistic Shocks

While finite element methods are used extensively by researchers solving computational fluid dynamics in fields other than astrophysics, their use in astrophysical fluid simulations has been predominantly overlooked. Current simulations using other methods such as finite difference and finite volume (based on finite difference) have shown remarkable results, but these methods are limited by their fundamental properties in aspects that are important for simulations with complex geometries and widely varying spatial and temporal scale differences. We have explored the use of finite element methods for astrophysical fluids in order to establish the validity of using such methods in astrophysical environments. We present our numerical technique applied to solving ultra-relativistic (Lorentz Factor Gamma >> 1) shocks which are prevalent in astrophysical studies including relativistic jets and gamma-ray burst studies. We show our finite element formulation applied to simulations where the Lorentz factor ranges up to 2236 and demonstrate its stability in solving ultra-relativistic flows. Our numerical method is based on the Flowfield Dependent Variation (FDV) Method, unique in that numerical diffusion is derived from physical parameters rather than traditional artificial viscosity methods. Numerical instabilities account for most of the difficulties when capturing shocks in this regime. Our method results in stable solutions and accurate results as compared with other methods.

Richardson, G. A.↗