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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Entropy-Based Approach To Nonlinear Stability

NASA technical memorandum suggests schemes for numerical solution of differential equations of flow made more accurate and robust by invoking second law of thermodynamics. Proposes instead of using artificial viscosity to suppress such unphysical solutions as spurious numerical oscillations and nonlinear instabilities, one should formulate equations so that rate of production of entropy within each cell of computational grid be nonnegative, as required by second law.

Merriam, Marshal L.↗
NASA NTRS1982

Theoretical analyses of Baroclinic flows

A stability analysis of a thin horizontal rotating fluid layer which is subjected to arbitrary horizontal and vertical temperature gradients is presented. The basic state is a nonlinear Hadley cell which contains both Ekman and thermal boundary layers; it is given in closed form. The stability analysis is based on the linearized Navier-Stokes equations, and zonally symmetric perturbations in the form of waves propagating in the meridional direction are considered. Numerical methods were used for the stability problem. It was found that the instability sets in when the Richardson number is close to unity and that the critical Richardson number is a non-monotonic function of the Prandtl number. Further, it was found that the critical Richardson number decreases with increasing Ekman number until a critical value of the Ekman number is reached beyond which the fluid is stable.

Antar, B.↗
NASA NTRS1998

Numerical Simulations of Unsteady Low-Reynolds-Number Flows Over the APEX Airfoil

Laminar and transitional separation bubbles are an important feature of low-Reynolds-number flows over airfoils. The separation bubbles are unsteady and have a significant impact on the aerodynamic properties of the airfoils. In this paper unsteady low-Reynolds-number separated flows over the APEX airfoil are calculated using a Navier-Stokes solver. The numerical results show the presence of unsteady separation bubbles in the flowfield. An analysis of the numerical results shows that flowfield disturbances are amplified significantly in the separation bubble, leading to periodic vortex shedding. A linear stability analysis of the separated boundary layer is performed and the results show that the dominant wavenumber and frequency in the numerical simulations agree with the most unstable wavenumber and frequency from the linear stability analysis. The numerical results also show the growth and interaction of disturbance waves in the separation bubble. For transonic flows over the APEX airfoil, the calculations show that the presence of shocks causes significant changes in the separation location and consequently, the overall flowfield.

Tatineni, Mahidhar↗
DOE OSTI2021

Elucidating the Redox Behavior during Atomic Layer Deposition on Lithium-Ion Battery Cathode Materials

Atomic layer deposition (ALD) has emerged as a promising technology for applying ultrathin protective coatings on lithium-ion battery (LIB) cathode surfaces to improve their cycling stability. While there have been numerous reports evaluating the electrochemical performance of these surface-modified cathode materials, the chemical changes induced on the surface of the cathode materials upon ALD coating are not fully studied. This paper describes a systematic investigation to understand the interfacial changes of 12 different cathode materials upon coating with aluminum oxide (Al 2 O 3 ) using trimethyl aluminum (TMA) and H 2 O, and aluminum fluoride (AlF 3 ) using TMA and hydrogen fluoride pyridine (HFPy). Here, the surface composition of these cathode materials, which range from simple transition metal oxides (e.g. NiO and MnO) to complex multi-element cathode materials (e.g. LiNi x Mn 1-x-y Co y O 2 , NMC), was studied via X-ray photoelectron spectroscopy (XPS). The XPS measurements reveal that the transition metals in the cathode materials undergo selective oxidation/reduction depending upon the nature of the precursor, the coating, and the cathode material. The ability to chemically modify the surface of cathode materials via vapor-phase precursor adsorption will open new avenues to systematically control the interface of the cathode materials in LIBs that are not possible by conventional coating methods.

25 ENERGY STORAGE↗
DOE OSTI2026

Boundary-layer receptivity to oblique freestream vorticity waves for a high-enthalpy hypersonic flow

The receptivity of a Mach 15 straight-cone boundary layer to oblique freestream vorticity waves is investigated using direct numerical simulation (DNS) alongside linear stability theory and the linear parabolized stability equations. A thermochemical nonequilibrium gas model is used. Oblique freestream vorticity waves at frequencies of 400, 800, and 1200 kHz are considered, with incident angles ranging from 0° to 29.4° at 400 kHz, 0° to 15.7° at 800 kHz, and 0° to 20.6° at 1200 kHz. The 400 kHz case is of primary interest due to the strong second-mode amplification at this frequency. Although the underlying base flow is axisymmetric, the oblique vorticity waves lead to a fully three-dimensional boundary-layer disturbance whose characteristics vary depending on the azimuthal ray relative to the freestream wave. Moving downstream within the second-mode instability region, some clear trends emerge in terms of the boundary-layer disturbance amplitudes; that is, disturbance amplitudes are highest at the leeward ray (relative to the freestream wave), but weakest about halfway between the windward and leeward rays. Increasing the incident angle causes the amplitudes to increase on the leeward ray and decrease on the windward ray. Moreover, the boundary-layer disturbance throughout contains a wide spectrum of azimuthal wavenumbers in which the disturbance energy falls off at higher wavenumbers. Increasing the incident angle causes the azimuthal spectrum of the boundary-layer disturbance to broaden overall. Qualitatively similar results are found for the two higher frequencies leading up to the peak-amplitude locations corresponding to the second-mode instability.

42 ENGINEERING↗
DOE OSTI2022

Universal Battery Supercharger

During the course of the UBS project a dc-fast charger was developed to be used in off board high power charging stations. High efficiency is targeted to reduce the charging price for the consumer while reducing the charge time. In the meantime, the size of the converter is minimized so that it could be used in small spaces and potentially urban areas. The topology is modular, and the power level can be expanded for charging higher power vehicles like trucks and buses. During the course of this research several key technologies was developed in LESES that improve the performance of the power electronics system used in EV charger. Integrated magnetics was developed for DMCR to improve the power density as well as efficiency of the magnetic components which is one of the main bottlenecks of EV development. Triangular mode EV current conduction mode is using to perform soft switching (increase efficiency) without adding any new components to the existing topology. Passive and active damping are developed to suppress the interaction between grid and converter, stabilizing the control system. Further numerous other technical problems have been solved to pave the way for even higher performance converters in the future.

25 ENERGY STORAGE↗
NASA NTRS1975

Second-order non-iterative ADI solution of non-linear partial differential equations

A new method for the solution of non-linear partial differential equations by an ADI procedure is described. Although the method is second order accurate in time, it does not require either iterations or predictor corrector methods to overcome the nonlinearity of the equations. Thus the computational effort required for the solution of the non-linear problem becomes similar to that required for the linear case. The method is applied to a two-dimensional 'extended Burgers equation'. Linear stability is studied, and some numerical solutions obtained. The improved accuracy obtained by the 2nd order truncation error is clearly manifested.

Wolfshtein, M.↗
NASA NTRS1976

The development of the Poincare-similar elements with eccentric anomaly as the independent variable

A new set of element differential equations for the perturbed two-body motion is derived. The elements are canonical and are similar to the classical canonical Poincare elements, which have time as the independent variable. The phase space is extended by introducing the total energy and time as canonically conjugated variables. The new independent variable is, to within an additive constant, the eccentric anomaly. These elements are compared to the Kustaanheimo-Stiefel (KS) element differential equations, which also have the eccentric anomaly as the independent variable. For several numerical examples, the accuracy and stability of the new set are equal to those of the KS solution. This comparable accuracy result can probably be attributed to the fact that both sets have the same time element and very similar energy elements. The new set has only 8 elements, compared to 10 elements for the KS set. Both sets are free from singularities due to vanishing eccentricity and inclination.

Bond, V. R.↗
NASA NTRS1983

Efficient linear and nonlinear heat conduction with a quadrilateral element

A method is presented for performing efficient and stable finite element calculations of heat conduction with quadrilaterals using one-point quadrature. The stability in space is obtained by using a stabilization matrix which is orthogonal to all linear fields and its magnitude is determined by a stabilization parameter. It is shown that the accuracy is almost independent of the value of the stabilization parameter over a wide range of values; in fact, the values 3, 2, and 1 for the normalized stabilization parameter lead to the 5-point, 9-point finite difference, and fully integrated finite element operators, respectively, for rectangular meshes and have identical rates of convergence in the L2 norm. Eigenvalues of the element matrices, which are needed for stability limits, are also given. Numerical applications are used to show that the method yields accurate solutions with large increases in efficiency, particularly in nonlinear problems.

Liu, W. K.↗
NASA NTRS1983

Convenient stability criteria for difference approximations of hyperbolic initial-boundary value problems

New convenient stability criteria are provided in this paper for a large class of finite difference approximations to initial-boundary value problems associated with the hyperbolic system u sub t = au sub x + Bu + f in the quarter plane x or = 0, t or = 0. Using the new criteria, stability is easily established for numerous combinations of well known basic schemes and boundary conditions, thus generalizing many special cases studied in recent literature.

Goldberg, M.↗
NASA NTRS1983

Barotropic instability of weakly non-parallel zonal flows

Integrations of the linear stability problem, in the present numerical investigation of weakly nonparallel zonal flows barotropic instability having localized intense shear regions, reveal the existence of unstable localized wave packets. The spatial structure and eigenfrequencies of these packets depend on two parameters measuring the degree of supercriticality and the zonal length scale of the shear region. It is found that instability structure is defined by conditions ensuring the decay of the wave packet at infinity, and the transition from long to short waves across a turning point region that is controlled by nonparallel effects.

Merkine, L.-O.↗
NASA NTRS1984

The propagation and stability of time-dependent galactodetonation waves

A new analytic approach is applied to the problem of star formation fronts in differentially rotating disks recently studied by Cowie and Rybicki. A quasi-linear equation for the inclination angle of the front is derived, and its characteristic equations are solved. The stability of these fronts found numerically in the previous study is revealed to be a direct consequence of the structure of the front characteristics. The time-dependent solutions found here are potentially useful in describing the evolution of local wave fronts shearing in a galactic velocity field, as observed in spirals lacking 'grand-design' structure.

Balbus, S. A.↗
NASA NTRS1985

Convenient stability criteria for difference approximations of hyperbolic initial-boundary value problems

New convenient stability criteria are provided in this paper for a large class of finite difference approximations to initial-boundary value problems associated with the hyperbolic system u sub t = au sub x + Bu + f in the quarter plane x or = 0, t or = 0. Using the new criteria, stability is easily established for numerous combinations of well known basic schemes and boundary conditins, thus generalizing many special cases studied in recent literature.

Goldberg, M.↗
NASA NTRS1985

Magnetohydrodynamic thermal instabilities in cool inhomogeneous atmospheres

The stability of magnetic loops to current-driven filamentation instabilities is investigated. The unperturbed atmosphere is assumed to be composed of an (upper) isothermal optically thin low-density portion and a (lower) higher-density portion which is in radiative equilibrium; in both cases, the atmosphere is in hydrostatic equilibrium, so that gravitational stratification is taken into account. In order to provide specific equilibrium conditions for evaluation of the dispersion relation, conditions appropriate for the surface of a solar-type star are adopted; i.e., a fairly low temperature (T = 5000 K) appropriate for a 'precoronal' state associated, for example, with magnetic flux emerging from photospheric levels under the action of magnetic buoyancy. A linear stability analysis is performed, and numerical results show that physically plausible current densities, which would be generated by typical loop-footpoint motions, are effective in driving MHD instabilities in such a plasma. The instability growth rates are strongly dependent on the assumed current density distribution and on the density scale height.

Bodo, G.↗
NASA NTRS1985

Dynamics of liquid sloshing in upright and inverted bladdered tanks

A system identification methodology was used to examine the dynamics of liquid sloshing in the upright and inverted bladdered hydrazine tanks of the Tracking and Data Relay Satellite, (TDRS) and to evaluate the effects of bladder stiffness on the sloshing parameters. Mechanical models of the two systems were developed using the numerical values derived from static stability tests and from slosh frequency response tests of a full-size model tank fitted with a prototype bladder. For the upright tank (liquid below the bladder) a modified conventional pendulum was used. In the inverted tank (liquid above the bladder) where sloshing is unconventional due to the highly nonsymmetrical orientation of the liquid held by the bladder, a mechanical model using an inverted pendulum which is able to undergo small oscillations as well as large reorientations was necessary. Both thrusting and low-gravity conditions are considered.

Dodge, F. T.↗
NASA NTRS1987

Aeroelastic analysis of launch vehicles in transonic flight

A numerical study of the aeroelastic stability of launch vehicles in transonic flight is performed. The unsteady aerodynanic flowfield about the vehicle is simulated using the three dimensional Reynolds-Averaged Navier-Stokes equations with the account for turbulent mixing done by a two-layer algebraic eddy viscosity model. The structural-dynamic formulation considers free-free flexural vibration of an elongated beam with variable properties and the equations are cast in modal form. The two sets of equations are integrated simultaneously in time, after a initial perturbation, to ascertain the aeroelastic stability of a given configuration. The method propossed is described in detail, and results that illustrate its application are presented.

Azevedo, Joao Luiz F.↗
NASA NTRS1988

A study of numerical methods for hyperbolic conservation laws with stiff source terms

The proper modeling of nonequilibrium gas dynamics is required in certain regimes of hypersonic flow. For inviscid flow this gives a system of conservation laws coupled with source terms representing the chemistry. Often a wide range of time scales is present in the problem, leading to numerical difficulties as in stiff systems of ordinary differential equations. Stability can be achieved by using implicit methods, but other numerical difficulties are observed. The behavior of typical numerical methods on a simple advection equation with a parameter-dependent source term was studied. Two approaches to incorporate the source term were utilized: MacCormack type predictor-corrector methods with flux limiters, and splitting methods in which the fluid dynamics and chemistry are handled in separate steps. Various comparisons over a wide range of parameter values were made. In the stiff case where the solution contains discontinuities, incorrect numerical propagation speeds are observed with all of the methods considered. This phenomenon is studied and explained.

Leveque, R. J.↗