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,675 records · Page 93

Development and application of linear and nonlinear methods for interpretation of lightning strikes to in-flight aircraft

Since 1980, NASA has been collecting direct strike lightning data by flying an instrumented F-106B aircraft into thunderstorms. The continuing effort to interpret the measured data is reported here. Both linear and nonlinear finite difference modeling techniques are applied to the problem of lightning triggered by an aircraft in a thunderstorm. Five different aircraft are analyzed to determine the effect of aircraft size and shape on lightning triggering. The effect of lightning channel impedance on aircraft response is investigated. The particle environment in thunderstorms and electric field enhancements by typical ice particles is also investigated.

Rudolph, Terence↗

Aeroacoustics: Fan noise and control; Duct acoustics; Rotor noise

The recent progress in duct acoustic research is reflected by fourteen papers which cover essentially the areas of linear duct acoustics, the characteristics of absorbent materials, and nonlinear duct acoustics. With respect to propagation in ducts, emphasis is placed on computational methods aimed at attaining effective economical modeling of propagation in variable-geometry and in hard-wall or soft-wall ducts. Computational methods developed to reduce computer processing and storage requirements are outlined, along with a wave envelope technique for reducing processing and storage needs. The latter technique is applied to the optimum segmented acoustic duct linear design. A finite difference technique is used in combination with conformal mapping to compute the optimum multisection duct linear design.

Schwartz, I. R.↗

Algorithm for Unsteady Potential Flow About Airfoils

Implicit finite-difference scheme efficiently computes unsteady potential flow about airfoils. Formulation uses density and velocity potential as dependent variables. Conservation form is retained to assure that shock wave location and speed are computed correctly. Scheme fills need for method to calculate efficiently unsteady potential flow about airfoils and to predict flutter and other unsteady aeroelastic phenomena in transonic flow regimes.

Chipman, R.↗

Numerical boundary condition procedures and multigrid methods; Proceedings of the Symposium, NASA Ames Research Center, Moffett Field, CA, October 19-22, 1981

Papers presented in this volume provide an overview of recent work on numerical boundary condition procedures and multigrid methods. The topics discussed include implicit boundary conditions for the solution of the parabolized Navier-Stokes equations for supersonic flows; far field boundary conditions for compressible flows; and influence of boundary approximations and conditions on finite-difference solutions. Papers are also presented on fully implicit shock tracking and on the stability of two-dimensional hyperbolic initial boundary value problems for explicit and implicit schemes.

Source record↗

Structural optimization by generalized, multilevel decomposition

The developments toward a general multilevel optimization capability and results for a three-level structural optimization are described. The method partitions a structure into a number of substructuring levels where each substructure corresponds to a subsystem in the general case of an engineering system. The method is illustrated by a portal framework that decomposes into individual beams. Each beam is a box that can be further decomposed into stiffened plates. Substructuring for this example spans three different levels: (1) the bottom level of finite elements representing the plates; (2) an intermediate level of beams treated as substructures; and (3) the top level for the assembled structure. The three-level case is now considered to be qualitatively complete.

Sobieszczanski-Sobieski, J.↗

A new approach to finite state modelling of unsteady aerodynamics

This paper presents a novel technique for formulating a high quality finite state unsteady aerodynamic model by applying Bode plot methods, used in control engineering. Indicial response functions for both fixed wing and rotary wing applications are obtained using these finite state unsteady aerodynamic models. It is shown that the rotary wing indicial response function has a fundamentally different characteristic when compared to fixed wing indicial response. The rotary wing indicial response function is oscillatory in nature while the fixed wing indicial response function is nonoscillatory. Furthermore it should be emphasized that this is the first that a rotary-wing indicial response function has been presented in the literature.

Venkatesan, C.↗

Characterization and modeling of an advanced flexible thermal protection material for space applications

This paper describes experimental and analytical characterization of a new flexible thermal protection material known as Tailorable Advanced Blanket Insulation (TABI). This material utilizes a three-dimensional ceramic fabric core structure and an insulation filler. TABI is the leading candidate for use in deployable aeroassisted vehicle designs. Such designs require extensive structural modeling, and the most significant in-plane material properties necessary for model development are measured and analytically verified in this study. Unique test methods are developed for damping measurements. Mathematical models are developed for verification of the experimental modulus and damping data, and finally, transverse properties are described in terms of the inplane properties through use of a 12-dof finite difference model of a simple TABI configuration.

Clayton, Joseph P.↗

Advanced Rotorcraft Transmission (ART) program status

Reported herein is work done on the Advanced Rotorcraft Transmission by McDonnell Douglas Helicopter Company under Army/NASA contract. The novel concept pursued includes the use of face gears for power transmission and a torque splitting arrangement. The design reduces the size and weight of the corner-turning hardware and the next reduction stage. New methods of analyzing face gears have increased confidence in their usefulness. Test gears have been designed and manufactured for power transmission testing on the NASA-Lewis spiral bevel test rig. Transmission design effort has included finite element modeling of the split torque paths to assure equal deflection under load. A finite element model of the Apache main transmission has been completed to substantiate noise prediction methods. A positive engagement overrunning clutch design is described. Test spur gears have been made by near-net-shape forging from five different materials. Three housing materials have been procured for evaluation testing.

Bossler, Robert↗

Active control of asymmetric vortical flows around cones using injection and heating

The effectiveness of certain active-control methods for asymmetric flows around circular cones is investigated by using computational solution of the unsteady, compressible full Navier-Stokes equations. Two main methods of active control which include flow injection and surface heating are used. For the flow-injection-control method, flow injection is used either in the normal direction to the surface or in the tangential direction to the surface. For the surface-heating-control method, the temperature of the cone surface is increased. The effectiveness of a hybrid method of flow control which combines normal injection with surface heating has also been studied. The Navier-Stokes equations, subjected to various surface boundary conditions, are solved by using an implicit, upwind, flux-difference splitting, finite-volume scheme for locally-conical flow solutions.

Kandil, Osama A.↗

On the decay of the magnetic fields of single radio pulsars

We investigate the statistical evidence for the decay of the magnetic field of single radio pulsars. We perform population syntheses for different assumed values for the time scale of field decay using a Monte Carlo method. We allow for the selection effects in pulsar surveys and compare the synthesized populations with the observed pulsars. We take account of the finite scale height of the distribution in the Galaxy of free electrons, which determine the dispersion measure and hence the apparent distance of radio pulsars. Our simulations give much better agreement with the observations if the time scale for the field decay is assumed to be longer than the typical active life time of a radio pulsar. This indicates that no significant field decay occurs in single radio pulsars.

Bhattacharya, Dipankar↗

A Parallel, Finite-Volume Algorithm for Large-Eddy Simulation of Turbulent Flows

A parallel, finite-volume algorithm has been developed for large-eddy simulation (LES) of compressible turbulent flows. This algorithm includes piecewise linear least-square reconstruction, trilinear finite-element interpolation, Roe flux-difference splitting, and second-order MacCormack time marching. Parallel implementation is done using the message-passing programming model. In this paper, the numerical algorithm is described. To validate the numerical method for turbulence simulation, LES of fully developed turbulent flow in a square duct is performed for a Reynolds number of 320 based on the average friction velocity and the hydraulic diameter of the duct. Direct numerical simulation (DNS) results are available for this test case, and the accuracy of this algorithm for turbulence simulations can be ascertained by comparing the LES solutions with the DNS results. The effects of grid resolution, upwind numerical dissipation, and subgrid-scale dissipation on the accuracy of the LES are examined. Comparison with DNS results shows that the standard Roe flux-difference splitting dissipation adversely affects the accuracy of the turbulence simulation. For accurate turbulence simulations, only 3-5 percent of the standard Roe flux-difference splitting dissipation is needed.

Bui, Trong T.↗

A Numerical and Experimental Study of Coflow Laminar Diffusion Flames: Effects of Gravity and Inlet Velocity

In this work, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting behavior of laminar coflow methane-air diffusion flames was investigated both computationally and experimentally. A series of flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) was assessed numerically under microgravity and normal gravity conditions with the fuel stream CH4 mole fraction ranging from 0.4 to 1.0. Computationally, the MC-Smooth vorticity-velocity formulation of the governing equations was employed to describe the reactive gaseous mixture; the soot evolution process was considered as a classical aerosol dynamics problem and was represented by the sectional aerosol equations. Since each flame is axisymmetric, a two-dimensional computational domain was employed, where the grid on the axisymmetric domain was a nonuniform tensor product mesh. The governing equations and boundary conditions were discretized on the mesh by a nine-point finite difference stencil, with the convective terms approximated by a monotonic upwind scheme and all other derivatives approximated by centered differences. The resulting set of fully coupled, strongly nonlinear equations was solved simultaneously using a damped, modified Newton's method and a nested Bi-CGSTAB linear algebra solver. Experimentally, the flame shape, size, lift-off height, and soot temperature were determined by flame emission images recorded by a digital camera, and the soot volume fraction was quantified through an absolute light calibration using a thermocouple. For a broad spectrum of flames in microgravity and normal gravity, the computed and measured flame quantities (e.g., temperature profile, flame shape, lift-off height, and soot volume fraction) were first compared to assess the accuracy of the numerical model. After its validity was established, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting tendency of laminar coflow methane-air diffusion flames was explored further by examining quantities derived from the computational results.

microgravity↗

Numerical Modeling of an Integrated Vehicle Fluids System Loop for Pressurizing a Cryogenic Tank

This paper presents a numerical model of the pressurization loop of the Integrated Vehicle Fluids (IVF) system using the Generalized Fluid System Simulation Program (GFSSP). The IVF propulsion system, being developed by United Launch Alliance to reduce system weight and enhance reliability, uses boiloff propellants to drive thrusters for the reaction control system as well as to run internal combustion engines to develop power and drive compressors to pressurize propellant tanks. NASA Marshall Space Flight Center (MSFC) conducted tests to verify the functioning of the IVF system using a flight-like tank. GFSSP, a finite volume based flow network analysis software developed at MSFC, has been used to support the test program. This paper presents the simulation of three different test series, comparison of numerical prediction and test data and a novel method of presenting data in a dimensionless form. The paper also presents a methodology of implementing a compressor map in a system level code.

LeClair, A. C.↗

Design of a Modular Monolithic Implicit Solver for Multi-Physics Applications

The design of a modular multi-physics high-order space-time finite-element framework is presented together with its extension to allow monolithic coupling of different physics. One of the main objectives of the framework is to perform efficient high- fidelity simulations of capsule/parachute systems. This problem requires simulating multiple physics including, but not limited to, the compressible Navier-Stokes equations, the dynamics of a moving body with mesh deformations and adaptation, the linear shell equations, non-re effective boundary conditions and wall modeling. The solver is based on high-order space-time - finite element methods. Continuous, discontinuous and C1-discontinuous Galerkin methods are implemented, allowing one to discretize various physical models. Tangent and adjoint sensitivity analysis are also targeted in order to conduct gradient-based optimization, error estimation, mesh adaptation, and flow control, adding another layer of complexity to the framework. The decisions made to tackle these challenges are presented. The discussion focuses first on the "single-physics" solver and later on its extension to the monolithic coupling of different physics. The implementation of different physics modules, relevant to the capsule/parachute system, are also presented. Finally, examples of coupled computations are presented, paving the way to the simulation of the full capsule/parachute system.

Carton De Wiart, Corentin↗

A fast Cauchy-Riemann solver

The inhomogeneous Cauchy-Riemann equations in a rectangle are discretized by a finite difference approximation. Several different boundary conditions are treated explicitly, leading to algorithms which have overall second-order accuracy. All boundary conditions with either u or v prescribed along a side of the rectangle can be treated by similar methods. The algorithms presented here have nearly minimal time and storage requirements and seem suitable for development into a general-purpose direct Cauchy-Riemann solver for arbitrary boundary conditions.

Ghil, M.↗

A fast hidden line algorithm for plotting finite element models

Effective plotting of finite element models requires the use of fast hidden line plot techniques that provide interactive response. A high speed hidden line technique was developed to facilitate the plotting of NASTRAN finite element models. Based on testing using 14 different models, the new hidden line algorithm (JONES-D) appears to be very fast: its speed equals that for normal (all lines visible) plotting and when compared to other existing methods it appears to be substantially faster. It also appears to be very reliable: no plot errors were observed using the new method to plot NASTRAN models. The new algorithm was made part of the NPLOT NASTRAN plot package and was used by structural analysts for normal production tasks.

Jones, G. K.↗

Stability analysis of flexible wind turbine blades using finite element method

Static vibration and flutter analysis of a straight elastic axis blade was performed based on a finite element method solution. The total potential energy functional was formulated according to linear beam theory. The inertia and aerodynamic loads were formulated according to the blade absolute acceleration and absolute velocity vectors. In vibration analysis, the direction of motion of the blade during the first out-of-lane and first in-plane modes was examined; numerical results involve NASA/DOE Mod-0, McCauley propeller, north wind turbine and flat plate behavior. In flutter analysis, comparison cases were examined involving several references. Vibration analysis of a nonstraight elastic axis blade based on a finite element method solution was performed in a similar manner with the straight elastic axis blade, since it was recognized that a curved blade can be approximated by an assembly of a sufficient number of straight blade elements at different inclinations with respect to common system of axes. Numerical results involve comparison between the behavior of a straight and a curved cantilever beam during the lowest two in-plane and out-of-plane modes.

Kamoulakos, A.↗

Theory of molecular formation by radiative association in interstellar clouds

A theory of molecular formation by radiative association is presented which is consistent with angular momentum being conserved during the encounter and which incorporates explicitly the long-range attraction between the reactants. It is pointed out that the activated complex would not have a Boltzmann energy distribution should the rotational and kinetic temperatures of the reactants differ, and it is shown how allowance for this may be made. Account is taken of the inaccessibility of a considerable fraction of the nuclear spin states of the complex. Methods are given for treating the effect which the finiteness of the dissociation frequency may have on the association rate. Calculations on some reactions of interest are performed. A very simple semiempirical formula for the rate coefficient for radiative association is also given.

Bates, D. R.↗