SEARCH · Engineering Papers

Results for “boundary element method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Transition to turbulence in plane channel flows

Results obtained from a numerical simulation of the final stages of transition to turbulence in plane channel flow are described. Three dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time evolution of two and three dimensional finite amplitude disturbances. Computations are performed on CYBER-203 vector processor for a 32x51x32 grid. Results are presented for no-slip boundary conditions at the solid walls as well as for periodic suction blowing to simulate active control of transition by mass transfer. Solutions indicate that the method is capable of simulating the complex character of vorticity dynamics during the various stages of transition and final breakdown. In particular, evidence points to the formation of a lambda-shape vortex and the subsequent system of horseshoe vortices inclined to the main flow direction as the main elements of transition. Calculations involving periodic suction-blowing indicate that interference with a wave of suitable phase and amplitude reduces the disturbance growth rates.

Biringen, S.↗

NASA NTRS1983

A finite element method for diffusion dominated unsteady viscous flows

A general conforming finite element scheme for computing viscous flows is presented which is of second-order accuracy in space and time. Viscous terms are treated implicitly and advection terms are treated explicitly in the time marching segment of the algorithm. A method for solving the algebraic equations at each time step is given. The method is demonstrated on two test problems, one of them being a plane vortex flow for which asymptotic methods are used to obtain suitable numerical boundary conditions at each time step.

Gunzburger, M. D.↗

NASA NTRS1982

Natural convection in melt crystal growth - The influence of flow pattern on solute segregation

The results of two lines of research aimed at calculating the structure of the flows driven by buoyancy in small-scale crystal growth systems and at understanding the coupling between these flows, the shape of the solidification interface, and dopant segregation in the crystal are reviewed. First, finite-element methods are combined with computer-aided methods for detecting multiple steady solutions to analyze the structure of the buoyancy-driven axisymmetric flows in a vertical cylinder heated from below. This system exhibits onset of convection, multiple steady flows, and loss of the primary stable flow beyond a critical value of the Rayleigh number. Second, results are presented for calculations of convection, melt/solid interface shape, and dopant segregation within a vertical ampoule with thermal boundary conditions that represent a prototype of the vertical Bridgman growth system.

Brown, R. A.↗

NASA NTRS1987

Robustness issues in boundary feedback of flexible structures

Transfer function models for basic structural elements with boundary control reveal certain inherent properties relevant to questions of control system realization, design, and analysis of robustness. The transfer functions involved are not strictly proper, may be nonminimum phase (except for the special case of colocated actuation and sensing), and often have large number of poles on the imaginary axis. In this paper these basic questions are considered in terms of modeling and control system design for robustness. The application of algebraic methods for computing stabilizing control, based on certain exact, irrational transfer functions is investigated. Boundary control of the wave equation is considered and extension of the method to more general problems is suggested. Special attention is given to implementation issues associated with a class of infinite dimensional control laws.

Bennett, W. H.↗

NASA NTRS1956

Aerodynamic Design of Axial-flow Compressors. Volume III

Chapters XI to XIII concern the unsteady compressor operation arising when compressor blade elements stall. The fields of compressor stall and surge are reviewed in Chapters XI and XII, respectively. The part-speed operating problem in high-pressure-ratio multistage axial-flow compressors is analyzed in Chapter XIII. Chapter XIV summarizes design methods and theories that extend beyond the simplified two-dimensional approach used previously in the report. Chapter XV extends this three-dimensional treatment by summarizing the literature on secondary flows and boundary layer effects. Charts for determining the effects of errors in design parameters and experimental measurements on compressor performance are given in Chapters XVI. Chapter XVII reviews existing literature on compressor and turbine matching techniques.

COMPRESSORS - PRESSURE RATIO↗

NASA NTRS2004

A New Approach to Identify Optimal Properties of Shunting Elements for Maximum Damping of Structural Vibration Using Piezoelectric Patches

The use of shunted piezoelectric patches in reducing vibration and sound radiation of structures has several advantages over passive viscoelastic elements, e.g., lower weight with increased controllability. The performance of the piezoelectric patches depends on the shunting electronics that are designed to dissipate vibration energy through a resistive element. In past efforts most of the proposed tuning methods were based on modal properties of the structure. In these cases, the tuning applies only to one mode of interest and maximum tuning is limited to invariant points when based on den Hartog's invariant points concept. In this study, a design method based on the wave propagation approach is proposed. Optimal tuning is investigated depending on the dynamic and geometric properties that include effects from boundary conditions and position of the shunted piezoelectric patch relative to the structure. Active filters are proposed as shunting electronics to implement the tuning criteria. The developed tuning methods resulted in superior capabilities in minimizing structural vibration and noise radiation compared to other tuning methods. The tuned circuits are relatively insensitive to changes in modal properties and boundary conditions, and can applied to frequency ranges in which multiple modes have effects.

Park, Junhong↗

NASA NTRS1989

The simulation of 2D compressible viscous high speed flow by the finite element method

Am implicit/explicit procedure for the solution of problems of two-dimensional steady compressible viscous high-speed flows is presented. In the vicinity of solid walls, a grid which need only exhibit structure in the normal direction is employed while, away from this region, the grid is totally unstructured. The implicit form of the algorithm is used near solid walls, with the grid structure being utilized in an equation solution approach, based upon line relaxation. The explicit form of the algorithm is used elsewhere. Grid adaptation is achieved by means of adaptive remeshing. To illustrate the performance of the proposed method, problems of shock-boundary layer interaction and flow over a simulated forebody at high Mach number are included.

Hassan, O.↗

NASA NTRS1983

Manipulator interactive design with interconnected flexible elements

This paper describes the development of an analysis tool for the interactive design of control systems for manipulators and similar electro-mechanical systems amenable to representation as structures in a topological chain. The chain consists of a series of elastic bodies subject to small deformations and arbitrary displacements. The bodies are connected by hinges which permit kinematic constraints, control, or relative motion with six degrees of freedom. The equations of motion for the chain configuration are derived via Kane's method, extended for application to interconnected flexible bodies with time-varying boundary conditions. A corresponding set of modal coordinates has been selected. The motion equations are imbedded within a simulation that transforms the vector-dyadic equations into scalar form for numerical integration. The simulation also includes a linear, time-invariant controler specified in transfer function format and a set of sensors and actuators that interface between the structure and controller. The simulation is driven by an interactive set-up program resulting in an easy-to-use analysis tool.

Singh, R. P.↗

NASA NTRS1975

Analytical study of mixing and reacting three-dimensional supersonic combustor flow fields

An analytical investigation is presented of mixing and reacting hydrogen jets injected from multiple orifices transverse and parallel to a supersonic airstream. The COMOC computer program, based upon a finite-element solution algorithm, was developed to solve the governing equations for three-dimensional, turbulent, reacting, boundary-region, and confined flow fields. The computational results provide a three-dimensional description of the velocity, temperature, and species-concentration fields downstream of hydrogen injection. Detailed comparisons between cold-flow data and results of the computational analysis have established validity of the turbulent-mixing model based on the elementary mixing-length hypothesis. A method is established to initiate computations for reacting flow fields based upon cold-flow correlations and the appropriate experimental parameters of Mach number, injector spacing, and pressure ratio. Key analytical observations on mixing and combustion efficiency for reacting flows are presented and discussed.

Baker, A. J.↗

NASA NTRS1983

An investigation of the accuracy of finite difference methods in the solution of linear elasticity problems

The accuracy of the finite difference method in the solution of linear elasticity problems that involve either a stress discontinuity or a stress singularity is considered. Solutions to three elasticity problems are discussed in detail: a semi-infinite plane subjected to a uniform load over a portion of its boundary; a bimetallic plate under uniform tensile stress; and a long, midplane symmetric, fiber reinforced laminate subjected to uniform axial strain. Finite difference solutions to the three problems are compared with finite element solutions to corresponding problems. For the first problem a comparison with the exact solution is also made. The finite difference formulations for the three problems are based on second order finite difference formulas that provide for variable spacings in two perpendicular directions. Forward and backward difference formulas are used near boundaries where their use eliminates the need for fictitious grid points.

Bauld, N. R., Jr.↗

NASA NTRS1983

Transition to turbulence in plane channel flow

A numerical simulation of the final stages of transition to turbulence in plane channel flow is reported. Three dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time-evolution of two and three dimensional finite amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32x51x32 grid. Results are presented for no-slip boundary conditions at the solid walls as well as for periodic suction-blowing to simulate active control of transition by mass transfer. Solutions indicate that the method is capable of simulating the complex character of vorticity dynamics during the various stages of transition and final breakdown. In particular, evidence points to the formation of a lambda-shape vortex and the subsequent system of horseshoe vortices inclined to the main flow direction as the main elements of transition. Calculations involving suction-blowing indicate that interference with a wave of suitable phase and amplitude reduces the disturbance growth rates.

Biringen, S.↗

NASA NTRS1989

Extension-torsion coupling behavior of advanced composite tilt-rotor blades

An analytic model was developed to study the extension-bend-twist coupling behavior of an advanced composite helicopter or tilt-rotor blade. The outer surface of the blade is defined by rotating an arbitrary cross section about an initial twist axis. The cross section can be nonhomogeneous and composed of generally anisotropic materials. The model is developed based upon a three dimensional elasticity approach that is recast as a coupled two-dimensional boundary value problem defined in a curvilinear coordinate system. Displacement solutions are written in terms of known functions that represent extension, bending, and twisting and unknown functions for local cross section deformations. The unknown local deformation functions are determined by applying the principle of minimum potential energy to the discretized two-dimensional cross section. This is an application of the Ritz method, where the trial function family is the displacement field associated with a finite element (8-node isoparametric quadrilaterals) representation of the section. A computer program was written where the cross section is discretized into 8-node quadrilateral subregions. Initially the program was verified using previously published results (both three-dimensional elasticity and technical beam theory) for pretwisted isotropic bars with an elliptical cross section. In addition, solid and thin-wall multi-cell NACA-0012 airfoil sections were analyzed to illustrate the pronounced effects that pretwist, initial twist axis location, and spar location has on coupled behavior. Currently, a series of advanced composite airfoils are being modeled in order to assess how the use of laminated composite materials interacts with pretwist to alter the coupling behavior of the blade. These studies will investigate the use of different ply angle orientations and the use of symmetric versus unsymmetric laminates.

Kosmatka, J. B.↗

NASA NTRS2007

Method and program product for determining a radiance field in an optical environment

A hybrid method is presented by which Monte Carlo techniques are combined with iterative relaxation techniques to solve the Radiative Transfer Equation in arbitrary one-, two- or three-dimensional optical environments. The optical environments are first divided into contiguous regions, or elements, with Monte Carlo techniques then being employed to determine the optical response function of each type of element. The elements are combined, and the iterative relaxation techniques are used to determine simultaneously the radiance field on the boundary and throughout the interior of the modeled environment. This hybrid model is capable of providing estimates of the under-water light field needed to expedite inspection of ship hulls and port facilities. It is also capable of providing estimates of the subaerial light field for structured, absorbing or non-absorbing environments such as shadows of mountain ranges within and without absorption spectral bands such as water vapor or CO.sub.2 bands.

Reinersman, Phillip N.↗

NASA NTRS1993

Data Systems Dynamic Simulator

The Data System Dynamic Simulator (DSDS) is a discrete event simulation tool. It was developed for NASA for the specific purpose of evaluating candidate architectures for data systems of the Space Station era. DSDS provides three methods for meeting this requirement. First, the user has access to a library of standard pre-programmed elements. These elements represent tailorable components of NASA data systems and can be connected in any logical manner. Secondly, DSDS supports the development of additional elements. This allows the more sophisticated DSDS user the option of extending the standard element set. Thirdly, DSDS supports the use of data streams simulation. Data streams is the name given to a technique that ignores packet boundaries, but is sensitive to rate changes. Because rate changes are rare compared to packet arrivals in a typical NASA data system, data stream simulations require a fraction of the CPU run time. Additionally, the data stream technique is considerably more accurate than another commonly-used optimization technique.

Rouff, Christopher↗

NASA NTRS1994

Non-Reflecting Regions for Finite Difference Methods in Modeling of Elastic Wave Propagation in Plates

Solution of the wave equation using techniques such as finite difference or finite element methods can model elastic wave propagation in solids. This requires mapping the physical geometry into a computational domain whose size is governed by the size of the physical domain of interest and by the required resolution. This computational domain, in turn, dictates the computer memory requirements as well as the calculation time. Quite often, the physical region of interest is only a part of the whole physical body, and does not necessarily include all the physical boundaries. Reduction of the calculation domain requires positioning an artificial boundary or region where a physical boundary does not exist. It is important however that such a boundary, or region, will not affect the internal domain, i.e., it should not cause reflections that propagate back into the material. This paper concentrates on the issue of constructing such a boundary region.

Kishoni, Doron↗

NASA NTRS1959

Source Distribution Method for Unsteady One-Dimensional Flows With Small Mass, Momentum, and Heat Addition and Small Area Variation

A source distribution method is presented for obtaining flow perturbations due to small unsteady area variations, mass, momentum, and heat additions in a basic uniform (or piecewise uniform) one-dimensional flow. First, the perturbations due to an elemental area variation, mass, momentum, and heat addition are found. The general solution is then represented by a spatial and temporal distribution of these elemental (source) solutions. Emphasis is placed on discussing the physical nature of the flow phenomena. The method is illustrated by several examples. These include the determination of perturbations in basic flows consisting of (1) a shock propagating through a nonuniform tube, (2) a constant-velocity piston driving a shock, (3) ideal shock-tube flows, and (4) deflagrations initiated at a closed end. The method is particularly applicable for finding the perturbations due to relatively thin wall boundary layers.

Mirels, Harold↗

NASA NTRS2007

One-Dimensional Ablation with Pyrolysis Gas Flow Using a Full Newton's Method and Finite Control Volume Procedure

The development and verification of a one-dimensional material thermal response code with ablation is presented. The implicit time integrator, control volume finite element spatial discretization, and Newton's method for nonlinear iteration on the entire system of residual equations have been implemented and verified for the thermochemical ablation of internally decomposing materials. This study is a continuation of the work presented in "One-Dimensional Ablation with Pyrolysis Gas Flow Using a Full Newton's Method and Finite Control Volume Procedure" (AIAA-2006-2910), which described the derivation, implementation, and verification of the constant density solid energy equation terms and boundary conditions. The present study extends the model to decomposing materials including decomposition kinetics, pyrolysis gas flow through the porous char layer, and a mixture (solid and gas) energy equation. Verification results are presented for the thermochemical ablation of a carbon-phenolic ablator which involves the solution of the entire system of governing equations.

Amar, Adam J.↗

NASA NTRS1974

Integral representation of field variables for the finite element solution of viscous flow problems

The kinematics aspects of compressible and incompressible viscous flow problems are recasted into an integral representation for the velocity vector. For the incompressible flow problem, the kinetics aspect is similarly recasted into an integral representation. The result is a system of integral equations ideally suited for the finite element method. The integral representations are shown to permit the confinement of the computation field to the region of nonnegligible vorticity and dilatation which, for the incompressible flow, is identified with the viscous region. This drastic reduction in computation field results in a superior solution speed which is further improved by the use of a flowfield segmentation technique. The new approach is further shown to remove the difficulties associated with previous methods in specifying the far field and extraneous boundary conditions.

Wu, J. C.↗