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 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Calculating Transient Vibrations Of Coupled Substructures

Finite-element, numerical-integration method for estimating transient vibrational response of structure composed of coupled substructures entails less computation. Responses of substructures to external forces and to forces of interaction between substructures computed. Equations of motion of each substructure solved independently. Applicable to any number and configuration of linearly responding substructures and to both determinate and indeterminate boundary conditions at interfaces between substructures. Also applicable to changing interface boundary conditions.

Admire, J. R.↗

NASA NTRS1992

Stress analysis and damage evaluation of flawed composite laminates by hybrid-numerical methods

Structural components in flight vehicles is often inherited flaws, such as microcracks, voids, holes, and delamination. These defects will degrade structures the same as that due to damages in service, such as impact, corrosion, and erosion. It is very important to know how a structural component can be useful and survive after these flaws and damages. To understand the behavior and limitation of these structural components researchers usually do experimental tests or theoretical analyses on structures with simulated flaws. However, neither approach has been completely successful. As Durelli states that 'Seldom does one method give a complete solution, with the most efficiency'. Examples of this principle is seen in photomechanics which additional strain-gage testing can only average stresses at locations of high concentration. On the other hand, theoretical analyses including numerical analyses are implemented with simplified assumptions which may not reflect actual boundary conditions. Hybrid-Numerical methods which combine photomechanics and numerical analysis have been used to correct this inefficiency since 1950's. But its application is limited until 1970's when modern computer codes became available. In recent years, researchers have enhanced the data obtained from photoelasticity, laser speckle, holography and moire' interferometry for input of finite element analysis on metals. Nevertheless, there is only few of literature being done on composite laminates. Therefore, this research is dedicated to this highly anisotropic material.

Yang, Yii-Ching↗

NASA NTRS1994

Dynamically consistent hydrography and absolute velocity in the eastern North Atlantic Ocean

The problem of mapping a dynamically consistent hydrographic field and associated absolute geostrophic flow in the eastern North Atlantic between 24 deg and 36 deg N is related directly to the solution of the so-called thermocline equations. A nonlinear optimization problem involving Needler's P equation is solved to find the hydrography and resulting flow that minimizes the vertical mixing above about 1500 m in the ocean and is simultaneously consistent with the observations. A sharp minimum (at least in some dimensions) is found, apparently corresponding to a solution nearly conserving potential vorticity and with vertical eddy coefficient less than about 10(exp -5) sq m/s. Estimates of `residual' quantities such as eddy coefficients are extremely sensitive to slight modifications to the observed fields. Boundary conditions, vertical velocities, etc., are a product of the optimization and produce estimates differing quantitatively from prior ones relying directly upon observed hydrography. The results are generally insensitive to particular elements of the solution methodology, but many questions remain concerning the extent to which different synoptic sections can be asserted to represent the same ocean. The method can be regarded as a practical generalization of the beta spiral and geostrophic balance inverses for the estimate of absolute geostrophic flows. Numerous improvements to the methodology used in this preliminary attempt are possible.

Wunsch, Carl↗

NASA NTRS1969

Interdisciplinary Approach to the Lubrication of Concentrated Contacts

In an effort to further the advances in the complex and multidisciplinary subject of lubrication, the National Aeronautics and Space Administration (NASA) inaugurated in 1967 an interdisciplinary lubrication symposium series. The first meeting of the series, held in San Antonio, Texas, on November 28-30, 1967, was devoted to a critical appraisal of the level of understanding and the needed research in the area of sliding friction and wear under unlubricated and boundary-lubrication conditions. The second meeting, held in Cleveland, Ohio, on November 19-21, 1968, again emphasized unlubricated and boundary-lubricated sliding friction and wear but attempted to focus attention on some selected theoretical and practical problems by discussing in small working group their implications and possible methods of attack. Having thus provided for an examination of sliding friction and wear under unlubricated and boundary-lubrication conditions in considerable breadth and depth, the Steering Committee recommended to ASA in September 196 that further attention be turned to an evaluation of the collective understanding and the needed research in the lubrication of highly loaded machine elements in a combined rolling and sliding situation. The recommendation was accepted by ASA and, accordingly, a symposium entitled Interdisciplinary Approach to the Lubrication of Concentrated Contacts was held on the p-rounds of Rensselaer Polytechnic Institute, Troy, New York, on July 15-17, 1969.

P M Ku↗

NASA NTRS1977

A three-dimensional viscous/potential flow interaction analysis method for multi-element wings

An analysis method and computer program were developed for the calculation of the viscosity dependent aerodynamic characteristics of multi-element, finite wings in incompressible flow. A fully-three dimensional potential flow program is used to determine the inviscid pressure distribution about the configuration. The potential flow program uses surface source and vortex singularities to represent the inviscid flow. The method is capable of analysing configurations having at most one slat, a main element, and two slotted flaps. Configurations are limited to full span slats or flaps. The configuration wake is allowed to relax as a force free wake, although roll up is not allowed at this time. Once the inviscid pressure distribution is calculated, a series of boundary layer computations are made along streamwise strips.

Dvorak, F. A.↗

NASA NTRS1979

An evaluation of four single element airfoil analytic methods

A comparison of four computer codes for the analysis of two-dimensional single element airfoil sections is presented for three classes of section geometries. Two of the computer codes utilize vortex singularities methods to obtain the potential flow solution. The other two codes solve the full inviscid potential flow equation using finite differencing techniques, allowing results to be obtained for transonic flow about an airfoil including weak shocks. Each program incorporates boundary layer routines for computing the boundary layer displacement thickness and boundary layer effects on aerodynamic coefficients. Computational results are given for a symmetrical section represented by an NACA 0012 profile, a conventional section illustrated by an NACA 65A413 profile, and a supercritical type section for general aviation applications typified by a NASA LS(1)-0413 section. The four codes are compared and contrasted in the areas of method of approach, range of applicability, agreement among each other and with experiment, individual advantages and disadvantages, computer run times and memory requirements, and operational idiosyncrasies.

Freuler, R. J.↗

NASA NTRS1980

Effect of grid system on finite element calculation

Detailed parametric studies of the effect of grid system on finite element calculation for potential flows were made. These studies led to the formulation of a design criteria for optimum mesh system and the development of two methods to generate the optimum mesh system. The guidelines for optimum mesh system are: (1) the mesh structure should be regular; (2) the element should be as regular and equilateral as possible; (3) the distribution of size of element should be consistent with that of flow variables to insure maximum uniformity in error distribution; (4) for non-Dirichlet boundary conditions, smaller boundary elements or higher order interpolation functions should be used; and (5) the mesh should accommodate the boundary geometry as accurately as possible. The results of the parametric studies are presented.

Lee, K. D.↗

NASA NTRS1984

Newton like: Minimal residual methods applied to transonic flow calculations

A computational technique for the solution of the full potential equation is presented. The method consists of outer and inner iterations. The outer iterate is based on a Newton like algorithm, and a preconditioned Minimal Residual method is used to seek an approximate solution of the system of linear equations arising at each inner iterate. The present iterative scheme is formulated so that the uncertainties and difficulties associated with many iterative techniques, namely the requirements of acceleration parameters and the treatment of additional boundary conditions for the intermediate variables, are eliminated. Numerical experiments based on the new method for transonic potential flows around the NACA 0012 airfoil at different Mach numbers and different angles of attack are presented, and these results are compared with those obtained by the Approximate Factorization technique. Extention to three dimensional flow calculations and application in finite element methods for fluid dynamics problems by the present method are also discussed. The Inexact Newton like method produces a smoother reduction in the residual norm, and the number of supersonic points and circulations are rapidly established as the number of iterations is increased.

Wong, Y. S.↗

NASA NTRS1985

Newton-like minimal residual methods applied to transonic flow calculations

Wong, Y. S.↗

NASA NTRS1991

Least-squares solution of incompressible Navier-Stokes equations with the p-version of finite elements

A p-version of the least squares finite element method, based on the velocity-pressure-vorticity formulation, is developed for solving steady state incompressible viscous flow problems. The resulting system of symmetric and positive definite linear equations can be solved satisfactorily with the conjugate gradient method. In conjunction with the use of rapid operator application which avoids the formation of either element of global matrices, it is possible to achieve a highly compact and efficient solution scheme for the incompressible Navier-Stokes equations. Numerical results are presented for two-dimensional flow over a backward facing step. The effectiveness of simple outflow boundary conditions is also demonstrated.

Jiang, Bo-Nan↗

NASA NTRS1993

Analytical evaluation of a surface integral expressing the coupling between interior and exterior volumes in a FE-IE approach

Recently an approach which combines the finite element technique and an integral equation to determine the fields scattered by inhomogeneous bodies of complicated shape has been proposed. Basically, a mathematical surface which encloses the scatterers is introduced, thus dividing the space into an interior and an exterior volume, in which the finite element technique and an integral equation for EM scattering, respectively, are applied. The integral equation is set up for the tangential components of the fields at the surface, while the interior volume the unknowns are the total fields. Continuity of the tangential fields at the boundary, as required by Maxwell's equations, is imposed, thus coupling the two methods to obtain a consistent solution. The coupling term is expressed by a surface integral formed by the dot product of a FE basis function and an IE testing function, or viceversa. By choosing the boundary to be a surface of revolution and by making a convenient selection of IE basis (testing) functions, it is possible to evaluate the integrals analytically on surfaces such as curved triangles, curved quadrilaterals and curved pentagons. We will illustrate the salient steps involved in setting up and carrying out these integrals and discuss what class of basis (testing) functions and analytic surfaces of revolution they are applicable to. Analytic calculations offer the advantage of better accuracy than purely numerical ones, and, when combined with them, often shed light on issues of numerical convergence and limiting values. Furthermore, they may reduce computation time and storage requirements.

Zuffada, C.↗

NASA NTRS1995

Sensitivity analysis of frictional contact response of axisymmetric composite structures

A computational procedure is presented for evaluating the sensitivity coefficients of the static frictional contact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of nodal displacements, and Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for the frictional contact of a composite spherical cap pressed against a rigid plate.

Karaoglan, I.↗

NASA NTRS1996

Three-Dimensional Simulations of Marangoni-Benard Convection in Small Containers by the Least-Squares Finite Element Method

This paper reports a numerical study of the Marangoni-Benard (MB) convection in a planar fluid layer. The least-squares finite element method (LSFEM) is employed to solve the three-dimensional Stokes equations and the energy equation. First, the governing equations are reduced to be first-order by introducing variables such as vorticity and heat fluxes. The resultant first-order system is then cast into a div-curl-grad formulation, and its ellipticity and permissible boundary conditions are readily proved. This numerical approach provides an equal-order discretization for velocity, pressure, vorticity, temperature, and heat conduction fluxes, and therefore can provide high fidelity solutions for the complex flow physics of the MB convection. Numerical results reported include the critical Marangoni numbers (M(sub ac)) for the onset of the convection in containers with various aspect ratios, and the planforms of supercritical MB flows. The numerical solutions compared favorably with the experimental results reported by Koschmieder et al..

Yu, Sheng-Tao↗

NASA NTRS1971

Borehole probes: Summary of previous work

Previous research explored the applicability of some devices to lunar borehold testing and, in particular, examined the use of such instruments to assess the strength of soils and rocks. Special attention was given to borehole jacks, which essentially conduct a plate bearing test across opposing wall areas. Analytical work employed finite element analysis, the theory of plasticity, and the theory of elasticity using the complex variable method. A number of solutions were reached for tractable subproblems in the set of problems posed when a jack is used to expand sectors of a borehole. In addition, an investigation with physical models was made to explore the modes of behavior for varying boundary conditions and materials, including both elastic and plastic media.

Goodman, R. E.↗

NASA NTRS1974

Sensitivity analysis of torsional vibration characteristics of helicopter rotor blades. Part 1: Structural dynamics analysis

A theoretical investigation of structural vibration characteristics of rotor blades was carried out. Coupled equations of motion for flapwise bending and torsion were formulated for rotor blades with noncollinear elastic and mass axes. The finite element method was applied for a detailed representation of blade structural properties. Coupled structural mass and stiffness coefficients were evaluated. The range of validity of a set of coupled equations of motion linearized with respect to eccentricity between elastic and mass axes was investigated. The sensitivity of blade vibration characteristics to torsion were evaluated by varying blade geometric properties, boundary conditions, and eccentricities between mass and elastic axes.

Bratanow, T.↗

NASA NTRS1979

Upgraded viscous flow analysis of multi-element airfoils

A description of an improved version of the NASA/Lockheed multi-element airfoil analysis computer program is presented. The improvements include several major modifications of the aerodynamic model as well as substantial changes of the computer code. The modifications of the aerodynamic model comprise the representation of the boundary layer and wake displacement effects with an equivalent source distribution, the prediction of wake parameters with Green's lag-entrainment method, the calculation of turbulent boundary layer separation with the method of Nash and Hicks, the estimation of the onset of confluent boundary layer separation with a modified form of Goradia's method, and the prediction of profile drag with the formula of Squire and Young. The modifications of the computer program for which the structured approach to computer software development was employed are also described. Important aspects of the structured program development such as the functional decomposition of the aerodynamic theory and its numerical implementation, the analysis of the data flow within the code, and the application of a pseudo code are discussed.

Brune, G. W.↗

NASA NTRS1985

Probabilistic Finite Element Development

The probabilistic finite element computer program known as Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) is being developed for the analysis of critical structural components for reusable space propulsion systems. First year efforts involve the formulation of the probabilistic analysis strategy and the development of a probabilistic linear analysis code. The ultimate goal of the 3-year program is the development of a finite element code capable of performing nonlinear dynamic analysis of structures having stochastic material properties, geometry, and boundary conditions and subjected to random loading. Three levels of sophistication are envisioned for the stochastic description of the structural problem, namely: (1) homogeneous random variable for stiffness, mass, damping, and external loading; (2) stochastic characterization of variables at the element level, with specified interelement correlations; and (3) stochastic interpolation of variables within a finite element. Two alternative probabilistic analysis methods will be developed, allowing for all three levels of modeling sophistication.

Nagtegaal, J.↗

NASA NTRS1987

Finite element techniques for the Navier-Stokes equations in the primitive variable formulation and the vorticity stream-function formulation

Finite element procedures for the Navier-Stokes equations in the primitive variable formulation and the vorticity stream-function formulation have been implemented. For both formulations, streamline-upwind/Petrov-Galerkin techniques are used for the discretization of the transport equations. The main problem associated with the vorticity stream-function formulation is the lack of boundary conditions for vorticity at solid surfaces. Here an implicit treatment of the vorticity at no-slip boundaries is incorporated in a predictor-multicorrector time integration scheme. For the primitive variable formulation, mixed finite-element approximations are used. A nine-node element and a four-node + bubble element have been implemented. The latter is shown to exhibit a checkerboard pressure mode and a numerical treatment for this spurious pressure mode is proposed. The two methods are compared from the points of view of simulating internal and external flows and the possibilities of extensions to three dimensions.

Glaisner, F.↗