SEARCH · Engineering Papers

Results for “discrete element modeling”

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 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Prediction of Acoustical Response of Three-dimensional Cavities Using an Indirect Boundary Element Method

A technique which characterizes the acoustics of generalized cavities with the minimum model possible is developed. All boundary element methods have two advantages over finite element methods: (1) the models are smaller; and (2) the assumed variable behavior, inherent in the method to allow discretization, is harmonic rather than polynomial. Further, IBEM often requires one rather than two numerical boundary integrals as required by DBEM. Thus, a quadratic, isoparametric IBEM program was developed. The source distribution in this solution is continuous and quadratically variable rather than continuous and constant. The program was also formulated to include the additional capability of interior point sources and impedance boundary conditions. To test the quadratic, isoparametric IBEM program, several simple cavity enclosure problems where studied. Results are shown.

Bernhard, R. J.↗

NASA NTRS1992

Iterative methods for large scale static analysis of structures on a scalable multiprocessor supercomputer

A parallel Preconditioned Conjugate Gradient (PCG) iterative solver has been developed and implemented on the iPSC-860 scalable hypercube. This new implementation makes use of the Parallel Automated Runtime Toolkit at ICASE (PARTI) primitives to efficiently program irregular communications patterns that exist in general sparse matrices and in particular in the finite element sparse stiffness matrices. The iterative PCG has been used to solve the finite element equations that result from discretizing large scale aerospace structures. In particular, the static response of the High Speed Civil Transport (HSCT) finite element model is solved on the iPSC-860.

Sobh, Nahil Atef↗

NASA NTRS1993

Current capabilities for simulating the extreme distortion of thin structures subjected to severe impacts

The explicit transient dynamics technology in use today for simulating the impact and subsequent transient dynamic response of a structure has its origins in the 'hydrocodes' dating back to the late 1940's. The growth in capability in explicit transient dynamics technology parallels the growth in speed and size of digital computers. Computer software for simulating the explicit transient dynamic response of a structure is characterized by algorithms that use a large number of small steps. In explicit transient dynamics software there is a significant emphasis on speed and simplicity. The finite element technology used to generate the spatial discretization of a structure is based on a compromise between completeness of the representation for the physical processes modelled and speed in execution. That is, since it is expected in every calculation that the deformation will be finite and the material will be strained beyond the elastic range, the geometry and the associated gradient operators must be reconstructed, as well as complex stress-strain models evaluated at every time step. As a result, finite elements derived for explicit transient dynamics software use the simplest and barest constructions possible for computational efficiency while retaining an essential representation of the physical behavior. The best example of this technology is the four-node bending quadrilateral derived by Belytschko, Lin and Tsay. Today, the speed, memory capacity and availability of computer hardware allows a number of the previously used algorithms to be 'improved.' That is, it is possible with today's computing hardware to modify many of the standard algorithms to improve their representation of the physical process at the expense of added complexity and computational effort. The purpose is to review a number of these algorithms and identify the improvements possible. In many instances, both the older, faster version of the algorithm and the improved and somewhat slower version of the algorithm are found implemented together in software. Specifically, the following seven algorithmic items are examined: the invariant time derivatives of stress used in material models expressed in rate form; incremental objectivity and strain used in the numerical integration of the material models; the use of one-point element integration versus mean quadrature; shell elements used to represent the behavior of thin structural components; beam elements based on stress-resultant plasticity versus cross-section integration; the fidelity of elastic-plastic material models in their representation of ductile metals; and the use of Courant subcycling to reduce computational effort.

Key, Samuel W.↗

NASA NTRS2002

A Multifunctional Interface Method for Coupling Finite Element and Finite Difference Methods: Two-Dimensional Scalar-Field Problems

A multifunctional interface method with capabilities for variable-fidelity modeling and multiple method analysis is presented. The methodology provides an effective capability by which domains with diverse idealizations can be modeled independently to exploit the advantages of one approach over another. The multifunctional method is used to couple independently discretized subdomains, and it is used to couple the finite element and the finite difference methods. The method is based on a weighted residual variational method and is presented for two-dimensional scalar-field problems. A verification test problem and a benchmark application are presented, and the computational implications are discussed.

Ransom, Jonathan B.↗

NASA NTRS1986

A comparison of the structureborne and airborne paths for propfan interior noise

A comparison is made between the relative levels of aircraft interior noise related to structureborne and airborne paths for the same propeller source. A simple, but physically meaningful, model of the structure treats the fuselage interior as a rectangular cavity with five rigid walls. The sixth wall, the fuselage sidewall, is a stiffened panel. The wing is modeled as a simple beam carried into the fuselage by a large discrete stiffener representing the carry-through structure. The fuselage interior is represented by analytically-derived acoustic cavity modes and the entire structure is represented by structural modes derived from a finite element model. The noise source for structureborne noise is the unsteady lift generation on the wing due to the rotating trailing vortex system of the propeller. The airborne noise source is the acoustic field created by a propeller model consistent with the vortex representation. Comparisons are made on the basis of interior noise over a range of propeller rotational frequencies at a fixed thrust.

Eversman, W.↗

NASA NTRS2000

An Assessment of Five Modeling Approaches for Thermo-Mechanical Stress Analysis of Laminated Composite Panels

A study is made of the effects of variation in the lamination and geometric parameters, and boundary conditions of multi-layered composite panels on the accuracy of the detailed response characteristics obtained by five different modeling approaches. The modeling approaches considered include four two-dimensional models, each with five parameters to characterize the deformation in the thickness direction, and a predictor-corrector approach with twelve displacement parameters. The two-dimensional models are first-order shear deformation theory, third-order theory; a theory based on trigonometric variation of the transverse shear stresses through the thickness, and a discrete layer theory. The combination of the following four key elements distinguishes the present study from previous studies reported in the literature: (1) the standard of comparison is taken to be the solutions obtained by using three-dimensional continuum models for each of the individual layers; (2) both mechanical and thermal loadings are considered; (3) boundary conditions other than simply supported edges are considered; and (4) quantities compared include detailed through-the-thickness distributions of transverse shear and transverse normal stresses. Based on the numerical studies conducted, the predictor-corrector approach appears to be the most effective technique for obtaining accurate transverse stresses, and for thermal loading, none of the two-dimensional models is adequate for calculating transverse normal stresses, even when used in conjunction with three-dimensional equilibrium equations.

Noor, A. K.↗

NASA NTRS1976

Transonic flow calculations

The development of relaxation methods to calculate transonic flows is discussed. Rather accurate predictions can be made for a number of flows of interest using the transonic potential flow equation, which may be derived from the Euler equations for inviscid compressible flow by introducing the assumption that the flow is irrotational. Once the choice of a mathematical model has been settled, the numerical procedure for actually computing a solution contains two main elements: the construction of a discrete approximation which converges to the solution of the continuous problem in the limit as the mesh width is reduced to zero, and the solution of the resulting set of nonlinear difference equations by a convergent iterative scheme. The choice of an appropriate coordinate system and its influence on the accuracy of the discrete approximation is discussed. Several applications of the general method are described.

Jameson, A.↗

NASA NTRS2000

Forced Boundary-Layer Transition on X-43 (Hyper-X) in NASA LaRC 20-Inch Mach 6 Air Tunnel

Aeroheating and boundary layer transition characteristics for the X-43 (Hyper-X) configuration have been experimentally examined in the Langley 20-Inch Mach 6 Air Tunnel. Global surface heat transfer distributions, and surface streamline patterns were measured on a 0.333-scale model of the Hyper-X forebody. Parametric variations include angles-of-attack of 0-deg, 2-deg, and 4-deg; Reynolds numbers based on model length of 1.2 to 15.4 million; and inlet cowl door both open and closed. The effects of discrete roughness elements on the forebody boundary layer, which included variations in trip configuration and height, were investigated. This document is intended to serve as a release of preliminary data to the Hyper-X program; analysis is limited to observations of the experimental trends in order to expedite dissemination.

Berry, Scott A.↗

NASA NTRS2000

Forced Boundary-Layer Transition on X-43 (Hyper-X) in NASA LaRC 31-Inch Mach 10 Air Tunnel

Aeroheating and boundary layer transition characteristics for the X-43 (Hyper-X) configuration have been experimentally examined in the Langley 31-Inch Mach 10 Air Tunnel. Global surface heat transfer distributions, and surface streamline patterns were measured on a 0.333-scale model of the Hyper-X forebody. Parametric variations include angles-of-attack of 0-deg, 2-deg, 3-deg, and 4-deg; Reynolds numbers based on model length of 1.2 to 5.1 million; and inlet cowl door both open and closed. The effects of discrete roughness elements on the forebody boundary layer, which included variations in trip configuration and height, were investigated. This document is intended to serve as a release of preliminary data to the Hyper-X program; analysis is limited to observations of the experimental trends in order to expedite dissemination.

Berry, Scott A.↗

NASA NTRS2014

Elastic Model Transitions: a Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes both quadratically constrained least squares (LSQI) and Direct Shape Mapping (DSM) algorithms to determine physical displacements. This approach is applicable to the simulation of the elastic behavior of launch vehicles and other structures that utilize multiple LTI finite element model (FEM) derived mode sets that are propagated throughout time. The time invariant nature of the elastic data for discrete segments of the launch vehicle trajectory presents a problem of how to properly transition between models while preserving motion across the transition. In addition, energy may vary between flex models when using a truncated mode set. The LSQI-DSM algorithm can accommodate significant changes in energy between FEM models and carries elastic motion across FEM model transitions. Compared with previous approaches, the LSQI-DSM algorithm shows improvements ranging from a significant reduction to a complete removal of transients across FEM model transitions as well as maintaining elastic motion from the prior state.

Jurenko, Robert J.↗

NASA NTRSDate not supplied

Wall-Modeled Large Eddy Simulation Method for Unstructured-Grid Navier-Stokes Solvers

This paper reports on the implementation and assessment of a Wall-Modeled Large-Eddy Simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-Eddy Simulations (LES) resolve large-scale turbulent-flow features and filter out small-scale effects using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge-Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multielement airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

CFD; turbulence modeling; High-Lift flow simulatio↗

NASA NTRS1978

The dynamics and control of large flexible space structures. Part A: Discrete model and modal control

Attitude control techniques for the pointing and stabilization of very large, inherently flexible spacecraft systems were investigated. The attitude dynamics and control of a long, homogeneous flexible beam whose center of mass is assumed to follow a circular orbit was analyzed. First order effects of gravity gradient were included. A mathematical model which describes the system rotations and deflections within the orbital plane was developed by treating the beam as a number of discretized mass particles connected by massless, elastic structural elements. The uncontrolled dynamics of the system are simulated and, in addition, the effects of the control devices were considered. The concept of distributed modal control, which provides a means for controlling a system mode independently of all other modes, was examined. The effect of varying the number of modes in the model as well as the number and location of the control devices were also considered.

Bainum, P. M.↗

NASA NTRS1991

Computation of rotor aerodynamic loads with a constant vorticity contour free wake model

An analytical method is presented which facilitates the study of isolated rotors with an improved approach to wake simulation. Vortex filaments are simulated along contours of constant sheet strength for the sheet of vorticity resulting from each rotor blade. Curved vortex elements comprise the filaments which can be distorted by the local velocity field. Called the Constant Vorticity Contour wake model, the approach permits the simulation of the blades' wakes corresponding to the full span of the rotor blade. The discretization of the wake of the rotor blade produces spacing and structure that are consistent with the spatial and temporal variations in the loading. A vortex-lattice aerodynamic model of the blade is also included which introduces a finite-element structural model of the blade and consideration of the force and moment trim analysis. Results of the present version of the simulation, called RotorCRAFT, are found to correlate well with H-34 flight-test data.

Quackenbush, Todd R.↗

NASA NTRS1991

An aeroelastic analysis with a generalized dynamic wake

An aeroelastic model with generalized dynamic wake is developed for application in the integration of aerodynamic, dynamic, and structural optimization of a rotor blade. The investigation is carried out with special attention to efficiency and accuracy of aeroelastic modeling. Each blade is assumed to be an elastic beam undergoing flap bending, lead-lag bending, elastic twist and axial deflections. The nonuniform blade is discretized into finite beam elements, each of which consists of twelve degrees of freedom. Such important blade design variables as pretwist, and chordwise offsets of the blade center of gravity and of the aerodynamic center from the elastic axis have been included in the analysis. Aerodynamic loads are computed from unsteady blade element theory where the rotor three-dimensional unsteady wake is modeled using a generalized dynamic wake theory. The noncirculatory loads based on unsteady thin airfoil theory are also included.

He, Cheng J.↗

NASA NTRS2020

Toward an In-Depth Material Model for Cermet Nuclear Thermal Rocket Fuel Elements

The development and qualification of nuclear thermal propulsion (NTP) fuel element technologies would be aided by an in-depth model of material response and failure modes at operating conditions. Integrated computational materials engineering techniques have the potential to provide such a model, as demonstrated here through three case studies focused on a tungsten-uranium mononitride cermet fuel. The first case focuses on the erosion of tungsten (W, also named wolfram), a nominal coating/cladding material, in hot hydrogen. Ab initio techniques are used to calculate erosion rates and thermal expansion at NTP operating conditions. The second focuses on the stability of uranium mononitride (UN) fuels at high temperature and in the presence of hydrogen. Phase diagram techniques reveal potential instabilities and decomposition pathways at high hydrogen concentrations. The third focuses on using microstructure information to predict high temperature mechanical response and failure of tungsten, used in refractory cermet materials. Combined finite element and discrete dislocation dynamics techniques provide mechanical properties in agreement with experimental methods. The integration of these techniques for an all-encompassing material model is discussed.

cermet↗

NASA NTRS1978

Modal control of the planar motion of a long flexible beam in orbit

Attitude control techniques for the pointing and stabilization of very large, inherently flexible spacecraft systems are investigated. The attitude dynamics and control of a long, homogeneous flexible beam whose center of mass is assumed to follow a circular orbit is analyzed. In this study, first order effects of gravity-gradient are included, whereas external perturbations and related orbital station keeping maneuvers are neglected. A mathematical model which describes the system deflections within the orbital plane has been developed by treating the beam as having a maximum of three discretized mass particles connected by massless, elastic structural elements. The uncontrolled dynamics of this system are simulated and, in addition, the effects of the control devices are considered. The concept of distributed modal control, which provides a means for controlling a system mode independently of all other modes, is examined. The effect of varying the number of modes in the model as well as the number and location of the control devices are also considered.

Sellappan, R.↗

NASA NTRS1991

Modal analysis of multistage gear systems coupled with gearbox vibrations

An analytical procedure to simulate vibrations in gear transmission systems is presented. This procedure couples the dynamics of the rotor-bearing gear system with the vibration in the gear box structure. The model synthesis method is used in solving the overall dynamics of the system, and a variable time-stepping integration scheme is used in evaluating the global transient vibration of the system. Locally each gear stage is modeled as a multimass rotor-bearing system using a discrete model. The modal characteristics are calculated using the matrix-transfer technique. The gearbox structure is represented by a finite element models, and modal parameters are solved by using NASTRAN. The rotor-gear stages are coupled through nonlinear compliance in the gear mesh while the gearbox structure is coupled through the bearing supports of the rotor system. Transient and steady state vibrations of the coupled system are examined in both time and frequency domains. A typical three-geared system is used as an example for demonstration of the developed procedure.

Choy, F. K.↗

NASA NTRS1991

Guided waves in a multilayered composite and ultrasonic nondestructive evaluation

Analysis of wave propagation in a laminated composite plate with many layers is made difficult by the anisotropy of the laminae, different lay-ups, many layers and defects in the fabricated plate. A good understanding of the wave field in such a plate without defects is a prerequisite for use of ultrasonic techniques to evaluate the various defects that can substantially reduce the strength and service life of this structure. Results of modeling studies of wave propagation and scattering in a multilayered composite plate are presented. The analysis technique combines the finite element discretization through the thickness with wave form representation along the plate in order to calculate the dispersive modes of propagation in the plate. To study scattering of these waves by a matrix crack that grows into delaminations, a hybrid finite element representation of the field near the crack and the modal representation of the scattered exterior field are used. Results are presented showing the dependence of the reflection and transmission coefficients on the size of the crack with delamination which can be used to size the defect using ultrasonic techniques.

Datta, S. K.↗