Search NASASearch

SEARCH · Search NASA

Results for “BEM”

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 19 records

Development of an integrated BEM approach for hot fluid structure interaction: BEST-FSI: Boundary Element Solution Technique for Fluid Structure Interaction

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.

Dargush, G. F.

Time-domain transient elastodynamic analysis of 3-D solids by BEM

The BEM algorithm developed by Banerjee and Ahmad (1985) and Banerjee et al. (1986) for the analysis of three-dimensional time-dependent elastodynamics problems is extended and refined to improve the accuracy of the internal stress predictions. A quadratic spatial shape function and a linear temporal shape function are employed, and the complete numerical implementation (as part of BEST3D) is outlined. Applications of the improved BEM to typical problems are presented graphically and briefly characterized. Included are spherical cavities subjected to sudden radial expansion, triangular or rectangular pulses of radial pressure, or ramp loading; a circular loaded area on a half space; and a flexible square plate foundation on an elastic half space.

Ahmad, S.

A new BEM formulation for two- and three-dimensional elastoplasticity using particular integrals

New two- and three-dimensional BEM formulations are developed for elastoplastic stress analysis. These new procedures differ from previous work in that volume integration is not required to incorporate the nonlinear effects in the analysis. Instead, initial stress rates are introduced in the BEM system via particular integrals. The present formulation is implemented in a general-purpose multi-region system, and examples are presented to demonstrate the accuracy and versatility of the method.

Henry, Donald P., Jr.

Development of an integrated BEM approach for hot fluid structure interaction

In the present work, the boundary element method (BEM) is chosen as the basic analysis tool, principally because the definition of temperature, flux, displacement and traction are very precise on a boundary-based discretization scheme. One fundamental difficulty is, of course, that a BEM formulation requires a considerable amount of analytical work, which is not needed in the other numerical methods. Progress made toward the development of a boundary element formulation for the study of hot fluid-structure interaction in Earth-to-Orbit engine hot section components is reported. The primary thrust of the program to date has been directed quite naturally toward the examination of fluid flow, since boundary element methods for fluids are at a much less developed state.

Dargush, Gary F.

Development of an integrated BEM for hot fluid-structure interaction

One of the most difficult problems in engine structural component durability analysis is the determination of the temperatures and fluxes in the structural components directly in contact with the hot gas flow path. Currently there exists no rational analytical or numerical technique which can effectively deal with this problem. Since the temperature distribution in the structural components are strongly influenced by both the fluid flow and the deformation as well as the cooling system in the structure, the only effective way to deal with this problem is to develop an integrated solid mechanics, fluid mechanics and heat transfer analysis for this problem. Herein, the Boundary Element Method (BEM) is chosen as the basic analysis tool principally because the definition of quantities like fluxes, temperatures, displacements, and velocities are very precise on a boundary based discretization scheme. One fundamental difficulty is that a BEM analysis requires a considerable amount of analytical work which is not present in other numerical methods. During the past year, all of this analytical work was completed and a two dimensional, general purpose code was written. A portion of the work is summarized.

Dargush, G. F.

Unsteady three-dimensional thermal field prediction in turbine blades using nonlinear BEM

A time-and-space accurate and computationally efficient fully three dimensional unsteady temperature field analysis computer code has been developed for truly arbitrary configurations. It uses boundary element method (BEM) formulation based on an unsteady Green's function approach, multi-point Gaussian quadrature spatial integration on each panel, and a highly clustered time-step integration. The code accepts either temperatures or heat fluxes as boundary conditions that can vary in time on a point-by-point basis. Comparisons of the BEM numerical results and known analytical unsteady results for simple shapes demonstrate very high accuracy and reliability of the algorithm. An example of computed three dimensional temperature and heat flux fields in a realistically shaped internally cooled turbine blade is also discussed.

Martin, Thomas J.

Electromagnetic Scattering Analysis of Arbitrarily Shaped Material Cylinder by FEM-BEM Method

A hybrid method that combines the finite element method (FEM) and the boundary element method (BEM) is developed to analyze electromagnetic scattering from arbitrarily shaped material cylinders. By this method, the material cylinder is first enclosed by a fictitious boundary. Maxwell's equations are then solved by FEM inside and by BEM outside the boundary. Electromagnetic scattering from several arbitrarily shaped material cylinders is computed and compared with results obtained by other numerical techniques.

Deshpande, M. D.

User's Manual for FEM-BEM Method

A user's manual for using FORTRAN code to perform electromagnetic analysis of arbitrarily shaped material cylinders using a hybrid method that combines the finite element method (FEM) and the boundary element method (BEM). In this method, the material cylinder is enclosed by a fictitious boundary and the Maxwell's equations are solved by FEM inside the boundary and by BEM outside the boundary. The electromagnetic scattering on several arbitrarily shaped material cylinders using this FORTRAN code is computed to as examples.

Butler, Theresa

Issues and Methods Concerning the Evaluation of Hypersingular and Near-Hypersingular Integrals in BEM Formulations

It is known that higher order modeling of the sources and the geometry in Boundary Element Modeling (BEM) formulations is essential to highly efficient computational electromagnetics. However, in order to achieve the benefits of hIgher order basis and geometry modeling, the singular and near-singular terms arising in BEM formulations must be integrated accurately. In particular, the accurate integration of near-singular terms, which occur when observation points are near but not on source regions of the scattering object, has been considered one of the remaining limitations on the computational efficiency of integral equation methods. The method of singularity subtraction has been used extensively for the evaluation of singular and near-singular terms. Piecewise integration of the source terms in this manner, while manageable for bases of constant and linear orders, becomes unwieldy and prone to error for bases of higher order. Furthermore, we find that the singularity subtraction method is not conducive to object-oriented programming practices, particularly in the context of multiple operators. To extend the capabilities, accuracy, and maintainability of general-purpose codes, the subtraction method is being replaced in favor of the purely numerical quadrature schemes. These schemes employ singularity cancellation methods in which a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. An example of the sin,oularity cancellation approach is the Duffy method, which has two major drawbacks: 1) In the resulting integrand, it produces an angular variation about the singular point that becomes nearly-singular for observation points close to an edge of the parent element, and 2) it appears not to work well when applied to nearly-singular integrals. Recently, the authors have introduced the transformation u(x(prime))= sinh (exp -1) x(prime)/Square root of ((y prime (exp 2))+ z(exp 2) for integrating functions of the form I = Integral of (lambda(r(prime))((e(exp -jkR))/(4 pi R) d D where A (r (prime)) is a vector or scalar basis function and R = Square root of( (x(prime)(exp2) + (y(prime)(exp2) + z(exp 2)) is the distance between source and observation points. This scheme has all of the advantages of the Duffy method while avoiding the disadvantages listed above. In this presentation we will survey similar approaches for handling singular and near-singular terms for kernels with 1/R(exp 2) type behavior, addressing potential pitfalls and offering techniques to efficiently handle special cases.

Fink, P. W.

Development of an integrated BEM for hot fluid-structure interaction

The Boundary Element Method (BEM) is chosen as a basic analysis tool principally because the definition of quantities like fluxes, temperature, displacements, and velocities is very precise on a boundary base discretization scheme. One fundamental difficulty is, of course, that the entire analysis requires a very considerable amount of analytical work which is not present in other numerical methods. During the last 18 months all of this analytical work was completed and a two-dimensional, general purpose code was written. Some of the early results are described. It is anticipated that within the next two to three months almost all two-dimensional idealizations will be examined. It should be noted that the analytical work for the three-dimensional case has also been done and numerical implementation will begin next year.

Banerjee, P. K.

BEM for wave equation with boundary in arbitrary motion and applications to compressible potential aerodynamics of airplanes and helicopters

The wave equation for an object in arbitrary motion is investigated analytically using a BEM approach, and practical applications to potential flows of compressible fluids around aircraft wings and helicopter rotors are considered. The treatment accounts for arbitrary combined rotational and translational motion of the reference frame and for the wake motion. The numerical implementation as a computer algorithm is demonstrated on problems with prescribed and free wakes, the former in compressible flows and the latter for incompressible flows; results are presented graphically and briefly characterized.

Morino, Luigi

COMGEN-BEM: Boundary element model generation for composite materials micromechanical analysis

Composite Model Generation-Boundary Element Method (COMGEN-BEM) is a program developed in PATRAN command language (PCL) which generates boundary element models of continuous fiber composites at the micromechanical (constituent) scale. Based on the entry of a few simple parameters such as fiber volume fraction and fiber diameter, the model geometry and boundary element model are generated. In addition, various mesh densities, material properties, fiber orientation angles, loads, and boundary conditions can be specified. The generated model can then be translated to a format consistent with a boundary element analysis code such as BEST-CMS.

Goldberg, Robert K.

Applications of FEM and BEM in two-dimensional fracture mechanics problems

A comparison of the finite element method (FEM) and boundary element method (BEM) for the solution of two-dimensional plane strain problems in fracture mechanics is presented in this paper. Stress intensity factors (SIF's) were calculated using both methods for elastic plates with either a single-edge crack or an inclined-edge crack. In particular, two currently available programs, ANSYS for finite element analysis and BEASY for boundary element analysis, were used.

Min, J. B.

Micromechanical studies of composites by BEM

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.

Banerjee, P. K.

Refinement of Methods for Evaluation of Near-Hypersingular Integrals in BEM Formulations

In this paper, we present advances in singularity cancellation techniques applied to integrals in BEM formulations that are nearly hypersingular. Significant advances have been made recently in singularity cancellation techniques applied to 1 R type kernels [M. Khayat, D. Wilton, IEEE Trans. Antennas and Prop., 53, pp. 3180-3190, 2005], as well as to the gradients of these kernels [P. Fink, D. Wilton, and M. Khayat, Proc. ICEAA, pp. 861-864, Torino, Italy, 2005] on curved subdomains. In these approaches, the source triangle is divided into three tangent subtriangles with a common vertex at the normal projection of the observation point onto the source element or the extended surface containing it. The geometry of a typical tangent subtriangle and its local rectangular coordinate system with origin at the projected observation point is shown in Fig. 1. Whereas singularity cancellation techniques for 1 R type kernels are now nearing maturity, the efficient handling of near-hypersingular kernels still needs attention. For example, in the gradient reference above, techniques are presented for computing the normal component of the gradient relative to the plane containing the tangent subtriangle. These techniques, summarized in the transformations in Table 1, are applied at the sub-triangle level and correspond particularly to the case in which the normal projection of the observation point lies within the boundary of the source element. They are found to be highly efficient as z approaches zero. Here, we extend the approach to cover two instances not previously addressed. First, we consider the case in which the normal projection of the observation point lies external to the source element. For such cases, we find that simple modifications to the transformations of Table 1 permit significant savings in computational cost. Second, we present techniques that permit accurate computation of the tangential components of the gradient; i.e., tangent to the plane containing the source element.

Fink, Patricia W.

Development of an integrated BEM (Boundary Element Mesh) for hot fluid-structure interaction

A boundary integral representation for a coupled approach to fluid flow and solid deformation problems associated with the design of hot-section components such as those in the Space Shuttle Main Engine is discussed. The formulation is based on the fundamental analytical solution of the Navier-Stokes equation for fluid velocity in an infinite domain. This fundamental solution was obtained by decomposing a Navier-Stokes equation into vorticity and dilation transport equations. A boundary integral involving convolutions in time was then constructed in which the convective terms appear in the volume integral.

Dargush, G. F.

Free vibration analysis by BEM using particular integrals

A new method for the free-vibration analysis using the boundary element technique is presented. The method utilizes a fictitious vector function to approximate the inertia forces and then uses the well-known concept of complementary functions and particular integrals to solve the resulting governing differential equations. The necessary particular integrals are defined for the two and three-dimensional analyses, and the present formulation is applied to a number of two-dimensional problems to show its accuracy and efficiency in the solution of realistic engineering problems.

Ahmad, S.

Development of an integrated BEM approach for hot fluid structure interaction

Significant progress was made toward the goal of developing a general purpose boundary element method for hot fluid-structure interaction. For the solid phase, a boundary-only formulation was developed and implemented for uncoupled transient thermoelasticity in two dimensions. The elimination of volume discretization not only drastically reduces required modeling effort, but also permits unconstrained variation of the through-the-thickness temperature distribution. Meanwhile, for the fluids, fundamental solutions were derived for transient incompressible and compressible flow in the absence of the convective terms. Boundary element formulations were developed and described. For the incompressible case, the necessary kernal functions, under transient and steady-state conditions, were derived and fully implemented into a general purpose, multi-region boundary element code. Several examples were examined to study the suitability and convergence characteristics of the various algorithms.

Dargush, Gary F.