Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical schemes”

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,261 records · Page 70

Stress and Fracture Analyses Under Elastic-plastic and Creep Conditions: Some Basic Developments and Computational Approaches

A new hybrid-stress finite element algorith, suitable for analyses of large quasi-static deformations of inelastic solids, is presented. Principal variables in the formulation are the nominal stress-rate and spin. A such, a consistent reformulation of the constitutive equation is necessary, and is discussed. The finite element equations give rise to an initial value problem. Time integration has been accomplished by Euler and Runge-Kutta schemes and the superior accuracy of the higher order schemes is noted. In the course of integration of stress in time, it has been demonstrated that classical schemes such as Euler's and Runge-Kutta may lead to strong frame-dependence. As a remedy, modified integration schemes are proposed and the potential of the new schemes for suppressing frame dependence of numerically integrated stress is demonstrated. The topic of the development of valid creep fracture criteria is also addressed.

Reed, K. W.↗

Convenient stability criteria for difference approximations of hyperbolic initial-boundary value problems

New convenient stability criteria are provided in this paper for a large class of finite difference approximations to initial-boundary value problems associated with the hyperbolic system u sub t = au sub x + Bu + f in the quarter plane x or = 0, t or = 0. Using the new criteria, stability is easily established for numerous combinations of well known basic schemes and boundary conditions, thus generalizing many special cases studied in recent literature.

Goldberg, M.↗

Computation of transonic viscous-inviscid interacting flow

Transonic viscous-inviscid interaction is considered using the Euler and inverse compressible turbulent boundary-layer equations. Certain improvements in the inverse boundary-layer method are mentioned, along with experiences in using various Runge-Kutta schemes to solve the Euler equations. Numerical conditions imposed on the Euler equations at a surface for viscous-inviscid interaction using the method of equivalent sources are developed, and numerical solutions are presented and compared with experimental data to illustrate essential points. Previously announced in STAR N83-17829

Whitfield, D. L.↗

Second-order explicit finite-difference methods for transient-flow analysis

Three second-order accurate numerical methods - MacCormack's method, Lambda scheme and Gabutti scheme - are introduced to solve the quasi-linear, hyperbolic partial differential equations describing transient flows in closed conduits. The details of these methods and the treatment of boundary conditions are presented and the results computed by using these methods for a typical piping system are compared. It is shown that for the same accuracy, second-order methods require considerably lesser number of computational nodes and computer time as compared to those required by the first-order methods.

Chaudhry, M. H.↗

Discrete approximation methods for parameter identification in delay systems

Approximation schemes for parameter identification problems in which the governing state equation is a linear functional differential equation of retarded type are constructed. The basis of the schemes is the replacement of the parameter identification problem having an infinite dimensional state equation by a sequence of approximating parameter identification problems in which the states are given by finite dimensional discrete difference equations. The difference equations are constructed using linear semigroup theory and rational function approximations to the exponential. Sufficient conditions are given for the convergence of solutions to the approximating problems, which can be obtained using conventional methods, to solutions to the original parameter identification problem. Finite difference and spline based schemes using Paderational function approximations to the exponential are constructed, and shown to satisfy the sufficient conditions for convergence. A discussion and analysis of numerical results obtained through the application of the schemes to several examples is included.

Rosen, I. G.↗

Difference equation state approximations for nonlinear hereditary control problems

Discrete approximation schemes for the solution of nonlinear hereditary control problems are constructed. The methods involve approximation by a sequence of optimal control problems in which the original infinite dimensional state equation has been approximated by a finite dimensional discrete difference equation. Convergence of the state approximations is argued using linear semigroup theory and is then used to demonstrate that solutions to the approximating optimal control problems in some sense approximate solutions to the original control problem. Two schemes, one based upon piecewise constant approximation, and the other involving spline functions are discussed. Numerical results are presented, analyzed and used to compare the schemes to other available approximation methods for the solution of hereditary control problems. Previously announced in STAR as N83-33589

Rosen, I. G.↗

A 3-dimensional mass conserving element for compressible flows

A variety of finite element schemes has been used in the numerical approximation of compressible flows particularly in underwater acoustics. In many instances instabilities have been generated due to the lack of mass conservation. Two- and three-dimensional elements are developed which avoid these problems.

Fix, G.↗

Convenient stability criteria for difference approximations of hyperbolic initial-boundary value problems

New convenient stability criteria are provided in this paper for a large class of finite difference approximations to initial-boundary value problems associated with the hyperbolic system u sub t = au sub x + Bu + f in the quarter plane x or = 0, t or = 0. Using the new criteria, stability is easily established for numerous combinations of well known basic schemes and boundary conditins, thus generalizing many special cases studied in recent literature.

Goldberg, M.↗

New stability criteria for difference approximations of hyperbolic initial-boundary value problems

New convenient stability criteria are provided for a large class of finite difference approximations to initial-boundary value problems associated with the hyperbolic system u sub t = Au sub x + Bu + f in the quarter plane x greater than or equal to 0, t greater than or equal to 0. The criteria are used to easily establish stability for numerous combinations of well known basic schemes and boundary conditions, thus generalizing many special cases studied in the recent literature. A number of examples are examined, including the unitary unconditionally stable Crank-Nicholson scheme and an almost-dissipative unconditionally stable backward Euler scheme.

Goldberg, M.↗

Synopsis of the Fifth Conference on the Meteorology of the Stratosphere and Mesosphere

The papers presented at the Fifth Conference on the Meteorology of the Stratosphere and Mesosphere held on April 23-26, 1985, are reviewed. The observational aspects of large-scale circulation, such as summer and winter circulation in the Southern Hemisphere, and analysis schemes, like the multivariate statistical analysis scheme, are discussed. The topics of numerical simulations of the general circulation and sudden-warming are examined. Papers concerning processes of O3, NO2, H2, and HNO3 are described. Research on large-scale mixing processes in the stratosphere is presented. The topic of equatorial dynamics and stability is analyzed. Papers focusing on the effect of gravity waves on the general circulation are studied.

Schoeberl, M. R.↗

A three-dimensional boundary-layer method for flow over delta wings with leading-edge separation

A three-dimensional, laminar boundary-layer method is applied to the incompressible flow over a slender delta wing at incidence. The predictor-corrector finite-difference scheme of Matsuno is used to difference the governing equations. The method has the advantages that no iterations are required to advance the solution and the cross-flow derivatives are formed independent of the cross-field direction. The difference scheme is demonstrated to yield accurate numerical results when compared to the exact solution of the three-dimensional boundary-layer equations for parabolic flow over a moving flat plate. The method is applied to delta wings of various sweep angles at angles of attack up to 20 deg., with the inviscid solution determined using a higher-order, three-dimensional panel method.

Woodson, S. H.↗

Computational fluid mechanics utilizing the variational principle of modeling damping seals

A computational fluid dynamics code for application to traditional incompressible flow problems has been developed. The method is actually a slight compressibility approach which takes advantage of the bulk modulus and finite sound speed of all real fluids. The finite element numerical analog uses a dynamic differencing scheme based, in part, on a variational principle for computational fluid dynamics. The code was developed in order to study the feasibility of damping seals for high speed turbomachinery. Preliminary seal analyses have been performed.

Abernathy, J. M.↗

The minimal residual QR-factorization algorithm for reliably solving subset regression problems

A new algorithm to solve test subset regression problems is described, called the minimal residual QR factorization algorithm (MRQR). This scheme performs a QR factorization with a new column pivoting strategy. Basically, this strategy is based on the change in the residual of the least squares problem. Furthermore, it is demonstrated that this basic scheme might be extended in a numerically efficient way to combine the advantages of existing numerical procedures, such as the singular value decomposition, with those of more classical statistical procedures, such as stepwise regression. This extension is presented as an advisory expert system that guides the user in solving the subset regression problem. The advantages of the new procedure are highlighted by a numerical example.

Verhaegen, M. H.↗

Wave envelope technique for multimode wave guide problems

A numerical solution technique for two-dimensional EM waveguide problems is developed analytically and demonstrated. The frequency-domain waveguide problem is transformed to a time-dependent problem and treated following the limiting-amplitude principle; the wave-envelope technique is applied; and the resulting problem is solved numerically using a finite-difference scheme based on a variation of the finite-area method (Baumeister, 1977 and 1984). Results for several sample problems, obtained in 10 min or less of execution time on PCs, are presented in tables.

Hariharan, S. I.↗

Evaluation and applicability of transfinite element formulations with particular reference to radiation

The evaluation and applicability of hybrid transfinite element formulations are described for transient nonlinear thermal models with particular reference to account for radiation effects. The formulations are developed from basic concepts and evaluated for both one- and two-dimensional thermal models. The methodology presented is a hybrid approach as it combines finite elements, classical Galerkin schemes, and transform methods. Basic features of the hybrid formulations and the associated solution scheme used are first described. Therein, numerical test models are presented to evaluate the applicability for predicting the nonlinear transient response in one- and two-dimensional thermal models influenced due to radiation effects.

Tamma, Kumar K.↗

The evolution and discharge of electric fields within a thunderstorm

An analysis of the present three-dimensional thunderstorm electrical model and its finite-difference approximations indicates unconditional stability for the discretization that results from the approximation of the spatial derivatives by a box-schemelike method and of the temporal derivative by either a backward-difference or Crank-Nicholson scheme. Lightning propagation is treated through numerical techniques based on the inverse-matrix modification formula and Cholesky updates. The model is applied to a storm observed at the Kennedy Space Center in 1978, and numerical comparisons are conducted between the model and the theoretical results obtained by Wilson (1920) and Holzer and Saxon (1952).

Hager, William W.↗

Combustion of liquid fuel droplets in supercritical conditions

A comprehensive analysis of liquid-fuel droplet combustion in both sub- and super-critical environments has been conducted. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates finite-rate chemical kinetics and a full treatment of liquid-vapor phase equilibrium at the droplet surface. The governing equations and the associated interface boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to the supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures of 5-140 atm. Results indicate that the ambient gas pressure exerts significant control of droplet gasification and burning processes through its influences on the fluid transport, gas/liquid interface thermodynamics, and chemical reactions. The droplet gasification rate increases progressively with pressure. However, the data for the overall burnout time exhibits a significant variation near the critical burning pressure, mainly as a result of reduced mass-diffusion rate and latent heat of vaporization with increased pressure. The influence of droplet size on the burning characteristics is also noted.

Shuen, J. S.↗

Vaporization and combustion of fuel droplets at supercritical conditions

Vaporization and combustion liquid-fuel droplets in both sub- and super-critical environments have been examined. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates finite-rate chemical kinetics and a full treatment of liquid-vapor phase equilibrium at the droplet surface. The governing equations and the associated interface boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to the supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures of 5-140 atm. In addition, the dynamic responses of droplet vaporization and combustion to ambient-pressure oscillations are investigated. Results indicate that the droplet gasification and burning mechanisms depend greatly on the ambient pressure. In particular, a rapid enlargement of the vaporization and combustion responses occurs when the droplet surface reaches its critical point, mainly due to the strong variations of latent heat of vaporization and thermophysical properties at the critical state.

Yang, Vigor↗