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 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS2018

Quantification of Error Associated with Using Misaligned Meshes in Continuum Damage Mechanics Material Models for Matrix Crack Growth Predictions in Composites

The ability of a material model to capture in-plane matrix mode I and mode II crack growth is an essential component for modeling ply level damage evolution in composite structures. Previous studies using a continuum damage mechanics (CDM) approach have shown success in satisfying benchmark solutions for mode I and II crack growth. However, success was shown using a fiber-aligned meshing strategy, which encourages matrix cracks to propagate in a single band of elements, along the fiber direction. Generating a fiber-aligned mesh becomes a highly involved process for laminates including off-axis (non 0° or 90°) plies. The objective of this study is to quantify the effect of non-fiber aligned mesh discretization on predictions of in-plane matrix crack propagation. The approach taken incrementally varies the mesh orientation angle relative to the fiber orientation; more specifically, misaligned meshes are used to quantify the effect of element angle orientation relative to the initial crack orientation on the energy released during matrix crack propagation simulations using a CDM method. CDM solutions obtained with the misaligned meshes are evaluated against known benchmarks for mode I and II matrix crack growth. The CDM solutions reveal a near-polynomial trend of increased predicted failure stress with increased mesh misalignment angle; hence implying a potential relationship between element orientation angle and apparent fracture toughness.

Justusson, Brian↗

NASA NTRS1984

Static shape determination and control for large space structures. I - The flexible beam. II - A large space antenna

A method for determining and controlling the shape of large, continuous space structures by means of discrete or pointwise observations and control devices is presented. The general linear boundary value problem satisfied by a one-dimensional shape function is defined, and the existence of solutions is studied. The static shape control problems for one-dimensional systems with and without rigid body modes and the static shape estimation problem are presented and solved. Eigenfunction expansions are presented which provide approximations to the algorithm terms when the associated Green's function is not known. An integral operator approach is applied to the multidimensional static problem, and the results are illustrated with a finite element model of the disk of a large space antenna which assumes no rigid body modes. It is shown that the shape control algorithm must be modified for systems with rigid body modes.

Weeks, C. J.↗

NASA NTRS1996

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

This research program dealt with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in January 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled three-component problem were developed during 1994 and 1995. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor tor parallel versions of ENG10 was developed. During 1995 and 1996 we developed the capability tor the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames. Benchmark results were presented at the 1196 Computational Aeroscience meeting.

Felippa, C. A.↗

NASA NTRS2001

Characterizing Aeroelastic Systems Using Eigenanalysis, Explicitly Retaining The Aerodynamic Degrees of Freedom

Discrete time aeroelastic models with explicitly retained aerodynamic modes have been generated employing a time marching vortex lattice aerodynamic model. This paper presents analytical results from eigenanalysis of these models. The potential of these models to calculate the behavior of modes that represent damped system motion (noncritical modes) in addition to the simple harmonic modes is explored. A typical section with only structural freedom in pitch is examined. The eigenvalues are examined and compared to experimental data. Issues regarding the convergence of the solution with regard to refining the aerodynamic discretization are investigated. Eigenvector behavior is examined; the eigenvector associated with a particular eigenvalue can be viewed as the set of modal participation factors for that particular mode. For the present formulation of the equations of motion, the vorticity for each aerodynamic element appears explicitly as an element of each eigenvector in addition to the structural dynamic generalized coordinates. Thus, modal participation of the aerodynamic degrees of freedom can be assessed in M addition to participation of structural degrees of freedom.

Heeg, Jennifer↗

NASA NTRS1972

Analysis of a laterally loaded ring with a hinged cross section.

A ring assembly constructed by lacing together three elements into a basic channel cross section is analyzed. The ring is supported at three equidistant points and loaded by a uniform distribution of radial and transverse loads. Bulkheads may be introduced at discrete cross sections to prevent distortion. The purpose of this analysis is to determine the behavior of the deflection as the number and location of the bulkheads and the degree of rigidity of the supports are varied. The method of analysis is an application of the Principle of Virtual Work within the framework of small displacement theory. Numerical results are presented for the geometrical parameters of a model. An important result is the observation that bulkheads have virtually no effect on the deflections of the web.

Williams, H. E.↗

NASA NTRS1993

Remote sensing x ray fluorescence spectrometry for future lunar exploration missions

Measurement of discrete line X-ray emission from space can be used to obtain both qualitative and quantitative elemental surface composition information. Remote orbital X-ray elemental analysis is the measurement of characteristic X-rays following the interaction of solar X-rays with the surface of a given solid body. Thus, X-ray emission from a surface is strongly dependent on the incident solar spectrum as well as on the chemical composition of the surface. In this paper, the relationship between the solar flux and the major lunar element fluorescence and scattered X-ray flux is characterized. A model was developed to calculate lunar surface X-ray fluorescence emission spectra under a variety of solar conditions and, in particular, conditions that might be expected for a Lunar Scout mission proposed for 1995 or 1996. This time period should be during solar minimum.

Clark, Pamela E.↗

NASA NTRS1996

Analysis of Sensory/Active Piezoelectric Composite Structures in Thermal Environments

Although there has been extensive development of analytical methods for modeling the behavior of piezoelectric structures, only a limited amount of research has been performed concerning the implications of thermal effects on both the active and sensory response of smart structures. Thermal effects become important when the piezoelectric structure has to operate in either extremely hot or cold temperature environments. Consequently, the purpose of this paper is to extend the previously developed discrete layer formulation of Saravanos and Heyliger to account for the coupled mechanical, electrical, and thermal response in modern smart composite beams. The mechanics accounts for thermal effects which may arise in the elastic and piezoelectric media at the material level through the constitutive equations. The displacements, electric potentials, and temperatures are introduced as state variables, allowing them to be modeled as variable fields through the laminate thickness. This unified representation leads to an inherent capability to model both the active compensation of thermal distortions in smart structures and the resultant sensory voltage when thermal loads are applied. The corresponding finite element formulation is developed and numerical results demonstrate the ability to model both the active and sensory modes of composite beams with heterogeneous plies with attached piezoelectric layers under thermal loadings.

Lee, Ho-Jun↗

NASA NTRS1998

Control System and Flexible Satellite Interaction During Orbit Transfer Maneuver

In this paper the interaction between the attitude control system and the flexible structure of an artificial satellite during orbit transfer maneuver has been investigated. The satellite was modeled by a rigid central body with one or more flexible appendages. The dynamics equations were obtained by Lagrangean approach. The flexible appendages were treated as clamped-free beam and its displacement was discretized by assumed- mode method. In order to transfer the satellite, a typical Hohmann transfer and a burn-coast-burn strategy were used and the attitude was controlled by an on-off controller. During transfer procedure a global analysis of satellite has been done, such as: performance of control system, influence of elastic response in control system, thruster firing frequency, fuel consumption and variation of orbital elements. In order to avoid the interaction with structure motion, a control system with bandwidth of one decade bellow the fundamental frequency was used. In the simulations the firing frequency was evaluated in an approximately way but kept below the fundamental frequency of the structure. The control system has kept the attitude below the specifications. As a result, the orbit transfer maneuvering has been done correctly without excessive excitation of flexible appendage.

Adenilson Roberto da Silva↗

NASA NTRS2015

Grid Convergence for Turbulent Flows(Invited)

A detailed grid convergence study has been conducted to establish accurate reference solutions corresponding to the one-equation linear eddy-viscosity Spalart-Allmaras turbulence model for two dimensional turbulent flows around the NACA 0012 airfoil and a flat plate. The study involved three widely used codes, CFL3D (NASA), FUN3D (NASA), and TAU (DLR), and families of uniformly refined structured grids that differ in the grid density patterns. Solutions computed by different codes on different grid families appear to converge to the same continuous limit, but exhibit different convergence characteristics. The grid resolution in the vicinity of geometric singularities, such as a sharp trailing edge, is found to be the major factor affecting accuracy and convergence of discrete solutions, more prominent than differences in discretization schemes and/or grid elements. The results reported for these relatively simple turbulent flows demonstrate that CFL3D, FUN3D, and TAU solutions are very accurate on the finest grids used in the study, but even those grids are not sufficient to conclusively establish an asymptotic convergence order.

Diskin, Boris↗

NASA NTRS1983

A tire contact solution technique

An efficient method for calculating the contact boundary and interfacial pressure distribution was developed. This solution technique utilizes the discrete Fourier transform to establish an influence coefficient matrix for the portion of the pressurized tire surface that may be in the contact region. This matrix is used in a linear algebra algorithm to determine the contact boundary and the array of forces within the boundary that are necessary to hold the tire in equilibrium against a specified contact surface. The algorithm also determines the normal and tangential displacements of those points on the tire surface that are included in the influence coefficient matrix. Displacements within and outside the contact region are calculated. The solution technique is implemented with a finite-element tire model that is based on orthotropic, nonlinear shell of revolution elements which can respond to nonaxisymmetric loads. A sample contact solution is presented.

Tielking, J. T.↗

NASA NTRSDate not supplied

Bistable Deployable Composite Booms With Parabolic Cross-Sections

This paper investigates how stable equilibrium states in the extended and coiled configurations can be predicted in thin-shelled composite booms with parabolic cross-sections. These conic shapes potentially offer greater stiffness properties when compared to circular cross-sections, which is critical for improving the load bearing performance of deployed booms. Inducing bistability through the choice of composite shell layups in parabolic booms would allow for controllable self-deployment due to a less energetic coiled state when compared to monostable booms. An inextensional analytical model is used to predict the stable coiled diameters of tape spring and Collapsible Tubular Mast (CTM) booms with parabolic cross-sections. The parabolic section is discretized into circular segments using biarc spline interpolation, which allows them to be integrated into the strain energy minimization procedure used to obtain the equilibrium states. When the parabolic booms are parametrically compared against circular booms with identical layups, flattened height, and mass, the former are found to generally have better stiffness performance while being less efficient in stowed volume as evidenced by larger coiled diameters. Analytical coiled diameters and their strain energy are verified with finite element simulations for optimal parabolic tape spring and CTM booms.

Andrew J. Lee↗

NASA NTRS1993

Using neural networks to assist in OPAD data analysis

The space shuttle main engine (SSME) became the subject of plume emission spectroscopy in 1986 when researchers from NASA-Marshall Space Flight Center (MSFC), Arnold Engineering Development Center (AEDC), and Rocketdyne went to the SSME test stands at the NASA-Stennis Space Center and at Rocketdyne's Santa Susan Field Laboratory to optically observe the plume. Since then, plume spectral acquisitions have recorded many nominal tests and the qualitative spectral features of the SSME plume are now well established. Significant discoveries made with both wide-band and narrow-band plume emission spectroscopy systems led MSFC to promote the Optical Plume Anomaly Detection (OPAD) program with a goal of instrumenting all SSME test stands with customized spectrometer systems. A prototype OPAD system is now installed on the SSME Technology Test Bed (TTB) at MSFC. The OPAD system instrumentation consists of a broad-band, optical multiple-channel analyzer (OMA) and a narrow-band device called a polychrometer. The OMA is a high-resolution (1.5-2.0 Angstroms) 'super-spectrometer' covering the near-ultraviolet to near-infrared waveband (2800-7400 Angstroms), providing two scans per second. The polychrometer consists of sixteen narrow-band radiometers: fourteen monitoring discrete wavelengths of health and condition monitoring elements and two dedicated to monitoring background emissions. All sixteen channels are capable of providing 500 samples per second. To date, the prototype OPAD system has been used during 43 SSME firings on the TTB, collecting well over 250 megabytes of plume spectral data. One goal of OPAD data analysis is to determine interatively with the help of a computer code, SPECTRA4, developed at AEDC. Experience has shown that iteration with SPECTRA4 is an incredibly labor intensive task and not one to be performed by band. What is really needed is the 'inverse' of SPECTRA4 but the mathematical model for this inverse mapping is tenuous at best. However, the robustness of PSECTRA4 run in the 'forward' direction means that accurate input/output mappings can be obtained. If the mappings were inverted (i.e., input becomes output and output becomes input) then an 'inverse' of SPECTRA4 would be at hand but the 'model' would be specific to the data utilized and would in no way be general. Building a generalized model based upon known input/output mappings while ignoring the details of the governing physical model is possible through the use of a neural network. The research investigation described involves the development of a neural network to provide a generalized 'inverse' of SPECTRA4. The objectives of the research were to design an appropriate neural network architecture, train the network, and then evaluate its performance.

Whitaker, Kevin W.↗

NASA NTRS2008

Tone-Based Command of Deep Space Probes using Ground Antennas

A document discusses a technique for enabling the reception of spacecraft commands at received signal levels as much as three orders of magnitude below those of current deep space systems. Tone-based commanding deals with the reception of commands that are sent in the form of precise frequency offsets using an open-loop receiver. The key elements of this technique are an ultrastable oscillator and open-loop receiver onboard the spacecraft, both of which are part of the existing New Horizons (Pluto flyby) communications system design. This enables possible flight experimentation for tone-based commanding during the long cruise of the spacecraft to Pluto. In this technique, it is also necessary to accurately remove Doppler shift from the uplink signal presented to the spacecraft. A signal processor in the spacecraft performs a discrete Fourier transform on the received signal to determine the frequency of the received signal. Due to the long-term drift in the oscillators and orbit prediction model, the system is likely to be implemented differentially, where changes in the uplink frequency convey the command information.

Bokulic, Robert S.↗

NASA NTRS2015

Improved Convergence and Robustness of USM3D Solutions on Mixed Element Grids (Invited)

Several improvements to the mixed-element USM3D discretization and defect-correction schemes have been made. A new methodology for nonlinear iterations, called the Hierarchical Adaptive Nonlinear Iteration Scheme (HANIS), has been developed and implemented. It provides two additional hierarchies around a simple and approximate preconditioner of USM3D. The hierarchies are a matrix-free linear solver for the exact linearization of Reynolds-averaged Navier Stokes (RANS) equations and a nonlinear control of the solution update. Two variants of the new methodology are assessed on four benchmark cases, namely, a zero-pressure gradient flat plate, a bump-in-channel configuration, the NACA 0012 airfoil, and a NASA Common Research Model configuration. The new methodology provides a convergence acceleration factor of 1.4 to 13 over the baseline solver technology.

Pandya, Mohagna J.↗

NASA NTRS1995

Parallel automated adaptive procedures for unstructured meshes

Consideration is given to the techniques required to support adaptive analysis of automatically generated unstructured meshes on distributed memory MIMD parallel computers. The key areas of new development are focused on the support of effective parallel computations when the structure of the numerical discretization, the mesh, is evolving, and in fact constructed, during the computation. All the procedures presented operate in parallel on already distributed mesh information. Starting from a mesh definition in terms of a topological hierarchy, techniques to support the distribution, redistribution and communication among the mesh entities over the processors is given, and algorithms to dynamically balance processor workload based on the migration of mesh entities are given. A procedure to automatically generate meshes in parallel, starting from CAD geometric models, is given. Parallel procedures to enrich the mesh through local mesh modifications are also given. Finally, the combination of these techniques to produce a parallel automated finite element analysis procedure for rotorcraft aerodynamics calculations is discussed and demonstrated.

Shephard, M. S.↗

NASA NTRS2013

Spectrally Consistent Scattering, Absorption, and Polarization Properties of Atmospheric Ice Crystals at Wavelengths from 0.2 to 100 um

A data library is developed containing the scattering, absorption, and polarization properties of ice particles in the spectral range from 0.2 to 100 microns. The properties are computed based on a combination of the Amsterdam discrete dipole approximation (ADDA), the T-matrix method, and the improved geometric optics method (IGOM). The electromagnetic edge effect is incorporated into the extinction and absorption efficiencies computed from the IGOM. A full set of single-scattering properties is provided by considering three-dimensional random orientations for 11 ice crystal habits: droxtals, prolate spheroids, oblate spheroids, solid and hollow columns, compact aggregates composed of eight solid columns, hexagonal plates, small spatial aggregates composed of 5 plates, large spatial aggregates composed of 10 plates, and solid and hollow bullet rosettes. The maximum dimension of each habit ranges from 2 to 10,000 microns in 189 discrete sizes. For each ice crystal habit, three surface roughness conditions (i.e., smooth, moderately roughened, and severely roughened) are considered to account for the surface texture of large particles in the IGOM applicable domain. The data library contains the extinction efficiency, single-scattering albedo, asymmetry parameter, six independent nonzero elements of the phase matrix (P11, P12, P22, P33, P43, and P44), particle projected area, and particle volume to provide the basic single-scattering properties for remote sensing applications and radiative transfer simulations involving ice clouds. Furthermore, a comparison of satellite observations and theoretical simulations for the polarization characteristics of ice clouds demonstrates that ice cloud optical models assuming severely roughened ice crystals significantly outperform their counterparts assuming smooth ice crystals.

Optical properties↗

NASA NTRS1977

Finite element analysis of aeroacoustic jet-flap flows

A computational analysis was performed on the steady, turbulent aerodynamic flowfields associated with a jet-blown flap. For regions devoid of flow separation, a parabolic approximation to the governing time-averaged Navier-Stokes equations was applied. Numerical results are presented for the symmetry plane flow of a slot-nozzle planar jet flap geometry, including prediction of flowfield evolution within the secondary mixing region immediately downstream of the trailing edge. Using a two equation turbulence kinetic energy closure model, rapid generation and decay of large spatial gradients in mean and correlated fluctuating velocity components within the immediate wake region were predicted. Modifications to the turbulent flow structure, as induced by porous surface treatment of the flap, were evaluated. The recirculating flow within a representative discrete slot in the surface was evaluated, using the two dimensional, time-averaged Navier-Stokes equations.

Baker, A. J.↗

NASA NTRS1992

Electromagnetic scattering and radiation from microstrip patch antennas and spirals residing in a cavity

A new hybrid method is presented for the analysis of the scattering and radiation by conformal antennas and arrays comprised of circular or rectangular elements. In addition, calculations for cavity-backed spiral antennas are given. The method employs a finite element formulation within the cavity and the boundary integral (exact boundary condition) for terminating the mesh. By virtue of the finite element discretization, the method has no restrictions on the geometry and composition of the cavity or its termination. Furthermore, because of the convolutional nature of the boundary integral and the inherent sparseness of the finite element matrix, the storage requirement is kept very low at O(n). These unique features of the method have already been exploited in other scattering applications and have permitted the analysis of large-size structures with remarkable efficiency. In this report, we describe the method's formulation and implementation for circular and rectangular patch antennas in different superstrate and substrate configurations which may also include the presence of lumped loads and resistive sheets/cards. Also, various modelling approaches are investigated and implemented for characterizing a variety of feed structures to permit the computation of the input impedance and radiation pattern. Many computational examples for rectangular and circular patch configurations are presented which demonstrate the method's versatility, modeling capability and accuracy.

Volakis, J. L.↗