Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite difference methods”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Time-marching solution of incompressible Navier-Stokes equations for internal flow

Primitive variables with central differencing on a staggered grid are used in the present, factored ADI finite-difference scheme for artificial compressibility method solution of the incompressible Navier-Stokes equations, leading to a close coupling between velocity and pressure that both enhances stability and eliminates the need for artificial damping. Computational efficiency is enhanced through the use of a spatially variable, fixed Courant number-based time-step. The numerical results obtained for a driven cavity at Re of 10,000, with local cell Re as high as 100, exhibits no flow variable spatial oscillations on a 40 x 40 stretched grid solution.

Soh, W. Y.↗
NASA NTRS1989

Stable explicit schemes for equations of Schroedinger type

A method for constructing explicit finite-difference schemes which can be used to solve Schroedinger-type partial-differential equations is presented. A forward Euler scheme that is conditionally stable is given by the procedure. The results presented are based on the analysis of the simplest Schroedinger type equation.

Mickens, Ronald E.↗
NASA NTRS1989

A CFD study of tilt rotor flowfields

The download on the wing produced by the rotor wake of a tilt rotor vehicle in hover is of major concern because of its severe impact on payload-carrying capability. In a concerted effort to understand the fundamental fluid dynamics that cause this download, and to help find ways to reduce it, computational fluid dynamics (CFD) is employed to study this problem. The thin-layer Navier-Stokes equations are used to describe the flow, and an implicit, finite difference numerical algorithm is the method of solution. The methodology is developed to analyze the tilt rotor flowfield. Included are discussions of computations of an airfoil and wing in freestream flows at -90 degrees, a rotor alone, and wing/rotor interaction in two and three dimensions. Preliminary results demonstrate the feasibility and great potential of the present approach. Recommendations are made for both near-term and far-term improvements to the method.

Fejtek, Ian↗
NASA NTRS1990

Evolution of midplate hotspot swells: Numerical solutions

The evolution of midplate hotspot swells on an oceanic plate moving over a hot, upwelling mantle plume is numerically simulated. The plume supplies a Gaussian-shaped thermal perturbation and thermally-induced dynamic support. The lithosphere is treated as a thermal boundary layer with a strongly temperature-dependent viscosity. The two fundamental mechanisms of transferring heat, conduction and convection, during the interaction of the lithosphere with the mantle plume are considered. The transient heat transfer equations, with boundary conditions varying in both time and space, are solved in cylindrical coordinates using the finite difference ADI (alternating direction implicit) method on a 100 x 100 grid. The topography, geoid anomaly, and heat flow anomaly of the Hawaiian swell and the Bermuda rise are used to constrain the models. Results confirm the conclusion of previous works that the Hawaiian swell can not be explained by conductive heating alone, even if extremely high thermal perturbation is allowed. On the other hand, the model of convective thinning predicts successfully the topography, geoid anomaly, and the heat flow anomaly around the Hawaiian islands, as well as the changes in the topography and anomalous heat flow along the Hawaiian volcanic chain.

Liu, Mian↗
NASA NTRS1992

Sensitivity analysis of dynamic aeroelastic responses

This paper summarizes ongoing research on the sensitivity analysis of dynamic aeroelastic response of wings. Two approaches are being used to express the unsteady aerodynamic loads: (1) the frequency-domain approach, and (2) the state-space approach. The frequency-domain approach is demonstrated on a three-dimensional box wing and the state-space domain approach is demonstrated on a simple two-dimensional sectional model. Three different methods are used to find the sensitivities: (1) a purely finite difference approach, (2) a semi-analytical approach in which an analytical expression is used for calculating the sensitivity of an eigenvalue of the complex valued aeroelastic matrix, however, the derivatives of the components of the matrix are obtained using finite difference, and (3) a semi-analytic approach that differs from (2) in the sense that the sensitivity of the aerodynamic matrix is now obtained analytically. A good agreement is seen between the three sets of results. For the two-dimensional sectional model the results for the sensitivities of the flutter speed with respect to the various parameters are obtained.

Kapania, Rakesh K.↗
NASA NTRS1992

Finite difference time domain modeling of spiral antennas

The objectives outlined in the original proposal for this project were to create a well-documented computer analysis model based on the finite-difference, time-domain (FDTD) method that would be capable of computing antenna impedance, far-zone radiation patterns, and radar cross-section (RCS). The ability to model a variety of penetrable materials in addition to conductors is also desired. The spiral antennas under study by this project meet these requirements since they are constructed of slots cut into conducting surfaces which are backed by dielectric materials.

Penney, Christopher W.↗
NASA NTRS1992

The eigenvalue spectrum of the Rayleigh equation for a plane shear layer

The eigenvalue spectrum of the Rayleigh equation is examined using three different solution techniques. In particular, a simple second-order finite difference scheme and two spectral methods, the Chebyshev tau and Chebyshev collocation methods, are used to discretize the equation. All of the approximation methods are shown to be capable of predicting the discrete spectrum as well as the continuous spectrum associated with the critical point singularity for the Rayleigh equation. The global eigenvalue methods considered here provide an efficient way of obtaining either an approximation to the complete eigenvalue spectrum or initial guesses for a local shooting procedure for the discrete part of the spectrum.

Liou, William W.-W.↗
NASA NTRS1993

Two pass method and radiation interchange processing when applied to thermal-structural analysis of large space truss structures

A method of efficient and automated thermal-structural processing of very large space structures is presented. The method interfaces the finite element and finite difference techniques. It also results in a pronounced reduction of the quantity of computations, computer resources and manpower required for the task, while assuring the desired accuracy of the results.

Warren, Andrew H.↗
NASA NTRS1993

Nonlinear evolution of resistive tearing mode instability with shear flow and viscosity

The effect of shear flow on the nonlinear evolution of the tearing mode is investigated via numerical solutions of the resistive MHD equations in slab geometry, using a finite-difference alternative-direction implicit method. It was found that, when the shear flow is small (V less than 0.3), the tearing mode saturates within one resistive time, whereas for larger flows the nonlinear saturation develops on longer time scales. The magnetic energy release decreases and the saturation time increases with increasing values of V for both small and large resistivity. Shear flow was found to decrease the saturated magnetic island width and to generate currents far from the tearing layer. Results suggest that equilibrium shear flow may improve the confinement of tokamak plasma.

Ofman, L.↗
NASA NTRS1994

Low-dissipation and -dispersion Runge-Kutta schemes for computational acoustics

In this paper, we investigate accurate and efficient time advancing methods for computational acoustics, where non-dissipative and non-dispersive properties are of critical importance. Our analysis pertains to the application of Runge-Kutta methods to high-order finite difference discretization. In many CFD applications multi-stage Runge-Kutta schemes have often been favored for their low storage requirements and relatively large stability limits. For computing acoustic waves, however, the stability consideration alone is not sufficient, since the Runge-Kutta schemes entail both dissipation and dispersion errors. The time step is now limited by the tolerable dissipation and dispersion errors in the computation. In the present paper, it is shown that if the traditional Runge-Kutta schemes are used for time advancing in acoustic problems, time steps greatly smaller than that allowed by the stability limit are necessary. Low-Dissipation and -Dispersion Runge-Kutta (LDDRE) schemes are proposed, based on an optimization that minimizes the dissipation and dispersion errors for wave propagation. Order optimizations of both single-step and two-step alternating schemes are considered. The proposed LDDRK schemes are remarkably more efficient than the classical Runge-Kutta schemes for acoustic computations. Moreover, low storage implementations of the optimized schemes are discussed. Special issues of implementing numerical boundary conditions in the LDDRK schemes are also addressed.

Hu, F. Q.↗
NASA NTRS1995

Application of low dissipation and dispersion Runge-Kutta schemes to benchmark problems in computational aeroacoustics

We investigate accurate and efficient time advancing methods for computational aeroacoustics, where non-dissipative and non-dispersive properties are of critical importance. Our analysis pertains to the application of Runge-Kutta methods to high-order finite difference discretization. In many CFD applications, multi-stage Runge-Kutta schemes have often been favored for their low storage requirements and relatively large stability limits. For computing acoustic waves, however, the stability consideration alone is not sufficient, since the Runge-Kutta schemes entail both dissipation and dispersion errors. The time step is now limited by the tolerable dissipation and dispersion errors in the computation. In the present paper, it is shown that if the traditional Runge-Kutta schemes are used for time advancing in acoustic problems, time steps greatly smaller than that allowed by the stability limit are necessary. Low Dissipation and Dispersion Runge-Kutta (LDDRK) schemes are proposed, based on an optimization that minimizes the dissipation and dispersion errors for wave propagation. Optimizations of both single-step and two-step alternating schemes are considered. The proposed LDDRK schemes are remarkably more efficient than the classical Runge-Kutta schemes for acoustic computations. Numerical results of each Category of the Benchmark Problems are presented. Moreover, low storage implementations of the optimized schemes are discussed. Special issues of implementing numerical boundary conditions in the LDDRK schemes are also addressed.

Hu, F. Q.↗
NASA NTRS1996

Nonstandard Finite Difference Schemes: Relations Between Time and Space Step-Sizes in Numerical Schemes for PDE's That Follow from Positivity Condition

A large class of physical phenomena can be modeled by evolution and wave type Partial Differential Equations (PDE). Few of these equations have known explicit exact solutions. Finite-difference techniques are a popular method for constructing discrete representations of these equations for the purpose of numerical integration. However, the solutions to the difference equations often contain so called numerical instabilities; these are solutions to the difference equations that do not correspond to any solution of the PDE's. For explicit schemes, the elimination of this behavior requires functional relations to exist between the time and space steps-sizes. We show that such functional relations can be obtained for certain PDE's by use of a positivity condition. The PDE's studied are the Burgers, Fisher, and linearized Euler equations.

Mickens, Ronald E.↗
NASA NTRS1998

A Virtual World of Visualization

In 1990, Sterling Software, Inc., developed the Flow Analysis Software Toolkit (FAST) for NASA Ames on contract. FAST is a workstation based modular analysis and visualization tool. It is used to visualize and animate grids and grid oriented data, typically generated by finite difference, finite element and other analytical methods. FAST is now available through COSMIC, NASA's software storehouse.

Source record↗
NASA NTRS2018

Does Orbital Angular Momentum Have Effect on Laser’s Scattering by Molecular Atmosphere?

Lasers with orbital angular momentum (OAM) have potential applications in communication technology, manipulation of particles, and remote sensing. Because of its unusual light-scattering properties, the OAM laser’s interaction with a molecular atmosphere must be studied to ensure that it is not lossy for communication or remote-sensing applications that involve its transmission through an atmospheric environment. In this study, the finite-difference time-domain (FDTD) method [21] is applied to calculate the light scattering of the purely azimuthal (the radial mode number is assumed to be zero) Laguerre-Gaussian (LG) beams with OAM by very small dielectric particles. Not like Lorentz-Mie solutions, the FDTD method can calculate for particles off the central axis of the LG beam. It is found that when the particles are very small, and the topological charge number of the OAM of a laser is not extremely large, the laser’s OAM has little effect on the scattering phase function. This suggests that Rayleigh theory can be applied directly to calculate the light scattering by atmospheric molecules. The transmission of a laser beam with OAM in a molecular atmosphere is not different from that of a regular Gaussian beam.

Molecular atmosphere↗
NASA NTRS1976

Advances in the application of fast semidirect computational methods in transonic flow

The paper uses finite-difference algorithms called 'fast direct elliptic solvers' within an iteration scheme for the rapid solution of the equations of inviscid transonic aerodynamics. The methods are called 'direct' (or semidirect) because the entire computational field is solved at once rather than in successive traverses over the field. These semidirect iterative methods have been limited here to the investigation of two-dimensional steady inviscid flow over airfoils in subsonic free stream.

Martin, E. D.↗
NASA NTRS1984

Structural sensitivity analysis: Methods, applications, and needs

Some innovative techniques applicable to sensitivity analysis of discretized structural systems are reviewed. These techniques include a finite-difference step-size selection algorithm, a method for derivatives of iterative solutions, a Green's function technique for derivatives of transient response, a simultaneous calculation of temperatures and their derivatives, derivatives with respect to shape, and derivatives of optimum designs with respect to problem parameters. Computerized implementations of sensitivity analysis and applications of sensitivity derivatives are also discussed. Finally, some of the critical needs in the structural sensitivity area are indicated along with Langley plans for dealing with some of these needs.

Adelman, H. M.↗
NASA NTRS1996

Algebraic Multigrid by Smoothed Aggregation for Second and Fourth Order Elliptic Problems

Multigrid methods are very efficient iterative solvers for system of algebraic equations arising from finite element and finite difference discretization of elliptic boundary value problems. The main principle of multigrid methods is to complement the local exchange of information in point-wise iterative methods by a global one utilizing several related systems, called coarse levels, with a smaller number of variables. The coarse levels are often obtained as a hierarchy of discretizations with different characteristic meshsizes, but this requires that the discretization is controlled by the iterative method. To solve linear systems produced by existing finite element software, one needs to create an artificial hierarchy of coarse problems. The principal issue is then to obtain computational complexity and approximation properties similar to those for nested meshes, using only information in the matrix of the system and as little extra information as possible. Such algebraic multigrid method that uses the system matrix only was developed by Ruge. The prolongations were based on the matrix of the system by partial solution from given values at selected coarse points. The coarse grid points were selected so that each point would be interpolated to via so-called strong connections. Our approach is based on smoothed aggregation introduced recently by Vanek. First the set of nodes is decomposed into small mutually disjoint subsets. A tentative piecewise constant interpolation (in the discrete sense) is then defined on those subsets as piecewise constant for second order problems, and piecewise linear for fourth order problems. The prolongation operator is then obtained by smoothing the output of the tentative prolongation and coarse level operators are defined variationally.

Vanek, Petr↗
NASA NTRS2015

Impedance Eduction in Sound Fields With Peripherally Varying Liners and Flow

A two-dimensional impedance eduction theory is extended to three-dimensional sound fields and peripherally varying duct liners. The approach is to first measure the acoustic pressure field at a series of flush-mounted wall microphones located around the periphery of the flow duct. The numerical solution for the acoustic pressure field at these microphones is also obtained by solving the three-dimensional convected Helmholtz equation using the finite element method. A quadratic objective function based on the difference between the measured and finite element solution is constructed and the unknown impedance function is obtained by minimizing this objective function. Impedance spectra educed for two uniform-structure liners (a wire-mesh and a conventional liner) and a hard-soft-hard peripherally varying liner (for which the soft segment is that of the conventional liner) are presented. Results are presented at three mean flow Mach numbers and fourteen sound source frequencies. The impedance spectra of the uniform-structure liners are also computed using a two-dimensional impedance eduction theory. The primary conclusions of the study are: 1) when measured data is used with the uniform-structure liners, the three-dimensional theory reproduces the same impedance spectra as the two-dimensional theory except for frequencies corresponding to very low or very high liner attenuation; and 2) good agreement between the educed impedance spectra of the uniform structure conventional liner and the soft segment of the peripherally varying liner is obtained.

Watson, W. R.↗