Engineering Papers⌕ Search

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.

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

At least 631 records · Page 35

The analysis and design of transonic two-element airfoil systems

The multiphase effort in the development of tools for the analysis and design of two-element airfoil systems, that is, airfoils with a slat or a flap at transonic speeds is described. The first phase involved the development of a method to compute the inviscid flow over such configurations. In the second phase the inviscid code was coupled to a boundary layer calculation program in order to compute the loss in performance due to viscous effects. An inverse code that constructs the airfoil system corresponding to a desired pressure distribution is described.

Volpe, G.↗

Navier-Stokes calculations for unsteady three-dimensional vortical flows in unbounded domains

Finite-difference Navier-Stokes calculations for unsteady, three-dimensional, incompressible, viscous flows induced by initial vorticity distributions are presented and discussed in this paper. The initial vorticity distributions are assumed to be embedded in a flow field of infinite extent that is quiescent at infinity. These vorticity distributions are typical of vortex rings and other closed vortical tubes or structures. Such structures are important elements in fluid flows such as jets, atmospheric convection and the far-field wakes of aircraft; studies of their interaction may aid in an understanding of complex fluid flows. The calculations employ a method recently proposed by Ting to approximate the infinite-domain boundary value problem with a finite boundary computational domain, and this method is shown to yield accurate three-dimensional results for reasonable expenditures of computer time. Because of the efficiency of the boundary condition technique and the resulting Navier-Stokes code, a 16-bit minicomputer with virtual memory was capable of performing the calculations for the unsteady motion of two obliquely colliding vortex rings. The results of these calculations are presented in the paper.

Chamberlain, J. P.↗

A first-order Green's function approach to supersonic oscillatory flow: A mixed analytic and numeric treatment

A frequency domain Green's Function Method for unsteady supersonic potential flow around complex aircraft configurations is presented. The focus is on the supersonic range wherein the linear potential flow assumption is valid. In this range the effects of the nonlinear terms in the unsteady supersonic compressible velocity potential equation are negligible and therefore these terms will be omitted. The Green's function method is employed in order to convert the potential flow differential equation into an integral one. This integral equation is then discretized, through standard finite element technique, to yield a linear algebraic system of equations relating the unknown potential to its prescribed co-normalwash (boundary condition) on the surface of the aircraft. The arbitrary complex aircraft configuration (e.g., finite-thickness wing, wing-body-tail) is discretized into hyperboloidal (twisted quadrilateral) panels. The potential and co-normalwash are assumed to vary linearly within each panel. The long range goal is to develop a comprehensive theory for unsteady supersonic potential aerodynamic which is capable of yielding accurate results even in the low supersonic (i.e., high transonic) range.

Freedman, M. I.↗

Scattering of Lamb waves in a composite plate

A combined analytical and finite element technique is developed to gain a better understanding of the scattering of elastic waves by defects. This hybrid method is capable of predicting scattered displacements from arbitrary shaped defects as well as inclusions of different material. The continuity of traction and displacements at the boundaries of the two areas provided the necessary equations to find the nodal displacements and expansion coefficients. Results clearly illustrate the influence of increasing crack depth on the scattered signal.

Bratton, Robert↗

An advancing front Delaunay triangulation algorithm designed for robustness

A new algorithm is described for generating an unstructured mesh about an arbitrary two-dimensional configuration. Mesh points are generated automatically by the algorithm in a manner which ensures a smooth variation of elements, and the resulting triangulation constitutes the Delaunay triangulation of these points. The algorithm combines the mathematical elegance and efficiency of Delaunay triangulation algorithms with the desirable point placement features, boundary integrity, and robustness traditionally associated with advancing-front-type mesh generation strategies. The method offers increased robustness over previous algorithms in that it cannot fail regardless of the initial boundary point distribution and the prescribed cell size distribution throughout the flow-field.

Mavriplis, D. J.↗

An advancing front Delaunay triangulation algorithm designed for robustness

A new algorithm is described for generating an unstructured mesh about an arbitrary two-dimensional configuration. Mesh points are generated automatically by the algorithm in a manner which ensures a smooth variation of elements, and the resulting triangulation constitutes the Delaunay triangulation of these points. The algorithm combines the mathematical elegance and efficiency of Delaunay triangulation algorithms with the desirable point placement features, boundary integrity, and robustness traditionally associated with advancing-front-type mesh generation strategies. The method offers increased robustness over previous algorithms in that it cannot fail regardless of the initial boundary point distribution and the prescribed cell size distribution throughout the flow-field.

Mavriplis, D. J.↗

Nonlinear vibrations of thin arbitrarily laminated composite plates subjected to harmonic excitations using DKT elements

A finite element formulation is presented for determining the large-amplitude free and steady-state forced vibration response of arbitrarily laminated anisotropic composite thin plates using the Discrete Kirchhoff Theory (DKT) triangular elements. The nonlinear stiffness and harmonic force matrices of an arbitrarily laminated composite triangular plate element are developed for nonlinear free and forced vibration analyses. The linearized updated-mode method with nonlinear time function approximation is employed for the solution of the system nonlinear eigenvalue equations. The amplitude-frequency relations for convergence with gridwork refinement, triangular plates, different boundary conditions, lamination angles, number of plies, and uniform versus concentrated loads are presented.

Chiang, C. K.↗

Models of space averaged energetics of plates

The analysis of high frequency vibrations in plates is of particular interest in the study of structure borne noise in aircrafts. The current methods of analysis are either too expensive (finite element method) or may have a confidence band wider than desirable (Statistical Energy Analysis). An alternative technique to model the space and time averaged response of structural acoustics problems with enough detail to include all significant mechanisms of energy generation, transmission, and absorption is highly desirable. The focus of this paper is the development of a set of equations which govern the space and time averaged energy density in plates. To solve this equation, a new type of boundary value problem must be treated in terms of energy density variables using energy and intensity boundary conditions. A computer simulation verification study of the energy governing equation is performed. A finite element formulation of the new equations is also implemented and several test cases are analyzed and compared to analytical solutions.

Bouthier, O. M.↗

Self-adaptive closed constrained solution algorithms for nonlinear conduction

Self-adaptive solution algorithms are developed for nonlinear heat conduction problems encountered in analyzing materials for use in high temperature or cryogenic conditions. The nonlinear effects are noted to occur due to convection and radiation effects, as well as temperature-dependent properties of the materials. Incremental successive substitution (ISS) and Newton-Raphson (NR) procedures are treated as extrapolation schemes which have solution projections bounded by a hyperline with an externally applied thermal load vector arising from internal heat generation and boundary conditions. Closed constraints are formulated which improve the efficiency and stability of the procedures by employing closed ellipsoidal surfaces to control the size of successive iterations. Governing equations are defined for nonlinear finite element models, and comparisons are made of results using the the new method and the ISS and NR schemes for epoxy, PVC, and CuGe.

Padovan, J.↗

Vectorizable algorithms for adaptive schemes for rapid analysis of SSME flows

An initial study into vectorizable algorithms for use in adaptive schemes for various types of boundary value problems is described. The focus is on two key aspects of adaptive computational methods which are crucial in the use of such methods (for complex flow simulations such as those in the Space Shuttle Main Engine): the adaptive scheme itself and the applicability of element-by-element matrix computations in a vectorizable format for rapid calculations in adaptive mesh procedures.

Oden, J. Tinsley↗

Computational considerations in application of the finite element method for analysis of unsteady flow around airfoils.

Difficulties related to the application of the finite element method for numerical solution of the Navier-Stokes equation are discussed. The unsteady flow around a stationary and oscillating NACA 0012 airfoil was analyzed. In an effort to establish the range of applicability of the developed numerical method, the accuracy and stability problems in representing unsteady flow patterns and determining pressure distribution around the airfoil were investigated. Error involved in determining the velocity field and in the numerical integration of the discretized equations were analyzed in terms of geometry of the finite element gridwork, boundary conditions for stationary and oscillating airfoils, angles of attack and flow conditions.

Bratanow, T.↗

High-Lift Separated Flow About Airfoils

TRANSEP Calculates flow field about low-speed single-element airfoil at high-angle-of-attack and high-lift conditions with massive boundary-layer separation. TRANSEP includes effects of weak viscous interactions and can be used for subsonic/transonic airfoil design and analysis. The approach used in TRANSEP is based on direct-inverse method and its ability to use either displacement surface or pressure as airfoil boundary condition.

Carlson, L. A.↗

A feasibility study on the numerical modeling of interior noise fields

New developments in electronics and computer control have made active methods of noise suppression feasible. To date, however, active methods of noise control have been applied principally in geometrically uniform spaces with localized sound sources. The purpose of this paper is to demonstrate how finite element techniques may be used to model sound fields and active control inside geometrically complex spaces with non-uniform acoustic boundaries and multiple spatially separated sources. A closed-form solution is derived for the case of a one-dimensional tube with a vibrating piston at each end. Results are compared with finite element solutions and the problem is then generalized to two-dimensions with non-uniform geometry and non-uniform acoustic boundaries.

Abrahamson, A. L.↗

Member vibration effects on LSS behavior

This paper evaluates the sensitivity of modal characteristics of large deployable lattice-type space structures to joint boundary conditions. The evaluation of joint rotational boundary conditions is accomplished by a review of modal strain energy among elements of the analytic dynamic model. This review shows that space structure designs having no bending energy in lattice members in the low frequency range are insensitive to joint rotational boundary conditions in that frequency range. Evaluation of joint translation boundary conditions is accomplished by the 'modal freeplay method' which is described in this paper. The modal freeplay method relates modal frequency to joint translational freeplay and to amplitude of applied forces.

White, Charles W.↗

Finite element developments for two dimensional multiple-interface phase change problems

Finite element developments for multiple phase change problems in two-dimensional models are presented for the first time. The enthalpy method is used to simulate latent heat release in conjunction with fixed grid techniques. An unconditionally stable implicit method is used for the time integration. The effects of boundary conditions and the different phase regions on the multiple phase front developments are examined for numerous examples. Discussions and conclusions are appropriately addressed.

Ouyang, Tianhong↗

Impedance Eduction in a Duct Using the Linearized Euler Equations

An impedance eduction technique that is based upon a numerical solution to the linearized Euler equations coupled with the minimization of an objective function is presented. The linearized Euler equations are solved numerically using a cubic finite element method and the minimum of the objective function is obtained via a gradient based optimizer. The mean flow is allowed to have a gradient in two coordinate directions and the impedance of the liner may vary arbitrarily over the liner surface. A new regularized wall impedance boundary condition presented by Rienstra which contains the effects of the mean boundary layer thickness is implemented. Impedances are educed for a perforate over honeycomb liner using test data acquired in the Langley Grazing Flow Impedance Tube. The impedance of the liner is educed both with an upstream source and a downstream source. The method is validated by comparing the educed impedance to a benchmark method. The primary conclusion of the study is that the impedances educed for upstream and downstream sources are not well matched and that the primary effect of the boundary thickness is to reduce the resistance of the liner.

Watson, Willie R.↗

Solidification of II-VI Compounds in a Rotating Magnetic Field

This project is aimed at using a rotating magnetic field (RMF) to control fluid flow and transport during directional solidification of elemental and compound melts. Microgravity experiments have demonstrated that small amounts of residual acceleration of less than a micro-g can initiate and prolong fluid flow, particularly when there is a static component of the field perpendicular to the liquid solid interface. Thus a true diffusion boundary layer is not formed, and it becomes difficult to verify theories of solidification or to achieve diffusion controlled solidification. The RMF superimposes a stirring effect on an electrically conducting liquid, and with appropriate field strengths and frequencies, controlled transport of material through a liquid column can be obtained. As diffusion conditions are precluded and complete mixing conditions prevail, the technique is appropriate for traveling solvent zone or float zone growth methods in which the overall composition of the liquid can be maintained throughout the growth experiment. Crystals grown by RMF techniques in microgravity in previous, unrelated missions have shown exceptional properties. The objective of the project is two-fold, namely (1) using numerical modeling to simulate the behavior of a solvent zone with applied thermal boundary conditions and demonstrate the effects of decreasing gravity levels, or an increasing applied RMF, or both, and (2) to grow elements and II-VI compounds from traveling solvent zones both with and without applied RMFs, and to determine objectively how well the modeling predicts solidification parameters. Numerical modeling has demonstrated that, in the growth of CdTe from a tellurium solution, a rotating magnetic field can advantageously modify the shape of the liquid solid interface such that the interface is convex as seen from the liquid. Under such circumstances, the defect structure is reduced as any defects which are formed tend to grow out and not propagate. The flow of liquid, however, is complex due to the competing flow induced by the rotating magnetic field and the buoyancy driven convection. When the acceleration forces are reduced to one thousandth of gravity, the flow pattern is much simplified and well controlled material transport through the solvent zone can be readily achieved. Triple axis diffractometry and x-ray synchrotron topography have demonstrated that there is no significant improvement in crystal quality for HgCdTe grown on earth from a tellurium solution when a rotating magnetic field is applied. However, modeling shows that the flow in microgravity with a rotating magnetic field would produce a superior product.

Gillies, D. C.↗

Flow prediction over a transport multi-element high-lift system and comparison with flight measurements

Accurate prediction of surface-pressure distributions, merging boundary-layers, and separated-flow regions over multi-element high-lift airfoils is required to design advanced high-lift systems for efficient subsonic transport aircraft. The availability of detailed measurements of pressure distributions and both averaged and time-dependent boundary-layer flow parameters at flight Reynolds numbers is critical to evaluate computational methods and to model the turbulence structure for closure of the flow equations. Several detailed wind-tunnel measurements at subscale Reynolds numbers were conducted to obtain detailed flow information including the Reynolds-stress component. As part of a subsonic-transport high-lift research program, flight experiments are conducted using the NASA-Langley B737-100 research aircraft to obtain detailed flow characteristics for support of computational and wind-tunnel efforts. Planned flight measurements include pressure distributions at several spanwise locations, boundary-layer transition and separation locations, surface skin friction, as well as boundary-layer profiles and Reynolds stresses in adverse pressure-gradient flow.

Vijgen, P. M. H. W.↗