SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Multigrid solution of the Navier-Stokes equations on triangular meshes

A Navier-Stokes algorithm for use on unstructured triangular meshes is presented. Spatial discretization of the governing equations is achieved using a finite element Galerkin approximation, which can be shown to be equivalent to a finite volume approximation for regular equilateral triangular meshes. Integration steady-state is performed using a multistage time-stepping scheme, and convergence is accelerated by means of implicit residual smoothing and an unstructured multigrid algorithm. Directional scaling of the artificial dissipation and the implicit residual smoothing operator is achieved for unstructured meshes by considering local mesh stretching vectors at each point. The accuracy of the scheme for highly stretched triangular meshes is validated by comparing computed flat-plate laminar boundary layer results with the well known similarity solution, and by comparing laminar airfoil results with those obtained from various well-established structured quadrilateral-mesh codes. The convergence efficiency of the present method is also shown to be competitive with those demonstrated by structured quadrilateral-mesh algorithms.

Mavriplis, Dimitri J.↗

NASA NTRS1990

Multigrid solution of the Navier-Stokes equations on triangular meshes

Mavriplis, D. J.↗

NASA NTRS1994

Correlation of finite element free vibration predictions using random vibration test data

Finite element analysis is regularly used during the engineering cycle of mechanical systems to predict the response to static, thermal, and dynamic loads. The finite element model (FEM) used to represent the system is often correlated with physical test results to determine the validity of analytical results provided. Results from dynamic testing provide one means for performing this correlation. One of the most common methods of measuring accuracy is by classical modal testing, whereby vibratory mode shapes are compared to mode shapes provided by finite element analysis. The degree of correlation between the test and analytical mode shapes can be shown mathematically using the cross orthogonality check. A great deal of time and effort can be exhausted in generating the set of test acquired mode shapes needed for the cross orthogonality check. In most situations response data from vibration tests are digitally processed to generate the mode shapes from a combination of modal parameters, forcing functions, and recorded response data. An alternate method is proposed in which the same correlation of analytical and test acquired mode shapes can be achieved without conducting the modal survey. Instead a procedure is detailed in which a minimum of test information, specifically the acceleration response data from a random vibration test, is used to generate a set of equivalent local accelerations to be applied to the reduced analytical model at discrete points corresponding to the test measurement locations. The static solution of the analytical model then produces a set of deformations that once normalized can be used to represent the test acquired mode shapes in the cross orthogonality relation. The method proposed has been shown to provide accurate results for both a simple analytical model as well as a complex space flight structure.

Chambers, Jeffrey A.↗

NASA NTRS2010

Modeling Optical Properties of Mineral Aerosol Particles by Using Nonsymmetric Hexahedra

We explore the use of nonsymmetric geometries to simulate the single-scattering properties of airborne dust particles with complicated morphologies. Specifically, the shapes of irregular dust particles are assumed to be nonsymmetric hexahedra defined by using the Monte Carlo method. A combination of the discrete dipole approximation method and an improved geometric optics method is employed to compute the single-scattering properties of dust particles for size parameters ranging from 0.5 to 3000. The primary optical effect of eliminating the geometric symmetry of regular hexahedra is to smooth the scattering features in the phase function and to decrease the backscatter. The optical properties of the nonsymmetric hexahedra are used to mimic the laboratory measurements. It is demonstrated that a relatively close agreement can be achieved by using only one shape of nonsymmetric hexahedra. The agreement between the theoretical results and their measurement counterparts can be further improved by using a mixture of nonsymmetric hexahedra. It is also shown that the hexahedron model is much more appropriate than the "equivalent sphere" model for simulating the optical properties of dust particles, particularly, in the case of the elements of the phase matrix that associated with the polarization state of scattered light.

Bi, Lei↗

DOE OSTI2024

A layered solid finite element formulation with interlaminar enhanced displacements for the modeling of laminated composite structures

Accurate modeling of layered composite structures often requires the use of detailed finite element models which can sufficiently resolve the kinematics and material behavior within each layer of the composite. However, individually discretizing each material layer into finite elements presents a prohibitive computational expensive given the large number of thin layers comprising some laminated composites. To address these challenges, an 8-node layered solid hexahedral finite element is formulated with the aim of striking an appropriate balance between efficiency and fidelity. The element is discretized into an arbitrary number of distinct material layers, and employs reduced in-plane integration within each layer. The chosen reduced integration scheme is supplemented by a novel physical stabilization approach which includes layerwise enhancements to mitigate various forms of locking phenomena. The proposed framework additionally supports the inclusion of interlaminar enhanced displacements to better represent the kinematics of general layered composite materials. Finally, the described element formulation has been implemented in the ParaDyn finite element code, and its efficacy for modeling laminated composite structures is demonstrated on a variety of verification problems.

42 ENGINEERING↗

DOE OSTI2023

High-Order Mesh Morphing for Boundary and Interface Fitting to Implicit Geometries

Here, we propose a method that morphs high-order meshes such that their boundaries and interfaces coincide/align with implicitly defined geometries. Our focus is particularly on the case when the target surface is prescribed as the zero isocontour of a smooth discrete function. Common examples of this scenario include using level set functions to represent material interfaces in multimaterial configurations, and evolving geometries in shape and topology optimization. The proposed method formulates the mesh optimization problem as a variational minimization of the sum of a chosen mesh-quality metric using the Target-Matrix Optimization Paradigm (TMOP) and a penalty term that weakly forces the selected faces of the mesh to align with the target surface. The distinct features of the method are use of a source mesh to represent the level set function with sufficient accuracy, and adaptive strategies for setting the penalization weight and selecting the faces of the mesh to be fit to the target isocontour of the level set field. We demonstrate that the proposed method is robust for generating boundary- and interface-fitted meshes for curvilinear domains using different element types in 2D and 3D.

97 MATHEMATICS AND COMPUTING↗

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↗

DOE OSTI2023

On a fully-implicit VMS-stabilized FE formulation for low Mach number compressible resistive MHD with application to MCF

This study presents the development and evaluation of a fully-implicit variational multiscale (VMS) stabilized unstructured finite element (FE) formulation for compressible magnetohydrodynamics (MHD) model, at low Mach number regime. The model describes the dynamics of a compressible conducting fluid in the low Mach number limit in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics phenomena. The specific applications that motivate this study are macroscopic simulations of the longer time-scale stability and disruptions of magnetic confinement fusion (MCF) devices, specifically the ITER tokamak. The discussion considers the development of the VMS FE representation, the structure of the stabilizing terms that deal with significant convective flows, the stabilization of the nearly incompressible response of the fluid flow, and the stabilization of the constraint that enforces the solenoidal involution on the magnetic field. The nonlinear discretized system is solved with scalable preconditioned Newton–Krylov iterative methods, which employs a multiphysics block preconditioning method based on approximate block factorizations and Schur complements. The study presents an evaluation of the VMS method on a 2D cartesian tearing mode instability, and illustrates the scalability of the solvers on MCF relevant problems. A set of results are also presented for longer time-scale stability and disruptions for the ITER tokamak. These include a vertical displacement event (VDE), and a (1,1) internal kink mode. Here, the formulation is demonstrated to be scalable and also reasonably robust with respect to the Lundquist number scaling.

42 ENGINEERING↗

NASA NTRS1991

The first ANDES elements: 9-DOF plate bending triangles

New elements are derived to validate and assess the assumed natural deviatoric strain (ANDES) formulation. This is a brand new variant of the assumed natural strain (ANS) formulation of finite elements, which has recently attracted attention as an effective method for constructing high-performance elements for linear and nonlinear analysis. The ANDES formulation is based on an extended parametrized variational principle developed in recent publications. The key concept is that only the deviatoric part of the strains is assumed over the element whereas the mean strain part is discarded in favor of a constant stress assumption. Unlike conventional ANS elements, ANDES elements satisfy the individual element test (a stringent form of the patch test) a priori while retaining the favorable distortion-insensitivity properties of ANS elements. The first application of this formulation is the development of several Kirchhoff plate bending triangular elements with the standard nine degrees of freedom. Linear curvature variations are sampled along the three sides with the corners as gage reading points. These sample values are interpolated over the triangle using three schemes. Two schemes merge back to conventional ANS elements, one being identical to the Discrete Kirchhoff Triangle (DKT), whereas the third one produces two new ANDES elements. Numerical experiments indicate that one of the ANDES element is relatively insensitive to distortion compared to previously derived high-performance plate-bending elements, while retaining accuracy for nondistorted elements.

Militello, Carmelo↗

DOE OSTI2022

Bringing Trimmed Serendipity Methods to Computational Practice in Firedrake

We present an implementation of the trimmed serendipity finite element family, using the open-source finite element package Firedrake. The new elements can be used seamlessly within the software suite for problems requiring H 1 , H (curl), or H (div)-conforming elements on meshes of squares or cubes. To test how well trimmed serendipity elements perform in comparison to traditional tensor product elements, we perform a sequence of numerical experiments including the primal Poisson, mixed Poisson, and Maxwell cavity eigenvalue problems. Overall, we find that the trimmed serendipity elements converge, as expected, at the same rate as the respective tensor product elements, while being able to offer significant savings in the time or memory required to solve certain problems.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2012

Derivation of Apollo 14 High-Al Basalts at Discrete Times: Rb-Sr Isotopic Constraints

Pristine Apollo 14 (A-14) high-Al basalts represent the oldest volcanic deposits returned from the Moon [1,2] and are relatively enriched in Al2O3 (>11 wt%) compared to other mare basalts (7-11 wt%). Literature Rb-Sr isotopic data suggest there are at least three different eruption episodes for the A-14 high-Al basalts spanning the age range approx.4.3 Ga to approx.3.95 Ga [1,3]. Therefore, the high-Al basalts may record lunar mantle evolution between the formation of lunar crust (approx.4.4 Ga) and the main basin-filling mare volcanism (<3.85 Ga) [4]. The high-Al basalts were originally classified into five compositional groups [5,6], and then regrouped into three with a possible fourth comprising 14072 based on the whole-rock incompatible trace element (ITE) ratios and Rb-Sr radiometric ages [7]. However, Rb-Sr ages of these basalts from different laboratories may not be consistent with each other because of the use of different 87Rb decay constants [8] and different isochron derivation methods over the last four decades. This study involved a literature search for Rb-Sr isotopic data previously reported for the high-Al basalts. With the re-calculated Rb-Sr radiometric ages, eruption episodes of A-14 high-Al basalts were determined, and their petrogenesis was investigated in light of the "new" Rb-Sr isotopic data and published trace element abundances of these basalts.

Hui. Hejiu↗

DOE OSTI2023

Conformational space exploration of cryo-EM structures by variability refinement

Cryo-EM observation of biological samples enables visualization of sample heterogeneity, in the form of discrete states that are separable, or continuous heterogeneity as a result of local protein motion before flash freezing. Variability analysis of this continuous heterogeneity describes the variance between a particle stack and a volume, and results in a map series describing the various steps undertaken by the sample in the particle stack. While this observation is absolutely stunning, it is very hard to pinpoint structural details to elements of the maps. Here, in order to bridge the gap between observation and explanation, we designed a tool that refines an ensemble of structures into all the maps from variability analysis. Using this bundle of structures, it is easy to spot variable parts of the structure, as well as the parts that are not moving. Comparison with molecular dynamics simulations highlights the fact that the movements follow the same directions, albeit with different amplitudes. Ligand can also be investigated using this method. Variability refinement is available in the Phenix software suite, accessible under the program name phenix.varref.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1985

Contouring Trivariate Data

In many applications, the data consist of discrete 3-D points at which one or more parameters are given. To display contours, the data are represented by a continuous function which is evaluated at any point as needed for contouring. Contouring results are presented which are applicable both to artitrarily spaced data and to data which lie on a topologically rectangular three dimensional grid. the contours are assumed to be described by 3-D display lists for viewing on a dynamic color graphics device; that is, they will not simply be projected into 2-D and viewed as a static image on a frame buffer. Dynamic viewing of color contours may be essential to the proper interpretation of results. The gridded data lie on a topologically rectangular grid although two or more nodes may be the same point. Parametric tensor product methods may be used to fit the gridded data and to generate the contours. Rectangular elements are convenient but not necessary. For example, there are other methods which are effective for contouring over tetrahedral elements.

Stead, S. E.↗

NASA NTRS2016

Simulation of Satellite, Airborne and Terrestrial LiDAR with DART (I):Waveform Simulation with Quasi-Monte Carlo Ray Tracing

Light Detection And Ranging (LiDAR) provides unique data on the 3-D structure of atmosphere constituents and the Earth's surface. Simulating LiDAR returns for different laser technologies and Earth scenes is fundamental for evaluating and interpreting signal and noise in LiDAR data. Different types of models are capable of simulating LiDAR waveforms of Earth surfaces. Semi-empirical and geometric models can be imprecise because they rely on simplified simulations of Earth surfaces and light interaction mechanisms. On the other hand, Monte Carlo ray tracing (MCRT) models are potentially accurate but require long computational time. Here, we present a new LiDAR waveform simulation tool that is based on the introduction of a quasi-Monte Carlo ray tracing approach in the Discrete Anisotropic Radiative Transfer (DART) model. Two new approaches, the so-called "box method" and "Ray Carlo method", are implemented to provide robust and accurate simulations of LiDAR waveforms for any landscape, atmosphere and LiDAR sensor configuration (view direction, footprint size, pulse characteristics, etc.). The box method accelerates the selection of the scattering direction of a photon in the presence of scatterers with non-invertible phase function. The Ray Carlo method brings traditional ray-tracking into MCRT simulation, which makes computational time independent of LiDAR field of view (FOV) and reception solid angle. Both methods are fast enough for simulating multi-pulse acquisition. Sensitivity studies with various landscapes and atmosphere constituents are presented, and the simulated LiDAR signals compare favorably with their associated reflectance images and Laser Vegetation Imaging Sensor (LVIS) waveforms. The LiDAR module is fully integrated into DART, enabling more detailed simulations of LiDAR sensitivity to specific scene elements (e.g., atmospheric aerosols, leaf area, branches, or topography) and sensor configuration for airborne or satellite LiDAR sensors.

Gastellu-Etchegorry, Jean-Philippe↗

DOE OSTI2021

Sparse Approximate Multifrontal Factorization with Butterfly Compression for High-Frequency Wave Equations

In this work, we present a fast and approximate multifrontal solver for large-scale sparse linear systems arising from finite-difference, finite-volume or finite-element discretization of high-frequency wave equations. The proposed solver leverages the butterfly algorithm and its hierarchical matrix extension for compressing and factorizing large frontal matrices via graph-distance guided entry evaluation or randomized matrix-vector multiplication-based schemes. Complexity analysis and numerical experiments demonstrate $\mathcal{O}(N\log^2 N)$ computation and $\mathcal{O}(N)$ memory complexity when applied to an $N\times N$ sparse system arising from 3D high-frequency Helmholtz and Maxwell problems.

97 MATHEMATICS AND COMPUTING↗

NASA NTRSDate not supplied

Variational Asymptotic Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites

The objective of this paper is to develop a constitutive model for finitely deformed viscoelastic-viscoplastic materials and a micromechanics approach to homogenizing composites consisting of such materials. The development of the constitutive model involves establishing a thermodynamic framework based on finite strain theory, developing a viscoelasticity and a viscoplasticity model based on the thermodynamic framework, developing a radial return algorithm based on a classic framework, and deriving a closed-from incremental constitutive relation in the spatial configuration. The development of the micromechanics approach involves pulling-back the above constitutive relation to the material configuration, formulating a variational statement with the resulting constitutive relation, discretizing variational statement in a finite-dimensional space, and solving the discretized variational statement using an Euler–Newton predictor–corrector method. The constitutive model is calibrated via monotonic uniaxial compression tests on a polymer, and the calibrated model is validated by comparing its predictions with the cyclic test data. It is shown capable of characterizing viscoelasticity, viscoplasticity, and complex loading paths. The micromechanics approach’s capabilities are demonstrated through homogenizing a unidirectional fiber-reinforced composite, subjected to uniaxial, biaxial, and shear loading, at different strain rates. It is demonstrated to be capable of handling rate dependence and complex loading paths. The present framework can be further improved by implementing more sophisticated viscoelasticity and viscoplasticity models in future work.

Finite element analysis↗

NASA NTRS1995

Large-eddy simulation of flow past a circular cylinder

Some of the most challenging applications of large-eddy simulation are those in complex geometries where spectral methods are of limited use. For such applications more conventional methods such as finite difference or finite element have to be used. However, it has become clear in recent years that dissipative numerical schemes which are routinely used in viscous flow simulations are not good candidates for use in LES of turbulent flows. Except in cases where the flow is extremely well resolved, it has been found that upwind schemes tend to damp out a significant portion of the small scales that can be resolved on the grid. Furthermore, it has been found that even specially designed higher-order upwind schemes that have been used successfully in the direct numerical simulation of turbulent flows produce too much dissipation when used in conjunction with large-eddy simulation. The objective of the current study is to perform a LES of incompressible flow past a circular cylinder at a Reynolds number of 3900 using a solver which employs an energy-conservative second-order central difference scheme for spatial discretization and compare the results obtained with those of Beaudan & Moin (1994) and with the experiments in order to assess the performance of the central scheme for this relatively complex geometry.

Mittal, R.↗

DOE OSTI2021

Scalable Approaches to Selecting Key Entities in Large Networked Infrastructure Systems

This work aims at bringing advances in discrete optimization algorithms to solving practical engineering problems at scale. Often times, in many engineering design problems, there is a need to select a small set of influential or representative elements from a large ground set of entities in an optimal fashion. Submodular optimization provides for a formal way to solve such problems. Common examples with infrastructure systems involve sensor placement and identification of key entities with certain objectives. However, scaling these approaches to large infrastructure systems can be challenging because of the high computational complexity of the overall framework that include the optimization algorithms as well as high-complexity compute-oracles that provide the necessary objective function values. In this work, we explore a well-studied and widely-applicable paradigm, namely leader-selection in a multi-agent networked setting in the context of scalable methodologies. We demonstrate novel frameworks that utilize variations of accelerated submodular optimization algorithms along with linear-algebraic methods that can help accelerate the oracle computations. We further explore this combination in conjunction with graph partitioning paradigms to take advantage of the accelerated algorithms in a distributed setting. Finally we demonstrate the key findings on a practical problem in an operational setting. For this, we leverage an example road network with approximately 18k nodes and 27k edges in a traffic control application, where we seek a limited number of k=200 key intersections. This problem can be solved in a serial setting in just under 5 hours providing more than 2 orders of magnitude speed-up over methods that do not consider acceleration techniques.

Visweswara Sathanur, Arun↗