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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Numerical study of an explosion in a non-homogeneous medium with and without magnetic fields

A version of the two step Lax-Wendroff difference method with second order accuracy is used to seek solutions of the unsteady magnetohydrodynamic (MHD) equations for the study of an explosion in a non-homogeneous medium with and without magnetic fields. The explosion is initiated by introducing a finite amount of energy in a small volume of gas which leads to an instantaneous increase of gas temperature. Two types of magnetic field configurations (i.e. open and closed) are considered to illustrate the dependence and differences of magnetic and a non-magnetic flow motion. Numerical results show the strong dependence of induced MHD flow field on the magnetic field configuration and strength. The development and decay of the fast and slow MHD shock waves, as well as the ordinary gasdynamic shock waves, are also represented well by the present computer simulations. In conclusion, we have demonstrated that the numerical scheme we have used is a reliable one for studying 2-dimensional MHD problems in non-homogeneous medium.

Han, S. M.↗
NASA NTRS1980

Iterative methods based upon residual averaging

Iterative methods for solving boundary value problems for systems of nonlinear partial differential equations are discussed. The methods involve subtracting an average of residuals from one approximation in order to arrive at a subsequent approximation. Two abstract methods in Hilbert space are given and application of these methods to quasilinear systems to give numerical schemes for such problems is demonstrated. Potential theoretic matters related to the iteration schemes are discussed.

Neuberger, J. W.↗
NASA NTRS1980

Numerical techniques for solving nonlinear instability problems in smokeless tactical solid rocket motors

The selection of a satisfactory numerical method for calculating the propagation of steep fronted shock life waveforms in a solid rocket motor combustion chamber is discussed. A number of different numerical schemes were evaluated by comparing the results obtained for three problems: the shock tube problems; the linear wave equation, and nonlinear wave propagation in a closed tube. The most promising method--a combination of the Lax-Wendroff, Hybrid and Artificial Compression techniques, was incorporated into an existing nonlinear instability program. The capability of the modified program to treat steep fronted wave instabilities in low smoke tactical motors was verified by solving a number of motor test cases with disturbance amplitudes as high as 80% of the mean pressure.

Baum, J. D.↗
NASA NTRS1981

A new payload integration method

This paper introduces a new payload integration method. The system equations of motion are derived in terms of interface-restrained booster and payload modes. These system equations are directly integrated, avoiding the solution of an expensive eigenvalue problem. A modified Newmark-Chan-Beta numerical scheme is used to perform this integration. Due to the unique form of the system equations it is possible to achieve significant savings. A similar savings is feasible for the computation of the internal loads. This new approach does not involve approximations. Although completely general, the presented technique is most effective when the number of interface degrees of freedom is relatively small. Both determinate and indeterminate interfaces are allowed.

Engels, R. C.↗
NASA NTRS1981

Analysis of cracks emanating from a circular hole in unidirectional fiber reinforced composites, part 2

An analytical method is developed for cracks emanating from a circular hole in an off-axis unidirectional fiber-reinforced composite. The method which is formulated by using conservation laws of elasticity and fundamental relationships in anisotropic fracture mechanics, provides a convenient and accurate means to examine the complicated crack behavior, when used in conjunction with a suitable numerical scheme such as the finite element method. The formulation is eventually reduced to a system of linear algebraic equations of mixed-mode stress intensity factors. Fracture parameters, describing crack-tip deformation and fracture in the composite, are obtained explicitly. Effects of material anisotropy and crack/hole geometry are examined also. Of particular interest are the energy release rates associated with crack extension; their values are evaluated for various cases. Results show that mixed-mode stress intensity factors and energy release rates associated with the cracks emanating from a hole change very appreciably with fiber orientation in the composite. K sub 1 and G increase monotonically with increasing theta; but K sub 2 reaches its maximum at theta = 45 deg, and then decreases gradually as theta increases further.

Wang, S. S.↗
NASA NTRS1982

Conjugated heat transfer in a laminar boundary layer with heat source at the wall

This work presents the solution of the temperature field in a two-dimensional laminar incompressible flow over a conducting solid plate with a line heat source located at the fluid-solid interface perpendicular to the flow direction. A numerical scheme was used to obtain the temperature profiles as a function of the source strength, and of the properties of the fluid and the solid. The heat conduction and forced convection in the fluid and the heat conduction in the solid were solved for the case of moderate temperature rise, where the assumption of constant properties applies. The model enables the improvement of an instrument for the detection of boundary layer separation. It was found that for the actual parameters of the separation detector, a distance of 4 to 24 mm between the sensors gives an indication of 70 percent of the maximum temperature difference.

Brosh, A.↗
NASA NTRS1982

An iterative radiative transfer code for ocean-atmosphere systems

The details of an iterative radiative transfer code for computing the intensity and degree of polarization of diffuse radiation in models of the ocean-atmosphere system are described. The present code neglects the upwelling radiation from below the ocean surface and as such can be applied to the part of the spectrum where the absorption by water is strong. To establish the reliability of the numerical scheme and the computer code, the results are compared with those of Fraser and Walker (1968), Dave (1972), and Mullamaa (1964); they are found to be in excellent agreement. The computations also show that both the intensity and the degree of polarization of the upwelling diffuse radiation at the top of the atmosphere vary significantly when the rough ocean at the base of the atmosphere is replaced by a Lambertian surface that reflects the same energy as the rough ocean.

Ahmad, Z.↗
NASA NTRS1981

Numerical analysis of turbulent coaxial flow with internal heat generation

A computational method with which to obtain a physical understanding of the turbulent field of two coaxial jets entering an axisymmetric chamber is developed. Even the laminar field of this flow is quite complicated. This is due to the many different domains which exist in the field especially in the entrance region. Physically, three regions may be identified: the wall region, the initial region near the axis of symmetry and the mixing region. Advancing downstream, these regions change relative size with the ratio of the two jets' mass fluxes as the main parameter. The turbulent field of these flows is much more complicated due to the difference in the effective transport coefficients and turbulence level from region to region. However, being aware beforehand of the complications and the different regions of this field, the appropriate turbulence model and numerical scheme can be adjusted to treat the problem.

Lin, A.↗
NASA NTRS1984

Spectral methods in time for hyperbolic equations

A pseudospectral numerical scheme for solving linear, periodic, hyperbolic problems is described. It has infinite accuracy both in time and in space. The high accuracy in time is achieved without increasing the computational work and memory space which is needed for a regular, one step explicit scheme. The algorithm is shown to be optimal in the sense that among all the explicit algorithms of a certain class it requires the least amount of work to achieve a certain given resolution. The class of algorithms referred to consists of all explicit schemes which may be represented as a polynomial in the spatial operator.

Tal-Ezer, H.↗
NASA NTRS1984

Contact stresses in pin-loaded orthotropic plates

The effects of pin elasticity, friction, and clearance on the stresses near the hole in a pin-loaded orthotropic plate are described. The problem is modeled as a contact elasticity problem using complex variable theory, the pin and the plate being two elastic bodies interacting through contact. This modeling is in contrast to previous works which assumed that the pin is rigid or that it exerts a known cosinusoidal radial traction on the hole boundary. Neither of these approaches explicitly involves a pin. A collocation procedure and iteration were used to obtain numerical results for a variety of plate and pin elastic properties and various levels of friction and clearance. Collocation was used to enforce the boundary and iteration was used to find the contact and no-slip regions on the boundary. Details of the numerical scheme are discussed.

Hyer, M. W.↗
NASA NTRS1984

A numerical study of the controlled flow tunnel for a high lift model

A controlled flow tunnel employs active control of flow through the walls of the wind tunnel so that the model is in approximately free air conditions during the test. This improves the wind tunnel test environment, enhancing the validity of the experimentally obtained test data. This concept is applied to a three dimensional jet flapped wing with full span jet flap. It is shown that a special treatment is required for the high energy wake associated with this and other V/STOL models. An iterative numerical scheme is developed to describe the working of an actual controlled flow tunnel and comparisons are shown with other available results. It is shown that control need be exerted over only part of the tunnel walls to closely approximate free air flow conditions. It is concluded that such a tunnel is able to produce a nearly interference free test environment even with a high lift model in the tunnel.

Parikh, P. C.↗
NASA NTRS1980

Solution of 3-dimensional time-dependent viscous flows. Part 2: Development of the computer code

There is considerable interest in developing a numerical scheme for solving the time dependent viscous compressible three dimensional flow equations to aid in the design of helicopter rotors. The development of a computer code to solve a three dimensional unsteady approximate form of the Navier-Stokes equations employing a linearized block emplicit technique in conjunction with a QR operator scheme is described. Results of calculations of several Cartesian test cases are presented. The computer code can be applied to more complex flow fields such as these encountered on rotating airfoils.

Weinberg, B. C.↗
NASA NTRS1984

Legendre-tau approximation for functional differential equations. Part 2: The linear quadratic optimal control problem

The numerical scheme based on the Legendre-tau approximation is proposed to approximate the feedback solution to the linear quadratic optimal control problem for hereditary differential systems. The convergence property is established using Trotter ideas. The method yields very good approximations at low orders and provides an approximation technique for computing closed-loop eigenvalues of the feedback system. A comparison with existing methods (based on averaging and spline approximations) is made.

Ito, K.↗
NASA NTRS1984

Boundary layer transition effects on flow separation around V/STOL engine inlets at high incidence

Numerical methods for calculating laminar and turbulent boundary layer development around vertical-short take off and landing engine inlets at high incidence angles are investigated. Various transition models were compared and evaluated in calculations off flow separation bound inside the inlet. Results of the transition effects on the boundary layer characteristics at onset of separation for two types of engine inlet geometries are presented. Some of the numerical results are compared with existing wind-tunnel test data for scaled inlet models to demonstrate the effects of transition models in the numerical scheme. The effects of transition modeling on the boundary layer development are illustrated for typical engine operating conditions.

Chou, D. C.↗
NASA NTRS1984

A two-dimensional analysis of laser heat addition in converging nozzles

The two-dimensional equations of motion describing the interaction between a laser beam and a flowing gas are considered. An implicit numerical scheme is used to solve these equations for unchoked flow through a converging-diverging nozzle. Separate grids are used for the fluid dynamics and the radiation equations. The effects of beam focusing and cross-beam intensity profiles are included. The calculations are based upon real gas properties for all quantities except the gas absorptivity, which is taken as a constant. The solutions contain the expected hot central core region with cool gas near the walls. This results in steep temperature gradients in both the streamwise and cross-stream directions. The absorption zone acts as a blockage in the nozzle causing a nonuniform velocity profile at the inlet and an overall decrease in mass flow. The absorption region also forces the streamlines to move away from the axis of symmetry, although this effect is not strong.

Molvik, G. A.↗
NASA NTRS1983

An eigenvalue method for solving transient heat conduction problems

The eigenvalue method, which has been used by researchers in structure mechanics, is applied to problems in heat conduction. Its formulation is decribed in terms of an examination of transient heat conduction in a square slab. Taking advantage of the availability of the exact solution, we compare the accuracy and other numerical properties of the eigenvalue method with those of existing numerical schemes. The comparsion shows that, overall, the eigenvalue method appears to be fairly attractive. Furthermore, only a few dominant eigenvalues and their corresponding eigenvectors need to be computed and retained to yield reasonably high accuracy. Greater savings are attained in the computation time for a transient problem with long time duration and a large computational domain.

Shih, T. M.↗
NASA NTRS1982

Numerical solution of viscous flow in unbounded fluid

Incompressible viscous flow fields induced by initial vorticity distributions with bounded support or exponential decay in the far field are investigated. A numerical scheme for the solution of the vorticity distribution and the velocity field is presented with special emphasis on the treatment of the boundary data. The efficiency of the scheme is demonstrated. The present method has been applied to the study of the merging and collision of vortex rings.

Liu, C. H.↗
NASA NTRS1984

Dynamic response of shock waves in transonic diffuser and supersonic inlet - An analysis with the Navier-Stokes equations and adaptive grid

An existing method which solves the multi-dimensional ensemble-averaged compressible time-dependent Navier-Stokes equations in conjunction with mixing length turbulence model and shock capturing technique has been extended to include the shock-tracking adaptive grid systems. The numerical scheme for solving the governing equations is based on a linearized block implicit approach. The effects of grid-motion and grid-distribution on the calculated flow solutions have been studied in relative detail and this is carried out in the context of physically steady, shocked flows computed with non-stationary grids. Subsequently, the unsteady dynamics of the flows occurring in a supercritically operated transonic diffuser and a mixed compression supersonic inlet have been investigated with the adaptive grid systems by solving the Navier-Stokes equations.

Liu, N.-S.↗