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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

The structure and stability of coronal magnetic fields

The high degree of symmetry often assumed in studies of the structure and stability of coronal magnetic field configurations is restrictive and can yield misleading results. Fully three-dimensional numerical methods for constructing force-free equilibria and for examining their stability properties are developed which make no assumptions about symmetry. A test of the stability analysis has been performed by applying it to the Gold-Hoyle twisted flux tube, which is known to be kink-unstable if the helical field makes more than about one turn between the line-tying end-plates. The preliminary result is that the critical number of turns is about 1.1, in good agreement with the previous best estimate. However, it is found that the growth rate, which has not been discussed previously, is orders of magnitude smaller than expected, even when the flux tube is twisted far beyond the stability limit.

Mcclymont, A. N.↗
NASA NTRS2005

Numerical Analysis of Flow Evolution in a Helium Jet Injected into Ambient Air

A computational model to study the stability characteristics of an evolving buoyant helium gas jet in ambient air environment is presented. Numerical formulation incorporates a segregated approach to solve for the transport equations of helium mass fraction coupled with the conservation equations of mixture mass and momentum using a staggered grid method. The operating parameters correspond to the Reynolds number varying from 30 to 300 to demarcate the flow dynamics in oscillating and non-oscillating regimes. Computed velocity and concentration fields were used to analyze the flow structure in the evolving jet. For Re=300 case, results showed that an instability mode that sets in during the evolution process in Earth gravity is absent in zero gravity, signifying the importance of buoyancy. Though buoyancy initiates the instability, below a certain jet exit velocity, diffusion dominates the entrainment process to make the jet non-oscillatory as observed for the Re=30 case. Initiation of the instability was found to be dependent on the interaction of buoyancy and momentum forces along the jet shear layer.

Satti, Rajani P.↗
NASA NTRS1998

Developments in Cylindrical Shell Stability Analysis

Today high-performance computing systems and new analytical and numerical techniques enable engineers to explore the use of advanced materials for shell design. This paper reviews some of the historical developments of shell buckling analysis and design. The paper concludes by identifying key research directions for reliable and robust methods development in shell stability analysis and design.

Knight, Norman F., Jr.↗
NASA NTRS1975

A numerical method for solving the Navier-Stokes equations with application to shock-boundary layer interactions

A numerical method for solving the compressible form of the unsteady Navier-Stokes equations is described. This method was originally presented in 1970 and has since been modified during the development of computer programs at Ames for implementing models that account for the effects of turbulence in shock-induced separated flows. Although this paper does not describe the turbulence models themselves, a complete description of the basic numerical method is given with emphasis on the choice of a computational mesh for high Reynolds number flows, finite-difference approximations for mixed partial derivatives, extension of the Courant-Friedrichs-Lewy stability condition for viscous flows, mesh boundary conditions, and numerical smoothing for strong shock-wave calculations.

Maccormack, R. W.↗
DOE OSTI2018

Analysis of the Spectral Stability of the Generalized Runge–Kutta Methods Applied to Initial-Boundary-Value Problems for Equations of the Parabolic Type. I. Explicit Methods

We develop a general algorithm for the analysis the spectral stability of generalized multistage Runge–Kutta methods of various orders of accuracy as applied to the numerical integration with respect to time of an initial-boundary-value problem for the second-order parabolic equation. The expression for the function of spectral stability is obtained in two alternative forms: on the basis of matrix relations and in the determinant form. We study specific realizations of various generalized explicit Runge–Kutta methods and their spectral stability. It is shown that all explicit generalized Runge–Kutta methods possess the property of conditional spectral stability and the property of conditional monotonicity of the numerical solution in time whose violation leads to the appearance of false oscillations of the approximate solution. The stability function for these methods is polynomial. It is shown that the application of two-stage generalized explicit Runge–Kutta methods leads to the appearance of predictor-corrector-type schemes. In the case of nonstationary one-dimensional heat-conduction problem, on the basis of a one-stage generalized Runge–Kutta method, we obtain the classical two-layer conditionally stable explicit finite-difference scheme on a four-point template. It is demonstrated that the five-stage generalized Runge–Kutta–Merson method is characterized by the weakest condition of spectral stability as compared with all investigated explicit generalized Runge–Kutta methods.

42 ENGINEERING↗
NASA NTRS1993

On the shear stabilization of capillary break-up of finite liquid bridges

In this paper we consider an isothermal finite liquid bridge under zero-gravity. A sinusoidal interfacial shear stress drives a nonparallel flow in the bridge. The linear stability of this flow to three-dimensional disturbances (which may deform the cylindrical gas-liquid interface) is determined numerically by solving an elliptic eigenvalue problem. Previous results on shear stabilization of capillary break-up of axially unbounded cylindrical interfaces containing a parallel flow are hereby extended to a nonparallel flow. The turning flow regions influence the stability significantly. However, for the particular cases considered, a small area in parameter space remains where the capillary instability is suppressed through interfacial shear. Second, non-axisymmetric oscillatory instabilities are found which originate from an interaction of the interface deformation and the basic flow. These instabilities may be the isothermal limit of the oscillatory instabilities observed in float-zone crystal growth.

Dijkstra, Henk A.↗
DOE OSTI2026

Neural entropy-stable conservative flux form neural networks for learning hyperbolic conservation laws

We propose a neural entropy-stable conservative flux form neural network (NESCFN) for learning hyperbolic conservation laws and their associated entropy functions directly from solution trajectories, without requiring any predefined numerical discretization. While recent neural network architectures have successfully integrated classical numerical principles into learned models, most rely on prior knowledge of the governing equations or assume a fixed discretization. Our approach removes this dependency by embedding entropy-stable design principles into the learning process itself, enabling the discovery of physically consistent dynamics in a fully data-driven setting. By jointly learning both the flux function and a corresponding entropy, NESCFN promotes conservation and entropy dissipation, which is critical for long-term stability and fidelity in the system of hyperbolic conservation laws. Furthermore, numerical results demonstrate that the method achieves stability and conservation over extended time horizons and accurately captures shock propagation speeds, even without oracle access to future-time solution profiles in the training data.

Conservative flux form↗
DOE OSTI2017

High- β equilibrium and ballooning stability of the low aspect ratio CNT stellarator

In the paper, the existence and ballooning-stability of low aspect ratio stellarator equilibria is predicted for the Columbia Neutral Torus (CNT) with the aid of 3D numerical tools. In addition to having a low aspect ratio, CNT is characterized by a low magnetic field and small plasma volume. Also, highly overdense plasmas were recently heated in CNT by means of microwaves. These characteristics suggest that CNT might attain relatively high values of plasma beta and thus be of use in the experimental study of stellarator stability to high-beta instabilities such as ballooning modes. As a first step in that direction, here the ballooning stability limit is found numerically. Depending on the particular magnetic configuration we expect volume-averaged β limits in the range 0.9%–3.0%, and possibly higher, and observe indications of a second region of ballooning stability. As the aspect ratio is reduced, stability is found to increase in some configurations and decrease in others. Energy-balance estimates using stellarator scaling laws indicate that the lower β limit may be attainable with overdense heating at powers of 40 to 100 kW. The present study serves the additional purpose of testing VMEC and other stellarator codes at high values of β and at low aspect ratios. For this reason, the study was carried out both for free boundary, for maximum fidelity to experiment, as well as with a fixed boundary, as a numerical test.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The arbitrary order mixed mimetic finite difference method for the diffusion equation

Here, we propose an arbitrary-order accurate mimetic finite difference (MFD) method for the approximation of diffusion problems in mixed form on unstructured polygonal and polyhedral meshes. As usual in the mimetic numerical technology, the method satisfies local consistency and stability conditions, which determines the accuracy and the well-posedness of the resulting approximation. The method also requires the definition of a high-order discrete divergence operator that is the discrete analog of the divergence operator and is acting on the degrees of freedom. The new family of mimetic methods is proved theoretically to be convergent and optimal error estimates for flux and scalar variable are derived from the convergence analysis. A numerical experiment confirms the high-order accuracy of the method in solving diffusion problems with variable diffusion tensor. It is worth mentioning that the approximation of the scalar variable presents a superconvergence effect.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A Homotopy-Based Method for Robust Computation of Controlling Unstable Equilibrium Points

This paper proposes a homotopy-based method for efficiently computing and precisely determining controlling unstable equilibrium points (CUEPs) in transient stability analysis and screening using direct methods. The proposed method is intended to overcome potential numerical convergence problems associated with computing CUEPs. These convergence problems are mostly due to irregularities in the region of convergence of CUEPs which make it difficult to determine suitable initial points. The most commonly used initial points in computing CUEPs are minimum gradient points (MGPs) on the stability boundary which in turn rely on the determination of exit points (EPs-points at which sustained-fault trajectories exit system stability boundaries). However, determination of MGPs and EPs are computationally involved and sophisticated numerical methods have been introduced such as stability-boundary-following procedure and shadowing method to correct the solution trajectory. The proposed method does not depend on MGPs and also does not require accurate EPs. The proposed method maps the solution from EPs to the CUEPs. The proposed method is applied on several known systems including the NE 39 bus system and the reduced WECC system and the results are provided. The results are compared with the traditional methods such as Boundary of stability region based CUEP method (BCU method).

42 ENGINEERING↗
DOE OSTI2017

An asymptotic-preserving stochastic Galerkin method for the radiative heat transfer equations with random inputs and diffusive scalings

In this paper, we develop an Asymptotic-Preserving (AP) stochastic Galerkin scheme for the radiative heat transfer equations with random inputs and diffusive scalings. In this problem the random inputs arise due to uncertainties in cross section, initial data or boundary data. We use the generalized polynomial chaos based stochastic Galerkin (gPC-SG) method, which is combined with the micro–macro decomposition based deterministic AP framework in order to handle efficiently the diffusive regime. For linearized problem we prove the regularity of the solution in the random space and consequently the spectral accuracy of the gPC-SG method. We also prove the uniform (in the mean free path) linear stability for the space-time discretizations. Several numerical tests are presented to show the efficiency and accuracy of proposed scheme, especially in the diffusive regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1984

On the stability of the decelerating laminar boundary layer

The stability of a decelerating boundary-layer flow is investigated experimentally and numerically. Experimentally, a flat plate having a Blasius boundary layer is decelerated in an 18 m towing tank. The boundary layer becomes unstable to two-dimensional waves, which break down into three-dimensional patterns, hairpin vortices, and finally turbulent bursts when the vortices lift off the wall. The unsteady boundary-layer equations are solved numerically to generate instantaneous velocity profiles for a range of boundary and initial conditions. A quasi-steady approximation is invoked and the stability of local velocity profiles is determined by solving the Orr-Sommerfeld equation using Chebyshev matrix methods. Comparisons are made between the numerical predictions and the experimentally observed instabilities.

Gad-El-hak, M.↗
NASA NTRS1992

Bleed-boundary conditions for numerically simulated mixed-compression supersonic inlet flow

A numerical study was performed to investigate boundary conditions for bleed openings that are needed in the numerical simulation of critical flow through turbojet inlets with bleed. Nine different boundary conditions (BCs) were devised and examined. Of the BCs examined, three were able to stabilize the terminal shock near the inlet's throat. This numerical study was based on the ensemble-averaged, unsteady, three-dimensional, compressible Navier-Stokes equations closed by the Baldwin-Lomax algebraic turbulence model. Solutions to the Navier-Stokes equations were obtained by using a hybrid implicit-explicit method (Beam and Warming/MacCormack). The multiblock grid system used in the solution procedure was generated by solving a set of elliptic partial differential equations. Solutions obtained were compared with available experimental data.

Chyu, W. J.↗
DOE OSTI2017

Development of Nanoparticle-Stabilized Foams to Improve Performance of Water-less Hydraulic Fracturing

We have successfully created ultra dry carbon-dioxide-in-water and nitrogen-in-water foams (with water content down to 2-5% range), that are remarkably stable at high temperatures (up to 120 deg, C) and pressures (up to 3000psi) and viscous enough (100-200 cP tunable range) to carry proppant. Two generations of these ultra-dry foams have been developed; they are stabilized either with a synergy of surfactants and nanoparticle, or just with viscoelastic surfactants that viscosify the aqueous phase. Not only does this reduce water utilization and disposal, but it minimizes fluid blocking of hydrocarbon production. Further, the most recent development shows successful use of environmentally friendly surfactants at high temperature and pressure. We pay special attention to the role of nanoparticles in stabilization of the foams, specifically for high salinity brines. The preliminary numerical simulation for which shows they open wider fractures with shorter half-length and require less clean-up due to minimal water use. We also tested the stability and sand carrying properties of these foams at high pressure, room temperature conditions in sapphire cell. We performed on a preliminary numerical investigation of applicability for improved oil recovery applications. The applicability was evaluated by running multiphase flow injection simulations in a case-study oil reservoir. The results of this research thus expand the options available to operators for hydraulic fracturing and can simplify the design and field implementation of foamed fracturing fluids.

36 MATERIALS SCIENCE↗
NASA NTRS2010

Guidelines for Computing Longitudinal Dynamic Stability Characteristics of a Subsonic Transport

A systematic study is presented to guide the selection of a numerical solution strategy for URANS computation of a subsonic transport configuration undergoing simulated forced oscillation about its pitch axis. Forced oscillation is central to the prevalent wind tunnel methodology for quantifying aircraft dynamic stability derivatives from force and moment coefficients, which is the ultimate goal for the computational simulations. Extensive computations are performed that lead in key insights of the critical numerical parameters affecting solution convergence. A preliminary linear harmonic analysis is included to demonstrate the potential of extracting dynamic stability derivatives from computational solutions.

Thompson, Joseph R.↗
DOE OSTI2014

Stability of Hall equilibria in neutron star crusts

In the solid crusts of neutron stars, the advection of the magnetic field by the current-carrying electrons, an effect known as Hall drift, should play a very important role as the ions remain essentially fixed (as long as the solid does not break). Although Hall drift preserves the magnetic field energy, it has been argued that it may drive a turbulent cascade to scales at which ohmic dissipation becomes effective, allowing a much faster decay in objects with very strong fields. On the other hand, it has been found that there are 'Hall equilibria', i.e., field configurations that are unaffected by Hall drift. Here we address the crucial question of the stability of these equilibria through axially symmetric (two-dimensional (2D)) numerical simulations of Hall drift and ohmic diffusion, with the simplifying assumption of uniform electron density and conductivity. We demonstrate the 2D stability of a purely poloidal equilibrium, for which ohmic dissipation makes the field evolve toward an attractor state through adjacent stable configurations, around which damped oscillations occur. For this field, the decay scales with the ohmic timescale. We also study the case of an unstable equilibrium consisting of both poloidal and toroidal field components that are confined within the crust. This field evolves into a stable configuration, which undergoes damped oscillations superimposed on a slow evolution toward an attractor, just as the purely poloidal one.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Numerical experiments with the extended virtual element method for the Laplace problem with strong discontinuities

The virtual element method (VEM) is a stabilized Galerkin formulation on arbitrary polytopal meshes. In the VEM, the basis functions are implicit (virtual) — they are not known nor do they need to be computed within the problem domain. Suitable projection operators are used to decompose the bilinear form at the element-level into two parts: a consistent term that reproduces a given polynomial space and a correction term that ensures stability. In this study, we consider a low-order extended virtual element method (X-VEM) that is in the spirit of the extended finite element method for crack problems. Herein, we focus on the two-dimensional Laplace crack problem. In the X-VEM, we enrich the standard virtual element space with additional discontinuous functions through the framework of partition-of-unity. The nodal basis functions in the VEM are chosen as the partition-of-unity functions, and we study means to stabilize the standard and enriched sub-matrices that constitute the element stiffness matrix. Numerical experiments are performed on the problem of a cracked membrane under mode III loading to affirm the accuracy, and to establish the optimal convergence in energy of the method.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Secondary instability in three-dimensional magnetic reconnection

We consider the transition to turbulence in three-dimensional reconnection of a magnetic neutral sheet. We find that the transition can occur via a three-step process. First, the sheet undergoes the usual tearing instability. Second, the tearing mode saturates to form a two-dimensional quasi-steady state. Third, this secondary equilibrium is itself unstable when it is perturbed by three-dimensional disturbances. Most of this paper is devoted to the analysis and simulation of the three-dimensional linear stability properties of the two-dimensional saturated tearing layer. The numerical simulations are performed with a semi-implicit, pseudospectral-Fourier collocation algorithm. We identify a three-dimensional secondary linear stability which grows on the ideal timescale. An examination of the modal energetics reveals that the largest energy transfer is from the mean field to the three-dimensional field, with the two-dimensional field acting as a catalyst.

Dahlburg, R. B.↗