Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discrete element models”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-J-149: Secondary Emission Detection for Improved Proton Relative Stopping Power Identification

Purpose: This research investigates application of secondary prompt gamma (PG) emission spectra, resulting from nuclear reactions induced by protons, to characterize tissue composition along the particle path. The objective of utilizing the intensity of discrete high-energy peaks of PG is to improve the accuracy of relative stopping power (RSP) values available for proton therapy treatment planning on a patient specific basis and to reduce uncertainty in dose depth calculations. Methods: In this research, MCNP6 was used to simulate PG emission spectra generated from proton induced nuclear reactions in medium of varying composition of carbon, oxygen, calcium and nitrogen, the predominant elements found in human tissue. The relative peak intensities at discrete energies predicted by MCNP6 were compared to the corresponding atomic composition of the medium. Results: The results have shown a good general agreement with experimentally measured values reported by other investigators. Unexpected divergence from experimental spectra was noted in the peak intensities for some cases depending on the source of the cross-section data when using compiled proton table libraries vs. physics models built into MCNP6. While the use of proton cross-section libraries is generally recommended when available, these libraries lack data for several less abundant isotopes. This limits the range of their applicability and forces the simulations to rely on physics models for reactions with natural atomic compositions. Conclusion: Current end-of-range proton imaging provides an average RSP for the total estimated track length. The accurate identification of tissue composition along the incident particle path using PG detection and characterization allows for improved determination of the tissue RSP on the local level. While this would allow for more accurate depth calculations resulting in tighter treatment margins, precise understanding of proton beam behavior in tissue of various compositions is necessary requiring detailed simulations with a high degree of accuracy.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1981

Stability of a diffuse linear pinch with axial boundaries

A formulation of the stability behavior of a finite-length pinch is presented. A general initial perturbation is expressed as a uniformly convergent sum over a complete discrete k set. A variational calculation is then performed, based on the energy principle, in which the end-boundary conditions appear as constraints. The requisite Lagrange multipliers mutually couple the elemental periodic excitations. The resulting extended form of delta-W still admits a proper second-variation treatment so that the minimization and stability considerations of Newcomb remain applicable. Comparison theorems are discussed as is the relevance of this end-effect model to the stability of solar coronal loops.

Einaudi, G.↗
NASA NTRS2022

Multi-Scale Modelling of the Bound Metal Deposition Manufacturing of Ti6Al4V

Nonlinear shrinkage of the metal part during manufacturing by bound metal deposition, both on the ground and under microgravity, is considered. A multi-scale physics-based approach is developed to address the problem. It spans timescales from atomistic dynamics on the order of nanoseconds to full-part shrinkage on the order of hours. This approach enables estimation of the key parameters of the problem, including the widths of grain boundaries, the coefficient of surface diffusion, the initial redistribution of particles during the debinding stage, the evolution of the microstructure from round particles to densely-packed grains, the corresponding changes in the total and chemical free energies, and the sintering stress. The method has been used to predict shrinkage at the levels of two particles, of the filament cross-section, of the sub-model, and of the whole green, brown, and metal parts.

Nonlinear shrinkage↗
DOE OSTI2013

Correlation measurements in nuclear {beta}-decay using traps and polarized low energy beams

Precision measurements in nuclear {beta}-decay provide sensitive means to test discrete symmetries in the weak interaction and to determine some of the fundamental constants in semi-leptonic decays, like the coupling of the lightest quarks to charged weak bosons. The main motivation of such measurements is to find deviations from Standard Model predictions as possible indications of new physics. In this contribution I will focus on two topics related to precision measurements in nuclear {beta}-decay: i) the determination of the V{sub ud} element of the Cabibbo-Kobayashi-Maskawa quark mixing matrix from nuclear mirror transitions and ii) the search for exotic scalar or tensor contributions from {beta}{nu} angular correlations. The purpose is to underline the role being played by experimental techniques based on the confinement of radioactive species with atom and ion traps as well as the plans to use low energy polarized beams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Design and Analysis of Multifidelity Finite Element Simulations

Abstract The numerical accuracy of finite element analysis (FEA) depends on the number of finite elements used in the discretization of the space, which can be varied using the mesh size. The larger the number of elements, the more accurate the results are. However, the computational cost increases with the number of elements. In current practice, the experimenter chooses a mesh size that is expected to produce a reasonably accurate result, and for which the computer simulation can be completed in a reasonable amount of time. Improvements to this approach have been proposed using multifidelity modeling by choosing two or three mesh sizes. However, mesh size is a continuous parameter, and therefore, multifidelity simulations can be performed easily by choosing a different value for the mesh size for each of the simulations. In this article, we develop a method to optimally find the mesh sizes for each simulation and satisfy the same time constraints as a single or a double mesh size experiment. A range of different mesh sizes used in the proposed method allows one to fit multifidelity models more reliably and predict the outcome when meshes approach infinitesimally small, which is impossible to achieve in actual simulations. We illustrate our approach using an analytical function and a cantilever beam finite element analysis experiment.

Engineering↗
NASA NTRS1995

Slot spiral antenna modeling using hybrid/mixed FE-BI technique

This report is concerned with the numerical simulation of the printed slot spiral antenna (SSA) using the finite element - boundary integral (FE-BI) analysis. It has been reported that the FE-BI technique is suited for modeling microstrip patch antennas of any shape, printed on layered planar structure or cylindrical platform and fed with a coaxial cable or a microstripline network underneath the radiating elements. However, direct application of the hybrid FE-BI technique to thin slot spirals requires excessive sampling rates to accurately simulate the geometry. To alleviate the meshing/modeling difficulties encountered with cavity-backed slot antennas, we describe a mixed finite element-boundary integral formulation. As in the past, the boundary integral is used to describe the radiation of the slot above the cavity. However the cavity is now modeled using a suitable mix of edge and node elements. The latter are used only at the aperture of the thin slot so that the nodes follow the center line of the slot. In this manner, regular size elements can be used regardless of the slot's width and any meshing restrictions are substantially relaxed. The proposed mixed element FE-BI formulation introduces three different computational regions and as expected this complicates the generation of the discrete system.

Gong, Jian↗
NASA NTRS2001

Equivalent-Continuum Modeling of Nano-Structured Materials

A method has been developed for modeling structure-property relationships of nano-structured materials. This method serves as a link between computational chemistry and solid mechanics by substituting discrete molecular structures with an equivalent-continuum model. It has been shown that this substitution may be accomplished by equating the vibrational potential energy of a nano-structured material with the strain energy of representative truss and continuum models. As an important example with direct application to the development and characterization of single-walled carbon nanotubes, the model has been applied to determine the effective continuum geometry of a graphene sheet. A representative volume element of the equivalent-continuum model has been developed with an effective thickness. This effective thickness has been shown to be similar to, but slightly smaller than, the interatomic spacing of graphite.

Odegard, Gregory M.↗
DOE OSTI2024

Hybridized Discontinuous Galerkin Methods for Computational Fluid Dynamics

Hybridizable Discontinuous Galerkin (HDG) methods hold promise for any applications with significant advection character, including thermal hydraulics in light water reactors and advanced reactor concepts and fluid models of plasmas in magnetic confinement fusion. Its features include natural upwinding, local element conservation, and extensibility to arbitrarily high order accuracy. In the last fiscal year we have implemented HDG in the Multiphysics Object-Oriented Simulation Environment (MOOSE). We developed a first-of-its-kind automatic static condensation system in MOOSE’s underlying finite element library libMesh which can condense out arbitrarily many internal variables. Finally, we developed the first preconditioner for HDG discretizations of the Navier-Stokes equations which shows robust performance across a wide range of problem sizes and Reynolds numbers. This preconditioner yields solution times that are equivalent to the fastest developed for industry standard finite volume methods. Moreover, the arbitrarily high-order nature of HDG makes it a prime candidate for acceleration via graphical processing units (GPUs). We believe these developments will hold significant importance for future DOE Nuclear Energy (NE) and Fusion Energy Science (FES) programs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Far-field Fano resonance in nanoring lattices modeled from extracted, point dipole polarizability

Coupling and extinction of light among particles representable as point dipoles can be characterized using the coupled dipole approximation (CDA). The analytic form for dipole polarizability of spheroidal particles supports rapid electrodynamic analysis of nanoparticle lattices using CDA. However, computational expense increases for complex shapes with non-analytical polarizabilities which require discrete dipole (DDA) or higher order approximations. This work shows fast CDA analysis of assembled nanorings is possible using a single dipole nanoring polarizability extrapolated from a DDA calculation by summing contributions from individual polarizable volume elements. Plasmon resonance wavelengths of nanorings obtained using extracted polarizabilities blueshift as wall dimensions-to-inner radius aspect ratio increases, consistent with published theory and experiment. Calculated far-field Fano resonance energy maximum and minimum wavelengths were within 1% of full volume element results. Considering polarizability allows a more complete physical picture of predicting plasmon resonance location than metal dielectric alone. This method reduces time required for calculation of diffractive coupling more than 40 000-fold in ordered nanoring systems for 400–1400 nm incident wavelengths. Extension of this technique beyond nanorings is possible for more complex shapes that exhibit dipolar or quadrupole radiation patterns.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Development of a discontinuous approach for modeling fluid flow in heterogeneous media using the numerical manifold method

Summary In the numerical modeling of fluid flow in heterogeneous geological media, large material contrasts associated with complexly intersected material interfaces are challenging, not only related to mesh discretization but also for the accurate realization of the corresponding boundary constraints. To address these challenges, we developed a discontinuous approach for modeling fluid flow in heterogeneous media using the numerical manifold method (NMM) and the Lagrange multiplier method (LMM) for modeling boundary constraints. The advantages of NMM include meshing efficiency with fixed mathematical grids (covers), the convenience of increasing the approximation precision, and the high integration precision provided by simplex integration. In this discontinuous approach, the elements intersected by material interfaces are divided into different elements and linked together using the LMM. We derive and compare different forms of LMMs and arrive at a new LMM that is efficient in terms of not requiring additional Lagrange multiplier topology, yet stringently derived by physical principles, and accurate in numerical performance. To demonstrate the accuracy and efficiency of the NMM with the developed LMM for boundary constraints, we simulate a number of verification and demonstration examples, involving a Dirichlet boundary condition and dense and intersected material interfaces. Last, we applied the developed model for modeling fluid flow in heterogeneous media with several material zones containing a fault and an opening. We show that the developed discontinuous approach is very suitable for modeling fluid flow in strongly heterogeneous media with good accuracy for large material contrasts, complex Dirichlet boundary conditions, or complexly intersected material interfaces. Copyright © 2015 John Wiley & Sons, Ltd.

Hu, Mengsu↗
DOE OSTI2019

OpenACC acceleration for the P N – P N - 2 algorithm in Nek5000

Due to its high performance and throughput capabilities, GPU-accelerated computing is becoming a popular technology in scientific computing, in particular using programming models such as CUDA and OpenACC. The main advantage with OpenACC is that it enables to simply port codes in their ‘‘original’’ form to GPU systems through compiler directives, thus allowing an incremental approach. An OpenACC implementation is applied to the CFD code Nek5000 for simulation of incompressible flows, based on the spectral-element method. The work follows up previous implementations and focuses now on the P N - P N-2 method for the spatial discretization of the Navier–Stokes equations. Performance results of the ported code show a speed-up of up to 3.1 on multi-GPU for a polynomial order N > 11.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1983

Stable feedback control of distributed parameter systems: Time and frequency domain conditions

Large space structures, or any mechanically flexible structures, are inherently distributed parameter systems (DPS) whose dynamics are modeled by partial, rather than ordinary, differential equations. Such DPS are described by operator equations on an infinite-dimensional Hilbert (or Banach) space. However, any feedback controller for such a DPS must be a finite-dimensional (and discrete-time) system in order to be implemented with on-line digital computers and a finite (small) number of actuators and sensors. There are many ways to synthesize such controllers; the Galerkin or finite-element approach is emphasized.

Balas, M. J.↗
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.↗
DOE OSTI2026

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↗
NASA NTRS1985

Geometric-optical Modeling of a Conifer Forest Canopy

The objective of this research is to explore how the geometry of trees in forest stands influences the reflectance of the forest as imaged from space. Most plant canopy modeling has viewed the canopy as an assemblage of plane-parallel layers on top of a soil surface. For these models, leaf angle distribution, leaf area index, and the angular transmittance and reflectance of leaves are the primary optical and geometric parameters. Such models are now sufficiently well developed to explain most of the variance in angular reflectance measurements observed from homogeneous plant canopies. However, forest canopies as imaged by airborne and spaceborne scanners exhibit considerable variance at quite a different scale. Brightness values vary strongly from one pixel to the next primarily as a function of the number of trees they contain. At this scale, the forest canopy is nonuniform and discontinuous. This research focuses on a discrete-element, geometric-optical view of the forest canopy.

Strahler, A. H.↗
NASA NTRS1992

Comparison of radiated noise from shrouded and unshrouded propellers

The ducted propeller in a free field is modeled using the finite element method. The generation, propagation, and radiation of sound from a ducted fan is described by the convened wave equation with volumetric body forces. Body forces are used to introduce the blade loading for rotating blades and stationary exit guide vanes. For an axisymmetric nacelle or shroud, the problem is formulated in cylindrical coordinates. For a specified angular harmonic, the angular coordinate is eliminated, resulting in a two-dimensional representation. A finite element discretization based on nine-node quadratic isoparametric elements is used.

Eversman, Walter↗
NASA NTRS1995

Resolution of the buoyancy in the 8-foot high temperature tunnel combustor

Currently, the 8-Foot High Temperature Tunnel (8-Ft. HTT) combustor produces a good profile at only one point (2000 psia and 3650 R with oxygen enrichment). Air is enriched with oxygen (liquid) so that the combustor product gas will contain the volumetric amount of oxygen normally found in air. The oxygen enriched air has a large fraction that is not reacted and flows through the outer periphery of the fuel injector. This ring of cold air in addition to the relatively cold walls of the combustor set up buoyancy forces that produce a segregation of relatively cool gases at the bottom of the combustor exit. The basic problem is to produce a test gas that has uniform properties at all combustor conditions. The combustor temperature may be as high as 3700 R or as low as 2000 R. Combustor pressures can be as high as 3500 psia (no oxygen enrichment) and as low as 600 psia. The segregation is most severe with oxygen enriched air, since its temperature is lower and its density is high. The combustor is lined with nickel 201 and can be operated at about 1600 R maximum. A global mixing process is desired that produces an acceptable profile of temperature, species, and velocity at the exit of the combustor. The ultimate goal is a temperature profile with about 100 R variance and about 2 percent variance in oxygen. The exit total temperature must not be lowered significantly by the mixing apparatus or mechanisms employed. If immersed bodies are used, they must also be kept very hot. All combustor wall modifications must be able to survive the heat and structural conditions of the varied operating conditions. Our approach to resolving this issue is being conducted in three stages: (1) Consider mixing exclusively, (2) Resolve the heat transfer concerns resulting from the chosen mixing strategy, and (3) Solve the material and structural problems resulting from stages (1) and (2). Since the 8-Ft. HTT is unavailable for experimentation, the study is conducted exclusively with computational fluid dynamic (CPD) codes (Fluent/Uns and Rampant 3.1) using unstructured grid through body fitted coordinates. Both CFD codes are general purpose Navier-Stokes solution packages that can solve integral conservation equations for conservation of mass, momentum and energy. The governing equations are discretized using a control-volume finite-element method on unstructured triangular 2-D grids. In the interest of time, a 3-D tetrahedral grid was used to check the 2-D results on one mixing strategy and the 2-D results were confirmed. Preliminary results indicate that excellent mixing can be achieved with a body placed in the center of the flowing hot fluid with a minor modification to the combustor wall similar to a model positioned in the test section of a wind tunnel. The concept here, is to create longitudinal vortices strong enough to bring both fluids into intimate contact with each other near the wall where the cold fluid resides. However, there may be a trade-off in the choice of mixing strategy, heat transfer and structural requirements. Currently, we are examining the more promising geometries for heat transfer concerns and developing strategies for the material of construction of the center body.

Loney, Norman W.↗
NASA NTRS2007

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↗