ALBERT

Description: As part of the continuing effort at NASA LeRC to improve both the durability and reliability of hot section Earth-to-orbit engine components, significant enhancements must be made in existing finite element and finite difference methods, and advanced techniques, such as the boundary element method (BEM), must be explored. The BEM was chosen as the basic analysis tool because the critical variables (temperature, flux, displacement, and traction) can be very precisely determined with a boundary-based discretization scheme. Additionally, model preparation is considerably simplified compared to the more familiar domain-based methods. Furthermore, the hyperbolic character of high speed flow is captured through the use of an analytical fundamental solution, eliminating the dependence of the solution on the discretization pattern. The price that must be paid in order to realize these advantages is that any BEM formulation requires a considerable amount of analytical work, which is typically absent in the other numerical methods. All of the research accomplishments of a multi-year program aimed toward the development of a boundary element formulation for the study of hot fluid-structure interaction in Earth-to-orbit engine hot section components are detailed. Most of the effort was directed toward the examination of fluid flow, since BEM's for fluids are at a much less developed state. However, significant strides were made, not only in the analysis of thermoviscous fluids, but also in the solution of the fluid-structure interaction problem.

Description: The development of a comprehensive fluid-structure interaction capability within a boundary element computer code is described. This new capability is implemented in a completely general manner, so that quite arbitrary geometry, material properties and boundary conditions may be specified. Thus, a single analysis code can be used to run structures-only problems, fluids-only problems, or the combined fluid-structure problem. In all three cases, steady or transient conditions can be selected, with or without thermal effects. Nonlinear analyses can be solved via direct iteration or by employing a modified Newton-Raphson approach. A number of detailed numerical examples are included at the end of these two sections to validate the formulations and to emphasize both the accuracy and generality of the computer code. A brief review of the recent applicable boundary element literature is included for completeness. The fluid-structure interaction facility is discussed. Once again, several examples are provided to highlight this unique capability. A collection of potential boundary element applications that have been uncovered as a result of work related to the present grant is given. For most of those problems, satisfactory analysis techniques do not currently exist.

Description: BEST-CMS (boundary element solution technology - composite modeling system) is an advanced engineering system for the micro-analysis of fiber composite structures. BEST-CMS is based upon the boundary element program BEST3D which was developed for NASA by Pratt and Whitney Aircraft and the State University of New York at Buffalo under contract NAS3-23697. BEST-CMS presently has the capabilities for elastostatic analysis, steady-state and transient heat transfer analysis, steady-state and transient concurrent thermoelastic analysis, and elastoplastic and creep analysis. The fibers are assumed to be perfectly bonded to the composite matrix, or in the case of static or steady-state analysis, the fibers may be assumed to have spring connections, thermal resistance, and/or frictional sliding between the fibers and the composite matrix. The primary objective of this user's manual is to provide an overview of all BEST-CMS capabilities, along with detailed descriptions of the input data requirements. In the next chapter, a brief review of the theoretical background is presented for each analysis category. Then, chapter three discusses the key aspects of the numerical implementation, while chapter four provides a tutorial for the beginning BEST-CMS user. The heart of the manual, however, is in chapter five, where a complete description of all data input items is provided. Within this chapter, the individual entries are grouped on a functional basis for a more coherent presentation. Chapter six includes sample problems and should be of considerable assistance to the novice. Chapter seven includes capsules of a number of fiber-composite analysis problems that have been solved using BEST-CMS. This chapter is primarily descriptive in nature and is intended merely to illustrate the level of analysis that is possible within the present BEST-CMS system. Chapter eight contains a detail description of the BEST-CMS Neutral File which is helpful in writing an interface between BEST-CMS and any graphic post-processor program. Finally, all pertinent references are listed in chapter nine.

Description: Over the last five years BEM has been developed into a very effective tool for micromechanical studies of composites. The analyses developed range from static elastic analysis, steady state heat transfer, steady state and transient thermoelastic analysis, periodic and transient dynamic analysis. Currently some of these analyses are being extended to nonlinear micromechanical analyses involving viscoplasticity, elastoplasticity, combined plasticity fracture and creep. The paper describes some of the results for elastic analysis, heat transfer and thermoelastic analyses.

Description: A comprehensive boundary element method is presented for transient thermoelastic analysis of hot section Earth-to-Orbit engine components. This time-domain formulation requires discretization of only the surface of the component, and thus provides an attractive alternative to finite element analysis for this class of problems. In addition, steep thermal gradients, which often occur near the surface, can be captured more readily since with a boundary element approach there are no shape functions to constrain the solution in the direction normal to the surface. For example, the circular disc analysis indicates the high level of accuracy that can be obtained. In fact, on the basis of reduced modeling effort and improved accuracy, it appears that the present boundary element method should be the preferred approach for general problems of transient thermoelasticity.

Description: Details on the progress made during the first three years of a five-year program towards the development of a boundary element code are presented. This code was designed for the micromechanical studies of advance ceramic composites. Additional effort was made in generalizing the implementation to allow the program to be applicable to real problems in the aerospace industry. The ceramic composite formulations developed were implemented in the three-dimensional boundary element computer code BEST3D. BEST3D was adopted as the base for the ceramic composite program, so that many of the enhanced features of this general purpose boundary element code could by utilized. Some of these facilities include sophisticated numerical integration, the capability of local definition of boundary conditions, and the use of quadratic shape functions for modeling geometry and field variables on the boundary. The multi-region implementation permits a body to be modeled in substructural parts; thus dramatically reducing the cost of the analysis. Furthermore, it allows a body consisting of regions of different ceramic matrices and inserts to be studied.

Description: Different algorithms for the geometric connected component labeling (GCCL) problem are defined each of which involves d stages of message passing, for a d-dimensional hypercube. The major idea is that in each stage a hypercube multiprocessor increases its knowledge of domain. The algorithms under consideration include the QUAD algorithm for small number of processors and the Overlap Quad algorithm for large number of processors, subject to the locality of the connected sets. These algorithms differ in their run time, memory requirements, and message complexity. They were implemented on an Intel iPSC2/D4/MX hypercube.

Description: It is evident that for proper micromechanical analysis of ceramic composites, one needs to use a numerical method that is capable of idealizing the individual fibers or individual bundles of fibers embedded within a three-dimensional ceramic matrix. The analysis must be able to account for high stress or temperature gradients from diffusion of stress or temperature from the fiber to the ceramic matrix and allow for interaction between the fibers through the ceramic matrix. The analysis must be sophisticated enough to deal with the failure of fibers described by a series of increasingly sophisticated constitutive models. Finally, the analysis must deal with micromechanical modeling of the composite under nonlinear thermal and dynamic loading. This report details progress made towards the development of a boundary element code designed for the micromechanical studies of an advanced ceramic composite. Additional effort has been made in generalizing the implementation to allow the program to be applicable to real problems in the aerospace industry.

Description: The boundary element method is applied to transient viscous incompressible flow. The time-domain formulation allows a boundary-only solution for linear Stokes flow. For higher speed flows in which the nonlinear convective effects cannot be ignored, a volume integral must be retained. However, the introduction of reference velocities often limits the nonlinear region to the vicinity of obstacles or boundary surfaces. Additionally, the volume terms are rewritten to eliminate the need for the calculation of velocity gradients. A general purpose numerical implementation of this new formulation then produces a very attractive tool for engineering analysis. This implementation includes a Newton-Raphson algorithm, permitting accurate solutions up to the moderate Reynolds number range. Several numerical examples are provided to validate the present approach.

Description: An advanced level development of the boundary element method is presented for the elastic stress analysis of a three-dimensional solid containing a large number of small diameter, tubular shaped holes. The formulation has been developed such that these holes can be modeled by a system of curvilinear line elements, resulting in substantial savings in both data preparation and computing costs. This is accomplished by assuming a variation in the displacement field along the circumference defined in terms of trigonometric functions together with a linear or quadratic variation of displacements along the longitudinal direction.