Search NASA⌕ Search

SEARCH · Search NASA

Results for “discrete 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Correction factory techniques for improving aerodynamic prediction methods

A method for correcting discrete element lifting surface theory to reflect given experimental data is presented. Theoretical pressures are modified such that imposed constraints are satisfied while minimizing the changes to the pressures. Several types of correction procedures are presented and correlated; (1) scaling of pressures; (2) scaling of downwash values; and (3) addition of an increment to the downwash that is proportioned to pressure. Some special features are included in these methods and they include: (1) consideration of experimental data from multiple deflection modes, (2) limitation of the amplitudes of the correction factors, and (3) the use of correction factor mode shapes. These methods are correlated for cases involving all three Mach Number ranges using a FORTRAN IV computer program. Subsonically, a wing with an oscillating partial span control surface and a wing with a leading edge droop are presented. Transonically a two-dimensional airfoil with an oscillating flap is considered. Supersonically an arrow wing with and without camber is analyzed. In addition to correction factor methods an investigation is presented dealing with a new simplified transonic modification of the two-dimensional subsonic lifting surface theory. Correlations are presented for an airfoil with an oscillating flap.

Giesing, J. P.↗
NASA NTRS1989

Nonlinear finite element formulation for the large displacement analysis in multibody system dynamics

A total Lagrangian finite element formulation for the deformable bodies in multibody mechanical systems that undergo finite relative rotations is developed. The deformable bodies are discretized using finite element methods. The shape functions that are used to describe the displacement field are required to include the rigid body modes that describe only large translational displacements. This does not impose any limitations on the technique because most commonly used shape functions satisfy this requirement. The configuration of an element is defined using four sets of coordinate systems: Body, Element, Intermediate element, Global. The body coordinate system serves as a unique standard for the assembly of the elements forming the deformable body. The element coordinate system is rigidly attached to the element and therefore it translates and rotates with the element. The intermediate element coordinate system, whose axes are initially parallel to the element axes, has an origin which is rigidly attached to the origin of the body coordinate system and is used to conveniently describe the configuration of the element in undeformed state with respect to the body coordinate system.

Rismantab-Sany, J.↗
NASA NTRS1989

Error analysis of finite element solutions for postbuckled cylinders

A general method of error analysis and correction is investigated for the discrete finite-element results for cylindrical shell structures. The method for error analysis is an adaptation of the method of successive approximation. When applied to the equilibrium equations of shell theory, successive approximations derive an approximate continuous solution from the discrete finite-element results. The advantage of this continuous solution is that it contains continuous partial derivatives of an order higher than the basis functions of the finite-element solution. Preliminary numerical results are presented in this paper for the error analysis of finite-element results for a postbuckled stiffened cylindrical panel modeled by a general purpose shell code. Numerical results from the method have previously been reported for postbuckled stiffened plates. A procedure for correcting the continuous approximate solution by Newton's method is outlined.

Sistla, Rajaram↗
NASA NTRS1984

Methods for analysis of cracks in three-dimensional solids

Analytical and numerical methods evaluating the stress-intensity factors for three-dimensional cracks in solids are presented, with reference to fatigue failure in aerospace structures. The exact solutions for embedded elliptical and circular cracks in infinite solids, and the approximate methods, including the finite-element, the boundary-integral equation, the line-spring models, and the mixed methods are discussed. Among the mixed methods, the superposition of analytical and finite element methods, the stress-difference, the discretization-error, the alternating, and the finite element-alternating methods are reviewed. Comparison of the stress-intensity factor solutions for some three-dimensional crack configurations showed good agreement. Thus, the choice of a particular method in evaluating the stress-intensity factor is limited only to the availability of resources and computer programs.

Raju, I. S.↗
NASA NTRS1990

Scattering of impact wave by a crack in composite plate

The surface responses due to impact load on an infinite uniaxial graphite/epoxy plate with and without delamination cracks are investigated both in time and frequency domain by using a hybrid method combining the finite element discretization of the near-field with boundary integral representation of the field outside a contour enclosing completely the crack. This combined method leads to a set of linear unsymmetric complex matrix equations which are solved to obtain the response in the frequency domain by biconjugate gradient method. The time domain response is then obtained by using an FFT. In order to capture the time-domain characteristics accurately, high order finite elements have been used. Also, both the six node singular elements and eight node transition elements are used around the crack tips to model the crack-tip singularity. It is shown that from the numerical results for surface responses both depth and length of this crack can be identified.

Ju, T. H.↗
NASA NTRS1989

Numerical algorithms for steady and unsteady incompressible Navier-Stokes equations

The numerical analysis of the incompressible Navier-Stokes equations are becoming important tools in the understanding of some fluid flow problems which are encountered in research as well as in industry. With the advent of the supercomputers, more realistic problems can be studied with a wider choice of numerical algorithms. An alternative formulation is presented for viscous incompressible flows. The incompressible Navier-Stokes equations are cast in a velocity/vorticity formulation. This formulation consists of solving the Poisson equations for the velocity components and the vorticity transport equation. Two numerical algorithms for the steady two-dimensional laminar flows are presented. The first method is based on the actual partial differential equations. This uses a finite-difference approximation of the governing equations on a staggered grid. The second method uses a finite element discretization with the vorticity transport equation approximated using a Galerkin approximation and the Poisson equations are obtained using a least squares method. The equations are solved efficiently using Newton's method and a banded direct matrix solver (LINPACK). The method is extended to steady three-dimensional laminar flows and applied to a cubic driven cavity using finite difference schemes and a staggered grid arrangement on a Cartesian mesh. The equations are solved iteratively using a plane zebra relaxation scheme. Currently, a two-dimensional, unsteady algorithm is being developed using a generalized coordinate system. The equations are discretized using a finite-volume approach. This work will then be extended to three-dimensional flows.

Hafez, Mohammed↗
NASA NTRS1981

A Galerkin type finite element method for rotary-wing aeroelasticity in hover and forward flight

A Galerkin finite element method for the spatial discretization of the nonlinear, nonselfadjoint, partial differential equations governing rotary-wing aeroelasticity is presented. This method reduces algebraic manipulative labor significantly when compared to the global Galerkin method based on assumed modes. Furthermore, the Galerkin finite element method is ideally suited to treat rotor blades with discontinuous mass and stiffness distribution and structurally redundant configurations as they appear in bearingless rotors. Implementation of the method is illustrated for the coupled flap-lag aeroelastic problem of hingeless rotor blades in hover and forward flight. Numerical results for stability and response illustrate the numerical properties and convergence behavior of the method. It is concluded that the Galerkin finite element method is a practical tool for solving rotary-wing aeroelastic stability and response problems.

Straub, F. K.↗
NASA NTRS1985

Probabilistic finite elements for transient analysis in nonlinear continua

The probabilistic finite element method (PFEM), which is a combination of finite element methods and second-moment analysis, is formulated for linear and nonlinear continua with inhomogeneous random fields. Analogous to the discretization of the displacement field in finite element methods, the random field is also discretized. The formulation is simplified by transforming the correlated variables to a set of uncorrelated variables through an eigenvalue orthogonalization. Furthermore, it is shown that a reduced set of the uncorrelated variables is sufficient for the second-moment analysis. Based on the linear formulation of the PFEM, the method is then extended to transient analysis in nonlinear continua. The accuracy and efficiency of the method is demonstrated by application to a one-dimensional, elastic/plastic wave propagation problem. The moments calculated compare favorably with those obtained by Monte Carlo simulation. Also, the procedure is amenable to implementation in deterministic FEM based computer programs.

Liu, W. K.↗
NASA NTRS1989

Inverse dynamics of a 3 degree of freedom spatial flexible manipulator

A technique is presented for solving the inverse dynamics and kinematics of 3 degree of freedom spatial flexible manipulator. The proposed method finds the joint torques necessary to produce a specified end effector motion. Since the inverse dynamic problem in elastic manipulators is closely coupled to the inverse kinematic problem, the solution of the first also renders the displacements and rotations at any point of the manipulator, including the joints. Furthermore the formulation is complete in the sense that it includes all the nonlinear terms due to the large rotation of the links. The Timoshenko beam theory is used to model the elastic characteristics, and the resulting equations of motion are discretized using the finite element method. An iterative solution scheme is proposed that relies on local linearization of the problem. The solution of each linearization is carried out in the frequency domain. The performance and capabilities of this technique are tested through simulation analysis. Results show the potential use of this method for the smooth motion control of space telerobots.

Bayo, Eduardo↗
NASA NTRS1984

The numerical solution of the Helmholtz equation for wave propagation problems in underwater acoustics

The Helmholtz Equation (-delta-K(2)n(2))u=0 with a variable index of refraction, n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. A numerical algorithm was developed and a computer code implemented that can effectively solve this equation in the intermediate frequency range. The equation is discretized using the finite element method, thus allowing for the modeling of complicated geometrices (including interfaces) and complicated boundary conditions. A global radiation boundary condition is imposed at the far field boundary that is exact for an arbitrary number of propagating modes. The resulting large, non-selfadjoint system of linear equations with indefinite symmetric part is solved using the preconditioned conjugate gradient method applied to the normal equations. A new preconditioner is developed based on the multigrid method. This preconditioner is vectorizable and is extremely effective over a wide range of frequencies provided the number of grid levels is reduced for large frequencies. A heuristic argument is given that indicates the superior convergence properties of this preconditioner.

Bayliss, A.↗
NASA NTRS1985

Coupled atmosphere/canopy model for remote sensing of plant reflectance features

Solar radiative transfer through a coupled system of atmosphere and plant canopy is modeled as a multiple-scattering problem through a layered medium of random scatterers. The radiative transfer equation is solved by the discrete-ordinates finite-element method. Analytic expressions are derived that allow the calculation of scattering and absorption cross sections for any plant canopy layer form measurable biophysical parameters such as the leaf area index, leaf angle distribution, and individual leaf reflectance and transmittance data. An expression for a canopy scattering phase function is also given. Computational results are in good agreement with spectral reflectance measurements directly above a soybean canopy, and the concept of greenness- and brightness-transforms of Landsat MSS data is reconfirmed with the computed results. A sensitivity analysis with the coupled atmosphere/canopy model quantifies how satellite-sensed spectral radiances are affected by increased atmospheric aerosols, by varying leaf area index, by anisotropic leaf scattering, and by non-Lambertian soil boundary conditions. Possible extensions to a 2-D model are also discussed.

Gerstl, S. A.↗
NASA NTRS1985

Multidimensional Modeling of Atmospheric Effects and Surface Heterogeneities on Remote Sensing

The overall goal of this project is to establish a modeling capability that allows a quantitative determination of atmospheric effects on remote sensing including the effects of surface heterogeneities. This includes an improved understanding of aerosol and haze effects in connection with structural, angular, and spatial surface heterogeneities. One important objective of the research is the possible identification of intrinsic surface or canopy characteristics that might be invariant to atmospheric perturbations so that they could be used for scene identification. Conversely, an equally important objective is to find a correction algorithm for atmospheric effects in satellite-sensed surface reflectances. The technical approach is centered around a systematic model and code development effort based on existing, highly advanced computer codes that were originally developed for nuclear radiation shielding applications. Computational techniques for the numerical solution of the radiative transfer equation are adapted on the basis of the discrete-ordinates finite-element method which proved highly successful for one and two-dimensional radiative transfer problems with fully resolved angular representation of the radiation field.

Gerstl, S. A. W.↗
NASA NTRS1987

TRANAIR - A computer code for transonic analyses of arbitrary configurations

Attention is given to a new approach to solving full potential equations about arbitrary configurations. Numerical algorithms from such fields as finite elements, preconditioned Krylov subspace methods, discrete Fourier analysis, and integral equations are combined to take advantage of the size and speed of current and emerging supercomputers. On the basis of this appraoch, a robust, efficient and easy to use computer code referred to as TRANAIR has been developed for transonic analysis of complex geometries.

Samant, S. S.↗
NASA NTRS1986

A new approach to the solution of boundary value problems involving complex configurations

A new approach for solving certain types of boundary value problems about complex configurations is presented. Numerical algorithms from such diverse fields as finite elements, preconditioned Krylov subspace methods, discrete Fourier analysis, and integral equations are combined to take advantage of the memory, speed and architecture of current and emerging supercomputers. Although the approach has application to many branches of computational physics, the present effort is concentrated in areas of Computational Fluid Dynamics (CFD) such as steady nonlinear aerodynamics, time harmonic unsteady aerodynamics, and aeroacoustics. The most significant attribute of the approach is that it can handle truly arbitrary boundary geometries and eliminates the difficult task of generating surface fitted grids.

Rubbert, P. E.↗
NASA NTRS1990

Numerical investigation of the spreading of self-excited stratified jets

The structure and evolution of self-excited subsonic periodic arrays of jets of constant and variable density are studied using spectral-element direct numerical simulations. The governing equation of motion is presented, and a method based on spectral element discretizations appropriate for simulating arbitrarily complex geometry jets and large density variations for subsonic flows is developed. Variable density fields are found to be more unstable than the corresponding uniform density fields with much higher rms values; as a result, their spreading is also considerably larger. There is a dramatic increase in spreading after a few pairings occur. Findings presented for low and high side-momentum flux reveal a shifting of the origin of instability from the near-field to the far-field, respectively, and suggest possible routes of stabilization.

Batcho, P. F.↗
NASA NTRS1980

Steady, Oscillatory, and Unsteady Subsonic and Supersonic Aerodynamics, production version (SOUSSA-P 1.1). Volume 1: Theoretical manual

Recent developments of the Green's function method and the computer program SOUSSA (Steady, Oscillatory, and Unsteady Subsonic and Supersonic Aerodynamics) are reviewed and summarized. Applying the Green's function method to the fully unsteady (transient) potential equation yields an integro-differential-delay equation. With spatial discretization by the finite-element method, this equation is approximated by a set of differential-delay equations in time. Time solution by Laplace transform yields a matrix relating the velocity potential to the normal wash. Premultiplying and postmultiplying by the matrices relating generalized forces to the potential and the normal wash to the generalized coordinates one obtains the matrix of the generalized aerodynamic forces. The frequency and mode-shape dependence of this matrix makes the program SOUSSA useful for multiple frequency and repeated mode-shape evaluations.

Morino, L.↗
NASA NTRS1988

Theoretical studies on flapped delta wings

The effects of leading edge flaps on the aerodynamic characteristics of a low aspect-ratio delta wing are studied theoretically. As an extension of the classical crossflow plane analysis and in order to include separated shear layers, an analogy between three dimensional steady conical and two dimensional unsteady self-similar flows is explored. This analogy provides a simple steady-unsteady relationship. The criteria for the validity of the steady-unsteady analogy are also examined. Two different theoretical techniques are used to represent the separated shear layers based on the steady-unsteady analogy, neglecting the trailing edge effect. In the first approach, each vortex system is represented by a pair of concentrated vortices connected to the separation points by straight feeding sheets. In the second approach, the vortex cloud method is adopted for simulating the flow field in the crossflow plane. The separated shear layers are replaced with a cloud of discrete vortices and the boundary element method is employed to represent the wing trace by a vorticity distribution. A simple merging scheme is used to model the core region of the vortical flow as a single vortex by imposing a restriction on the shear layer rotation angle. The results are compared with experiments and with results from 3-D panel calculations.

Oh, S.↗
NASA NTRS1987

Discrete element weld model, phase 2

A numerical method was developed for analyzing the tungsten inert gas (TIG) welding process. The phenomena being modeled include melting under the arc and the flow in the melt under the action of buoyancy, surface tension, and electromagnetic forces. The latter entails the calculation of the electric potential and the computation of electric current and magnetic field therefrom. Melting may occur at a single temperature or over a temperature range, and the electrical and thermal conductivities can be a function of temperature. Results of sample calculations are presented and discussed at length. A major research contribution has been the development of numerical methodology for the calculation of phase change problems in a fixed grid framework. The model has been implemented on CHAM's general purpose computer code PHOENICS. The inputs to the computer model include: geometric parameters, material properties, and weld process parameters.

Prakash, C.↗