Search NASA⌕ Search

SEARCH · Search NASA

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

An assessment and application of turbulence models for hypersonic flows

The current approach to the Accurate Computation of Complex high-speed flows is to solve the Reynolds averaged Navier-Stokes equations using finite difference methods. An integral part of this approach consists of development and applications of mathematical turbulence models which are necessary in predicting the aerothermodynamic loads on the vehicle and the performance of the propulsion plant. Computations of several high speed turbulent flows using various turbulence models are described and the models are evaluated by comparing computations with the results of experimental measurements. The cases investigated include flows over insulated and cooled flat plates with Mach numbers ranging from 2 to 8 and wall temperature ratios ranging from 0.2 to 1.0. The turbulence models investigated include zero-equation, two-equation, and Reynolds-stress transport models.

Coakley, T. J.↗
NASA NTRS1991

Effect of leading-edge geometry on boundary-layer receptivity to freestream sound

The receptivity to freestream sound of the laminar boundary layer over a semi-infinite flat plate with an elliptic leading edge is simulated numerically. The incompressible flow past the flat plate is computed by solving the full Navier-Stokes equations in general curvilinear coordinates. A finite-difference method which is second-order accurate in space and time is used. Spatial and temporal developments of the Tollmien-Schlichting wave in the boundary layer, due to small-amplitude time-harmonic oscillations of the freestream velocity that closely simulate a sound wave travelling parallel to the plate, are observed. The effect of leading-edge curvature is studied by varying the aspect ratio of the ellipse. The boundary layer over the flat plate with a sharper leading edge is found to be less receptive. The relative contribution of the discontinuity in curvature at the ellipse-flat-plate juncture to receptivity is investigated by smoothing the juncture with a polynomial. Continuous curvature leads to less receptivity. A new geometry of the leading edge, a modified super ellipse, which provides continuous curvature at the juncture with the flat plate, is used to study the effect of continuous curvature and inherent pressure gradient on receptivity.

Lin, Nay↗
NASA NTRS1992

Development of an integrated BEM approach for hot fluid structure interaction: BEST-FSI: Boundary Element Solution Technique for Fluid Structure Interaction

As part of the continuing effort at NASA LeRC to improve both the durability and reliability of hot section Earth-to-orbit engine components, significant enhancements must be made in existing finite element and finite difference methods, and advanced techniques, such as the boundary element method (BEM), must be explored. The BEM was chosen as the basic analysis tool because the critical variables (temperature, flux, displacement, and traction) can be very precisely determined with a boundary-based discretization scheme. Additionally, model preparation is considerably simplified compared to the more familiar domain-based methods. Furthermore, the hyperbolic character of high speed flow is captured through the use of an analytical fundamental solution, eliminating the dependence of the solution on the discretization pattern. The price that must be paid in order to realize these advantages is that any BEM formulation requires a considerable amount of analytical work, which is typically absent in the other numerical methods. All of the research accomplishments of a multi-year program aimed toward the development of a boundary element formulation for the study of hot fluid-structure interaction in Earth-to-orbit engine hot section components are detailed. Most of the effort was directed toward the examination of fluid flow, since BEM's for fluids are at a much less developed state. However, significant strides were made, not only in the analysis of thermoviscous fluids, but also in the solution of the fluid-structure interaction problem.

Dargush, G. F.↗
NASA NTRS1991

Thermal calibration and analysis of a Bridgman flight furnace

A combination of computer modeling and experimental measurements is proposed which permits the calibration of the temperature profile for a spacecraft flight furnace. The Fluent code based on the finite difference method uses the furnace set-point temperatures as boundary conditions to model a Bridgman directional solidification furnace. Experimental measurements are employed as a comparison for the numerical analysis, and the radial gradients and wall temperatures are selected when the calculated values are sufficiently close to those from measurements. The calculations are found to agree with experimental data, and some results of the modeling are given which relate to specific furnace parameters. The CFD approach facilitates the transposition of furnace calibration to microgravitational operation in addition to providing precise initial radial temperature gradients and furnace-wall temperatures.

Knuteson, D. J.↗
NASA NTRS1991

Waveform effects of a metastable olivine tongue in subducting slabs

Velocity models of subducting slabs with a kinetically-depressed olivine to beta- and gamma-spinel transition are constructed, and the effect that such structures would have on teleseismic P waveforms are examined using a full-wave finite-difference method. These 2D calculations yielded waveforms at a range of distances in the downdip direction. The slab models included a wedge-shaped, low-velocity metastable olivine tongue (MOTO) to a depth of 670 km, as well as a plausible thermal anomaly; one model further included a 10-km-thick fast layer on the surface of the slab. The principal effect of MOTO is to produce grazing reflections at wide angles off the phase boundary, generating a secondary arrival 0 to 4 seconds after the initial arrival depending on the take-off angle. The amplitude and timing of this feature vary with the lateral location of the seismic source within the slab cross-section.

Vidale, John E.↗
NASA NTRS1992

Numerical models for the collapse and fragmentation of centrally condensed molecular cloud cores

The gravitational collapse and fragmentation of centrally condensed molecular cloud cores are investigated using a new hydrodynamical code. The numerical scheme is second-order accurate and uses explicit finite difference methods to advance the fluid variables on a 3D Cartesian grid. Two initial power-law density profiles, rho varies as r exp -1 and rho varies as r exp -2, are considered, as well as two initial density perturbations in the azimuthal coordinate theta, rho-i - rho(1 + a cos 2theta) where a = 0.1 and 0.5. Fragmentation is found to be possible in these centrally condensed cores if the initial conditions also include differential rotation. Models which collapse with initial uniform rotation do not produce fragments. If molecular cloud cores are indeed centrally condensed, as suggested by observations of star-forming regions and by studies of ambipolar diffusion, then differential rotation may be a mechanism for producing binary protostars during gravitational collapse.

Myhill, Elizabeth A.↗
NASA NTRS1992

Convection in superposed fluid and porous layers

Thermal convection due to heating from below in a porous layer underlying a fluid layer has been analyzed using the Navier-Stokes equations for the fluid layers and the extended Darcy equation (including Brinkman and Forchheimer terms) for the porous layer. The flow is assumed to be two-dimensional and periodic in the horizontal direction. The numerical scheme used is a combined Galerkin and finite-difference method, and appropriate boundary conditions are applied at the interface. Results have been obtained for depth ratios of 0, 0.1, 0.2, 0.5, and 1.0, where this ratio is defined as the ratio of the thickness of the fluid layer to that of the porous layer. For the depth ratio of 0.1, the convection is dominated by the porous layer, similar to the situation at onset, even though the Rayleigh number for the fluid layer is well into the supercritical regime.

Chen, Falin↗
NASA NTRS1992

Boundary-layer receptivity due to a wall suction and control of Tollmien-Schlichting waves

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 suction slot on an otherwise flat surface was carried out using finite difference 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 T-S waves generated by the interaction between the free-stream disturbance and the suction slot, 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 and the demonstration of the possible active control of the growth of T-S waves.

Bodonyi, R. J.↗
NASA NTRS1992

Effect of spatial resolution on apparent sensitivity to initial conditions of a decaying flow as it becomes turbulent

In order to check for spurious chaos and obtain superior solutions for decaying Navier-Stokes flows, an investigation is conducted of the effect of spatial resolution on numerical results. The fourth-order finite difference method results obtained with grids of 32-cubed and 64-cubed points, and those of a pseudospectral method for 128-cubed points, indicate that the sensitivity of initially neighboring solutions to small changes in initial conditions increases with improving spatial resolution.

Deissler, Robert G.↗
NASA NTRS1992

Finite-difference algorithms for the time-domain Maxwell's equations - A numerical approach to RCS analysis

The applications of two CFD-based finite-difference methods to computational electromagnetics are investigated. In the first method, the time-domain Maxwell's equations are solved using the explicit Lax-Wendroff scheme and in the second method, the second-order wave equations satisfying the Maxwell's equations are solved using the implicit Crank-Nicolson scheme. The governing equations are transformed to a generalized curvilinear coordinate system and solved on a body-conforming mesh using the scattered-field formulation. The induced surface current and the bistatic radar cross section are computed and the results are validated for several two-dimensional test cases involving perfectly-conducting scatterers submerged in transverse-magnetic plane waves.

Vinh, Hoang↗
NASA NTRS1991

Heat transfer in a discrete-velocity gas

Heat transfer of a discrete-velocity gas composed of identical hard sphere molecules between two parallel plates is studied numerically. A finite-difference method recently developed by the authors is used for solving the discrete Boltzmann equations. The temperature distributions between the plates of temperature ratios 2.0 and 4.0, and the heat transfer of temperature ratio 4.0 are presented and are compared with results of other methods.

Inamuro, Takaji↗
NASA NTRS1993

Numerical Simulation of a High Mach Number Jet Flow

The recent efforts to develop accurate numerical schemes for transition and turbulent flows are motivated, among other factors, by the need for accurate prediction of flow noise. The success of developing high speed civil transport plane (HSCT) is contingent upon our understanding and suppression of the jet exhaust noise. The radiated sound can be directly obtained by solving the full (time-dependent) compressible Navier-Stokes equations. However, this requires computational storage that is beyond currently available machines. This difficulty can be overcome by limiting the solution domain to the near field where the jet is nonlinear and then use acoustic analogy (e.g., Lighthill) to relate the far-field noise to the near-field sources. The later requires obtaining the time-dependent flow field. The other difficulty in aeroacoustics computations is that at high Reynolds numbers the turbulent flow has a large range of scales. Direct numerical simulations (DNS) cannot obtain all the scales of motion at high Reynolds number of technological interest. However, it is believed that the large scale structure is more efficient than the small-scale structure in radiating noise. Thus, one can model the small scales and calculate the acoustically active scales. The large scale structure in the noise-producing initial region of the jet can be viewed as a wavelike nature, the net radiated sound is the net cancellation after integration over space. As such, aeroacoustics computations are highly sensitive to errors in computing the sound sources. It is therefore essential to use a high-order numerical scheme to predict the flow field. The present paper presents the first step in a ongoing effort to predict jet noise. The emphasis here is in accurate prediction of the unsteady flow field. We solve the full time-dependent Navier-Stokes equations by a high order finite difference method. Time accurate spatial simulations of both plane and axisymmetric jet are presented. Jet Mach numbers of 1.5 and 2.1 are considered. Reynolds number in the simulations was about a million. Our numerical model is based on the 2-4 scheme by Gottlieb & Turkel. Bayliss et al. applied the 2-4 scheme in boundary layer computations. This scheme was also used by Ragab and Sheen to study the nonlinear development of supersonic instability waves in a mixing layer. In this study, we present two dimensional direct simulation results for both plane and axisymmetric jets. These results are compared with linear theory predictions. These computations were made for near nozzle exit region and velocity in spanwise/azimuthal direction was assumed to be zero.

Hayder, M. Ehtesham↗
NASA NTRS1992

Some Aspects of Bifurcation Structure of Laminar Flow in Curved Ducts

A bifurcation study is made of laminar flow in curved ducts. The problem is formulated in a curvilinear coordinate system, and the governing equations, after orthogonal mapping is applied, are solved numerically by an iterative finite-difference method. Many computer runs were made with various duct cross-sections ranging from a circle to a square, to learn the transition of bifurcation structure with this change in cross-section and to reconcile the differences between them. In addition, a simpler technique is proposed to generate symmetric four-cell solutions in a circular pipe and a means is put forward to stabilize four-vortex structures in a complete cross-section.

Kao, Hsiao C.↗
NASA NTRS1992

Finite element-finite difference thermal/structural analysis of large space truss structures

A technique of automated and efficient thermal-structural processing of truss structures that interfaces the finite element and finite difference method was developed. The thermal-structural analysis tasks include development of the thermal and structural math models, thermal analysis, development of an interface and data transfer between the models, and finally an evaluation of the thermal stresses and displacements in the structure. Consequently, the objective of the developed technique was to minimize the model development time, in order to assure an automatic transfer of data between the thermal and structural models as well as to minimize the computer resources needed for the analysis itself. The method and techniques described are illustrated on the thermal/structural analysis of the Space Station Freedom main truss.

Warren, Andrew H.↗
NASA NTRS1992

The Burnett shock structures in low density hypersonic flows

One-dimensional shock structures in low density hypersonic flows have been studied numerically by solving the Burnett and Navier-Stokes equations. A new lower-upper scheme has been developed to obtain the steady and stationary normal shock solutions by the time-marching technique. The new scheme directly employs the implicit finite difference method without computing the eigenvalues, and it is first-order accurate in time and second-order accurate in space. The computed shock structures have been compared with available experimental data, and they are in good agreements. Numerical results show that the Burnett calculations has less numerical diffusion than the Navier-Stokes calculations does. As to the shock thickness, the Burnett solutions are apparently superior than the Navier-Stokes solutions.

Liaw, Goang-Shin↗
NASA NTRS1993

Efficient development and processing of thermal math models of very large space truss structures

As the spacecraft moves along the orbit, the truss members are subjected to direct and reflected solar, albedo and planetary infra-red (IR) heating rates, as well as IR heating and shadowing from other spacecraft components. This is a transient process with continuously changing heating loads and the shadowing effects. The resulting nonuniform temperature distribution may cause nonuniform thermal expansion, deflection and stress in the truss elements, truss warping and thermal distortions. There are three challenges in the thermal-structural analysis of the large truss structures. The first is the development of the thermal and structural math models, the second - model processing, and the third - the data transfer between the models. All three tasks require considerable time and computer resources to be done because of a very large number of components involved. To address these challenges a series of techniques of automated thermal math modeling and efficient processing of very large space truss structures were developed. In the process the finite element and finite difference methods are interfaced. A very substantial reduction of the quantity of computations was achieved while assuring a desired accuracy of the results. The techniques are illustrated on the thermal analysis of a segment of the Space Station main truss.

Warren, Andrew H.↗
NASA NTRS1993

Sensitivity analysis of flutter response of a typical section and a wing in transonic flow

A sensitivity analysis of flutter response of a two-degree of freedom airfoil with plunging and pitching degrees of freedom in transonic flow was performed using a state-space representation of the unsteady aerodynamic behavior. The structural equations of motion of the airfoil with bending and torsional degrees of freedom are coupled to the unsteady airloads, and the aeroelastic system so modeled is solved as an eigenvalue problem to determine the stability. The results of eigenanalysis showed good agreement with flutter calculations performed using a time-integration of the aeroelastic equations. The sensitivities of the flutter speed with respect to the mass and stiffness parameters wee computed by both the analytical and finite difference methods, showing excellent agreement.

Kapania, Rakesh K.↗
NASA NTRS1993

When clusters collide - A numerical Hydro/N-body simulation of merging galaxy clusters

A 3D numerical simulation of two merging clusters of galaxies, using a hybrid Hydro/N-body code, is presented. The hydrodynamics of the code is solved by an Eulerian finite difference method. Initial results disclose that the X-ray emission of the dominant cluster becomes elongated and broadened; heating occurs at the core of the dominant cluster as a result of multiple shocks, and high velocity gas motions within the intracluster medium. It is predicted that clusters which have undergone recent mergers and do not have cooling flows will have high peculiar gas velocities and that the shocks and turbulence generated during the merger may power cluster-wide radio halos. Prolonged high-velocity gas motions through the dominant cluster core possibly play a major role in the formation and shaping of wide-angle tailed radio sources associated with central dominant galaxies. The N-body component of the simulation reveals the subcluster to be dispersed as it passes through the dominant cluster.

Roettiger, Kurt↗