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.

SourceResearch connectionFrom yearThrough year

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

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Current research in cavitating fluid films

A review of the current research of cavitation in fluid films is presented. Phenomena and experimental observations include gaseous cavitation, vapor cavitation, and gas entrainment. Cavitation in flooded, starved, and dynamically loaded journal bearings, as well as squeeze films are reviewed. Observations of cavitation damage in bearings and the possibility of cavitation between parallel plates with microasperities were discussed. The transcavity fluid transport process, meniscus motion and geometry or form of the film during rupture, and reformation were summarized. Performance effects were related to heat transfer models in the cavitated region and hysteresis influence on rotor dynamics coefficients. A number of cavitation algorithms was presented together with solution procedures using the finite difference and finite element methods. Although Newtonian fluids were assumed in most of the discussions, the effect of non-Newtonian fluids on cavitation was also discussed.

Brewe, D. E.↗

NASA NTRS1989

Time-dependent viscous incompressible Navier-Stokes equations - The finite difference Galerkin formulation and streamfunction algorithms

Numerical techniques are developed to solve the Navier-Stokes equations for unsteady incompressible flow. The extension of the finite-difference Galerkin (FDG) method of Stephens et al. (1984) to the continuous-time case in two or three space dimensions is explained, and the numerical implementation of the method is discussed with particular attention to the staggered-MAC-grid primitive-variable discretization, the application of discrete mass balance to avoid problems inherent in FDG schemes, the direct interpretation of the FDG expansion variables as a discrete streamfunction, and a mass-balance approach to two-dimensional problems with throughflow or obstacles. Numerical results are presented graphically for the evolution of asymptotic steady flow in a driven cavity at Reynolds number 400, 1000, or 3200; good agreement with published experimental data is demonstrated, with accurate predictions of secondary-vortex formation from wall bubble recirculations at Reynolds number 1000.

Goodrich, John W.↗

NASA NTRS1989

A numerical study of the interaction between unsteady free-stream disturbances and localized variations in surface geometry

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small free-stream disturbance and a small localized variation of the surface geometry has been carried out using both finite-difference and spectral methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T-S waves generated by the interaction between the free-stream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T-S waves.

Bodonyi, R. J.↗

NASA NTRS1991

Finite difference time domain calculation of transients in antennas with nonlinear loads

In this paper transient fields for antennas with more general geometries are calculated directly using Finite Difference Time Domain methods. In each FDTD cell which contains a nonlinear load, a nonlinear equation is solved at each time step. As a test case the transient current in a long dipole antenna with a nonlinear load excited by a pulsed plane wave is computed using this approach. The results agree well with both calculated and measured results previously published. The approach given here extends the applicability of the FDTD method to problems involving scattering from targets including nonlinear loads and materials, and to coupling between antennas containing nonlinear loads. It may also be extended to propagation through nonlinear materials.

Luebbers, Raymond J.↗

NASA NTRS1993

Pre- and postprocessing techniques for determining goodness of computational meshes

Research in error estimation, mesh conditioning, and solution enhancement for finite element, finite difference, and finite volume methods has been incorporated into AUDITOR, a modern, user-friendly code, which operates on 2D and 3D unstructured neutral files to improve the accuracy and reliability of computational results. Residual error estimation capabilities provide local and global estimates of solution error in the energy norm. Higher order results for derived quantities may be extracted from initial solutions. Within the X-MOTIF graphical user interface, extensive visualization capabilities support critical evaluation of results in linear elasticity, steady state heat transfer, and both compressible and incompressible fluid dynamics.

Oden, J. Tinsley↗

NASA NTRS1993

Numerical simulation of laminar-turbulent transition in a spatially-developing flat plate wake

Laminar-turbulent transition of an incompressible flat-plate wake is investigated by direct numerical integration of the Navier-Stokes equations. For the numerical integration, a combination of finite-difference and spectral methods along with an ADI/Crank-Nicolson/Adams-Bashforth time integration scheme is employed. Subject to 2D forcing, the wake exhibited a rapidly-growing fundamental disturbance that quickly saturated. This saturation was due partly to the stabilizing effect of the mean flow distortion. Downstream of the saturation point, disturbance energy was concentrated in the fundamental disturbance, the second harmonic, and the mean flow distortion component. At large amplitude levels, a Karman vortex street formed. Variations in the 2D forcing level did not alter the qualitative behavior of the disturbances. Simulations of 3D breakdown indicates that the presence of large-amplitude, 2D disturbances tends to initially suppress small-amplitude 3D disturbance growth. Following this initial suppression, a resumption of 3D growth is observed that may have been due to a secondary instability mechanism. For high levels of 3D disturbance energy, lambda-shaped vortical structures formed between adjacent Karman vortices.

Dratler, D. I.↗

NASA NTRS1994

Theoretical and Experimental Study of Microstrip-to-Slot Line Uniplanar Transition

Recent advances in MMCI technology make it possible to construct transitions from CPW-to-microstrip with via hole, microstrip-to-slot line and microshield line-to-CPW all of which have potential applications in the feed network of antennas. In this study we investigate the characteristics of the microstrip-to-slot line uniplanar transition using the finite element methods (FEM) and finite difference time domain (FDTD) techniques, and compared the theoretical results with the measurements. In both cases, the results agree with the measurements within a few percent.

Yook, Jong-Gwan↗

NASA NTRS1994

Steady potential solver for unsteady aerodynamic analyses

Development of a steady flow solver for use with LINFLO was the objective of this report. The solver must be compatible with LINFLO, be composed of composite mesh, and have transonic capability. The approaches used were: (1) steady flow potential equations written in nonconservative form; (2) Newton's Method; (3) implicit, least-squares, interpolation method to obtain finite difference equations; and (4) matrix inversion routines from LINFLO. This report was given during the NASA LeRC Workshop on Forced Response in Turbomachinery in August of 1993.

Hoyniak, Dan↗

NASA NTRS1993

Direct time integration of Maxwell's equations in two-dimensional dielectric waveguides for propagation and scattering of femtosecond electromagnetic solitons

Solution of 2D vector nonlinear Maxwell's equations is described for material media with linear and nonlinear instantaneous and Lorentz dispersion effects in the electric polarization. Finite difference time domain method is used in the analysis. The method discussed here achieves robustness by enforcing vector-field boundary conditions at all interfaces of dissimilar media.

VECTOR FIELDS↗

NASA NTRS1995

Large-eddy simulation of flow past a circular cylinder

Some of the most challenging applications of large-eddy simulation are those in complex geometries where spectral methods are of limited use. For such applications more conventional methods such as finite difference or finite element have to be used. However, it has become clear in recent years that dissipative numerical schemes which are routinely used in viscous flow simulations are not good candidates for use in LES of turbulent flows. Except in cases where the flow is extremely well resolved, it has been found that upwind schemes tend to damp out a significant portion of the small scales that can be resolved on the grid. Furthermore, it has been found that even specially designed higher-order upwind schemes that have been used successfully in the direct numerical simulation of turbulent flows produce too much dissipation when used in conjunction with large-eddy simulation. The objective of the current study is to perform a LES of incompressible flow past a circular cylinder at a Reynolds number of 3900 using a solver which employs an energy-conservative second-order central difference scheme for spatial discretization and compare the results obtained with those of Beaudan & Moin (1994) and with the experiments in order to assess the performance of the central scheme for this relatively complex geometry.

Mittal, R.↗

NASA NTRS1996

Transient Heat Transfer in a Semitransparent Radiating Layer with Boundary Convection and Surface Reflections

Surface convection and refractive index are examined during transient radiative heating or cooling of a grey semitransparent layer with internal absorption, emission and conduction. Each side of the layer is exposed to hot or cold radiative surroundings, while each boundary is heated or cooled by convection. Emission within the layer and internal reflections depend on the layer refractive index. The reflected energy and heat conduction distribute energy across the layer and partially equalize the transient temperature distributions. Solutions are given to demonstrate the effect of radiative heating for layers with various optical thicknesses, the behavior of the layer heated by radiation on one side and convectively cooled on the other, and a layer heated by convection while being cooled by radiation. The numerical method is an implicit finite difference procedure with non-uniform space and time increments. The basic method developed in earlier work is expanded to include external convection and incident radiation.

Siegel, Robert↗

NASA NTRS1997

Lightning Effects in the Payload Changeout Room

Analytical and empirical studies have been performed to provide better understanding of the electromagnetic environment inside the Payload Changeout Room and Orbiter payload bay resulting from lightning strikes to the launch pad lightning protection system. The analytical studies consisted of physical and mathematical modeling of the pad structure and the Payload Changeout Room. Empirical testing was performed using a lightning simulator to simulate controlled (8 kA) lightning strikes to the catenary wire lightning protection system. In addition to the analyses and testing listed above, an analysis of the configuration with the vehicle present was conducted, in lieu of testing, by the Finite Difference, Time Domain method.

Thomas, Garland L.↗

NASA NTRS2000

A Technique of Treating Negative Weights in WENO Schemes

High order accurate weighted essentially non-oscillatory (WENO) schemes have recently been developed for finite difference and finite volume methods both in structural and in unstructured meshes. A key idea in WENO scheme is a linear combination of lower order fluxes or reconstructions to obtain a high order approximation. The combination coefficients, also called linear weights, are determined by local geometry of the mesh and order of accuracy and may become negative. WENO procedures cannot be applied directly to obtain a stable scheme if negative linear weights are present. Previous strategy for handling this difficulty is by either regrouping of stencils or reducing the order of accuracy to get rid of the negative linear weights. In this paper we present a simple and effective technique for handling negative linear weights without a need to get rid of them.

Shi, Jing↗

NASA NTRS2008

Numerical Studies of Scattering Properties of Leaves and Leaf Moisture Influences on the Scattering at Microwave Wavelengths

This study uses 3-dimensional finite difference time domain method to accurately calculate single-scattering properties of randomly orientated leaves and evaluate the influences of vegetation water content (VWC) on these properties at 19 and 37 GHz frequencies. The studied leaves are assumed to be thin elliptic disks with two different sizes and have various VWC values. Although the leaf moisture produces considerable absorption during scattering processes, the effective efficiencies of extinction and scattering of leaves still near-linearly increase with VWC. Calculated asymmetry factors and phase functions indicate that there are significant amounts of scattering at large scattering angles in microwave wavelengths, which provides good opportunities for off-nadir microwave remote sensing of forests. This study lays a basic foundation in future quantifications of the relations between satellite measurements and physical properties of vegetation canopies.

Lin, Bing↗

NASA NTRS1994

The Application of the FDTD Method to Millimeter-Wave Filter Circuits Including the Design and Analysis of a Compact Coplanar

The finite difference time domain (FDTD) method is applied to the analysis of microwave, millimeter-wave and submillimeter-wave filter circuits. In each case, the validity of this method is confirmed by comparison with measured data. In addition, the FDTD calculations are used to design a new ultra-thin coplanar-strip filter for feeding a THz planar-antenna mixer.

Millimeter-wave filter Coplanar-Strip filter↗

NASA NTRS1985

Higher order methods for convection-diffusion problems

This paper applies C1 cubic Hermite polynomials embedded in an orthogonal collocation scheme to the spatial discretization of the unsteady nonlinear Burgers equation as a model of the equations of fluid mechanics. The temporal discretization is carried out by means of either a noniterative finite difference or an iterative finite difference procedure. Results of this method are compared with those of a second-order finite difference scheme and a splined-cubic Taylor's series scheme. Stability limits are derived and the matrix structure of the several schemes are compared.

Murphy, J. D.↗

NASA NTRS1987

Progress in the prediction of unsteady heat transfer on turbines blades

Progress toward developing a general method for predicting unsteady heat transfer on turbine blades subject to blade-passing frequencies and Reynolds numbers relevant to the Space Shuttle Main Engine (SSME) is discussed. The method employs an invisid/viscous interactive procedure which has been tested extensively for steady subsonic and transonic external airfoil problems. One such example is shown. The agreement with experimental data and with Navier-Stokes calculations yields confidence in the method. The technique is extended to account for wake generated unsteadiness. The flow reversals around the stagnation point caused by the nonuniform onset velocity are accounted for by using the Characteristic Box scheme developed by Cebeci and Stewartson. The coupling between the inviscid and viscous methods is achieved by using a special procedure, which, with a novel inverse finite-difference boundary-layer method, allows the calculations to be performed for a wide range of flow conditions, including separation. Preliminary results are presented for the stagnation region of turbine blades for both laminar and turbulent flows. A laminar model problem corresponding to a flow on a circular cylinder which experiences the periodic passing of wakes from turbine blades is presented to demonstrate the ability of the method to calculate flow reversals around the stagnation region.

Cebeci, T.↗

NASA NTRS1978

Calculation of supersonic viscous flow over delta wings with sharp subsonic leading edges

Two complementary procedures were developed to calculate the viscous supersonic flow over conical shapes at large angles of attack, with application to cones and delta wings. In the first approach the flow is assumed to be conical and the governing equations are solved at a given Reynolds number with a time-marching explicit finite-difference algorithm. In the second method the parabolized Navier-Stokes equations are solved with a space-marching implicit noniterative finite-difference algorithm. This latter approach is not restricted to conical shapes and provides a large improvement in computational efficiency over published methods. Results from the two procedures agree very well with each other and with available experimental data.

Vigneron, Y. C.↗