Analysis of the Runge-Kutta-Fehlberg method
Stability and error analysis of Runge-Kutta- Fehlberg numerical integration method
NUMERICAL INTEGRATION↗
SEARCH · Engineering Papers
Results for “stabilization method”
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.
Stabilization of computational procedures for constrained dynamical systems
A new stabilization method of treating constraints in multibody dynamical systems is presented. By tailoring a penalty form of the constraint equations, the method achieves stabilization without artificial damping and yields a companion matrix differential equation for the constraint forces; hence, the constraint forces are obtained by integrating the companion differential equation for the constraint forces in time. A principal feature of the method is that the errors committed in each constraint condition decay with its corresponding characteristic time scale associated with its constraint force. Numerical experiments indicate that the method yields a marked improvement over existing techniques.
Park, K. C.↗
Algebraic Nonoverlapping Domain Decomposition Methods for Stabilized FEM and FV Discretizations
We consider preconditioning methods for convection dominated fluid flow problems based on a nonoverlapping Schur complement domain decomposition procedure for arbitrary triangulated domains. The triangulation is first partitioned into a number of subdomains and interfaces which induce a natural 2 x 2 partitioning of the p.d.e. discretization matrix. We view the Schur complement induced by this partitioning as an algebraically derived coarse space approximation. This avoids the known difficulties associated with the direct formation of an effective coarse discretization for advection dominated equations. By considering various approximations of the block factorization of the 2 x 2 system, we have developed a family of robust preconditioning techniques. A computer code based on these ideas has been developed and tested on the IBM SP2 using MPI message passing protocol. A number of 2-D CFD calculations will be presented for both scalar advection-diffusion equations and the Euler equations discretized using stabilized finite element and finite volume methods. These results show very good scalability of the preconditioner for various discretizations as the number of processors is increased while the number of degrees of freedom per processor is fixed.
Barth, Timothy J.↗
A method to stabilize linear systems using eigenvalue gradient information
Formal optimization methods and eigenvalue gradient information are used to develop a stabilizing control law for a closed loop linear system that is initially unstable. The method was originally formulated by using direct, constrained optimization methods with the constraints being the real parts of the eigenvalues. However, because of problems in trying to achieve stabilizing control laws, the problem was reformulated to be solved differently. The method described uses the Davidon-Fletcher-Powell minimization technique to solve an indirect, constrained minimization problem in which the performance index is the Kreisselmeier-Steinhauser function of the real parts of all the eigenvalues. The method is applied successfully to solve two different problems: the determination of a fourth-order control law stabilizes a single-input single-output active flutter suppression system and the determination of a second-order control law for a multi-input multi-output lateral-directional flight control system. Various sets of design variables and initial starting points were chosen to show the robustness of the method.
Wieseman, C. D.↗
Analysis of Instabilities in Non-Axisymmetric Hypersonic Boundary Layers Over Cones
Hypersonic flows over circular cones constitute one of the most important generic configurations for fundamental aerodynamic and aerothermodynamic studies. In this paper, numerical computations are carried out for Mach 6 flows over a 7-degree half-angle cone with two different flow incidence angles and a compression cone with a large concave curvature. Instability wave and transition-related flow physics are investigated using a series of advanced stability methods ranging from conventional linear stability theory (LST) and a higher-fidelity linear and nonlinear parabolized stability equations (PSE), to the 2D eigenvalue analysis based on partial differential equations. Computed N factor distribution pertinent to various instability mechanisms over the cone surface provides initial assessments of possible transition fronts and a guide to corresponding disturbance characteristics such as frequency and azimuthal wave numbers. It is also shown that strong secondary instability that eventually leads to transition to turbulence can be simulated very efficiently using a combination of advanced stability methods described above.
Li, Fei↗
Liapunov's second method
Stability theorem based on Liapunov function, and using invariance property of limit sets of solutions to differential equations
THEOREM PROVING↗
Analysis of random signal combinations for spacecraft pointing stability
Methods for obtaining the probability density function of random signal combustions are discussed. These methods provide a realistic criteria for the design of control systems subjected to external noise with several important applications for aerospace problems.
Howell, L.↗
Apparatus and method for stabilized phase detection for binary signal tracking loops
Apparatus and method is presented for phase detection in binary signal tracking loops wherein two bandpass detectors are alternately interchanged between electrical connection with two local code reference tracking signals in order to cancel any adverse effect of gain imbalance in the bandpass detectors and direct current offset or drift. The detectors are time shared in multiplex fashion between the two local reference signals.
Hopkins, P. M.↗
Modeling Metallic Fuel using Peridynamics
Based on available modeling and simulation capabilities of peridynamics module in MOOSE framework for oxide fuel, the overall goal of this project is to further develop the peridynamics capabilities for modeling metallic fuel. It includes two major tasks: 1) develop validated scheme to handle the shape tensor singularity due to insufficient active neighbors of a material particle in the peridynamic correspondence model for fracture problems, and 2) develop failure modeling scheme including failure criterion for metallic fuels. Before the peridynamics can be applied to model metallic fuel, the formulation instability of the peridynamic correspondence model should be addressed. The PI first worked on developing new stabilization method to improve the performance of the peridynamic correspondence model and reduce the possibility of getting a singular shape tensor while applying the model for fracture problems. The new stabilization scheme uses bond-associated weight function rather than bond-associated horizon. Compared to bondassociated horizon stabilized method, this new stabilization scheme has better performance with improved prediction accuracy and reduced free surface effect. Using this newly developed stabilization, materials models from BISON can be directly used in peridynamics for metallic fuels, such as fission rate and burnup dependent creep and swell models. Publication of this work is under preparation.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Computational Methods for Stability and Control (COMSAC): The Time Has Come
Powerful computational fluid dynamics (CFD) tools have emerged that appear to offer significant benefits as an adjunct to the experimental methods used by the stability and control community to predict aerodynamic parameters. The decreasing costs for and increasing availability of computing hours are making these applications increasingly viable as time goes on and the cost of computing continues to drop. This paper summarizes the efforts of four organizations to utilize high-end computational fluid dynamics (CFD) tools to address the challenges of the stability and control arena. General motivation and the backdrop for these efforts will be summarized as well as examples of current applications.
Hall, Robert M.↗
On the Application of the Energy Method to Stability Problems
Since stability problems have come into the field of vision of engineers, energy methods have proved to be one of the most powerful aids in mastering them. For finding the especially interesting critical loads special procedures have evolved that depart somewhat from those customary in the usual elasticity theory. A clarification of the connections seemed desirable,especially with regard to the post-critical region, for the treatment of which these special methods are not suited as they are. The present investigation discusses this question-complex (made important by shell construction in aircraft) especially in the classical example of the Euler strut, because in this case - since the basic features are not hidden by difficulties of a mathematical nature - the problem is especially clear. The present treatment differs from that appearing in the Z.f.a.M.M. (1938) under the title "Uber die Behandlung von Stabilittatsproblemen mit Hilfe der energetischen Methode" in that, in order to work out the basic ideas still more clearly,it dispenses with the investigation of behavior at large deflections and of the elastic foundation;in its place the present version gives an elaboration of the 6th section and (in its 7 th and 8th secs.)a new example that shows the applicability of the general criterion to a stability problem that differs from that of Euler in many respects.
Marguerre, Karl↗
Development of a transfer function method for dynamic stability measurement
Flutter testing method based on transfer function measurements is developed. The error statistics of several dynamic stability measurement methods are reviewed. It is shown that the transfer function measurement controls the error level by averaging the data and correlating the input and output. The method also gives a direct estimate of the error in the response measurement. An algorithm is developed for obtaining the natural frequency and damping ratio of low damped modes of the system, using integrals of the transfer function in the vicinity of a resonant peak. Guidelines are given for selecting the parameters in the transfer function measurement. Finally, the dynamic stability measurement technique is applied to data from a wind tunnel test of a proprotor and wing model.
Johnson, W.↗
Time Domain Stability Margin Assessment Method
The baseline stability margins for NASA's Space Launch System (SLS) launch vehicle were generated via the classical approach of linearizing the system equations of motion and determining the gain and phase margins from the resulting frequency domain model. To improve the fidelity of the classical methods, the linear frequency domain approach can be extended by replacing static, memoryless nonlinearities with describing functions. This technique, however, does not address the time varying nature of the dynamics of a launch vehicle in flight. An alternative technique for the evaluation of the stability of the nonlinear launch vehicle dynamics along its trajectory is to incrementally adjust the gain and/or time delay in the time domain simulation until the system exhibits unstable behavior. This technique has the added benefit of providing a direct comparison between the time domain and frequency domain tools in support of simulation validation.
Clements, Keith↗
Methods of obtaining stability derivatives
Various methods of obtaining aerodynamic derivatives are reviewed. These methods include wind-tunnel tests, aerodynamic theory, and extraction from flight tests, with the emphasis being on flight tests. Two techniques for extracting derivatives from flight are covered in sufficient detail to indicate the procedures involved. These methods are a strict least-squares approach and an iteration technique. References are listed for additional background and details of various techniques.
M. J. Queijo↗
Stability of multistep methods in numerical integration
Stability of multistep methods in numerical integration
STABILITY↗
Multigrid method for stability problems
The problem of calculating the stability of steady state solutions of differential equations is addressed. Leading eigenvalues of large matrices that arise from discretization are calculated, and an efficient multigrid method for solving these problems is presented. The resulting grid functions are used as initial approximations for appropriate eigenvalue problems. The method employs local relaxation on all levels together with a global change on the coarsest level only, which is designed to separate the different eigenfunctions as well as to update their corresponding eigenvalues. Coarsening is done using the FAS formulation in a nonstandard way in which the right-hand side of the coarse grid equations involves unknown parameters to be solved on the coarse grid. This leads to a new multigrid method for calculating the eigenvalues of symmetric problems. Numerical experiments with a model problem are presented which demonstrate the effectiveness of the method.
Ta'asan, Shlomo↗
A unified framework of stabilized finite element method for solving the Boltzmann transport equation
This paper presents a unified framework of stabilized finite element method for solving the Boltzmann transport equation. Unlike the traditional Petrov-Galerkin finite element method which modifies the test function to construct the stabilization term, we derive the stabilization methods from the standard Galerkin weak form with Sub-grid scale model. The basic idea of this method is to decompose the unknowns into its numerical solution and residual, with an approximation for the residual and embeds it in the Galerkin weak form to yield a stabilized variational formula. Different approximations of the residual lead to different stabilization methods, all the frequently used stabilized methods, including the Streamline Upwinding Petrov-Galerkin (SUPG) method, Galerkin/Least-Square (GLS) method, and Algebraic Sub-Grid Scale (ASGS) method can be obtained from this framework. The similarities and differences of the different approximations are compared in this paper. The numerical results show that the behaviors of the different methods area similar with the same stabilization parameter, and all these stabilized techniques can obtain a correct and stable solution. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Multigrid method for stability problems
The problem of calculating the stability of steady state solutions of differential equations is treated. Leading eigenvalues (i.e., having maximal real part) of large matrices that arise from discretization are to be calculated. An efficient multigrid method for solving these problems is presented. The method begins by obtaining an initial approximation for the dominant subspace on a coarse level using a damped Jacobi relaxation. This proceeds until enough accuracy for the dominant subspace has been obtained. The resulting grid functions are then used as an initial approximation for appropriate eigenvalue problems. These problems are being solved first on coarse levels, followed by refinement until a desired accuracy for the eigenvalues has been achieved. The method employs local relaxation on all levels together with a global change on the coarsest level only, which is designed to separate the different eigenfunctions as well as to update their corresponding eigenvalues. Coarsening is done using the FAS formulation in a non-standard way in which the right hand side of the coarse grid equations involves unknown parameters to be solved for on the coarse grid. This in particular leads to a new multigrid method for calculating the eigenvalues of symmetric problems. Numerical experiments with a model problem demonstrate the effectiveness of the method proposed. Using an FMG algorithm a solution to the level of discretization errors is obtained in just a few work units (less than 10), where a work unit is the work involved in one Jacobi relization on the finest level.
Taasan, Shlomo↗