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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

On the use and the performance of software reliability growth models

We address the problem of predicting future failures for a piece of software. The number of failures occurring during a finite future time interval is predicted from the number failures observed during an initial period of usage by using software reliability growth models. Two different methods for using the models are considered: straightforward use of individual models, and dynamic selection among models based on goodness-of-fit and quality-of-prediction criteria. Performance is judged by the relative error of the predicted number of failures over future finite time intervals relative to the number of failures eventually observed during the intervals. Six of the former models and eight of the latter are evaluated, based on their performance on twenty data sets. Many open questions remain regarding the use and the performance of software reliability growth models.

Keiller, Peter A.↗
NASA NTRS1994

Eulerian-Lagrangian Simulations of Transonic Flutter Instabilities

This paper presents an overview of recent applications of Eulerian-Lagrangian computational schemes in simulating transonic flutter instabilities. This approach, the fluid-structure system is treated as a single continuum dynamics problem, by switching from an Eulerian to a Lagrangian formulation at the fluid-structure boundary. This computational approach effectively eliminates the phase integration errors associated with previous methods, where the fluid and structure are integrated sequentially using different schemes. The formulation is based on Hamilton's Principle in mixed coordinates, and both finite volume and finite element discretization schemes are considered. Results from numerical simulations of transonic flutter instabilities are presented for isolated wings, thin panels, and turbomachinery blades. The results suggest that the method is capable of reproducing the energy exchange between the fluid and the structure with significantly less error than existing methods. Localized flutter modes and panel flutter modes involving traveling waves can also be simulated effectively with no a priori knowledge of the type of instability involved.

Bendiksen, Oddvar O.↗
DOE OSTI2022

Neutronics and thermal-hydraulics simulation of generic pebble-bed fluoride-salt-cooled high-temperature reactor (gFHR)

The fluoride-salt-cooled high-temperature reactor (FHR) is one type of the advanced reactors and has been attracting great interest from the research institutes and commercial companies in the recent years. However, currently the technology is relatively immature. To facilitate the design and safety analysis of FHRs, Kairos Power published the generic FHR (gFHR) benchmark. In this paper, a hybrid method combining the stochastic code and deterministic code is developed to simulate the gFHR benchmark. Serpent 2 is employed as a few-group cross section generator and the cross sections are applied to the finite difference neutron diffusion code AGREE. A good consistency with the gFHR benchmark is achieved. The agreement between the Serpent 2 results and the AGREE results shows that the hybrid method is applicable to FHRs. The thermal-hydraulics are coupled to neutronics for the steady-state calculations and good agreement between AGREE and SAM is achieved. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1987

Aerodynamic analysis for aircraft with nacelles, pylons, and winglets at transonic speeds

A computational method has been developed to provide an analysis for complex realistic aircraft configurations at transonic speeds. Wing-fuselage configurations with various combinations of pods, pylons, nacelles, and winglets can be analyzed along with simpler shapes such as airfoils, isolated wings, and isolated bodies. The flexibility required for the treatment of such diverse geometries is obtained by using a multiple nested grid approach in the finite-difference relaxation scheme. Aircraft components (and their grid systems) can be added or removed as required. As a result, the computational method can be used in the same manner as a wind tunnel to study high-speed aerodynamic interference effects. The multiple grid approach also provides high boundary point density/cost ratio. High resolution pressure distributions can be obtained. Computed results are correlated with wind tunnel and flight data using four different transport configurations. Experimental/computational component interference effects are included for cases where data are available. The computer code used for these comparisons is described in the appendices.

Boppe, Charles W.↗
NASA NTRS1975

Calculation of three-dimensional compressible laminar and turbulent boundary layers. An implicit finite-difference procedure for solving the three-dimensional compressible laminar, transitional, and turbulent boundary-layer equations

An implicit finite-difference procedure is presented for solving the compressible three-dimensional boundary-layer equations. The method is second-order accurate, unconditionally stable (conditional stability for reverse cross flow), and efficient from the viewpoint of computer storage and processing time. The Reynolds stress terms are modeled by (1) a single-layer mixing length model and (2) a two-layer eddy viscosity model. These models, although simple in concept, accurately predicted the equilibrium turbulent flow for the conditions considered. Numerical results are compared with experimental wall and profile data for a cone at an angle of attack larger than the cone semiapex angle. These comparisons clearly indicate that the numerical procedure and turbulence models accurately predict the experimental data with as few as 21 nodal points in the plane normal to the wall boundary.

Harris, J. E.↗
NASA NTRS1981

Time and frequency domain analysis of sampled data controllers via mixed operation equations

Specification of the mathematical equations required to define the dynamic response of a linear continuous plant, subject to sampled data control, is complicated by the fact that the digital components of the control system cannot be modeled via linear ordinary differential equations. This complication can be overcome by introducing two new mathematical operations; namely, the operation of zero order hold and digial delay. It is shown that by direct utilization of these operations, a set of linear mixed operation equations can be written and used to define the dynamic response characteristics of the controlled system. It also is shown how these linear mixed operation equations lead, in an automatable manner, directly to a set of finite difference equations which are in a format compatible with follow on time and frequency domain analysis methods.

Frisch, H. P.↗
NASA NTRS1981

Linear/nonlinear behavior in unsteady transonic aerodynamics

The accurate calculation of the aerodynamic forces in unsteady transonic flow requires the solution of the nonlinear flow equations. The aeroelastician, on the other hand, seeks to treat his problems (flutter, for example) by means of linear equations whenever possible. He may do this, even when the underlying flow is nonlinear, if the perturbation forces are linear over some (perhaps small) range of unsteady amplitude of motion. This paper assesses the range of parameters over which linear behavior occurs. In particular calculations are made for an NACA 64A006 airfoil oscillating in pitch over a range of amplitudes, frequencies, and Mach numbers. The primary aerodynamic method used is the well known LTRAN2 code of Ballhaus and Goorjian that provides a finite-difference solution to the low frequency, small disturbance, two-dimensional potential flow equation. Comparisons are made with linear subsonic theory, local linearization, and, for steady flow, with the full potential equation code of Bauer, Garabedian, and Korn.

Dowell, E. H.↗
NASA NTRS1984

A solution procedure for three-dimensional incompressible Navier-Stokes equation and its application

An implicit, finite-difference procedure is presented for numerically solving viscous incompressible flows. For convenience of applying the present method to three-dimensional problems, primitive variables, namely the pressure and velocities, are used. One of the major difficulties in solving incompressible flows that use primitive variables is caused by the pressure field solution method which is used as a mapping procedure to obtain a divergence-free velocity field. The present method is designed to accelerate the pressure-field solution procedure. This is achieved by the method of pseudocompressibility in which the time derivative pressure term is introduced into the mass conservation equation. The pressure wave propagation and the spreading of the viscous effect is investigated using simple test problems. The present study clarifies physical and numerical characteristics of the pseudo-compressible approach in simulating incompressible flows. Computed results for external and internal flows are presented to verify the present procedure. The present algorithm has been shown to be very robust and accurate if the selection of the pseudo-compressibility parameter has been made according to the guidelines given.

Kwak, D.↗
NASA NTRS1985

Vibration and buckling of rotating, pretwisted, preconed beams including Cooriolis effects

The effects of pretwist, precone, setting angle and Coriolis forces on the vibration and buckling behavior of rotating, torsionally rigid, cantilevered beams were studied. The beam is considered to be clamped on the axis of rotation in one case, and off the axis of rotation in the other. Two methods are employed for the solution of the vibration problem: (1) one based upon a finite-difference approach using second order central differences for solution of the equations of motion, and (2) based upon the minimum of the total potential energy functional with a Ritz type of solution procedure making use of complex forms of shape functions for the dependent variables. The individual and collective effects of pretwist, precone, setting angle, thickness ratio and Coriolis forces on the natural frequencies and the buckling boundaries are presented. It is shown that the inclusion of Coriolis effects is necessary for blades of moderate to large thickness ratios while these effects are not so important for small thickness ratio blades. The possibility of buckling due to centrifugal softening terms for large values of precone and rotation is shown.

Subrahmanyam, K. B.↗
NASA NTRS1985

Vibration and buckling of rotating, pretwisted, preconed beams including Coriolis effects

The effects of pretwist, precone, setting angle and Coriolis forces on the vibration and buckling behavior of rotating, torsionally rigid, cantilevered beams were studied. The beam is considered to be clamped on the axis of rotation in one case, and off the axis of rotation in the other. Two methods are employed for the solution of the vibration problem: (1) one based upon a finite-difference approach using second order central differences for solution of the equations of motion, and (2) based upon the minimum of the total potential energy functional with a Ritz type of solution procedure making use of complex forms of shape functions for the dependent variables. The individual and collective effects of pretwist, precone, setting angle, thickness ratio and Coriolis forces on the natural frequencies and the buckling boundaries are presented. It is shown that the inclusion of Coriolis effects is necessary for blades of moderate to large thickness ratios while these effects are not so important for small thickness ratio blades. The possibility of buckling due to centrifugal softening terms for large values of precone and rotation is shown.

Subrahmanyam, K. B.↗
NASA NTRS1987

A simplified computer solution for the flexibility matrix of contacting teeth for spiral bevel gears

A computer code, FLEXM, was developed to calculate the flexibility matrices of contacting teeth for spiral bevel gears using a simplified analysis based on the elementary beam theory for the deformation of gear and shaft. The simplified theory requires a computer time at least one order of magnitude less than that needed for the complete finite element method analysis reported earlier by H. Chao, and it is much easier to apply for different gear and shaft geometries. Results were obtained for a set of spiral bevel gears. The teeth deflections due to torsion, bending moment, shearing strain and axial force were found to be in the order 10(-5), 10(-6), 10(-7), and 10(-8) respectively. Thus, the torsional deformation was the most predominant factor. In the analysis of dynamic load, response frequencies were found to be larger when the mass or moment of inertia was smaller or the stiffness was larger. The change in damping coefficient had little influence on the resonance frequency, but has a marked influence on the dynamic load at the resonant frequencies.

Hsu, C. Y.↗
NASA NTRS1987

Numerical study of laminar and turbulent flows inside a turnaround duct with and without guide vanes

The purpose of this study is to examine in detail incompressible laminar and turbulent flows inside a turnaround duct with and without guide vanes and to investigate the effects of vanes on the flow characteristics. To perform this study, an implicit finite difference code cast in general curvilinear coordinates is further developed. The code is based on the method of pseudo-compressibility and utilize ADI or implicit approximate factorization algorithm to achieve computational efficiency. Method of segmental sweeping is developed to overcome the multiple-zone problem due to imposition of guide vanes. In the present study, several test cases have been computed. These test cases include laminar and turbulent flows inside a turnaround duct without and with two or three guide vanes. The study reveals that: (1) there exists large recirculation zones inside the duct if no vanes are present; (2) properly shaped and positional guide vanes are effective in eliminating flow separation; and (3) laminar and turbulent flows have similar flow features. But turbulent flow has less total pressure drop.

Lin, S.-J.↗
NASA NTRS1987

Numerical simulation of flow path in the oxidizer side hot gas manifold of the Space Shuttle main engine

The purpose of this study is to examine in detail incompressible laminar and turbulent flows inside the oxidizer side Hot Gas Manifold of the Space Shuttle Main Engine. To perform this study, an implicit finite difference code cast in general curvilinear coordinates is further developed. The code is based on the method of pseudo-compressibility and utilize ADI or implicit approximate factorization algorithm to achieve computational efficiency. A multiple-zone method is developed to overcome the complexity of the geometry. In the present study, the laminar and turbulent flows in the oxidizer side Hot Gas Manifold have been computed. The study reveals that: (1) there exists large recirculation zones inside the bowl if no vanes are present; (2) strong secondary flows are observed in the transfer tube; and (3) properly shaped and positioned guide vanes are effective in eliminating flow separation.

Lin, S. J.↗
NASA NTRS1987

Aerodynamic interaction between vortical wakes and lifting two-dimensional bodies

Unsteady rotor wake interactions with the empenage, tail boom, and other aerodynamic surfaces of a helicopter have a significant influence on its aerodynamic performance, the ride quality, and amount of vibration. A numerical method for computing the aerodynamic interaction between an interacting vortex wake and the viscous flow about arbitrary two-dimensional bodies has been developed to address this helicopter problem. The method solves for the flow field velocities on a body-fitted computational mesh using finite-difference techniques. The interaction of a rotor wake with the flow about a 4:1 elliptic cylinder at 45-deg incidence was calculated for a Reynolds number of 3000.

Stremel, Paul M.↗
NASA NTRS1989

Numerical computation of shock wave-turbulent boundary layer interaction in transonic flow over an axisymmetric curved hill

A control-volume based finite difference computation of a turbulent transonic flow over an axisymmetric curved hill is presented. The numerical method is based on the SIMPLE algorithm, and hence the conservation of mass equation is replaced by a pressure correction equation for compressible flows. The turbulence is described by a k-epsilon turbulence model supplemented by a near-wall turbulence model. In the method, the dissipation rate in the region very close to the wall is obtained from an algebraic equation and that for the rest of the flow domain is obtained by solving a partial differential equation for the dissipation rate. The other flow equations are integrated up to the wall. It is shown that the present turbulence model yields the correct location of the compression shock. The other computational results are also in good agreement with experimental data.

Kim, S.-W.↗
NASA NTRS1989

Aerodynamic interaction between vortical wakes and lifting two-dimensional bodies

Unsteady rotor wake interactions with the empennage, tail boom, and other aerodynamic surfaces of a helicopter have a significant influence on its aerodynamic performance, the ride quality, and vibration. A numerical method for computing the aerodynamic interaction between an interacting vortex wake and the viscous flow about arbitrary two-dimensional bodies was developed to address this helicopter problem. The method solves for the flow field velocities on a body-fitted computational mesh using finite-difference techniques. The interacting vortex wake is represented by an array of discrete vortices which, in turn, are represented by a finite-core model. The evolution of the interacting vortex wake is calculated by Lagrangian techniques. The viscous flow field of the two-dimensional body is calculated on an Eulerian grid. The flow around circular and elliptic cylinders in the absence of an interacting vortex wake was calculated. These results compare very well with other numerical results and with results obtained from experiment and thereby demonstrate the accuracy of the viscous solution. The interaction of a rotor wake with the flow about a 4 to 1 elliptic cylinder at 45 degree incidence was calculated for a Reynolds number of 3000. The results demonstrate the significant variations in the lift and drag on the elliptic cylinder in the presence of the interacting rotor wake.

Stremel, Paul M.↗
NASA NTRS1991

A discourse on sensitivity analysis for discretely-modeled structures

A descriptive review is presented of the most recent methods for performing sensitivity analysis of the structural behavior of discretely-modeled systems. The methods are generally but not exclusively aimed at finite element modeled structures. Topics included are: selections of finite difference step sizes; special consideration for finite difference sensitivity of iteratively-solved response problems; first and second derivatives of static structural response; sensitivity of stresses; nonlinear static response sensitivity; eigenvalue and eigenvector sensitivities for both distinct and repeated eigenvalues; and sensitivity of transient response for both linear and nonlinear structural response.

Adelman, Howard M.↗