SEARCH · Search NASA

Results for “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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Methods for Measuring and Computing the Reference Temperature in Newton’s Law of Cooling for External Flows

Newton’s law of cooling requires a reference temperature (𝑇 𝑟𝑒𝑓 ) to define the heat-transfer coefficient (ℎ). For external flows with multiple temperatures in the freestream, obtaining 𝑇 𝑟𝑒𝑓 is a challenge. One widely used method, referred to as the adiabatic-wall (AW) method, obtains 𝑇 𝑟𝑒𝑓 by requiring the surface of the solid exposed to convective heat transfer to be adiabatic. Another widely used method, referred to as the linear-extrapolation (LE) method, obtains 𝑇 𝑟𝑒𝑓 by measuring/computing the heat flux ($𝑞^{′′}_𝑠$) on the solid surface at two different surface temperatures (𝑇 𝑠 ) and then linearly extrapolating to $𝑞^{′′}_𝑠$ = 0. A third recently developed method, referred to as the state-space (SS) method, obtains 𝑇 𝑟𝑒𝑓 by probing the temperature space between the highest and lowest in the flow to account for the effects of 𝑇 𝑠 or $𝑞^{′′}_𝑠$ on 𝑇 𝑟𝑒𝑓 . This study examines the foundation and accuracy of these methods via a test problem involving film cooling of a flat plate where $𝑞^{′′}_𝑠$ switches signs on the plate’s surface. Results obtained show that only the SS method could guarantee a unique and physically meaningful 𝑇 𝑟𝑒𝑓 where 𝑇 𝑠 =𝑇 𝑟𝑒𝑓 on a nonadiabatic surface $𝑞^{′′}_𝑠$ = 0. The AW and LE methods both assume 𝑇 𝑟𝑒𝑓 to be independent of 𝑇 𝑠 , which the SS method shows to be incorrect. Though this study also showed the adiabatic-wall temperature, 𝑇 𝐴𝑊 , to be a good approximation of 𝑇 𝑟𝑒𝑓 (<10% relative error), huge errors can occur in ℎ about the solid surface where |𝑇 𝑠 −𝑇 𝐴𝑊 | is near zero because where 𝑇 𝑠 =𝑇 𝐴𝑊 , $𝑞^{′′}_𝑠$ ≠ 0.

Newton's law of cooling↗

NASA NTRS1972

An evaluation of random analysis methods for the determination of panel damping

An analysis is made of steady-state and non-steady-state methods for the measurement of panel damping. Particular emphasis is placed on the use of random process techniques in conjunction with digital data reduction methods. The steady-state methods considered use the response power spectral density, response autocorrelation, excitation-response crosspower spectral density, or single-sided Fourier transform (SSFT) of the response autocorrelation function. Non-steady-state methods are associated mainly with the use of rapid frequency sweep excitation. Problems associated with the practical application of each method are evaluated with specific reference to the case of a panel exposed to a turbulent airflow, and two methods, the power spectral density and the single-sided Fourier transform methods, are selected as being the most suitable. These two methods are demonstrated experimentally, and it is shown that the power spectral density method is satisfactory under most conditions, provided that appropriate corrections are applied to account for filter bandwidth and background noise errors. Thus, the response power spectral density method is recommended for the measurement of the damping of panels exposed to a moving airflow.

Bhat, W. V.↗

NASA NTRS1972

Interferometric spectropolarimetry - Alternate experimental methods.

Three alternate methods of obtaining spectra of the intensity and state of polarization of light are proposed. The methods make use of a two-beam amplitude division interferometer using the technique of Fourier spectroscopy. They do not require the presence of polarizers in the arms of the instrument. In one method a single analyzer is used in front of the detector with three successive orientations of its transmission axis azimuth. In another method, a (linear) polarizer assuming the same set of orientations is placed in the incident beam. A third method, a hybrid of the former two methods, makes use of both a polarizer and analyzer in the locations indicated. The first method permits the determination of all four Stokes parameters of polarization, whereas the other two methods cannot yield the ellipticity parameter. All methods require the recording of three interferograms. However, two interferograms can provide the intensity and degree of polarization in any of the methods described.

Fymat, A. L.↗

NASA NTRS1976

The Relation of Finite Element and Finite Difference Methods

Finite element and finite difference methods are examined in order to bring out their relationship. It is shown that both methods use two types of discrete representations of continuous functions. They differ in that finite difference methods emphasize the discretization of independent variable, while finite element methods emphasize the discretization of dependent variable (referred to as functional approximations). An important point is that finite element methods use global piecewise functional approximations, while finite difference methods normally use local functional approximations. A general conclusion is that finite element methods are best designed to handle complex boundaries, while finite difference methods are superior for complex equations. It is also shown that finite volume difference methods possess many of the advantages attributed to finite element methods.

Vinokur, M.↗

NASA NTRS1986

Methods for assessing wall interference in the 2- by 2-foot adaptive-wall wind tunnel

Discussed are two methods for assessing two-dimensional wall interference in the adaptive-wall test section of the NASA Ames 2 x 2-Foot Transonic Wind Tunnel: (1) a method for predicting free-air conditions near the walls of the test section (adaptive-wall methods); and (2) a method for estimating wall-induced velocities near the model (correction methods), both of which methods are based on measurements of either one or two components of flow velocity near the walls of the test section. Each method is demonstrated using simulated wind tunnel data and is compared with other methods of the same type. The two-component adaptive-wall and correction methods were found to be preferable to the corresponding one-component methods because: (1) they are more sensitive to, and give a more complete description of, wall interference; (2) they require measurements at fewer locations; (3) they can be used to establish free-stream conditions; and (4) they are independent of a description of the model and constants of integration.

Schairer, E. T.↗

NASA NTRS1987

Accuracies of southwell and force/stiffness methods in the prediction of buckling strength of hypersonic aircraft wing tubular panels

Accuracies of the Southwell method and the force/stiffness (F/S) method are examined when the methods were used in the prediction of buckling loads of hypersonic aircraft wing tubular panels, based on nondestructive buckling test data. Various factors affecting the accuracies of the two methods were discussed. Effects of load cutoff point in the nondestructive buckling tests on the accuracies of the two methods were discussed in great detail. For the tubular panels under pure compression, the F/S method was found to give more accurate buckling load predictions than the Southwell method, which excessively overpredicts the buckling load. It was found that the Southwell method required a higher load cutoff point, as compared with the F/S method. In using the F/S method for predicting the buckling load of tubular panels under pure compression, the load cutoff point of approximately 50 percent of the critical load could give reasonably accurate predictions.

Ko, William L.↗

NASA NTRS1988

An investigation of new methods for estimating parameter sensitivities

Parameter sensitivity is defined as the estimation of changes in the modeling functions and the design variables due to small changes in the fixed parameters of the formulation. There are currently several methods for estimating parameter sensitivities requiring either difficult to obtain second order information, or do not return reliable estimates for the derivatives. Additionally, all the methods assume that the set of active constraints does not change in a neighborhood of the estimation point. If the active set does in fact change, than any extrapolations based on these derivatives may be in error. The objective here is to investigate more efficient new methods for estimating parameter sensitivities when the active set changes. The new method is based on the recursive quadratic programming (RQP) method and in conjunction a differencing formula to produce estimates of the sensitivities. This is compared to existing methods and is shown to be very competitive in terms of the number of function evaluations required. In terms of accuracy, the method is shown to be equivalent to a modified version of the Kuhn-Tucker method, where the Hessian of the Lagrangian is estimated using the BFS method employed by the RPQ algorithm. Inital testing on a test set with known sensitivities demonstrates that the method can accurately calculate the parameter sensitivity. To handle changes in the active set, a deflection algorithm is proposed for those cases where the new set of active constraints remains linearly independent. For those cases where dependencies occur, a directional derivative is proposed. A few simple examples are included for the algorithm, but extensive testing has not yet been performed.

Beltracchi, Todd J.↗

NASA NTRS1991

Development of advanced modal methods for calculating transient thermal and structural response

Higher-order modal methods for predicting thermal and structural response are evaluated. More accurate methods or ones which can significantly reduce the size of complex, transient thermal and structural problems are desirable for analysis and are required for synthesis of real structures subjected to thermal and mechanical loading. A unified method is presented for deriving successively higher-order modal solutions related to previously-developed, lower-order methods such as the mode displacement and mode-acceleration methods. A new method, called the force-derivative method, is used to obtain higher-order modal solutions for both uncoupled (proportionally-damped) structural problems as well as thermal problems and coupled (non-proportionally damped) structural problems. The new method is called the force-derivative method because, analogous to the mode-acceleration method, it produces a term that depends on the forcing function and additional terms that depend on the time derivatives of the forcing function.

Camarda, Charles J.↗

NASA NTRS1997

Development of Improved Surface Integral Methods for Jet Aeroacoustic Predictions

The accurate prediction of aerodynamically generated noise has become an important goal over the past decade. Aeroacoustics must now be an integral part of the aircraft design process. The direct calculation of aerodynamically generated noise with CFD-like algorithms is plausible. However, large computer time and memory requirements often make these predictions impractical. It is therefore necessary to separate the aeroacoustics problem into two parts, one in which aerodynamic sound sources are determined, and another in which the propagating sound is calculated. This idea is applied in acoustic analogy methods. However, in the acoustic analogy, the determination of far-field sound requires the solution of a volume integral. This volume integration again leads to impractical computer requirements. An alternative to the volume integrations can be found in the Kirchhoff method. In this method, Green's theorem for the linear wave equation is used to determine sound propagation based on quantities on a surface surrounding the source region. The change from volume to surface integrals represents a tremendous savings in the computer resources required for an accurate prediction. This work is concerned with the development of enhancements of the Kirchhoff method for use in a wide variety of aeroacoustics problems. This enhanced method, the modified Kirchhoff method, is shown to be a Green's function solution of Lighthill's equation. It is also shown rigorously to be identical to the methods of Ffowcs Williams and Hawkings. This allows for development of versatile computer codes which can easily alternate between the different Kirchhoff and Ffowcs Williams-Hawkings formulations, using the most appropriate method for the problem at hand. The modified Kirchhoff method is developed primarily for use in jet aeroacoustics predictions. Applications of the method are shown for two dimensional and three dimensional jet flows. Additionally, the enhancements are generalized so that they may be used in any aeroacoustics problem.

Pilon, Anthony R.↗

NASA NTRS2000

Evaluation of Non-Ozone-Depleting-Chemical Cleaning Methods for Space Mechanisms Using a Vacuum Spiral Orbit Rolling Contact Tribometer

Because CFC 113, an ozone depleting chemical (ODC), can no longer be produced, alternative bearing cleaning methods must be studied. The objective of this work was to study the effect of the new cleaning methods on lubricant lifetime using a vacuum bearing simulator (spiral orbit rolling contact tribometer). Four alternative cleaning methods were studied: ultra-violet (UV) ozone, aqueous levigated alumina slurry (ALAS), super critical fluid (SCF) CO2 and aqueous Brulin 815GD. Baseline tests were done using CFC 113. Test conditions were the following: a vacuum of at least 1.3 x 10(exp -6) Pa, 440C steel components, a rotational speed of 10 RPM, a lubricant charge of between 60-75 micrograms, a perfluoropolyalkylether lubricant (Z-25), and a load of 200N (44.6 lbs., a mean Hertzian stress of 1.5 GPa). Normalized lubricant lifetime was determined by dividing the total number of ball orbits by the amount of lubricant. The failure condition was a friction coefficient of 0.38. Post-test XPS analysis was also performed, showing slight variations in post-cleaning surface chemistry. Statistical analysis of the resultant data was conducted and it was determined that the data sets were most directly comparable when subjected to a natural log transformation. The natural log life (NL-Life) data for each cleaning method were reasonably normally (statistically) distributed and yielded standard deviations that were not significantly different among the five cleaning methods investigated. This made comparison of their NL-Life means very straightforward using a Bonferroni multiple comparison of means procedure. This procedure showed that the ALAS, UV-ozone and CFC 113 methods were not statistically significantly different from one another with respect to mean NL-Life. It also found that the SCF CO2 method yielded a significantly higher mean NL-Life than the mean NL-Lives of the ALAS, UV-ozone and CFC 113 methods. It also determined that the aqueous Brulin 815GD method yielded a mean NL-Life that was statistically significantly higher than the mean NL-Lives of each of the other four methods. Baseline tests using CFC 113 cleaned parts yielded a mean NL-Life 3.62 orbits/micro-g. ALAS and UV-ozone yielded similar mean NL-Life (3.31 orbits/mg and 3.33 orbits/micro-g, respectively). SCF CO2, gave a mean NL-Life of 4.08 orbits/mg and aqueous Brulin 8l5GD data yielded the longest mean NL-Life (4.66 orbits/micro-g).

Jansen, Mark J.↗

NASA NTRS2001

A Unified Development of Basis Reduction Methods for Rotor Blade Analysis

The axial foreshortening effect plays a key role in rotor blade dynamics, but approximating it accurately in reduced basis models has long posed a difficult problem for analysts. Recently, though, several methods have been shown to be effective in obtaining accurate,reduced basis models for rotor blades. These methods are the axial elongation method,the mixed finite element method, and the nonlinear normal mode method. The main objective of this paper is to demonstrate the close relationships among these methods, which are seemingly disparate at first glance. First, the difficulties inherent in obtaining reduced basis models of rotor blades are illustrated by examining the modal reduction accuracy of several blade analysis formulations. It is shown that classical, displacement-based finite elements are ill-suited for rotor blade analysis because they can't accurately represent the axial strain in modal space, and that this problem may be solved by employing the axial force as a variable in the analysis. It is shown that the mixed finite element method is a convenient means for accomplishing this, and the derivation of a mixed finite element for rotor blade analysis is outlined. A shortcoming of the mixed finite element method is that is that it increases the number of variables in the analysis. It is demonstrated that this problem may be rectified by solving for the axial displacements in terms of the axial forces and the bending displacements. Effectively, this procedure constitutes a generalization of the widely used axial elongation method to blades of arbitrary topology. The procedure is developed first for a single element, and then extended to an arbitrary assemblage of elements of arbitrary type. Finally, it is shown that the generalized axial elongation method is essentially an approximate solution for an invariant manifold that can be used as the basis for a nonlinear normal mode.

Ruzicka, Gene C.↗

NASA NTRS2001

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗

NASA NTRS2002

The Fractional Step Method Applied to Simulations of Natural Convective Flows

This paper describes research done to apply the Fractional Step Method to finite-element simulations of natural convective flows in pure liquids, permeable media, and in a directionally solidified metal alloy casting. The Fractional Step Method has been applied commonly to high Reynold's number flow simulations, but is less common for low Reynold's number flows, such as natural convection in liquids and in permeable media. The Fractional Step Method offers increased speed and reduced memory requirements by allowing non-coupled solution of the pressure and the velocity components. The Fractional Step Method has particular benefits for predicting flows in a directionally solidified alloy, since other methods presently employed are not very efficient. Previously, the most suitable method for predicting flows in a directionally solidified binary alloy was the penalty method. The penalty method requires direct matrix solvers, due to the penalty term. The Fractional Step Method allows iterative solution of the finite element stiffness matrices, thereby allowing more efficient solution of the matrices. The Fractional Step Method also lends itself to parallel processing, since the velocity component stiffness matrices can be built and solved independently of each other. The finite-element simulations of a directionally solidified casting are used to predict macrosegregation in directionally solidified castings. In particular, the finite-element simulations predict the existence of 'channels' within the processing mushy zone and subsequently 'freckles' within the fully processed solid, which are known to result from macrosegregation, or what is often referred to as thermo-solutal convection. These freckles cause material property non-uniformities in directionally solidified castings; therefore many of these castings are scrapped. The phenomenon of natural convection in an alloy under-going directional solidification, or thermo-solutal convection, will be explained. The development of the momentum and continuity equations for natural convection in a fluid, a permeable medium, and in a binary alloy undergoing directional solidification will be presented. Finally, results for natural convection in a pure liquid, natural convection in a medium with a constant permeability, and for directional solidification will be presented.

Westra, Douglas G.↗

NASA NTRS2003

Finite Element Method for Capturing Ultra-relativistic Shocks

While finite element methods are used extensively by researchers solving computational fluid dynamics in fields other than astrophysics, their use in astrophysical fluid simulations has been predominantly overlooked. Current simulations using other methods such as finite difference and finite volume (based on finite difference) have shown remarkable results, but these methods are limited by their fundamental properties in aspects that are important for simulations with complex geometries and widely varying spatial and temporal scale differences. We have explored the use of finite element methods for astrophysical fluids in order to establish the validity of using such methods in astrophysical environments. We present our numerical technique applied to solving ultra-relativistic (Lorentz Factor Gamma >> 1) shocks which are prevalent in astrophysical studies including relativistic jets and gamma-ray burst studies. We show our finite element formulation applied to simulations where the Lorentz factor ranges up to 2236 and demonstrate its stability in solving ultra-relativistic flows. Our numerical method is based on the Flowfield Dependent Variation (FDV) Method, unique in that numerical diffusion is derived from physical parameters rather than traditional artificial viscosity methods. Numerical instabilities account for most of the difficulties when capturing shocks in this regime. Our method results in stable solutions and accurate results as compared with other methods.

Richardson, G. A.↗

NASA NTRS2004

Discontinuous Spectral Difference Method for Conservation Laws on Unstructured Grids

A new, high-order, conservative, and efficient discontinuous spectral finite difference (SD) method for conservation laws on unstructured grids is developed. The concept of discontinuous and high-order local representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) and the Spectral Volume (SV) methods, but while these methods are based on the integrated forms of the equations, the new method is based on the differential form to attain a simpler formulation and higher efficiency. Conventional unstructured finite-difference and finite-volume methods require data reconstruction based on the least-squares formulation using neighboring point or cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every point or cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In addition, the finite-difference method does not satisfy the integral conservation in general. By contrast, the DG and SV methods employ a local, universal reconstruction of a given order of accuracy in each cell in terms of internally defined conservative unknowns. Since the solution is discontinuous across cell boundaries, a Riemann solver is necessary to evaluate boundary flux terms and maintain conservation. In the DG method, a Galerkin finite-element method is employed to update the nodal unknowns within each cell. This requires the inversion of a mass matrix, and the use of quadratures of twice the order of accuracy of the reconstruction to evaluate the surface integrals and additional volume integrals for nonlinear flux functions. In the SV method, the integral conservation law is used to update volume averages over subcells defined by a geometrically similar partition of each grid cell. As the order of accuracy increases, the partitioning for 3D requires the introduction of a large number of parameters, whose optimization to achieve convergence becomes increasingly more difficult. Also, the number of interior facets required to subdivide non-planar faces, and the additional increase in the number of quadrature points for each facet, increases the computational cost greatly.

Liu, Yen↗

NASA NTRS2010

The ICARE Method

The ICARE method is a flexible, widely applicable method for systems engineers to solve problems and resolve issues in a complete and comprehensive manner. The method can be tailored by diverse users for direct application to their function (e.g. system integrators, design engineers, technical discipline leads, analysts, etc.). The clever acronym, ICARE, instills the attitude of accountability, safety, technical rigor and engagement in the problem resolution: Identify, Communicate, Assess, Report, Execute (ICARE). This method was developed through observation of Space Shuttle Propulsion Systems Engineering and Integration (PSE&I) office personnel approach in an attempt to succinctly describe the actions of an effective systems engineer. Additionally it evolved from an effort to make a broadly-defined checklist for a PSE&I worker to perform their responsibilities in an iterative and recursive manner. The National Aeronautics and Space Administration (NASA) Systems Engineering Handbook states, engineering of NASA systems requires a systematic and disciplined set of processes that are applied recursively and iteratively for the design, development, operation, maintenance, and closeout of systems throughout the life cycle of the programs and projects. ICARE is a method that can be applied within the boundaries and requirements of NASA s systems engineering set of processes to provide an elevated sense of duty and responsibility to crew and vehicle safety. The importance of a disciplined set of processes and a safety-conscious mindset increases with the complexity of the system. Moreover, the larger the system and the larger the workforce, the more important it is to encourage the usage of the ICARE method as widely as possible. According to the NASA Systems Engineering Handbook, elements of a system can include people, hardware, software, facilities, policies and documents; all things required to produce system-level results, qualities, properties, characteristics, functions, behavior and performance. The ICARE method can be used to improve all elements of a system and, consequently, the system-level functional, physical and operational performance. Even though ICARE was specifically designed for a systems engineer, any person whose job is to examine another person, product, or process can use the ICARE method to improve effectiveness, implementation, usefulness, value, capability, efficiency, integration, design, and/or marketability. This paper provides the details of the ICARE method, emphasizing the method s application to systems engineering. In addition, a sample of other, non-systems engineering applications are briefly discussed to demonstrate how ICARE can be tailored to a variety of diverse jobs (from project management to parenting).

Henke, Luke↗

NASA NTRS2004

Improved Method of Locating Defects in Wiring Insulation

An improved method of locating small breaches in insulation on electrical wires combines aspects of the prior dielectric withstand voltage (DWV) and time-domain reflectometry (TDR) methods. The method was invented to satisfy a need for reliably and quickly locating insulation defects in spacecraft, aircraft, ships, and other complex systems that contain large amounts of wiring, much of it enclosed in structures that make it difficult to inspect. In the DWV method, one applies a predetermined potential (usually 1.5 kV DC) to the wiring and notes whether the voltage causes any arcing between the wiring and ground. The DWV method does not provide an indication of the location of the defect (unless, in an exceptional case, the arc happens to be visible). In addition, if there is no electrically conductive component at ground potential within about 0.010 in. (approximately equal to 0.254 mm) of the wire at the location of an insulation defect, then the DWV method does not provide an indication of the defect. Moreover, one does not have the option to raise the potential in an effort to increase the detectability of such a defect because doing so can harm previously undamaged insulation. In the TDR method as practiced heretofore, one applies a pulse of electricity having an amplitude of less than 25 V to a wire and measures the round-trip travel time for the reflection of the pulse from a defect. The distance along the wire from the point of application of the pulse to the defect is then calculated as the product of half the round-trip travel time and the characteristic speed of a propagation of an electromagnetic signal in the wire. While the TDR method as practiced heretofore can be used to locate a short or open circuit, it does not ordinarily enable one to locate a small breach in insulation because the pulse voltage is too low to cause arcing and thus too low to induce an impedance discontinuity large enough to generate a measurable reflection. The present improved method overcomes the weaknesses of both the prior DWV and the prior TDR method.

Greulich, Owen R.↗

NASA NTRS2016

Diagonally Implicit Runge-Kutta Methods for Ordinary Differential Equations. A Review

A review of diagonally implicit Runge-Kutta (DIRK) methods applied to rst-order ordinary di erential equations (ODEs) is undertaken. The goal of this review is to summarize the characteristics, assess the potential, and then design several nearly optimal, general purpose, DIRK-type methods. Over 20 important aspects of DIRKtype methods are reviewed. A design study is then conducted on DIRK-type methods having from two to seven implicit stages. From this, 15 schemes are selected for general purpose application. Testing of the 15 chosen methods is done on three singular perturbation problems. Based on the review of method characteristics, these methods focus on having a stage order of two, sti accuracy, L-stability, high quality embedded and dense-output methods, small magnitudes of the algebraic stability matrix eigenvalues, small values of aii, and small or vanishing values of the internal stability function for large eigenvalues of the Jacobian. Among the 15 new methods, ESDIRK4(3)6L[2]SA is recommended as a good default method for solving sti problems at moderate error tolerances.

Kennedy, Christopher A.↗