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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

An efficient procedure for cascade aeroelastic stability determination using nonlinear, time-marching aerodynamic solvers

A numerical eigenvalue problem formulation and a practical calculation procedure for exact eigenvalues and corresponding eigenvectors are developed and applied to a nonlinear, two-dimensional, time-marching full potential solver for cascade aeroelastic stability analysis. This procedure is based on the Lanczos recursive method and it directly calculates stability information about a nonlinear steady state. It is compared to conventional approaches in the frequency and time domains developed earlier and is found to be 100-10.000 times more computationally efficient. Eigenvalue constellations and the flutter results for flow through a cascade SR5 propfan airfoil are presented.

Mahajan, Aparajit J.↗
DOE OSTI2015

Solitary Waves of a $\mathcal {P}$ $\mathcal {T}$ -Symmetric Nonlinear Dirac Equation

In our study we consider we consider a prototypical example of a mathcalP mathcalT-symmetric Dirac model. We discuss the underlying linear limit of the model and identify the threshold of the mathcalP mathcalT -phase transition in an analytical form. We then focus on the examination of the nonlinear model. We consider the continuation in the mathcalP mathcalT -symmetric model of the solutions of the corresponding Hamiltonian model and find that the solutions can be continued robustly as stable ones all the way up to the mathcalP mathcalT-transition threshold. In the latter, they degenerate into linear waves. We also examine the dynamics of the model. Given the stability of the waveforms in the mathcalP mathcalT-exact phase, we consider them as initial conditions for parameters outside of that phase. We also find that both oscillatory dynamics and exponential growth may arise, depending on the size of the corresponding “quench”. The former can be characterized by an interesting form of bifrequency solutions that have been predicted on the basis of the SU symmetry. Finally, we explore some special, analytically tractable, but not mathcalP mathcalT-symmetric solutions in the massless limit of t- e model.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

Stability and transition in supersonic boundary layers

The full three-dimensional time-dependent compressible Navier-Stokes equations are numerically solved by a Fourier-Chebyshev collocation algorithm to study the stability of supersonic flows over a flat plate. Several non-linear numerical experiments suggest the existence of a secondary instability which might provide a possible route to transition. The interaction of the modes involved in this secondary instability is possibly amenable to a Floquet theory. Pertinent differences between this instability and the more common incompressible K-type instabilities are pointed out.

Erlebacher, Gordon↗
NASA NTRS1994

The stabilizing effect of compressibility in turbulent shear flow

Direct numerical simulation of turbulent homogeneous shear flow is performed in order to clarify compressibility effects on the turbulence growth in the flow. The two Mach numbers relevant to homogeneous shear flow are the turbulent Mach number M(t) and the gradient Mach number M(g). Two series of simulations are performed where the initial values of M(g) and M(t) are increased separately. The growth rate of turbulent kinetic energy is observed to decrease in both series of simulations. This 'stabilizing' effect of compressibility on the turbulent energy growth rate is observed to be substantially larger in the DNS series where the initial value of M(g) is changed. A systematic companion of the different DNS cues shows that the compressibility effect of reduced turbulent energy growth rate is primarily due to the reduced level of turbulence production and not due to explicit dilatational effects. The reduced turbulence production is not a mean density effect since the mean density remains constant in compressible homogeneous shear flow. The stabilizing effect of compressibility on the turbulence growth is observed to increase with the gradient Mach number M(g) in the homogeneous shear flow DNS. Estimates of M(g) for the mixing and the boundary layer are obtained. These estimates show that the parameter M(g) becomes much larger in the high-speed mixing layer relative to the high-speed boundary layer even though the mean flow Mach numbers are the same in the two flows. Therefore, the inhibition of turbulent energy production and consequent 'stabilizing' effect of compressibility on the turbulence (over and above that due to the mean density variation) is expected to be larger in the mixing layer relative to the boundary layer in agreement with experimental observations.

Sarkar, S.↗
NASA NTRS1994

Coupled Marangoni-Benard/Rayleigh-Benard Instability with Temperature Dependent Viscosity

The onset of convection induced by coupled surface tension gradient and buoyancy forces is investigated with temperature dependent viscosity. Both surface tension and viscosity are assumed to vary linearly with temperature. The limiting case, Ma=0, is the buoyancy driven convection problem typically referred to as the Rayleigh-Benard stability problem. The other limiting case, Ra=0, is the surface tension (gradient) driven flow problem referred to as the Marangoni-Benard problem. The equations and boundary conditions obtained from the linear analysis are solved numerically as a generalized eigenvalue problem. Neutral stability curves for different viscosity slopes have been generated for the Marangoni-Benard and Rayleigh-Benard problems. It is shown that the curves can be collapsed to a single curve by appropriately scaling the results for each of the limiting cases. The critical Marangoni number is determined as a function of the slope of the viscosity temperature variation, epsilon, for different values of the Rayleigh number. When the viscosity decreases linearly with temperature, the coupled buoyancy-surface tension problem, including the limiting cases of Ra=0 and Ma=0, is found to be more stable than the constant viscosity case.

Skarda, J. Raymond Lee↗
NASA NTRS1971

Nutational stability of a dual-spin satellite under the influence of applied reaction torques

It is assumed that since the solar paddle attachments to the hub of the spacecraft are not rigidly locked, the effect of the solar panels can be replaced by a constant reaction torque acting on the hub of the spacecraft. This could result in the satellite having an equilibrium motion about an axis displaced from the nominal axis of symmetry. The variational equations of motion are developed about such an equilibrium position using the SAS-A spacecraft as a model. Energy dissipation on the rotor as well as the main body is included. This nonautonomous set of differential equations are linearized and transformed to an autonomous set using the Liapunov reducibility theorem. The stability of the kinematically similar system is examined numerically using representative SAS-A parameters for the case when either pair of solar panels is assumed to be loosely attached. Stability is verified for small system nutation angles (0.1 degree).

Bainum, P. M.↗
NASA NTRS1972

Nutational stability of a dual-spin satellite under the influence of applied reaction torques.

The nutational stability of a dual-spin satellite under the influence of an applied reaction torque is analyzed. It is assumed that since the solar paddle attachments to the hub of the spacecraft are not rigidly locked, the effect of the solar panels can be replaced by a constant reaction torque acting on the hub of the spacecraft. This could result in the satellite having an equilibrium motion about an axis displaced from the nominal axis of symmetry. The variational equations of motion are developed about such an equilibrium position using the SAS-A spacecraft as a model. Energy dissipation on the rotor as well as the main body is included. This nonautonomous set of differential equations are linearized and transformed to an autonomous set using the Lyapunov Reducibility Theorem. The stability of the kinematically similar system is examined numerically using representative SAS-A parameters for the case when either pair of solar panels is assumed to be loosely attached. Stability is verified for small system nutation angles (0.1 degree) although the time constant associated with the least damped mode is approximately one order of magnitude larger than that for the SAS-A system without the influence of the reaction torques.

Bainum, P. M.↗
NASA NTRS1994

A Numerical Comparison of Barrier and Modified Barrier Methods for Large-Scale Bound-Constrained Optimization

When a classical barrier method is applied to the solution of a nonlinear programming problem with inequality constraints, the Hessian matrix of the barrier function becomes increasingly ill-conditioned as the solution is approached. As a result, it may be desirable to consider alternative numerical algorithms. We compare the performance of two methods motivated by barrier functions. The first is a stabilized form of the classical barrier method, where a numerically stable approximation to the Newton direction is used when the barrier parameter is small. The second is a modified barrier method where a barrier function is applied to a shifted form of the problem, and the resulting barrier terms are scaled by estimates of the optimal Lagrange multipliers. The condition number of the Hessian matrix of the resulting modified barrier function remains bounded as the solution to the constrained optimization problem is approached. Both of these techniques can be used in the context of a truncated-Newton method, and hence can be applied to large problems, as well as on parallel computers. In this paper, both techniques are applied to problems with bound constraints and we compare their practical behavior.

Nash, Stephen G.↗
DOE OSTI2021

Ideal internal kink stability in presence of plasma flow and neoclassical toroidal viscosity due to energetic particles

Abstract The influence of energetic particles (EPs) on the ideal internal kink mode, in rotating tokamak plasmas, is numerically investigated by simultaneously solving MHD-kinetic hybrid equations together with a toroidal momentum balance equation utilizing the MARS-Q code (Liu et al 2013 Phys. Plasmas 20 042503). The neoclassical toroidal viscous (NTV) torque, induced by precessional drift resonances of trapped energetic particles, acts as the momentum sink term to damp the plasma flow. Quasi-linear initial value simulations show local reduction of the flow amplitude and enhancement of the flow shear near the q = 1 rational surface ( q is the safety factor) due to EP induced NTV. Both effects in turn destabilize the internal kink mode. These numerical findings are robust against the initial linear stability of internal kink, the initial plasma flow profile, as well as the equilibrium distribution model for EPs.

Physics↗
DOE OSTI2018

Stability and bifurcation analysis of Japanese encephalitis model with/without effects of some control parameters

In this paper, a mathematical model on transmission of Japanese Encephalitis disease has been developed considering some control parameters and time-dependent environmental-carrying capacity. Here, the total vector population is divided into two subpopulations such as susceptible mosquito and infected mosquito. Here also, total reservoir population (i.e., the population in which the encephalitis virus grows) such as pig, horse, etc., has been considered which is divided into three subpopulations such as susceptible reservoir, infected reservoir and recovered reservoir. Total human population is also divided into three subpopulations such as susceptible human, infected human and recovered human. The dynamical behaviors of the system have been investigated. Here, the basic reproduction number associated with the system has been analyzed with respect to control parameters both theoretically and numerically. The biological feasible equilibria and their stability properties have been discussed and the existence condition of the disease has been illustrated numerically. For a certain set of parametric values, effectiveness of control parameters of our proposed model has been checked numerically. At last, Hopf bifurcations have been made numerically without considering control parameters for the case of constant environmental-carrying capacity of mosquito.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2021

Wake Instability Behind Isolated Trip Near the Leading Edge of the BOLT-II Configuration

The BOLT-II (Holden Mission) configuration is an extended version of the BOLT flight article and depicts hypersonic boundary-layer transition in the presence of multiple and potentially interacting instability mechanisms. Several numerical studies of the boundary-layer instabilities over these configurations have been reported in the recent literature, including our previous studies of the modal instability characteristics of boundary-layer streaks adjacent to the minor-axis symmetry plane of the BOLT configuration and the wake instabilities behind a diamond planform (“pizza-box”) trip along the symmetry plane on the secondary side of the BOLT-II configuration. The present work extends the latter study to a scaled version of the same trip that is located in the region of nonzero crossflow in the vicinity of the leading edge at X/L = 0.5. Collectively, the trips along the symmetry plane and near the leading edge are the focus of the NASA roughness experiment on the secondary side of the BOLT-II configuration. The laminar basic state computation at the nominal flight design condition of Re ∞ = 5.44 x 10 6 /m and Re ∞ = 2.5 x 10 6 /m shows that the leading-edge trip with k/δ ≈ 0.70 and planform-side-length-to-height ratio of b/k = 3.0 induces multiple asymmetric, longitudinal streaks within the trip wake. The most prominent streak among these resembles a finite amplitude crossflow vortex, and it supports the amplification of multiple families of unstable modes. The application of multi-dimensional instability analysis to the wake flow reveals that the most amplified unstable mode can achieve a peak N-factor of up to 20 by X/L = 0.80, indicating that the onset of transition is more than likely to occur within eighty percent of the model length. To the best of our knowledge, the present study represents the first analysis including nonparallel and curvature effects on hypersonic tripwake instabilities in the presence of boundary-layer crossflow over a three-dimensional configuration.

Boundary layer transition↗
NASA NTRS1980

On the dynamics of large orbiting flexible beams and platforms oriented along the local horizontal

The dynamics and stability of large orbiting flexible beams and platforms oriented along the local horizontal are treated both analytically and numerically. It is assumed that such structures could be gravitationally stabilized by attaching a rigid lightweight dumbbell at the center of mass by a spring loaded hinge which also could provide viscous damping. For the beam it is seen that the small amplitude inplane pitch motion, dumbbell librational motion, and the anti-symmetric elastic modes are all coupled. The three dimensional equations of motion for a circular flat plate in orbit with a two-degree-of-freedom gimballed dumbbell are also developed and show that only those elastic modes described by a single nodal diameter line are influenced by the dumbbell motion. Stability criteria are developed for both examples and a parametric study of the least damped mode characteristics together with numerically simulated transient responses are carried out.

Bainum, P. M.↗
NASA NTRS1978

Investigation of ODE integrators using interactive graphics

Two FORTRAN programs using an interactive graphic terminal to generate accuracy and stability plots for given multistep ordinary differential equation (ODE) integrators are described. The first treats the fixed stepsize linear case with complex variable solutions, and generates plots to show accuracy and error response to step driving function of a numerical solution, as well as the linear stability region. The second generates an analog to the stability region for classes of non-linear ODE's as well as accuracy plots. Both systems can compute method coefficients from a simple specification of the method. Example plots are given.

Brown, R. L.↗
NASA NTRS1992

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems. Task 4: Advanced fan section aerodynamic analysis

The purpose of this study is the development of a three-dimensional Euler/Navier-Stokes flow analysis for fan section/engine geometries containing multiple blade rows and multiple spanwise flow splitters. An existing procedure developed by Dr. J. J. Adamczyk and associates and the NASA Lewis Research Center was modified to accept multiple spanwise splitter geometries and simulate engine core conditions. The procedure was also modified to allow coarse parallelization of the solution algorithm. This document is a final report outlining the development and techniques used in the procedure. The numerical solution is based upon a finite volume technique with a four stage Runge-Kutta time marching procedure. Numerical dissipation is used to gain solution stability but is reduced in viscous dominated flow regions. Local time stepping and implicit residual smoothing are used to increase the rate of convergence. Multiple blade row solutions are based upon the average-passage system of equations. The numerical solutions are performed on an H-type grid system, with meshes being generated by the system (TIGG3D) developed earlier under this contract. The grid generation scheme meets the average-passage requirement of maintaining a common axisymmetric mesh for each blade row grid. The analysis was run on several geometry configurations ranging from one to five blade rows and from one to four radial flow splitters. Pure internal flow solutions were obtained as well as solutions with flow about the cowl/nacelle and various engine core flow conditions. The efficiency of the solution procedure was shown to be the same as the original analysis.

Crook, Andrew J.↗
NASA NTRS1978

Investigation of ODE integrators using interactive graphics

In the present paper, two FORTRAN programs are described which use an interactive graphic terminal to generate accuracy and stability plots for a given multistep integrator of ordinary differential equations (ODE). One program treats the fixed-stepsize linear case with complex variable solutions, and generates plots to show accuracy and error response to step driving function of a numerical solution. It also generates the linear stability region. The other program generates an analog with an accuracy plot. Both systems are capable of computing method coefficients from a simple specification of the method. Example plots are given.

Brown, R. L.↗
DOE OSTI2019

A MUSCL-SCNI approach for meshfree modeling of shock waves in fluids

Here, a stable and nodally integrated meshfree formulation for modeling shock waves in fluids is developed. The reproducing kernel approximation is employed to discretize the conservation equations for compressible flow, and a flux vector splitting approach is applied to allow proper numerical treatments for the advection and pressure parts, respectively, based on the characteristics of each flux term. To capture the essential shock physics in fluids, including the Rankine–Hugoniot jump conditions and the entropy condition, local Riemann enrichment is introduced under the stabilized conforming nodal integration (SCNI) framework. Meanwhile, numerical instabilities associated with the advection flux are eliminated by adopting a modified upwind scheme. To further enhance accuracy, a MUSCL-type method is introduced in conjunction with an oscillation limiter to avoid Gibbs phenomenon and ensure monotonic piecewise linear reconstruction in the smooth region. The present meshfree formulation is free from tunable artificial parameters and is capable of capturing shock and rarefaction waves without over/undershoots. Finally, several numerical examples are analyzed to demonstrate the effectiveness of the proposed MUSCL-SCNI approach in meshfree modeling of complex shock phenomena, including shock diffraction, shock–vortex interaction, and high energy explosion processes.

42 ENGINEERING↗
DOE OSTI2024

Pressure-stabilized fixed-stress iterative solutions of compositional poromechanics

We consider the numerical behavior of the fixed-stress splitting method for coupled poromechanics as undrained regimes are approached. We explain that pressure stability is related to the splitting error of the scheme, not the fact that the discrete saddle point matrix never appears in the fixed-stress approach. This observation reconciles previous results regarding the pressure stability of the splitting method. Using examples of compositional poromechanics with application to geological CO sequestration, we see that solutions obtained using the fixed-stress scheme with a low order finite element-finite volume discretization which is not inherently inf-sup stable can exhibit the same pressure oscillations obtained with the corresponding fully implicit scheme. Moreover, pressure jump stabilization can effectively remove these spurious oscillations in the fixed-stress setting, while also improving the efficiency of the scheme in terms of the number of iterations required at every time step to reach convergence.

42 ENGINEERING↗
DOE OSTI2015

A numerical study of the diffusive-thermal instability of opposed nonpremixed tubular flames

The diffusive-thermal (D-T) instability of opposed nonpremixed tubular flames near extinction is investigated here using two-dimensional (2-D) direct numerical simulations together with the linear stability analysis. Two different initial conditions (IC), i.e. the perturbed IC and the C-shaped IC are adopted to elucidate the effects of small and large amplitude disturbances on the formation of flame cells, similar to conditions found in linear stability analysis and experiments, respectively. The characteristics of the D-T instability of tubular flames are identified by a critical Damköhler number, Da C , at which the D-T instability first occurs and the corresponding number of flame cells for three different tubular flames with different flame radii. It is found that Da C predicted through linear stability analysis shows good agreement with that obtained from the 2-D simulations performed with two different ICs. The flame cell number, N cell , from the 2-D simulations with the perturbed IC is also found to be equal to an integer close to the maximum wavenumber, k max , obtained from the linear stability analysis. However, N cell from the 2-D simulations with the C-shaped IC is smaller than k max and N cell found from the simulations with the perturbed IC. This is primarily because the strong reaction at the edges of the horseshoe-shaped cellular flame developed from the C-shaped IC is more likely to produce larger flame cells and reduce N cell . It is also found that for cases with the C-shaped IC, once the cellular instability occurs, the number of flame cells remains constant until global extinction occurs by incomplete reaction manifested by small Da. It is also verified through the displacement speed, S d , analysis that the two edges of the horseshoe-shaped cellular flame are stationary and therefore do not merge due to the diffusion–reaction balance at the edges. Moreover, large negative S d is observed at the local extinction points while small positive or negative S d features in the movement of flame cells as they adjust their location and size towards steady state.

42 ENGINEERING↗