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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Incompressible viscous flow computations for the pump components and the artificial heart

A finite-difference, three-dimensional incompressible Navier-Stokes formulation to calculate the flow through turbopump components is utilized. The solution method is based on the pseudocompressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Both steady and unsteady flow calculations can be performed using the current algorithm. In this work, the equations are solved in steadily rotating reference frames by using the steady-state formulation in order to simulate the flow through a turbopump inducer. Eddy viscosity is computed by using an algebraic mixing-length turbulence model. Numerical results are compared with experimental measurements and a good agreement is found between the two. Included in the appendix is a paper on incompressible viscous flow through artificial heart devices with moving boundaries. Time-accurate calculations, such as impeller and diffusor interaction, will be reported in future work.

Kiris, Cetin↗
NASA NTRS1992

Computation of incompressible viscous flows through turbopump components

A finite-difference, three-dimensional, incompressible Navier-Stokes formulation for calculating the flow through turbopump components is presented. The solution method is based on the pseudocompressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line-relaxation method. Both steady and unsteady flow calculations can be performed using the presented algorithm. In this paper, the equations are solved in steadily rotating reference frames by using the steady-state formulation in order to simulate the flow through a turbopump inducer. Eddy viscosity is computed by using the Baldwin-Lomax model. Numerical results are compared with experimental measurements and good agreement is found between the two. Time-accurate calculations will be reported in future publications.

Kiris, Cetin↗
NASA NTRS1992

The prediction of noise radiation from supersonic elliptic jets

This paper describes the prediction of noise radiation from supersonic elliptic jets. The noise is associated with the large scale structures in the jet mixing layer. These structures are described as instability waves. The local characteristics of the instability waves are determined from a compressible, linear, analysis. The jet mean velocity and density are described in elliptic cylindrical coordinates. The local eigensolution for the instability waves is determined from a finite difference solution of the non-separable boundary value problem. This inner solution which is formulated in terms of the method of multiple scales is matched with the radiated field using the method of matched asymptotic expansions. The form of the far-field directivity is derived. Predictions are presented for the noise radiation by the several modes of instability in the elliptic jet. The radiated field is not axisymmetric and certain modes radiate strongly in the directions of the major and minor axes of the jet. The extension of the present work to other geometries and flow fields is discussed.

Morris, Philip J.↗
NASA NTRS1995

Moving and adaptive grid methods for compressible flows

This paper describes adaptive grid methods developed specifically for compressible flow computations. The basic flow solver is a finite-volume implementation of Roe's flux difference splitting scheme or arbitrarily moving unstructured triangular meshes. The grid adaptation is performed according to geometric and flow requirements. Some results are included to illustrate the potential of the methodology.

Trepanier, Jean-Yves↗
NASA NTRS1998

Analysis of Linear Aerospike Plume Induced X-33 Base Heating Environment

Computational analysis is conducted to study the effect of an linear aerospike engine plume on the X-33 base-heating environment during ascent flight. To properly account for the freestream-body interaction and to allow for potential plume-induced flow-separation, the thermo-flowfield of the entire vehicle at several trajectory points is computed. A sequential grid-refinement technique is used in conjunction with solution-adaptive, patched, and embedded grid methods to limit the model to a manageable size. The computational methodology is based on a three-dimensional, finite-difference, viscous flow, chemically reacting, pressure-based computational fluid dynamics formulation, and a three-dimensional, finite-volume, spectral-line based weighted-sum-of-gray-gases absorption, computational radiation heat transfer formulation. The computed forebody and afterbody surface pressure coefficients and base pressure characteristic curves are compared with those of a cold-flow test. The predicted convective and radiative base-heat fluxes, the effect of base-bleed, and the potential of plume-induced flow separation are presented.

Wang, Ten-See↗
NASA NTRS1999

The Inclusion of In-Plane Stresses in Delamination Criteria

A study of delamination failure was conducted with emphasis on delamination criteria. Evidence is presented which supports the inclusion of the in-plane stresses in addition to the interlaminar stress terms in delamination criteria. The delamination is characterized as the failure of a resin rich region in between ply sets. The entire six component stress state in this resin layer is calculated through a finite element analysis, averaged over a dimension of 1.75 ply thicknesses, and used in a Modified von Mises Delamination Criterion. This criterion builds onto previous criteria by including all six stress components in the interply resin layer. The MVMDC shows good correlation to experimental data. The results show that the treatment of delamination as the failure of a finite interply resin layer is a valid method and that the MVMDC, considering the full stress state, accurately indicates delamination for different laminate families.

Fenske, Matthew T.↗
NASA NTRS2000

Magnetic Field Suppression of Flow in Semiconductor Melt

One of the most promising approaches for the reduction of convection during the crystal growth of conductive melts (semiconductor crystals) is the application of magnetic fields. Current technology allows the experimentation with very intense static fields (up to 80 KGauss) for which nearly convection free results are expected from simple scaling analysis in stabilized systems (vertical Bridgman method with axial magnetic field). However, controversial experimental results were obtained. The computational methods are, therefore, a fundamental tool in the understanding of the phenomena accounting during the solidification of semiconductor materials. Moreover, effects like the bending of the isomagnetic lines, different aspect ratios and misalignments between the direction of the gravity and magnetic field vectors can not be analyzed with analytical methods. The earliest numerical results showed controversial conclusions and are not able to explain the experimental results. Although the generated flows are extremely low, the computational task is a complicated because of the thin boundary layers. That is one of the reasons for the discrepancy in the results that numerical studies reported. Modeling of these magnetically damped crystal growth experiments requires advanced numerical methods. We used, for comparison, three different approaches to obtain the solution of the problem of thermal convection flows: (1) Spectral method in spectral superelement implementation, (2) Finite element method with regularization for boundary layers, (3) Multiquadric method, a novel method with global radial basis functions, that is proven to have exponential convergence. The results obtained by these three methods are presented for a wide region of Rayleigh and Hartman numbers. Comparison and discussion of accuracy, efficiency, reliability and agreement with experimental results will be presented as well.

Fedoseyev, A. I.↗
NASA NTRS1997

Dynamic Electrothermal Model of a Sputtered Thermopile Thermal Radiation Detector for Earth Radiation Budget Applications

The Clouds and the Earth's Radiant Energy System (CERES) is a program sponsored by the National Aeronautics and Space Administration (NASA) aimed at evaluating the global energy balance. Current scanning radiometers used for CERES consist of thin-film thermistor bolometers viewing the Earth through a Cassegrain telescope. The Thermal Radiation Group, a laboratory in the Department of Mechanical Engineering at Virginia Polytechnic Institute and State University, is currently studying a new sensor concept to replace the current bolometer: a thermopile thermal radiation detector. This next-generation detector would consist of a thermal sensor array made of thermocouple junction pairs, or thermopiles. The objective of the current research is to perform a thermal analysis of the thermopile. Numerical thermal models are particularly suited to solve problems for which temperature is the dominant mechanism of the operation of the device (through the thermoelectric effect), as well as for complex geometries composed of numerous different materials. Feasibility and design specifications are studied by developing a dynamic electrothermal model of the thermopile using the finite element method. A commercial finite element-modeling package, ALGOR, is used.

Weckmann, Stephanie↗
NASA NTRS1989

Prediction of the thermal environment and thermal response of simple panels exposed to radiant heat

A method of predicting the radiant heat flux distribution produced by a bank of tubular quartz heaters was applied to a radiant system consisting of a single unreflected lamp irradiating a flat metallic incident surface. In this manner, the method was experimentally verified for various radiant system parameter settings and used as a source of input for a finite element thermal analysis. Two finite element thermal analyses were applied to a thermal system consisting of a thin metallic panel exposed to radiant surface heating. A two-dimensional steady-state finite element thermal analysis algorithm, based on Galerkin's Method of Weighted Residuals (GFE), was formulated specifically for this problem and was used in comparison to the thermal analyzers of the Engineering Analysis Language (EAL). Both analyses allow conduction, convection, and radiation boundary conditions. Differences in the respective finite element formulation are discussed in terms of their accuracy and resulting comparison discrepancies. The thermal analyses are shown to perform well for the comparisons presented here with some important precautions about the various boundary condition models. A description of the experiment, corresponding analytical modeling, and resulting comparisons are presented.

Turner, Travis L.↗
NASA NTRS1987

Prediction of SA 349/2 GV blade loads in high speed flight using several rotor analyses

The influence of blade dynamics, dynamic stall, and transonic aerodynamics on the predictions of rotor loads in high-speed flight are presented. Data were obtained from an Aerospatiale Gazelle SA 349/2 helicopter with three Grande Vitesse blades. Several analyses are used for this investigation. First, blade dynamics effects on the correlation are studied using three rotor analyses which differ mainly in the method of calculating the blade elastic response. Next, an ONERA dynamic stall model is used to predict retreating blade stall. Finally, advancing blade aerodynamic loads are calculated using a NASA-developed rotorcraft analysis coupled with two transonic finite-difference analyses.

Gaubert, Michel↗
NASA NTRS1991

Basic research on design analysis methods for rotorcraft vibrations

The objective of the present work was to develop a method for identifying physically plausible finite element system models of airframe structures from test data. The assumed models were based on linear elastic behavior with general (nonproportional) damping. Physical plausibility of the identified system matrices was insured by restricting the identification process to designated physical parameters only and not simply to the elements of the system matrices themselves. For example, in a large finite element model the identified parameters might be restricted to the moduli for each of the different materials used in the structure. In the case of damping, a restricted set of damping values might be assigned to finite elements based on the material type and on the fabrication processes used. In this case, different damping values might be associated with riveted, bolted and bonded elements. The method itself is developed first, and several approaches are outlined for computing the identified parameter values. The method is applied first to a simple structure for which the 'measured' response is actually synthesized from an assumed model. Both stiffness and damping parameter values are accurately identified. The true test, however, is the application to a full-scale airframe structure. In this case, a NASTRAN model and actual measured modal parameters formed the basis for the identification of a restricted set of physically plausible stiffness and damping parameters.

Hanagud, S.↗
NASA NTRS1993

Effective thermoelastic and thermal properties of unidirectional fiber-reinforced composites and their sensitivity coefficients

Three-dimensional finite element models are used to assess the accuracy of the thermoelastic and thermal properties of unidirectional fiber-reinforced composites predicted by six different micromechanical models. The six models are: simple mechanics of materials type equations, fiber substructuring model, vanishing fiber diameter model, self-consistent model, Mori-Tanaka model, and method of cells. In addition, the finite element models are used to assess the accuracy of derivatives of the effective properties, with respect to each of the constituent material properties and fiber-volume ratio, computed using the six micromechanical models. The predictions of the micromechanical and finite element models for four advanced composite material systems are compared with experimental data. The results obtained in the present study show that the predictions of the Mori-Tanaka model and the method of cells are closer to those of the finite element models than those of all the other micromechanical models.

Noor, Ahmed K.↗
NASA NTRS1940

Flow around wings accompanied by separation of vortices

The flow around wings computed by the usual method leads in the case of a finite trailing edge to a stagnation point in the trailing edge due to the Kutta-Joukowsky condition of flow governing this region. As a result, the theoretical pressure distribution differs substantially from the experimental values in the vicinity of the trailing edge. The present report describes an alternative method of calculation in which the rear stagnation point no longer appears. The stream leaves the trailing edge tangentially on the pressure side and a similar tangential separation occurs on the suction side of the profile at a point slightly in front of the trailing edge.

Schmieden, C↗
NASA NTRS2014

Modified Displacement Transfer Functions for Deformed Shape Predictions of Slender Curved Structures with Varying Curvatives

To eliminate the need to use finite-element modeling for structure shape predictions, a new method was invented. This method is to use the Displacement Transfer Functions to transform the measured surface strains into deflections for mapping out overall structural deformed shapes. The Displacement Transfer Functions are expressed in terms of rectilinearly distributed surface strains, and contain no material properties. This report is to apply the patented method to the shape predictions of non-symmetrically loaded slender curved structures with different curvatures up to a full circle. Because the measured surface strains are not available, finite-element analysis had to be used to analytically generate the surface strains. Previously formulated straight-beam Displacement Transfer Functions were modified by introducing the curvature-effect correction terms. Through single-point or dual-point collocations with finite-elementgenerated deflection curves, functional forms of the curvature-effect correction terms were empirically established. The resulting modified Displacement Transfer Functions can then provide quite accurate shape predictions. Also, the uniform straight-beam Displacement Transfer Function was applied to the shape predictions of a section-cut of a generic capsule (GC) outer curved sandwich wall. The resulting GC shape predictions are quite accurate in partial regions where the radius of curvature does not change sharply.

displacement transfer functions↗
DOE OSTI2023

Effective sensor location for detection of change in structural dynamic response

This paper investigates the sensitivity of structural system response to the sensor location by investigating consequences of small changes in the location to the structural system response. The paper discusses how maximum observability (based on mode shape and the participation of that mode in the input provided) drives optimal location. The structural responses were investigated in terms of the g-rms response for various low-frequency inputs (pure sinusoids and real-life inputs such as an earthquake and trains). Results were then analyzed in the context of Modal Contributions Factors (MCF) and changes to the Force-to-response Transfer Functions (TRFs). A modal-matching process is first presented using a Matlab TM -based Finite Element Method (FEM) model of a cantilever beam and instrumentation to determine the location of a small mass based on three different criteria. Subsequently, the structural response is investigated using experiments and the FEM model. The accelerometer of small mass (at 1/3 height) was moved up or down to obtain changes in the structural response (TRF) to various realistic low-frequency inputs. Modal Contribution Factor (MCF) and derivative (slope) of the associated mode-shapes were correlated to the observed changes in TRFs. Results show how optimal sensor locations for detecting change in structural response can be based on the MCFs and the associated mode-shapes.

Woodall, J.↗
NASA NTRS1983

Gear Lubrication and Cooling Experiment and Analysis

A gear tooth temperature analysis was performed using a finite element method combined with a calculated heat input, a calculated oil jet impingement depth, and estimated heat transfer coefficients for the different parts of the gear tooth that are oil cooled and air cooled. Experimental measurements of gear tooth average surface temperature and gear tooth instantaneous surface temperature were made with a fast response, infrared, radiometric microscope. Increasing oil pressure has a significant effect on both average surface temperature and peak surface temperature at loads above 1895 N/cm(1083 lb/in) and speeds of 10,000 and 7500 rpm. Both increasing speed (from 5000 to 10,000 rpm) at constant speed cause a significant rise in the average surface temperature and in the instantaneous peak surface temperatures on the gear teeth. The oil jet pressure required to provide the best cooling for gears is the pressure required to obtain full gear tooth impingement. Calculated results for gear tooth temperatures were close to experimental results for various oil jet impingement depths for identical operating conditions.

Townsend, D. P.↗
NASA NTRS1989

A review of high-speed, convective, heat-transfer computation methods

The objective of this report is to provide useful engineering formulations and to instill a modest degree of physical understanding of the phenomena governing convective aerodynamic heating at high flight speeds. Some physical insight is not only essential to the application of the information presented here, but also to the effective use of computer codes which may be available to the reader. A discussion is given of cold-wall, laminar boundary layer heating. A brief presentation of the complex boundary layer transition phenomenon follows. Next, cold-wall turbulent boundary layer heating is discussed. This topic is followed by a brief coverage of separated flow-region and shock-interaction heating. A review of heat protection methods follows, including the influence of mass addition on laminar and turbulent boundary layers. Also discussed are a discussion of finite-difference computer codes and a comparison of some results from these codes. An extensive list of references is also provided from sources such as the various AIAA journals and NASA reports which are available in the open literature.

Tauber, Michael E.↗