Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite difference scheme”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Advanced rotor computations with a corrected potential method

An unsteady Full-Potential Rotor code (FPR) has been enhanced with modifications directed at improving its drag prediction capability. The potential code has been rewritten with modifications to increase the code accuracy. Also, the shock generated entropy has been included to provide solutions comparable to the Euler equations. Two different weakly interacted boundary layer models have also been coupled to FPR in order to estimate skin-friction drag. One is a two-dimensional integral method and the other is a three-dimensional finite-difference scheme. The new flow solver is able to find accurate inviscid drags without recourse to numerical error tares. This permits the resolution of drag distributions resulting from rotor geometric variations. Good comparisons have been obtained between computed and measured torque for a rectangular and a highly swept model rotor.

Bridgeman, John O.↗
NASA NTRS1990

Numerical simulation of solidification in a horizontal cylindrical annulus charged with an aqueous salt solution

A computational study of solidification of a binary Na2CO3 solution in a horizontal cylindrical annulus is performed using a continuum formulation with a control-volume based, finite-difference scheme. The initial conditions were selected to facilitate the study of counter thermal and solutal convection, accompanied by extensive mushy region growth. Numerical results are compared with experimental data with mixed success. Qualitative agreement is obtained for the overall solidification process and associated physical phenomena. However, the plume thickness calculated for the solutally-driven convective upflow is substantially smaller than the observed value. Evolution of double-diffusive layers is predicted, but over a time scale much smaller than that observed experimentally. Good agreement is obtained between predicted and measured results for solid growth, but the mushy region thickness is significantly overpredicted.

Neilson, D. G.↗
NASA NTRS1990

Navier-Stokes computations on swept-tapered wings, including flexibility

A procedure to couple the Navier-Stokes solutions with modal structural equations of motion is presented for computing aeroelastic responses of flexible fighter wings. The Navier-Stokes flow equations are solved by a finite-difference scheme with dynamic grids. The coupled aeroelastic equations of motion are solved using the linear-acceleration method. The configuration-adaptive dynamic grids are time-accurately generated using the aeroelastically deformed shape of the wing. The coupled calculations are compared with experiments when available. Effects of flexibility and pitch rate are demonstrated for flows with vortices. Turbulent flow computations are also compared with laminar flow computations.

Guruswamy, Guru P.↗
NASA NTRS1989

Basic numerical methods

Some of the basic finite-difference schemes that can be used to solve the nonlinear equations that describe unsteady inviscid and viscous transonic flow are reviewed. Numerical schemes for solving the unsteady Euler and Navier-Stokes, boundary-layer, and nonlinear potential equations are described. Emphasis is given to the elementary ideas used in constructing various numerical procedures, not specific details of any one procedure.

Steger, Joseph L.↗
NASA NTRS1989

Two-dimensional numerical simulation of a Stirling engine heat exchanger

The first phase of an effort to develop multidimensional models of Stirling engine components is described. The ultimate goal is to model an entire engine working space. Parallel plate and tubular heat exchanger models are described, with emphasis on the central part of the channel (i.e., ignoring hydrodynamic and thermal end effects). The model assumes laminar, incompressible flow with constant thermophysical properties. In addition, a constant axial temperature gradient is imposed. The governing equations describing the model have been solved using the Crack-Nicloson finite-difference scheme. Model predictions are compared with analytical solutions for oscillating/reversing flow and heat transfer in order to check numerical accuracy. Excellent agreement is obtained for flow both in circular tubes and between parallel plates. The computational heat transfer results are in good agreement with the analytical heat transfer results for parallel plates.

Ibrahim, Mounir↗
NASA NTRS1990

Effect of convection on the microstructure of a lamellar eutectic growing with a stepped interface

A two dimensional model developed to study the influence of convective flow on the concentration field ahead of a growing lamellar eutectic is described when one phase projects out into the melt creating a stepped interface. The two dimensional convective flow field, which is periodic in the horizontal direction, was computed numerically using the software FLUENT. The velocity field generated due to the flow of melt over the steps was then incorporated into a finite difference scheme employed to solve the concentration field. The average interfacial composition was calculated and converted to lamellar spacing using the Jackson and Hunt minimum supercooling criterion. It was found that a stepped interface is more sensitive to convection than a planar interface.

Seth, Jayshree↗
NASA NTRS1991

Radiative interactions in chemically reacting supersonic internal flows

The two-dimensional, elliptic Navier-Stokes equations are used to investigate supersonic flows with finite-rate chemistry and radiation for hydrogen-air systems. The chemistry source terms in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The explicit, unsplit MacCormack finite-difference scheme is used to advance the governing equations in time, until convergence is achieved. The specific problem considered is the premixed flow in a channel with a ten-degree compression ramp. Three different chemistry models are used, accounting for increasing number of reactions and participating species. Two chemistry models assure nitrogen as inert, while the third model accounts for nitrogen reactions and NO(x) formation. The tangent slab approximation is used in the radiative flux formulation. A pseudo-gray model is used to represent the absorption-emission characteristics of the participating species. Results obtained for specific conditions indicate that the radiative interactions vary substantially, depending on reactions involving HO2 and NO species and that this can have a significant influence on the flowfield.

Tiwari, S. N.↗
NASA NTRS1991

Asynchronous and corrected-asynchronous numerical solutions of parabolic PDES on MIMD multiprocessors

A major problem in achieving significant speed-up on parallel machines is the overhead involved with synchronizing the concurrent process. Removing the synchronization constraint has the potential of speeding up the computation. The authors present asynchronous (AS) and corrected-asynchronous (CA) finite difference schemes for the multi-dimensional heat equation. Although the discussion concentrates on the Euler scheme for the solution of the heat equation, it has the potential for being extended to other schemes and other parabolic partial differential equations (PDEs). These schemes are analyzed and implemented on the shared memory multi-user Sequent Balance machine. Numerical results for one and two dimensional problems are presented. It is shown experimentally that the synchronization penalty can be about 50 percent of run time: in most cases, the asynchronous scheme runs twice as fast as the parallel synchronous scheme. In general, the efficiency of the parallel schemes increases with processor load, with the time level, and with the problem dimension. The efficiency of the AS may reach 90 percent and over, but it provides accurate results only for steady-state values. The CA, on the other hand, is less efficient, but provides more accurate results for intermediate (non steady-state) values.

Amitai, Dganit↗
NASA NTRS1990

A discrete model of a modified Burgers' partial differential equation

A new finite-difference scheme is constructed for a modified Burger's equation. Three special cases of the equation are considered, and the 'exact' difference schemes for the space- and time-independent forms of the equation are presented, along with the diffusion-free case of Burger's equation modeled by a difference equation. The desired difference scheme is then obtained by imposing on any difference model of the initial equation the requirement that, in the appropriate limits, its difference scheme must reduce the results of the obtained equations.

Mickens, R. E.↗
NASA NTRS1990

Prediction of blade-vortex interaction noise using airloads generated by a finite-difference technique

The present numerical finite-difference scheme for helicopter blade-load prediction during realistic, self-generated three-dimensional blade-vortex interactions (BVI) derives the velocity field through a nonlinear superposition of the rotor flow-field yielded by the full potential rotor flow solver RFS2 for BVI, on the one hand, over the rotational vortex flow field computed with the Biot-Savart law. Despite the accurate prediction of the acoustic waveforms, peak amplitudes are found to have been persistently underpredicted. The inclusion of BVI noise source in the acoustic analysis significantly improved the perceived noise level-corrected tone prediction.

Tadghighi, Hormoz↗
NASA NTRS1991

Radiative interactions in a hydrogen-fueled supersonic combustor

The two-dimensional, elliptic Navier-Stokes equations are used to investigate supersonic flows with finite-rate chemistry and radiation, for hydrogen-air systems. The chemistry source term in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The explicit, unsplit MacCormack finite-difference scheme is used to advance the governing equations in time, until convergence is achieved. The specific problem considered is the premixed flow in a channel with a 10-deg compression ramp. Three different chemistry models are used, accounting for increasing number of reactions and participating species. Two chemistry models assume nitrogen as inert, while the third model accounts for nitrogen reactions and NO(x) formation. The tangent slab approximation is used in the radiative flux formulation. A pseudogray model is used to represent the absorption-emission characteristics of the participating species. Results otained for specific conditions indicate that the radiative interactions vary substantially, depending on reactions involving HO2 and NO species, and that this can have a significant influence on the flowfield.

Chandrasekhar, R.↗
NASA NTRS1991

Navier-Stokes computation of wing/rotor interaction for a tilt rotor in hover

A method has been developed to analyze the wing/rotor interaction of tilt rotor aircraft in hover. The unsteady, thin-layer compressible Navier-Stokes equations are solved using an implicit, finite difference scheme that employs LU-ADI factorization. The rotor is modeled as an actuator disk which imparts a radial and azimuthal distribution of pressure rise and swirl to the flowfield. The 'chimera' approach of grid point blanking is used to update the rotor boundary conditions. Results are presented for both a rotor alone and for wing/rotor interaction where the thrust coefficient is 0.0164 and wing flap deflection is 67 degrees. Many of the complex flow features are captured including the fountain effect, leading and trailing edge separation, and the unsteady wake beneath the wing. Wing surface pressures compare fairly well with experimental data although the time-averaged download is about twenty percent higher than the measured value. This discrepancy is due to a combination of factors that are discussed.

Fejtek, Ian↗
NASA NTRS1990

Asynchronous numerical solution of PDEs on parallel computers

Asynchronous and partially-synchronous finite difference schemes are investigated for the one-dimensional steady-state and nonsteady-state heat equation, respectively. It is shown experimentally that the asynchronous scheme runs twice as fast as the parallel synchronous scheme. For relatively small values of time levels, however, the asynchronous scheme produces the least accurate results, while the partially-asynchronous scheme is almost as accurate as its synchronous counterpart.

Amitai, Dganit↗
NASA NTRS1991

A Fourier-Chebyshev pseudospectral method for solving steady 3-D Navier-Stokes and heat equations in cylindrical cavities

A Fourier-Chebyshev pseudospectral method for solving steady 3D Navier-Stokes equations in cylindrical cavities is presented and discussed. The general method is pseudo-unsteady and uses a semi-implicit finite difference scheme for the time discretization. The generalized ADI (Alternating Direction Implicit) procedure is then applied to reduce the problem to a successive solution of one-dimensional problems. The spatial approximation uses a Fourier-Galerkin approximation in the periodic direction and a Chebyshev-collocation approximation in the other directions. Difficulties related to the pressure are surmounted by using the artifical compressibility method. A suitable variable change was chosen to avoid the problem of singularity at the axis generated by cylindrical coordinates. The method is first tested on an advection-diffusion equation and then on the Navier-Stokes equations. Finally, the method is illustrated by a convection problem of a differentially heated fluid.

Pulicani, J. P.↗
NASA NTRS1991

Leading-edge receptivity to a vortical freestream disturbance: A numerical analysis

The receptivity to freestream vorticity of the boundary layer over a flat plate with an elliptic leading edge is investigated numerically. The flow is simulated by solving the incompressible Navier-Stokes system in general curvilinear coordinates with the vorticity and stream function as dependent variables. A finite-difference scheme which is second-order accurate in both space and time is used. As a first step, the steady basic-state solution is computed. Then a small amplitude vortical disturbance is introduced at the upstream boundary and the governing equations are solved time-accurately to evaluate the spatial and temporal growth of the perturbations leading to instability waves (Tollmien-Schlichting waves) inside the boundary layer. Preliminary results for a symmetric, 2-D disturbance reveal the presence of Tollmien-Schlichting waves aft of the flat-plate/ellipse juncture.

Buter, Thomas A.↗
NASA NTRS1992

Sunspot: A program to model the behavior of hypervelocity impact damaged multilayer insulation in the Sunspot thermal vacuum chamber of Marshall Space Flight Center

The development of a computer program to predict the degradation of the insulating capabilities of the multilayer insulation (MLI) blanket of Space Station Freedom due to a hypervelocity impact with a space debris particle is described. A finite difference scheme is used for the calculations. The computer program was written in Microsoft BASIC. Also described is a test program that was undertaken to validate the numerical model. Twelve MLI specimens were impacted at hypervelocities with simulated debris particles using a light gas gun at Marshall Space Flight Center. The impact-damaged MLI specimens were then tested for insulating capability in the space environment of the Sunspot thermal vacuum chamber at MSFC. Two undamaged MLI specimens were also tested for comparison with the test results of the damaged specimens. The numerical model was found to adequately predict behavior of the MLI specimens in the Sunspot chamber. A parameter, called diameter ratio, was developed to relate the nominal MLI impact damage to the apparent (for thermal analysis purposes) impact damage based on the hypervelocity impact conditions of a specimen.

Rule, W. K.↗
NASA NTRS1991

The Development of High Order Numerical Techniques for Reentry Simulation of Hypersonic Spacecraft

The primary difficulty encountered when simulating hypersonic flow is that the flow normally includes strong nonlinear discontinuities. These discontinuities fall into three broad classes: shocks, slip-lines, and rarefaction waves. Moreover, in the hypersonic flow regime, the chemistry of hot gases plays a vital role and can not be neglected. These facts combine to make the numerical treatment of spacecraft reentry a most challenging problem. In this work, we develop a class of finite difference schemes that accurately resolve discontinuous solutions to spacecraft reentry flow and are simple to incorporate into existing spacecraft reentry codes.

Sanders, Richard↗
NASA NTRS1991

The study of PDF turbulence models in combustion

In combustion computations, it is known that the predictions of chemical reaction rates are poor if conventional turbulence models are used. The probability density function (pdf) method seems to be the only alternative that uses local instantaneous values of the temperature, density, etc., in predicting chemical reaction rates, and thus is the only viable approach for more accurate turbulent combustion calculations. The fact that the pdf equation has a very large dimensionality renders finite difference schemes extremely demanding on computer memories and thus impractical. A logical alternative is the Monte Carlo scheme. Since CFD has a certain maturity as well as acceptance, it seems that the use of a combined CFD and Monte Carlo scheme is more beneficial. Therefore, a scheme is chosen that uses a conventional CFD flow solver in calculating the flow field properties such as velocity, pressure, etc., while the chemical reaction part is solved using a Monte Carlo scheme. The discharge of a heated turbulent plane jet into quiescent air was studied. Experimental data for this problem shows that when the temperature difference between the jet and the surrounding air is small, buoyancy effect can be neglected and the temperature can be treated as a passive scalar. The fact that jet flows have a self-similar solution lends convenience in the modeling study. Futhermore, the existence of experimental data for turbulent shear stress and temperature variance make the case ideal for the testing of pdf models wherein these values can be directly evaluated.

Hsu, Andrew T.↗