Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability and convergence”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The arbitrary‐order virtual element method for linear elastodynamics models: convergence, stability and dispersion‐dissipation analysis

Abstract We design the conforming virtual element method for the numerical approximation of the two‐dimensional elastodynamics problem. We prove stability and convergence of the semidiscrete approximation and derive optimal error estimates under h ‐ and p ‐refinement in both the energy and the L 2 norms. The performance of the proposed virtual element method is assessed on a set of different computational meshes, including nonconvex cells up to order four in the h ‐refinement setting. Exponential convergence is also experimentally observed under p ‐refinement. Finally, we present a dispersion‐dissipation analysis for both the semidiscrete and fully discrete schemes, showing that polygonal meshes behave as classical simplicial/quadrilateral grids in terms of dispersion‐dissipation properties.

Antonietti, Paola F.↗
DOE OSTI2020

The arbitrary-order virtual element method for linear elastodynamics models. Convergence, stability and dispersion-dissipation analysis.

We design the conforming virtual element method for the numerical approximation of the two dimensional elastodynamics problem. We prove stability and convergence of the semi-discrete approximation and derive optimal error estimates under $\textit{h}$-refinement in both the energy and the $L^2$ norms, and optimal error estimates under $\textit{p}$-refinement in the energy norm. The performance of the proposed virtual element method is assessed on a set of different computational meshes, including non-convex cells up to order four in the h-refinement setting. Exponential convergence is also experimentally observed under p-refinement. Finally, we present a dispersion-dissipation analysis for both the semi-discrete and fully-discrete schemes, showing that polygonal meshes behave as classical simplicial/quadrilateral grids in terms of dispersion-dissipation properties.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Adaptive control: Stability, convergence, and robustness

The deterministic theory of adaptive control (AC) is presented in an introduction for graduate students and practicing engineers. Chapters are devoted to basic AC approaches, notation and fundamental theorems, the identification problem, model-reference AC, parameter convergence using averaging techniques, and AC robustness. Consideration is given to the use of prior information, the global stability of indirect AC schemes, multivariable AC, linearizing AC for a class of nonlinear systems, AC of linearizable minimum-phase systems, and MIMO systems decouplable by static state feedback.

Sastry, Shankar↗
DOE OSTI2024

Stability and Convergence of Solutions to Stochastic Inverse Problems Using Approximate Probability Densities

Data-consistent inversion is designed to solve a class of stochastic inverse problems where the solution is a pullback of a probability measure specified on the outputs of a quantities of interest (QoI) map. Here, this work presents stability and convergence results for the case where finite QoI data result in an approximation of the solution as a density. Given their popularity in the literature, separate results are proven for three different approaches to measuring discrepancies between probability measures: f-divergences, integral probability metrics, and L p metrics. In the context of integral probability metrics, we also introduce a pullback probability metric that is well-suited for data-consistent inversion. This fills a theoretical gap in the convergence and stability results for data-consistent inversion that have mostly focused on convergence of solutions associated with approximate maps. Numerical results are included to illustrate key theoretical results with intuitive and reproducible test problems that include a demonstration of convergence in the measure-theoretic "almost" sense.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1984

Stability and Convergence of Underintegrated Finite Element Approximations

The effects of underintegration on the numerical stability and convergence characteristics of certain classes of finite element approximations were analyzed. Particular attention is given to hourglassing instabilities that arise from underintegrating the stiffness matrix entries and checkerboard instabilities that arise from underintegrating constrain terms such as those arising from incompressibility conditions. A fundamental result reported here is the proof that the fully integrated stiffness is restored in some cases through a post-processing operation.

Oden, J. T.↗
NASA NTRS1983

Convergence and stability properties of minimal polynomial and reduced rank extrapolation algorithms

The minimal polynomial and reduced rank extrapolation algorithms are two acceleration of convergence methods for sequences of vectors. In a recent survey these methods were tested and compared with the scalar, vector, topological epsilon algorithms, and were observed to be more efficient than the latter. It was also observed that the two methods have similar convergence properties. The convergence and stability properties of these methods are analyzed and the performance of the acceleration methods when applied to a class of vector sequences that includes those sequences obtained from systems of linear equations by using matrix iterative methods is discussed.

Sidi, A.↗
DOE OSTI2025

Stability and convergence of nuclear detonations in white dwarf collisions

We investigate the numerical stability of thermonuclear detonations in 1D accelerated reactive shocks and 2D binary collisions of equal-mass magnetized and unmagnetized white dwarf stars. To achieve high resolution at initiation sites, we devised geometric gridding and mesh velocity strategies specially adapted to the unique requirements of head-on collisional geometries, scenarios in which one expects maximum production of iron-group products. We study the effects of grid resolution and the limiting of temperature, energy, and reactants for different stellar masses, separations, magnetic fields, inigenerationtial compositions, detonation mechanisms, and limiter parameters across a range of cell sizes from 1 to 100 km. Our results set bounds on the parameter space of limiter amplitudes for which both temperature- and energy-limiting procedures yield consistent and monotonically convergent solutions. Within these bounds, we find that grid resolutions of 5 km or better are necessary for uncertainties in total released energy and iron-group products to drop below 10%. Intermediate-mass products (e.g., calcium) exhibit similar convergence trends but with somewhat greater uncertainty. These conclusions apply equally to pure C/O white dwarfs, multispecies compositions (including helium shells), magnetized and unmagnetized cores, and either single or multiple detonation scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1982

Numerical analysis of confined turbulent flow

The considered investigation is concerned with the development of an efficient computational method for obtaining a physical understanding of an internal turbulent field. The employed approach makes use of a 'two equation' type model for the turbulence to obtain the numerical solution of a two-dimensional confined turbulent flow. The mean flow governing equations are considered along with the governing equation of the mean temperature and concentrations, and the boundary conditions. The numerical procedure for solving the turbulent flow is discussed, taking into account an approximation to the nonlinear terms, and the inner and outer coupling. Attention is given to a stability convergence analysis, the stability characteristics, and computational examples.

Lin, A.↗
NASA NTRS1986

On compactness of admissible parameter sets: Convergence and stability in inverse problems for distributed parameter systems

A series of numerical examples is reported and several algorithms compared for estimation of coefficients in differential equation models. Unconstrained, constrained and Tikhonov regularization methods are tested for their behavior with regard to both convergence (of approximation methods for the states and parameters) and stability (continuity of the estimates with respect to perturbations in the data or observed states).

Banks, H. T.↗
NASA NTRS1988

Semi-discrete approximations to nonlinear systems of conservation laws; consistency and L(infinity)-stability imply convergence

A convergence theory for semi-discrete approximations to nonlinear systems of conservation laws is developed. It is shown, by a series of scalar counter-examples, that consistency with the conservation law alone does not guarantee convergence. Instead, a notion of consistency which takes into account both the conservation law and its augmenting entropy condition is introduced. In this context it is concluded that consistency and L(infinity)-stability guarantee for a relevant class of admissible entropy functions, that their entropy production rate belongs to a compact subset of H(loc)sup -1 (x,t). One can now use compensated compactness arguments in order to turn this conclusion into a convergence proof. The current state of the art for these arguments includes the scalar and a wide class of 2 x 2 systems of conservation laws. The general framework of the vanishing viscosity method is studied as an effective way to meet the consistency and L(infinity)-stability requirements. How this method is utilized to enforce consistency and stability for scalar conservation laws is shown. In this context we prove, under the appropriate assumptions, the convergence of finite difference approximations (e.g., the high resolution TVD and UNO methods), finite element approximations (e.g., the Streamline-Diffusion methods) and spectral and pseudospectral approximations (e.g., the Spectral Viscosity methods).

Tadmor, Eitan↗
NASA NTRS2001

Interfacial Stability of Converging Plasma Jets for Magnetized Target Fusion

The merging of a spherical distribution of plasma jets to dynamically form a gaseous liner has been proposed for use in magnetized target fusion propulsion. In this paper, a study is made of the interfacial stability of the interaction of these jets. Specifically, the Orr-Sommerfeld equation is integrated to obtain the growth rate of a perturbation to the primary flow at the interface between the colliding jets. The analysis lead to an estimate on the tolerances on the relative flow velocities of the merging plasma jets to form a stable, imploding liner. The results show that during the merging of the jets to form a liner and before contact with the target plasma the growth of the perturbed flow at the jet interface is not likely to destabilize the liner. These data suggest that, as far as the stability of the interface between the merging jets is concerned, the formation of liner can withstand velocity variation up to 50% between the neighboring jets over the density and temperature ranges investigated.

Cassibry, J. T.↗
DOE OSTI2023

Accelerating and stabilizing the convergence of vibrational self-consistent field calculations via the direct inversion of the iterative subspace (vDIIS) algorithm

The vibrational self-consistent field (VSCF) method yields anharmonic states and spectra for molecular vibrations, and it serves as the starting point for more sophisticated correlated-vibration methods. Convergence of the iterative, non-linear optimization in VSCF calculations can be erratic or altogether unsuccessful, particularly for chemical systems involving low-frequency motions. In this work, a vibrational formulation of the Direct Inversion of the Iterative Subspace method of Pulay is presented and investigated. This formulation accounts for distinct attributes of the vibrational and electronic cases, including the expansion of each single-mode vibrational wavefunction in its own basis set. The resulting Direct Inversion of the Iterative Subspace method is shown to substantially accelerate VSCF convergence in all convergent cases as well as rectify many cases where Roothaan-based methods fail. Performance across systems ranging from small, rigid molecules to weakly bound molecular clusters is investigated in this analysis.

Chemistry↗
NASA NTRS2020

An Overview of the Null-Field Method. II: Convergence and Numerical Stability

In this paper we provide an analysis of the convergence and numerical stability of the null-field method with discrete sources. We show that (i) if the null-field scheme is numerically stable then we can decide whether or not convergence can be achieved; (ii) if the null-field scheme is numerically unstable then we cannot draw any conclusion about the convergence issue; and (iii) the numerical stability is closely related to the property of a tangential system of radiating discrete sources to form a Riesz basis. Our numerical analysis indicates that for prolate spheroids and localized vector spherical wave functions, the null-field scheme is numerically unstable (this system of vector functions does not form a Riesz basis), while for distributed vector spherical wave functions, the numerical instability is not so pronounced (this system of discrete sources almost possesses the property of being a Riesz basis). We also describe an analytical method for computing the surface integrals in the framework of the conventional null-field method with localized vector spherical wave functions which increases the stability of the numerical scheme.

Null-field method↗
DOE OSTI2022

Native Natural circulation in VERA

This work seeks to introduce a native natural circulation solver to VERA's thermohydraulics capabilities. This work begins with initial investigations into the viability of such a solver by determining the fundamental equations as well as necessary inputs. The method was then implemented in the form of a Python script coupled to CTF which was compared to TRACE. Once viability of the natural circulation solver was established, it was then implemented natively in MPACT for use with previously created simplified thermohydraulics solvers. This new implementation of the native natural circulation solver in MPACT is then tested for single physics problems for both a BWR and PWR single channel model and compared to TRACE where it shows an error of about 30% for BWRs and about 15% for PWRs. The solver is additionally tested for NuScale based models where it demonstrates an error to the reported flow rate of about 10%, resulting in an eigenvalue error of only 138 pcm when comparing fixed flow rate multiphysics runs. Finally, the solver is tested in multiphysics systems for a small PWR assembly, a small BWR assembly, and a full core approximate model of the NuScale design. In these multiphysics runs, it demonstrates no impact on multiphysics convergence stability, as well as negligible impact on runtimes. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗