Search results for "Numerical stability"
showing 9 items of 29 documents
Combining spectral and shock-capturing methods: A new numerical approach for 3D relativistic core collapse simulations
2005
We present a new three-dimensional general relativistic hydrodynamics code which is intended for simulations of stellar core collapse to a neutron star, as well as pulsations and instabilities of rotating relativistic stars. Contrary to the common approach followed in most existing three-dimensional numerical relativity codes which are based in Cartesian coordinates, in this code both the metric and the hydrodynamics equations are formulated and solved numerically using spherical polar coordinates. A distinctive feature of this new code is the combination of two types of accurate numerical schemes specifically designed to solve each system of equations. More precisely, the code uses spectra…
Stability of gyrotron operation in very high-order modes
2012
This study was motivated by the desire to increase the power, which can be delivered by gyrotrons in long pulse and continuous regimes. Since the admissible power level is determined by the density of ohmic losses in resonator walls, to increase the radiated power a gyrotron should operate in higher order modes. Using an existing gyrotron developed for plasma experiments in the International Thermonuclear Experimental Reactor as a base model, the stability of operation of such a gyrotron in modes with larger number of radial variations was studied. It is shown that the power level achievable in such gyrotrons in stable single mode regimes is close to 1.5 MW. The power level 1.7–1.8 MW can b…
On specific stability bounds for linear multiresolution schemes based on piecewise polynomial Lagrange interpolation
2009
Abstract The Deslauriers–Dubuc symmetric interpolation process can be considered as an interpolatory prediction scheme within Harten's framework. In this paper we express the Deslauriers–Dubuc prediction operator as a combination of either second order or first order differences. Through a detailed analysis of certain contractivity properties, we arrive to specific l ∞ -stability bounds for the multiresolution transform. A variety of tests indicate that these l ∞ bounds are closer to numerical estimates than those obtained with other approaches.
Efficient computation of root mean square deviations under rigid transformations
2013
The computation of root mean square deviations (RMSD) is an important step in many bioinformatics applications. If approached naively, each RMSD computation takes time linear in the number of atoms. In addition, a careful implementation is required to achieve numerical stability, which further increases runtimes. In practice, the structural variations under consideration are often induced by rigid transformations of the protein, or are at least dominated by a rigid component. In this work, we show how RMSD values resulting from rigid transformations can be computed in constant time from the protein's covariance matrix, which can be precomputed in linear time. As a typical application scenar…
On a class of compactly epi-Lipschitzian sets
2003
The paper is devoted to the study of the so-called compactly epi-Lipschitzian sets. These sets are needed for many aspects of generalized differentiation, particulary for necessary optimality conditions, stability of mathematical programming problems and calculus rules for subdifferentials and normal cones. We present general conditions under which sets defined by general constraints are compactly epi-Lipschitzian. This allows us to show how the compact epi-Lipschitzness properties behave under set intersections.
Applications and numerical convergence of the partial inverse method
2006
In 1983, J.E. Spingarn introduced what he called the Partial Inverse Method in the framework of Mathematical Programming. Since his initial articles, numerous applications have been given in various fields including Lagrangian multipliers methods, location theory, convex feasibility problems, analysis of data, economic equilibrium problems. In a first part of this paper we give a survey of these applications. Then by means of optimization problems relevant to location theory such as single and multifacility minimisum or minimax location problems, we examine the main advantages of the algorithm and we point out its drawbacks mainly concerning the rate of convergence. We study how different p…
A Positive Definite Advection Scheme Obtained by Nonlinear Renormalization of the Advective Fluxes
1989
Abstract A new method is developed to obtain a conservative and positive definite advection scheme that produces only small numerical diffusion. Advective fluxes are computed utilizing the integrated flux form of Tremback et al. These fluxes are normalized and then limited by upper and lower values. The resulting advection equation is numerically solved by means of the usual upstream procedure. The proposed treatment is not restricted to the integrated flux form but may also be applied to other known advection algorithms which are formulated in terms of advective fluxes. Different numerical tests are presented illustrating that the proposed scheme strongly reduces numerical and diffusion an…
Some efficient algorithms for the solution of a single nonlinear equation
1981
High order methods for the numerical solution of nonlinear scalar equations are proposed which are more efficient than known procedures, and a unified approach to various methods suggested in literature is given.
Approximate analytic and numerical solutions to Lane-Emden equation via fuzzy modeling method
2012
Published version in the journal: Mathematical Problems in Engineering. Also available from the publisher: http://dx.doi.org/10.1155/2012/259494 A novel algorithm, called variable weight fuzzy marginal linearization VWFML method, is proposed. Thismethod can supply approximate analytic and numerical solutions to Lane-Emden equations. And it is easy to be implemented and extended for solving other nonlinear differential equations. Numerical examples are included to demonstrate the validity and applicability of the developed technique.