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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Elastic Model Transitions Using Quadratic Inequality Constrained Least Squares

A technique is presented for initializing multiple discrete finite element model (FEM) mode sets for certain types of flight dynamics formulations that rely on superposition of orthogonal modes for modeling the elastic response. Such approaches are commonly used for modeling launch vehicle dynamics, and challenges arise due to the rapidly time-varying nature of the rigid-body and elastic characteristics. By way of an energy argument, a quadratic inequality constrained least squares (LSQI) algorithm is employed to e ect a smooth transition from one set of FEM eigenvectors to another with no requirement that the models be of similar dimension or that the eigenvectors be correlated in any particular way. The physically unrealistic and controversial method of eigenvector interpolation is completely avoided, and the discrete solution approximates that of the continuously varying system. The real-time computational burden is shown to be negligible due to convenient features of the solution method. Simulation results are presented, and applications to staging and other discontinuous mass changes are discussed

Orr, Jeb S.↗
NASA NTRS1991

Solution of geometrically nonlinear statics problems by the p-version of the finite element method

This project is concerned with the possibility of using computers for the simulation of structural systems with the same degree of reliability as full scale physical experiments. Reliable numerical simulation will make it possible to reduce the costs of engineering and improve the quality of engineering decisions based on computed information. An error of idealization is an error between the actual physical quantities on which engineering decisions are based (e.g., maximum principal stress, first natural frequency, etc.) and the same data corresponding to the exact solution of the mathematical model. An error of discretization is an error between the quantities of interest corresponding to the exact and approximate solutions of a mathematical model. A high degree of reliability can be achieved in numerical simulation only if both the errors of idealization and errors of discretization can be shown to be small.

Szabo, Barna A.↗
NASA NTRS1974

Nonlinear initial-boundary value solutions by the finite element method

The finite-element method is used to establish a numerical solution algorithm for the Navier-Stokes equations for two-dimensional flows of a viscous compressible fluid. Numerical experiments confirm the advection property for the finite-element equivalent of the nonlinear convection term for both unidirectional and recirculating flowfields. For linear functionals, the algorithm demonstrates good accuracy using coarse discretizations and h squared convergence with discretization refinement.

Baker, A. J.↗
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 NTRS2006

Nonlinear Spring Finite Elements for Predicting Mode I-Dominated Delamination Growth in Laminated Structure with Through-Thickness reinforcement

One particular concern of polymer matrix composite laminates is the relatively low resistance to delamination cracking, in particular when the dominant type of failure is mode I opening. One method proposed for alleviating this problem involves the insertion pultruded carbon pins through the laminate thickness. The pins, known as z-pins, are inserted into the prepreg laminate using an ultrasonic hammer prior to the curing process, resulting in a field of pins embedded normal to the laminate plane as illustrated in Figure. 1. Pin diameters range between 0.28-mm to 0.5-mm and standard areal densities range from 0.5% to 4%. The z-pins are provided by the manufacturer, Aztex(Registered TradeMark) , in a low-density foam preform, which acts to stabilize orientation of the pins during the insertion process [1-3]. Typical pin materials include boron and carbon fibers embedded in a polymer matrix. A number of methods have been developed for predicting delamination growth in laminates reinforced with z-pins. During a study on the effect of z-pin reinforcement on mode I delamination resistance, finite element analyses of z-pin reinforced double cantilever beam (DCB) specimens were performed by Cartie and Partridge [4]. The z-pin bridging stresses were modeled by applying equivalent forces at the pin locations. Single z-pin pull-out tests were performed to characterize the traction law of the pins under mode I loading conditions. Analytical solutions for delamination growth in z-pin reinforced DCB specimens were independently derived by Robinson and Das [5] and Ratcliffe and O'Brien [6]. In the former case, pin bridging stresses were modeled using a distributed load and in the latter example the bridging stresses were discretely modeled by way of grounded springs. Additionally, Robinson and Das developed a data reduction strategy for calculating mode I fracture toughness, G(sub Ic), from a z-pin reinforced DCB specimen test [5]. In both cases a traction law similar to that adopted by Cartie and Partridge was used to represent z-pin failure under mode I loading conditions. In the current work spring elements available in most commercial finite element codes were used to model z-pins. The traction law used in previous analyses [4-6] was employed to represent z-pin damage. This method is intended for and is limited to simulating z-pins in composite laminate structure containing mode I-dominated delamination cracking. The current technique differs from previous analyses in that spring finite elements (available in commercial codes) are employed for simulating zpins, reducing the complexity of the analysis construction process. Furthermore, the analysis method can be applied to general structure that experiences mode I-dominated delamination cracking, in contrast to existing analytical solutions that are only applicable to coupon DCB specimens.

Ratcliffe, James G.↗
NASA NTRS1989

Adaptation of boundary elements to the dynamics of elastic solids

The method of boundary elements is adapted to the dynamics of elastic axisymmetric solids. Three scalar potentials are used, each satisfying the Helmholtz equation and each represented by its own source distribution. The sources are distributed over a surface enclosing the solid boundary. Stress or displacement conditions on the boundary are realized by third-order derivatives of the potentials. Discrete ring elements on the source enclosure combine their influences at discrete ring elements on the solid boundary.

El-Raheb, Michael↗
NASA NTRS2004

Distributed antenna system and method

System and methods are disclosed for employing one or more radiators having non-unique phase centers mounted to a body with respect to a plurality of transmitters to determine location characteristics of the body such as the position and/or attitude of the body. The one or more radiators may consist of a single, continuous element or of two or more discrete radiation elements whose received signals are combined. In a preferred embodiment, the location characteristics are determined using carrier phase measurements whereby phase center information may be determined or estimated. A distributed antenna having a wide angle view may be mounted to a moveable body in accord with the present invention. The distributed antenna may be utilized for maintaining signal contact with multiple spaced apart transmitters, such as a GPS constellation, as the body rotates without the need for RF switches to thereby provide continuous attitude and position determination of the body.

Fink, Patrick W.↗
NASA NTRS1976

Extended applications of the vortex lattice method

The application of the vortex lattice method to problems not usually dealt with by this technique is considered. It is shown that if the discrete vortex lattice is considered as an approximation to surface-distributed vorticity, then the concept of the generalized principal part of an integral yields a residual term to the vortex-induced velocity that renders the vortex lattice method valid for supersonic flow. Special schemes for simulating non-zero thickness lifting surfaces and fusiform bodies with vortex lattice elements are presented. Thickness effects of wing-like components are simulated by a double vortex lattice layer, and fusiform bodies are represented by a vortex grid arranged on a series of concentric cylindrical surfaces. Numerical considerations peculiar to the application of these techniques are briefly discussed.

Miranda, L. R.↗
NASA NTRS1983

New discretization and solution techniques for incompressible viscous flow problems

This paper considers several topics arising in the finite element solution of the incompressible Navier-Stokes equations. Specifically, the question of choosing finite element velocity/pressure spaces is addressed, particularly from the viewpoint of achieving stable discretizations leading to convergent pressure approximations. Following this, the role of artificial viscosity in viscous flow calculations is studied, emphasizing recent work by several researchers for the anisotropic case. The last section treats the problem of solving the nonlinear systems of equations which arise from the discretization. Time marching methods and classical iterative techniques, as well as some recent modifications are mentioned.

Gunzburger, M. D.↗
NASA NTRS1995

Dynamic sensitivity analysis of frictional contact/impact response of axisymmetric composite structures

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic frictional contact/impact 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 to 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 the nodal displacements, and the 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. Temporal integration is performed by using Newmark method. 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 from the frictional contact/impact response of a composite spherical cap impacting on a rigid plate.

Karaoglan, Levent↗
NASA NTRS2004

A Re-programmable Platform for Dynamic Burn-in Test of Xilinx Virtexll 3000 FPGA for Military and Aerospace Applications

Field Programmable Gate Arrays (FPGA) have played increasingly important roles in military and aerospace applications. Xilinx SRAM-based FPGAs have been extensively used in commercial applications. They have been used less frequently in space flight applications due to their susceptibility to single-event upsets. Reliability of these devices in space applications is a concern that has not been addressed. The objective of this project is to design a fully programmable hardware/software platform that allows (but is not limited to) comprehensive static/dynamic burn-in test of Virtex-II 3000 FPGAs, at speed test and SEU test. Conventional methods test very few discrete AC parameters (primarily switching) of a given integrated circuit. This approach will test any possible configuration of the FPGA and any associated performance parameters. It allows complete or partial re-programming of the FPGA and verification of the program by using read back followed by dynamic test. Designers have full control over which functional elements of the FPGA to stress. They can completely simulate all possible types of configurations/functions. Another benefit of this platform is that it allows collecting information on elevation of the junction temperature as a function of gate utilization, operating frequency and functionality. A software tool has been implemented to demonstrate the various features of the system. The software consists of three major parts: the parallel interface driver, main system procedure and a graphical user interface (GUI).

Field Programmable Gate Arrays (FPGA)↗
NASA NTRS1995

Stellar coronal abundances. 3: The solar first ionization potential effect determined from full-disk observation

In this paper we reanalayze the full-disk quiet-sun spectrum of Mallinovsky & Heroux (1973) with modern atomic data. The purposes of this are to check our atomic data and methods in other investigations using data from nearby stars obtained with the NASA Extreme Ultraviolet Explorer (EUVE) satellite, and to confirm that the solar first ionization potential (FIP) effect investigated by previous authors studying discrete solar regions is the same as that found in full-disk spectra. We recover the usual solar FIP effect of a coronal abundance enhancement of elements with a low FIP of a factor approximately 3-4 for lines formed at temperatures greater than approximately 10(exp 6) K. For lower temperatures, the FIP effect seems to be substantially smaller, in qualitative agreement with other data. Comparing our full-disk result with those from discrete solar structures suggest that the FIP effect is a function of altitude, with the lower temperature full-disk emission being dominated by the super-granulation network. We also compare the recent ionization balance of Arnaud & Raymond (1992) with that of Arnaud & Rothenflug (1985).

Laming, J. Martin↗
DOE OSTI2025

Is tokenization needed for masked particle modeling?

In this work, we significantly enhance masked particle modeling (MPM), a self-supervised learning scheme for constructing highly expressive representations of unordered sets relevant to developing foundation models for high-energy physics. In MPM, a model is trained to recover the missing elements of a set, a learning objective that requires no labels and can be applied directly to experimental data. We achieve significant performance improvements over previous work on MPM by addressing inefficiencies in the implementation and incorporating a more powerful decoder. We compare several pre-training tasks and introduce new reconstruction methods that utilize conditional generative models without data tokenization or discretization. We show that these new methods outperform the tokenized learning objective from the original MPM on a new test bed for foundation models for jets, which includes using a wide variety of downstream tasks relevant to jet physics, such as classification, secondary vertex finding, and track identification.

conditional generative models↗
NASA NTRS1981

Method for alleviating thermal stress damage in laminates

The method is for metallic matrix composites, such as laminated sheet or foil composites. Non-intersecting discrete discontinuities are positively introduced into the interface between the layers so as to reduce the thermal stress produced by unequal expansion of the materials making up the composite. The discontinuities are preferably produced by drilling holes in the metallic matrix layer. However, a plurality of discrete elements may be used between the layers to carry out this purpose.

Hoffman, C. A.↗
NASA NTRS2018

Implementation of a Matrix Crack Spacing Parameter in a Continuum Damage Mechanics Finite Element Model

Continuum Damage Mechanics (CDM) based progressive damage and failure analysis (PDFA) methods have demonstrated success in a variety of finite element analysis (FEA) implementations. However, the technical maturity of CDM codes has not yet been proven for the full design space of composite materials in aerospace applications. CDM-based approaches represent the presence of damage by changing the local material stiffness definitions and without updating the original mesh or element integration schemes. Without discretely representing cracks and their paths through the mesh, damage in models with CDM-based materials is often distributed in a region of partially damaged elements ahead of stress concentrations. Having a series of discrete matrix cracks represented by a softened region may affect predictions of damage propagation and, thus, structural failure. This issue can be mitigated by restricting matrix damage development to discrete, fiber-aligned rows of elements; hence CDM-based matrix cracks can be implemented to be more representative of discrete matrix cracks. This paper evaluates the effect of restricting CDM matrix crack development to discrete, fiber-aligned rows where the spacing of these rows is controlled by a user-defined crack spacing parameter. Initially, the effect of incrementally increasing matrix crack spacing in a unidirectional center notch coupon is evaluated. Then, the lessons learned from the center notch specimen are applied to open-hole compression finite element models. Results are compared to test data, and the limitations, successes, and potential of the matrix crack spacing approach are discussed.

Hyder, Imran↗
DOE OSTI2025

10-th order of accuracy for numerical solution of 3-D elasticity equations for heterogeneous materials on unfitted Cartesian meshes

We have developed the Optimal Local Truncation Error Method (OLTEM) with 10-th order of accuracy on unfitted Cartesian meshes for a system of 3-D elasticity equations with smooth irregular interfaces. 5 x 5 x 5 = 125-point stencils (similar to those for quadratic finite elements) for elastic heterogeneous materials are used for OLTEM. There are no unknowns at the interface points between different materials; the structure of the global discrete equations is the same for homogeneous and heterogeneous materials. The calculation of unknown stencil coefficients is based on the minimization of the local truncation error of the stencil equations and yields the optimal 10-th order of accuracy for OLTEM on unfitted Cartesian meshes, i.e., the increase by 7 orders in accuracy compared to quadratic finite elements on conformal meshes. A new post-processing procedure provides the 9-th order of accuracy for stresses in the 3-D case. Similar to basic computations it uses OLTEM with the 125-point stencils, the interface conditions and the elasticity equations. It was shown that the use of the elasticity equations for post-processing improves the accuracy of 0.1% stresses by 6 orders compared to post-processing without the use of PDEs. At an accuracy of for stresses, OLTEM with the new post-processing procedure reduces the number of degrees of freedom by 360 - 8000 times compared to quadratic finite elements with similar stencils. OLTEM with the 125-point stencils yields even more accurate results than high-order finite elements with much wider stencils. OLTEM provides accurate numerical results for compressible and nearly incompressible materials.

elasticity equations↗
NASA NTRS1984

The solution of non-linear hyperbolic equation systems by the finite element method

A finite-element method for the solution of nonlinear hyperbolic systems of equations, such as those encountered in non-self-adjoint problems of transient phenomena in convection-diffusion or in the mixed representation of wave problems, is developed and demonstrated. The problem is rewritten in moving coordinates and reinterpolated to the original mesh by a Taylor expansion prior to a standard Galerkin spatial discretization, and it is shown that this procedure is equivalent to the time-discretization approach of Donea (1984). Numerical results for sample problems are presented graphically, including such shallow-water problems as the breaking of a dam, the shoaling of a wave, and the outflow of a river; compressible flows such as the isothermal flow in a nozzle and the Riemann shock-tube problem; and the two-dimensional scalar-advection, nonlinear-shallow-water, and Euler equations.

Loehner, R.↗
NASA NTRS1997

Isolated Roughness Induced Boundary-Layer Transition: Shuttle Orbiter Ground Tests and Flight Experience

Recent Orbiter wind tunnel data and flight data have been used to evaluate boundary-layer transition induced by discrete-roughness elements on the Orbiter windward surface. Orbiter flow field calculations have been used to compute transition parameters and disturbance parameters for correlating the results and comparing the trends. Existing transition correlations have been modified and applied to the Orbiter. These correlations provide a means to predict transition on the Orbiter given a known isolated roughness element. Furthermore, these data and data reduction methods provide information and guidance for the prediction of transition due to &screte-roughness elements on future winged reentry vehicles.

Bouslog, Stanley A.↗