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 577 records · Page 32

Modeling Deformable Linear Objects for Autonomous Robotic Outfitting of Lunar Surface Systems

This paper presents structural models of deformable linear objects (DLOs). DLOs are a subclass of deformable objects that encompasses common outfitting elements such as cables and ropes. Models are validated through hardware experiments, and integration in a robotic autonomy architecture for space environments is discussed. A persistent human presence on the lunar surface is one of the next major milestones in space exploration. This requires the development of robust extraplanetary construction technologies including structures and materials modeling and robotic systems. Previous robotic construction technology development has primarily focused on structural assembly, with significantly less focus on robotically performed outfitting tasks to instantiate subsystems providing power, data, life support, etc. These tasks involve manipulation of highly flexible elements, which are difficult to model, such as cable harnesses, ropes, and hoses. Robotic manipulation of DLOs, especially cable harnesses, is an active area of research as cable harnesses are essential for providing power and data to space assets. DLO models that can be used for robot manipulator trajectory generation are necessary for autonomous operation of lunar infrastructure. There are many proposed methods for modeling DLOs, and they primarily fall into three types: 1) discrete model-based, 2) continuum model-based, and 3) Neural Network-based. These types each have pros and cons, and the tradeoff between model accuracy and computational speed informs which type should be used. An understanding of this trade-off is imperative for real-time control of autonomous systems. High computational requirements reduce the speed of the model, making real-time control difficult, while accuracy is critical to preventing collisions. Discrete models, such as a mass-spring multibody representation, require relatively few calculations, and accuracy is directly tied to the step size of the discretization. Continuum models, such as a B-spline representation or a Cosserat rod model (a mix of continuous and discrete), are more informed of the structural properties of the cable and are much more accurate than a rigid body mass-spring model, but at significant computational cost. A Neural Network approach can provide an online solution with very few computational steps, but properly generating training data can be difficult and validation for an in-space application is not trivial. This paper explores the trade-off between different modeling approaches and compares accuracy and computational speed/complexity of the three types mentioned above. Model accuracy is evaluated using a cable in a static configuration. True cable shape is obtained using a depth camera for RGB images and point-cloud segmentation. The purpose of this experiment is to evaluate the trade-offs of different approaches to the DLO modeling problem. Understanding the tradeoffs between different cable modeling techniques paves the way for developing robotic control and planning architectures necessary for real-time manipulation of DLOs for lunar infrastructure outfitting. Real-time control is required for robotic systems to be able to actively manipulate a cable in a harsh environment where model and sensor errors compound, and environmental conditions can cause significant disturbances. Cable routing must be performed in areas with high density of objects/obstacles: through truss structures, near solar panels or mirror arrays, next to bundles of electrical equipment. Understanding the best way to plan and manipulate a cable without disrupting the environment or damaging the cable is imperative to robotic outfitting operations on the lunar surface.

Amy M Quartaro↗

Curvilinear Displacement Transfer Functions for Deformed Shape Predictions of Curved Structures Using Distributed Surface Strains

Curvilinear Displacement Transfer Functions were formulated for deformed shape predictions of different curved structures using surface strains. The embedded curved beam (depth-wise cross section of a curved structure along a surface strain-sensing line) was discretized into multiple small domains, with domain junctures matching the strain-sensing stations. Thus, the surface strain distribution can be described with a piecewise linear or a piecewise nonlinear function. The discrete approach enabled piecewise integrations of a curvature-strain differential equation for the embedded curved beam to yield closed-form Curvilinear Displacement Transfer Functions, which are written in terms of embedded curved-beam geometrical parameters and surface strains. By inputting the surface strain data, the Curvilinear Displacement Transfer Functions can transform surface strains into deflections along each embedded curved beam for mapping out the overall structural deformed shapes. The finite-element method was used to analytically generate the surface strains of the curved beams. The deformed shape prediction accuracies were then determined by comparing the theoretical deflections with the finite-element-generated deflections, which were used as yardsticks. By introducing the correction factors in simple mathematical forms, the Curvilinear Displacement Transfer Functions can be quite accurate for shape predictions of different curved-beam structures ranging from limit case of straight beam up to semicircular curved beam.

Ko, William L.↗

Development and application of a method for predicting rotor free wake positions and resulting rotor blade air loads. Volume 1: Model and results

Rotor wake geometries are predicted by a process similar to the startup of a rotor in a free stream. An array of discrete trailing and shed vortices is generated with vortex strengths corresponding to stepwise radial and azimuthal blade circulations. The array of shed and trailing vortices is limited to an arbitrary number of azimuthal steps behind each blade. The remainder of the wake model of each blade is an arbitrary number of trailing vortices. Vortex element end points were allowed to be transported by the resultant velocity of the free stream and vortex-induced velocities. Wake geometry, wake flow, and wake-induced velocity influence coefficients are generated by this program for use in the blade loads portion of the calculations. Blade loads computations include the effects of nonuniform inflow due to a free wake, nonlinear airfoil characteristics, and response of flexible blades to the applied loads. Computed wake flows and blade loads are compared with experimentally measured data. Predicted blade loads, response and shears and moments are obtained for a model rotor system having two independent rotors. The effects of advance ratio, vertical separation of rotors, different blade radius ratios, and different azimuthal spacing of the blades of one rotor with respect to the other are investigated.

Sadler, S. G.↗

Level-set topology optimization with PDE generated conformal meshes

This paper presents a level-set topology optimization approach that uses conformal meshes for the analysis of the displacement field. The structure’s boundary is represented by the iso-contour of a level-set field discretized on a fixed background design mesh. The conformal mesh is updated for each design iteration via a PDE based mesh morphing process that identifies the set of facets in the background mesh that are homeomorphic to the boundary and relaxes the homeomorphic mesh to conform to the structure’s boundary and ensure high element quality. The conformal mesh allows for a more accurate computation of the response versus density and some level-set based methods which interpolate material properties using the volume fraction. Numerical examples illustrate the proposed approach by optimizing linear-elastic two- and three-dimensional structures, wherein insight into the performance of the mesh morphing process is provided. The examples also highlight the scalability of the approach.

42 ENGINEERING↗

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗

Organic Analysis in Miller Range 090657 and Buckley Island 10933 CR2 Chondrites: Part 1 In-Situ Observation of Carbonaceous Material

Primitive carbonaceous chondrites contain a wide variety of organic material, ranging from soluble discrete molecules to insoluble unstructured kerogen-like component as well as structured nano-globules of macromolecular carbon. The relationship between the soluble organic molecules, macromolecular organic material, and host minerals are poorly understood. Due to the differences in extractability of soluble and insoluble organic materials, the analysis methods for each differ and are often performed independently. The combination of soluble and insoluble analyses, when performed concurrently, can provide a wider understanding on spatial distribution, and elemental, structural and isotopic composition of organic material in primitive meteorites. Furthermore, they can provide broader perspective on how extraterrestrial organic ma-terials potentially contributed to the synthesis of life's essential compounds such as amino acids, sugar acids, activated phosphates and nucleobases.

Cao, T.↗

Organic Analysis in the Miller Range 090657 CR2 Chondrite: Part 2 Amino Acid Analyses

Primitive carbonaceous chondrites contain a wide variety of organic material, ranging from soluble discrete molecules to insoluble, unstructured kerogen-like components, as well as structured nano-globules of macromolecular carbon. The relationship between the soluble organic molecules, macromolecular organic material, and host minerals are poorly understood. Due to the differences in extractability of soluble and insoluble organic materials, the analysis methods for each differ and are often performed independently. The combination of soluble and insoluble analyses, when performed concurrently, can provide a wider understanding of spatial distribution, and elemental, structural and isotopic composition of organic material in primitive meteorites. Using macroscale extraction and analysis techniques in combination with in situ microscale observation, we have been studying both insoluble and soluble organic material in the primitive CR2 chondrite Miller Range (MIL) 090657. In accompanying abstracts (Cao et al. and Messenger et al.) we discuss insoluble organic material in the samples. By performing the consortium studies, we aim to improve our understanding of the relationship between the meteorite minerals and the soluble and insoluble organic phases and to delineate which species formed within the meteorite and those that formed in nebular or presolar environments. In this abstract, we present the results of amino acid analyses of MIL 090657 by ultra performance liquid chromatography with fluorescence detection and quadrupole-time of flight mass spectrometry. Amino acids are of interest because they are essential to life on Earth, and because they are present in sufficient structural, enantiomeric and isotopic diversity to allow insights into early solar system chemical processes. Furthermore, these are among the most isotopically anomalous species, yet at least some fraction are thought to have formed by aqueously-mediated processes during parent body alteration.

Burton, A. S.↗

Interference Lattice-based Loop Nest Tilings for Stencil Computations

A common method for improving performance of stencil operations on structured multi-dimensional discretization grids is loop tiling. Tile shapes and sizes are usually determined heuristically, based on the size of the primary data cache. We provide a lower bound on the numbers of cache misses that must be incurred by any tiling, and a close achievable bound using a particular tiling based on the grid interference lattice. The latter tiling is used to derive highly efficient loop orderings. The total number of cache misses of a code is the sum of (necessary) cold misses and misses caused by elements being dropped from the cache between successive loads (replacement misses). Maximizing temporal locality is equivalent to minimizing replacement misses. Temporal locality of loop nests implementing stencil operations is optimized by tilings that avoid data conflicts. We divide the loop nest iteration space into conflict-free tiles, derived from the cache miss equation. The tiling involves the definition of the grid interference lattice an equivalence class of grid points whose images in main memory map to the same location in the cache-and the construction of a special basis for the lattice. Conflicts only occur on the boundaries of the tiles, unless the tiles are too thin. We show that the surface area of the tiles is bounded for grids of any dimensionality, and for caches of any associativity, provided the eccentricity of the fundamental parallelepiped (the tile spanned by the basis) of the lattice is bounded. Eccentricity is determined by two factors, aspect ratio and skewness. The aspect ratio of the parallelepiped can be bounded by appropriate array padding. The skewness can be bounded by the choice of a proper basis. Combining these two strategies ensures that pathologically thin tiles are avoided. They do not, however, minimize replacement misses per se. The reason is that tile visitation order influences the number of data conflicts on the tile boundaries. If two adjacent tiles are visited successively, there will be no replacement misses on the shared boundary. The iteration space may be covered with pencils larger than the size of the cache while avoiding data conflicts if the pencils are traversed by a scanning-face method. Replacement misses are incurred only on the boundaries of the pencils, and the number of misses is minimized by maximizing the volume of the scanning face, not the volume of the tile. We present an algorithm for constructing the most efficient scanning face for a given grid and stencil operator. In two dimensions it is based on a continued fraction algorithm. In three dimensions it follows Voronoi's successive minima algorithm. We show experimental results of using the scanning face, and compare with canonical loop orderings.

VanderWijngaart, Rob F.↗

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.↗

Modeling Malaria Transmission in Thailand and Indonesia

Malaria Modeling and Surveillance is a project in the NASA Applied Sciences Public Health Applications Program. The main objectives of this project are: 1) identification of the potential breeding sites for major vector species: 2) implementation of a malaria transmission model to identify they key factors that sustain or intensify malaria transmission; and 3) implementation of a risk algorithm to predict the occurrence of malaria and its transmission intensity. Remote sensing and GIs are the essential elements of this project. The NASA Earth science data sets used in this project include AVHRR Pathfinder, TRMM, MODIS, NSIPP and SIESIP. Textural-contextual classifications are used to identify small larval habitats. Neural network methods are used to model malaria cases as a function of precipitation, temperatures, humidity and vegetation. Hindcastings based on these environmental parameters have shown good agreement to epidemiological records. Examples for spatio-temporal modeling of malaria transmissions in Southeast Asia are given. Discrete event simulations were used for modeling the detailed interactions among the vector life cycle, sporogonic cycle and human infection cycle, under the explicit influences of selected extrinsic and intrinsic factors. The output of the model includes the individual infection status and the quantities normally observed in field studies, such as mosquito biting rates, sporozoite infection rates, gametocyte prevalence and incidence. Results are in good agreement with mosquito vector and human malaria data acquired by Coleman et al. over 4.5 years in Kong Mong Tha, a remote village in western Thailand. Application of our models is not restricted to Southeast Asia. The model and techniques are equally applicable to other regions of the world, when appropriate epidemiological and vector ecological parameters are used as input.

Kiang, Richard↗

Exploring constituent redistribution in irradiated U-19Pu-14Zr fuel via electron probe microanalysis

Here, the phenomena of constituent redistribution, wherein a previously homogeneous metallic fuel forms discrete, radially concentric compositional zones upon irradiation was investigated by examining an irradiated U-19Pu-14Zr fuel (where numbers represent wt. %) with a burnup of 11.5 at.% with electron probe microanalysis (EPMA) and quadruple inductively coupled plasma mass spectroscopy (Q-ICP-MS). EPMA-generated U, Pu, and Zr compositional data obtained from a diameter traverse of the sample was converted to mass and was used to: 1) compare the overall fuel element analysis results between the two methods, 2) determine the number of compositionally distinct zones forming as a result of constituent redistribution; and 3) quantify the post-irradiation loss or gain of U, Pu, and Zr atoms in each distinct compositional zone. Weight percent concentrations of U, Pu, and Zr for the overall cross section compare favorably between the two analytical methods, suggesting that the spatially resolved EPMA analysis complements bulk chemical analysis. Among the four identified compositional zones, post-irradiation quantification of U, Pu, and Zr elemental atom content changes shows that the quantity of U atoms lost from the innermost zone is slightly less than the quantity of U atoms gained by the middle two zones, and the quantity of Zr atoms lost from the high-U third zone is slightly less than is gained by the two innermost zones. Pu is lost from all four zones, although the innermost zone and the high-U third zone lose a significantly higher percentage (> 22 %) of their initial Pu atoms than the other two zones. For all three elements, EPMA cannot distinguish between atoms lost due to transport to a different zone from atoms lost due to nuclear processes; however, the insight gained from using this process can be used to experiment with new modeling techniques to predict constituent redistribution in U-Pu-Zr fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Leveraging a Neural Network-Enhanced Reproducing Kernel Particle Method for Multiphysics Degradation Modeling of Energy Storage Materials

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and resulting in reduced performance and service life. A coupled electro-chemo-mechanical reproducing kernel particle method (RKPM) formulation has been developed to analyze this system. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), pixel-based model construction by RKPM is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. Here, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1, 2] is introduced to effectively model damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RKPM is additionally used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022. [2] Baek, J., Chen, J. S., "A Neural Network-Based Enrichment of Reproducing Kernel Approximation for Modeling Brittle Fracture", Computer Methods in Applied Mechanics and Engineering Vol. 410, 116590, 2024.

degradation↗

A microstructural lattice model for strain oriented problems: A combined Monte Carlo finite element technique

A specialized, microstructural lattice model, termed MCFET for combined Monte Carlo Finite Element Technique, was developed which simulates microstructural evolution in material systems where modulated phases occur and the directionality of the modulation is influenced by internal and external stresses. In this approach, the microstructure is discretized onto a fine lattice. Each element in the lattice is labelled in accordance with its microstructural identity. Diffusion of material at elevated temperatures is simulated by allowing exchanges of neighboring elements if the exchange lowers the total energy of the system. A Monte Carlo approach is used to select the exchange site while the change in energy associated with stress fields is computed using a finite element technique. The MCFET analysis was validated by comparing this approach with a closed form, analytical method for stress assisted, shape changes of a single particle in an infinite matrix. Sample MCFET analytical for multiparticle problems were also run and in general the resulting microstructural changes associated with the application of an external stress are similar to that observed in Ni-Al-Cr alloys at elevated temperature.

Gayda, J.↗

A Posteriori Bounds for Linear-Functional Outputs of Crouzeix-Raviart Finite Element Discretizations of the Incompressible Stokes Problem

We present a finite element technique for the efficient generation of lower and upper bounds to outputs which are linear functionals of the solutions to the incompressible Stokes equations in two space dimensions; the finite element discretization is effected by Crouzeix-Raviart elements, the discontinuous pressure approximation of which is central to our approach. The bounds are based upon the construction of an augmented Lagrangian: the objective is a quadratic "energy" reformulation of the desired output; the constraints are the finite element equilibrium equations (including the incompressibility constraint), and the intersubdomain continuity conditions on velocity. Appeal to the dual max-min problem for appropriately chosen candidate Lagrange multipliers then yields inexpensive bounds for the output associated with a fine-mesh discretization; the Lagrange multipliers are generated by exploiting an associated coarse-mesh approximation. In addition to the requisite coarse-mesh calculations, the bound technique requires solution only of local subdomain Stokes problems on the fine-mesh. The method is illustrated for the Stokes equations, in which the outputs of interest are the flowrate past, and the lift force on, a body immersed in a channel.

Patera, Anthony T.↗

Progressive Hedging Decomposition for Solutions of Large-Scale Process Family Design Problems

Rapid, wide-scale deployment of green process systems, such as carbon capture or water desalination systems, is essential for combatting climate change. Methods relying on traditional design or modularity fail to capture the benefits of both economies of numbers and economies of scale. We have proposed process family design, which designs a family of processes simultaneously exploiting opportunities for common elements. In previous work, we explored different optimization formulations to solve this problem. In this work, we develop a decomposition approach to tackle larger problems efficiently. We solve a water desalination case study, which is too large to solve within a reasonable timeframe with the discretization formulation. We exploit the block angular structure of the discretization problem to decompose and solve using Progressive Hedging (PH). We use the open-source Python package mpi-sppy to execute PH which allows us to leverage parallelization and a HPC cluster to further improve solution time.

Stinchfield, Georgia↗

A Moving Discontinuous Galerkin Method with Interface Condition Enforcement for Reacting Hypersonic Flows

The necessity of enforcing conservation in computational elements or cells (element conservation) for discontinuous solutions is well understood and respected for solving conservation laws in computational fluid dynamics (CFD). In contrast, interface conservation, where the conservation across cell interfaces is enforced, is long ignored, and yet is also ruled and required by the underlying physics just like element conservation. Violation of the interface conservation across discontinuities is the root cause why an exact discontinuous solution can never be achieved in shock capturing methods. The interface conservation is examined and explored in this talk. Moving discontinuous Galerkin (MDG) finite element method with interface conservation enforcement (MDG-ICE) [1],[2] are then presented for solving compressible flow problems with discontinuities based on the observation that the interface conservation can only be satisfied, only when mesh interfaces are aligned with discontinuities. In the MDG-ICE formulation, both conservative quantities and grid geometry are considered as independent variables. A space-time DG formulation is used to solve the multi-material compressible Euler equations in the standard discontinuous solution space and the discrete grid geometry is solved using a variational formulation in a continuous space. A self-adaptive Levenberg-Marquardt method is utilized to solve the resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation. A number of numerical experiments for a variety of flow problems are conducted to assess the accuracy and performance of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to deliver the designed order of both h- and p-convergence even for discontinuous solutions, and detect all types of interfaces, via interface condition enforcement and satisfy, via grid movement, the compressible Euler equations and the associated interface condition.

Discontinuous Galekin↗

An upwind multigrid method for solving viscous flows on unstructured triangular meshes

A multigrid algorithm is combined with an upwind scheme for solving the two dimensional Reynolds averaged Navier-Stokes equations on triangular meshes resulting in an efficient, accurate code for solving complex flows around multiple bodies. The relaxation scheme uses a backward-Euler time difference and relaxes the resulting linear system using a red-black procedure. Roe's flux-splitting scheme is used to discretize convective and pressure terms, while a central difference is used for the diffusive terms. The multigrid scheme is demonstrated for several flows around single and multi-element airfoils, including inviscid, laminar, and turbulent flows. The results show an appreciable speed up of the scheme for inviscid and laminar flows, and dramatic increases in efficiency for turbulent cases, especially those on increasingly refined grids.

Bonhaus, Daryl Lawrence↗

Multilevel well modeling in aggregation-based nonlinear multigrid for multiphase flow in porous media

A full approximation scheme (FAS) nonlinear multigrid solver for two-phase flow and transport problems driven by wells with multiple perforations is developed here. It is an extension to our previous work on FAS solvers for diffusion and transport problems. The solver is applicable to discrete problems defined on unstructured grids as the coarsening algorithm is aggregation-based and algebraic. To construct coarse basis that can better capture the radial flow near wells, coarse grids in which perforated well cells are not near the coarse-element interface are desired. This is achieved by an aggregation algorithm proposed in this paper that makes use of the location of well cells in the cell-connectivity graph. Numerical examples in which the FAS solver is compared against Newton's method on benchmark problems are given. In particular, for a refined version of the SAIGUP model, the FAS solver is at least 35% faster than Newton's method for time steps with a CFL number greater than 10.

58 GEOSCIENCES↗