SEARCH · Engineering Papers
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.
Numerical methods and computer programs for two-dimensional and axisymmetric supersonic flow over blunt-nosed and flared bodies
A numerical method is presented for calculating the boundary-induced interference at subsonic speeds. The wind tunnel slot width or wall porosity can vary throughout the test section. The interference can be computed at any point in the test section. The model can be any configuration and can be located anywhere in the test section. Several examples are given, and comparison is made with other methods where available.
A numerical method, based on linearized theory, for designing minimum-drag supersonic wing camber surfaces of arbitrary planform for a given lift, with options for constraining the pitching moment and/or the surface deformation at the trailing edge of the root chord and for selecting any desired combination of eight specified wing-loading distributions to be employed in the optimization procedure is presented. Two examples are given to illustrate applications of the method. The results indicate that relatively small drag penalties are incurred in designing wings to be self-trimming and to have a reasonable camber surface.
The long-standing challenge of creating a Milky Way- (MW-) like disk galaxy from cosmological simulations has motivated significant developments in both numerical methods and physical models. We investigate these two fundamental aspects in a new comparison project using a set of cosmological hydrodynamic simulations of an MW-sized galaxy. In this study, we focus on the comparison of two particle-based hydrodynamics methods: an improved smoothed particle hydrodynamics (SPH) code Gadget, and a Lagrangian Meshless Finite-Mass (MFM) code Gizmo. All the simulations in this paper use the same initial conditions and physical models, which include star formation, “energy-driven” outflows, metal-dependent cooling, stellar evolution, and metal enrichment. We find that both numerical schemes produce a late-type galaxy with extended gaseous and stellar disks. However, notable differences are present in a wide range of galaxy properties and their evolution, including star-formation history, gas content, disk structure, and kinematics. Compared to Gizmo, the Gadget simulation produced a larger fraction of cold, dense gas at high redshift which fuels rapid star formation and results in a higher stellar mass by 20% and a lower gas fraction by 10% at z = 0, and the resulting gas disk is smoother and more coherent in rotation due to damping of turbulent motion by the numerical viscosity in SPH, in contrast to the Gizmo simulation, which shows a more prominent spiral structure. Given its better convergence properties and lower computational cost, we argue that the MFM method is a promising alternative to SPH in cosmological hydrodynamic simulations.
Numerical method for calculating Doppler effect of radio waves coherently emitted from satellite
Flexible conductive textiles have attracted an increasing amount of attention due to being flexible and stretchable with multi-functions, in addition to the strain sensing function. In this paper, an integrated numerical method is established to predict the strain sensing behavior of PPy-coated flexible conductive fibers. Based on our previous investigation on the strain sensing behavior of the PPy-coated elastic conductive fibers, it has been known that the main strain sensing mechanisms of the conductive fibers are micro-cracks opening under stretching and closing under unloading. Therefore, in the numerical approach, extended finite element method is used to simulate the multi-crack initiation and propagation. Then, by establishing a representative element model including a single crack, the resistance variation with the increase in crack width/depth can be obtained using the electric–thermal package in ABAQUS software. Combining the relationship between crack width and strain, the relationship between fractional increment in resistance (ΔR/R{sub 0}) and the strain can be obtained. The prediction of the fractional increment in resistance (ΔR/R{sub 0}) versus strain is consistent with the experimental results of PPy-coated Lycra fibers. The new integrated method is further verified by experimental results of PPy-coated XLA fibers, showing that the integrated numerical approach can capture the complicated strain sensing mechanisms of PPy-coated elastic fibers and predict the strain sensing behavior of flexible conductive fibers. This approach will make it possible to shorten the period, and reduce the cost in R&D of the PPy-coated electrically conductive fibers, comparing with the corresponding experimental work.
The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.
The following papers are presented: (1) numerical methods for the simulation of complex multi-body flows with applications for the Integrated Space Shuttle vehicle; (2) a generalized scheme for 3-D hyperbolic grid generation; (3) collar grids for intersecting geometric components within the Chimera overlapped grid scheme; and (4) application of the Chimera overlapped grid scheme to simulation of Space Shuttle ascent flows.
Differencing operators of arbitrarily high order can be constructed by interpolating a polynomial through a set of data followed by differentiation of this polynomial and finally evaluation of the polynomial at the point where a derivative approximation is desired. Furthermore, the interpolating polynomial can be constructed from algebraic, trigonometric, or, perhaps exponential polynomials. This paper begins with a comparison of such differencing operator construction. Next, the issue of proper grids for high order polynomials is addressed. Finally, an adaptive numerical method is introduced which adapts the numerical grid and the order of the differencing operator depending on the data. The numerical grid adaptation is performed on a Chebyshev grid. That is, at each level of refinement the grid is a Chebvshev grid and this grid is refined locally based on wavelet analysis.
These are slides for a lecture for the Parallel Computing Summer Research Internship at the National Security Education Center. This gives an introduction to numerical methods. Repetitive algorithms are used to obtain approximate solutions to mathematical problems, using sorting, searching, root finding, optimization, interpolation, extrapolation, least squares regression, Eigenvalue problems, ordinary differential equations, and partial differential equations. Many equations are shown. Discretizations allow us to approximate solutions to mathematical models of physical systems using a repetitive algorithm and introduce errors that can lead to numerical instabilities if we are not careful.
This paper examines some basic considerations underlying dynamic shell response analysis and the impact of these considerations upon the practical aspects of solution by numerical methods. Emphasis is placed on the solution of linear problems. The present states of development of the finite difference and finite element methods are reviewed, and techniques for the treatment of temporal variation are discussed. An examination is made of the frequency parameters characteristic of thin shell theory, applied excitations, and spatial mesh geometries, and the significance of these parameters with respect to computational convergence is illustrated.
Numerical optimization methods compared for nonlinear two point boundary value problem associated with optimum spacecraft trajectory
A fourth-order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method is the natural extension of the second-order Keller Box Scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary-layer equations. The efficiency of the present method is compared with other two-point and three-point higher-order methods; namely, the Keller Box Scheme with Richardson extrapolation, the method of deferred corrections, and the three-point spline methods. For equivalent accuracy, numerical results show the present method to be more efficient than the other higher-order methods for both laminar and turbulent flows.
A fourth-order method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method is the natural extension of the second-order Keller Box Scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary-layer equations for both attached and separated flows. The efficiency of the present method is compared with other higher-order methods; namely, the Keller Box Scheme with Richardson extrapolation, the method of deferred corrections, the three-point spline methods, and a modified finite-element method. For equivalent accuracy, numerical results show the present method to be more efficient than the other higher-order methods for both laminar and turbulent flows.
A fourth order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method, which is the natural extension of the second order box scheme to fourth order, was demonstrated with application to the incompressible, laminar and turbulent, boundary layer equations. The efficiency of the present method is compared with two point and three point higher order methods, namely, the Keller box scheme with Richardson extrapolation, the method of deferred corrections, a three point spline method, and a modified finite element method. For equivalent accuracy, numerical results show the present method to be more efficient than higher order methods for both laminar and turbulent flows.
The calculation of flow fields past aircraft configuration at flight Reynolds numbers is considered. Progress in devising accurate and efficient numerical methods, in understanding and modeling the physics of turbulence, and in developing reliable and powerful computer hardware is discussed. Emphasis is placed on efficient solutions to the Navier-Stokes equations.
The mathematical structure of hyperbolic systems and the scalar equation case of conservation laws are discussed. Linear, nonlinear systems and the Riemann problem for the Euler equations are also studied. The numerical methods for conservation laws are presented in a nonstandard manner which leads to large time steps generalizations and computations on irregular grids. The solution of conservation laws with stiff source terms is examined.
Sonic boom wave pressure history prediction on arbitrarily oriented plane walls by explicit fixed mesh time dependent numerical method