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.

SourceResearch connectionFrom yearThrough year

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

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Development of UMARC (University of Maryland Advanced Rotorcraft Code)

The University of Maryland Advanced Rotorcraft Code (UMARC) is a user-friendly, FEM-based comprehensive helicopter rotor simulation code of high numerical robustness and computational efficiency. UMARC formulates the rotor-fuselage equations using Hamilton's principle, and are discretized using finite elements in space and time. The FEM formulation allows the code to analyze a wide variety of rotor designs. Dynamic inflow modeling is used for unsteady wake inflow computations. Predicted stability, response, and blade-load data are validated with experimental data for several configurations, including representative articulated, hingeless, and bearingless rotors.

Bir, Gunjit↗

NASA NTRS1991

Turbulent flow calculations using unstructured and adaptive meshes

A method of efficiently computing turbulent compressible flow over complex two dimensional configurations is presented. The method makes use of fully unstructured meshes throughout the entire flow field, thus enabling the treatment of arbitrarily complex geometries and the use of adaptive meshing techniques throughout both viscous and inviscid regions of the flow field. Mesh generation is based on a locally mapped Delaunay technique in order to generate unstructured meshes with highly-stretched elements in the viscous regions. The flow equations are discretized using a finite element Navier-Stokes solver, and rapid convergence to steady-state is achieved using an unstructured multigrid algorithm. Turbulence modeling is performed using an inexpensive algebraic model, implemented for use on unstructured and adaptaive meshes. Compressible turbulent flow solutions about multiple-element airfoil geometries are computed and compared with expeirmental data.

Mavriplis, Dimitri J.↗

NASA NTRS1996

Simulation of Two-Fluid Flows by the Least-Squares Finite Element Method Using a Continuum Surface Tension Model

In this paper a numerical procedure for simulating two-fluid flows is presented. This procedure is based on the Volume of Fluid (VOF) method proposed by Hirt and Nichols and the continuum surface force (CSF) model developed by Brackbill, et al. In the VOF method fluids of different properties are identified through the use of a continuous field variable (color function). The color function assigns a unique constant (color) to each fluid. The interfaces between different fluids are distinct due to sharp gradients of the color function. The evolution of the interfaces is captured by solving the convective equation of the color function. The CSF model is used as a means to treat surface tension effect at the interfaces. Here a modified version of the CSF model, proposed by Jacqmin, is used to calculate the tension force. In the modified version, the force term is obtained by calculating the divergence of a stress tensor defined by the gradient of the color function. In its analytical form, this stress formulation is equivalent to the original CSF model. Numerically, however, the use of the stress formulation has some advantages over the original CSF model, as it bypasses the difficulty in approximating the curvatures of the interfaces. The least-squares finite element method (LSFEM) is used to discretize the governing equation systems. The LSFEM has proven to be effective in solving incompressible Navier-Stokes equations and pure convection equations, making it an ideal candidate for the present applications. The LSFEM handles all the equations in a unified manner without any additional special treatment such as upwinding or artificial dissipation. Various bench mark tests have been carried out for both two dimensional planar and axisymmetric flows, including a dam breaking, oscillating and stationary bubbles and a conical liquid sheet in a pressure swirl atomizer.

Wu, Jie↗

DOE OSTI2016

An energy-stable time-integrator for phase-field models

Here, we introduce a provably energy-stable time-integration method for general classes of phase-field models with polynomial potentials. We demonstrate how Taylor series expansions of the nonlinear terms present in the partial differential equations of these models can lead to expressions that guarantee energy-stability implicitly, which are second-order accurate in time. The spatial discretization relies on a mixed finite element formulation and isogeometric analysis. We also propose an adaptive time-stepping discretization that relies on a first-order backward approximation to give an error-estimator. This error estimator is accurate, robust, and does not require the computation of extra solutions to estimate the error. This methodology can be applied to any second-order accurate time-integration scheme. We present numerical examples in two and three spatial dimensions, which confirm the stability and robustness of the method. The implementation of the numerical schemes is done in PetIGA, a high-performance isogeometric analysis framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Modeling fission product diffusion in TRISO fuel particles with BISON

Diffusion of fission products in intact TRISO particles depends on particle geometry, fission product source rates, time, temperature, and temperature-dependent diffusion coefficients. Simulating this diffusion process requires models for source rates and diffusion coefficients, plus computation of the temperature field if not prescribed. In addition, simulation quality depends on discretization of the geometry, appropriate time stepping, and the accuracy of the solution method. In this paper, we explore the simulation of fission product diffusion in TRISO fuel particles using the finite element method via the fuel performance code Bison. Recent material model development has occurred in Bison for each material present in tri-structural isotropic (TRISO) fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. Also, new mesh generation and fission product release fraction capabilities have been added. Diffusion capabilities are shown to converge to the correct solution via formal verification tests. A large number of code benchmarking problems are also given, with good results, showing that Bison’s computed release fractions closely match those of other software tools. Finally, a significant validation effort is detailed in which fission product release, measured as part of the AGR-1 capsule experiments, is compared to Bison outputs. Bison outputs compare very well to the experimental data and to PARFUME results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Simulating Atmospheric Boundary Layer Turbulence with Nek5000/RS

We present large-eddy-simulation (LES) modeling approaches for the simulation of atmospheric boundary layer turbulence that are of direct relevance to wind energy production. In this report, we study a GABLS benchmark problem using high-order spectral element code Nek5000/RS, which is supported under the DOE’s Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms [1, 2]. We demonstrate our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. The model fidelity and scaling performance of Nek5000/RS on DOE’s leadership computing platforms in comparison to those of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, are discussed in [3].

17 WIND ENERGY↗

DOE OSTI2019

Surface complexation model of rare earth element adsorption onto bacterial surfaces with lanthanide binding tags

© 2019 Elsevier Ltd Lanthanide binding tags (LBTs) have been engineered onto the cell surface of E. coli to enhance biosorption and recovery of rare earth elements (REEs). The protonation behavior of the bacterial surfaces before and after LBT-display was compared by modeling acid-base titration data. A multiple discrete site, constant capacitance surface complexation model was constructed to examine rare earth (Tb) binding to cell surface functional groups, comparing wild type and LBT-engineered surfaces. Our acid-base titrations show similar pKa values between the two strains, suggesting induction of LBTs does not significantly alter cell surface protonation behavior. Tb sorption onto the wild type cell surface can be captured by a one-site carboxyl model. The LBT strain exhibited a higher metal loading that can be explained by the increase of sorption sites in the form of lanthanide binding tags. Furthermore, carboxyl site concentrations between wild type and LBT-induced cells were statistically indistinguishable. We thus attribute the engineered strains’ increase in Tb adsorption capacity and affinity to the addition of lanthanide binding tags to the cell surface. The Tb stability constant with the LBT site is two orders of magnitude higher than that with the carboxyl functional group. As a result, at low metal loading <10 μM, the Tb binding to the cell surface of the LBT-strain is controlled by the presence of high-affinity, but lower capacity, LBT sites. At higher metal loadings >10 μM, a more abundant but low affinity functional group becomes the main source of adsorption that results in an overall higher sorption capacity. This work demonstrates how surface complexation modeling can be implemented for bacterial surfaces engineered with a known protein tag to optimize REE recovery from fluids with variable pH and metal loadings.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2019

Benchmarking Mixed Mode Failure in Progressive Damage and Failure Analysis Methods

The verification and validation of progressive-damage-analysis finite element methods are difficult but critical tasks to undertake during their development. Verification exercises assess whether a predictive analysis tool produces results that are consistent with the fundamental concepts and assumptions of the tool under evaluation. Ideally, closed-form analytical solutions can be derived for which method verification results can be compared. Problems selected for computational tool verification are often simple and isolate individual features of the tool. In the case of progressive damage finite element methods, verifications should be performed to evaluate the ability of the model to predict the initiation of damage and its growth through the finite element mesh under a variety of conditions. Mabson et al. proposed a test case of a unidirectional, fiber-reinforced plate with a center crack subjected to tensile loads to evaluate matrix crack propagation predictions. The problem was modeled using the Abaqus Hashin continuum damage mechanics (CDM) model for fiber-reinforced composites. Different combinations of matrix strengths and element sizes were used in the simulations, and the results were compared to a closed-form solution based on linear elastic fracture mechanics (LEFM). It was determined that the Abaqus CDM model could predict the LEFM solution of Mode I cracks only when the finite element mesh density met specific requirements based on the material properties. This paper presents closed-form LEFM solutions for a center notch mixed mode (CNMM) verification problem. Parametric finite element analyses were developed using progressive damage analysis methods of both the Discrete Damage Mechanics (DDM) and CDM classes. The progressive damage analysis methods applied in the analyses of the CNMM problem include CompDam and the Floating Node Method. Analyses were conducted with various mode mixities and element sizes to verify that the damage models were working as intended and to identify any limits of applicability.

Leone, Frank↗

DOE OSTI2017

Experimental study and discrete element method simulation of Geldart Group A particles in a small-scale fluidized bed

Geldart Group A particles are of great importance in various chemical processes because of advantages such as ease of fluidization, large surface area, and many other unique properties. It is very challenging to model the fluidization behavior of such particles as widely reported in the literature. In this study, a pseudo-2D experimental column with a width of 5 cm, a height of 45 cm, and a depth of 0.32 cm was developed for detailed measurements of fluidized bed hydrodynamics of fine particles to facilitate the validation of computational fluid dynamic (CFD) modeling. The hydrodynamics of sieved FCC particles (Sauter mean diameter of 148 µm and density of 1300 kg/m3) and NETL-32D sorbents (Sauter mean diameter of 100 µm and density of 480 kg/m3) were investigated mainly through the visualization by a high-speed camera. Numerical simulations were then conducted by using NETL’s open source code MFIX-DEM. Both qualitative and quantitative information including bed expansion, bubble characteristics, and solid movement were compared between the numerical simulations and the experimental measurement. Furthermore, the cohesive van der Waals force was incorporated in the MFIX-DEM simulations and its influences on the flow hydrodynamics were studied.

42 ENGINEERING↗

NASA NTRS1982

Computation of the steady viscous flow over a tri-element 'augmentor wing' airfoil

The augmentor wing consists of a main airfoil with a slotted trailing edge for blowing, and two smaller aft airfoils which shroud the jet. This configuration has been modeled for numerical simulation by a novel discretization procedure which generates four separate grids: three surface-oriented airfoil grids and one outer free-stream grid. Grid lines and slopes are continuous across boundaries, so grid overlap at common boundaries provides boundary information without interpolation. A two-dimensional unsteady thin-layer Navier-Stokes code is used to calculate the flow for the no-blowing case at freestream Mach number = 0.7, Re = 12,600.000, and angles-of-incidence = 1.05 deg. Qualitative agreement with experimental data indicates the utility of this procedure in the analysis of multi-element configurations.

Lasinski, T. A.↗

DOE OSTI2022

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2019

Verification and Validation of a Three-Dimensional Composite Impact Model with Tabulated Input

A material model which incorporates several key capabilities which have been identified by the aerospace community as lacking in the composite impact models currently available in the commercial transient dynamic finite element code LS-DYNA® has been developed. The material model utilizes experimentally based tabulated input to define the evolution of plasticity and damage as opposed to specifying discrete input parameters (such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. The capability to account for the rate and temperature dependent deformation response of composites has also been incorporated into the material model. For the damage model, a strain equivalent formulation is utilized to allow for the uncoupling of the deformation and damage analyses. In the damage model, a diagonal damage tensor is defined to account for the directionally dependent variation of damage. However, the terms in the damage matrix are semi-coupled such that the damage in a particular coordinate direction is a function of the stresses and plastic strains in all of the coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure, which allows an arbitrarily shaped failure surface to be defined. A systematic series of validation and verification studies, at a variety of length scales ranging from single element simulations to simulations of a flat panel impact test, have been performed to fully exercise and evaluate the capabilities of the developed model.

Impact↗

DOE OSTI2024

A Simplified Method for Predicting Shaker Voltage in IMMATs

Impedance Matched Multi-Axis Tests (IMMATs) can replicate in-service vibration induced stress more accurately than single axis shaker table tests as they can better match a part’s operational boundary conditions and excite it in multiple degrees of freedom simultaneously. Here, the shakers used in IMMATs are less powerful than shaker tables, so shaker force limits can be exceeded during tests if they are not placed adequately for the desired environment. The ability to predict shaker voltage and force before performing a test is, therefore, helpful in selecting shaker locations so that their limits are not exceeded. In this study, electrodynamic shakers were modeled as discrete electromechanical systems, and the shaker parameters were chosen to match experimentally obtained acceleration/voltage frequency response functions (FRFs). These models were coupled to a finite element model of the device under test (DUT) via dynamic substructuring, and the substructured model was demonstrated to accurately predict shaker voltage as well as the error in reproducing the environment at multiple accelerometer locations. A simple method called the FRF Multiplication method, in which the FRF of the substructured system is approximated as the product of two separate FRFs of the shaker and DUT respectively, was proposed and applied to the same system, yielding similar voltage and error predictions to those obtained using substructuring. Simple case studies were presented to explore the applicability of the proposed method, and it was demonstrated to have similar accuracy to the substructuring method in a range of cases. Additionally, we showed that while it was not possible to derive a unique model of the shakers from acceleration/voltage FRFs alone, the models that could be obtained were sufficient to predict test error almost perfectly and shaker voltage with less than 40 percent error.

42 ENGINEERING↗

DOE OSTI2021

Atomic-scale mechanisms of void strengthening in tungsten

Neutron and heavy-ion irradiation of tungsten produces nanometer-size vacancy voids, gas-filled bubbles and dislocation loops. These defect features can affect mechanical properties and the impact can be quite significant because of their high density. Understanding the basic mechanisms of mechanical properties degradation is necessary for predicting radiation effects. Predictions can be made using discrete dislocation dynamics or/and finite element approaches which, however, need local interaction mechanisms as inputs. Such knowledge can be provided only by atomic-scale modeling. This paper reports the results of an extensive atomic-scale modeling study of the interactions between moving edge dislocations and voids in tungsten. Furthermore, the main focus is on the effects of the void size and ambient temperature. Critical resolved shear stress was calculated for voids up to 9 nm in diameter. Atomistic results are compared with the theoretical approach and with those obtained earlier for voids in body centered cubic (bcc) iron. An important role of the void surface has been revealed.

36 MATERIALS SCIENCE↗

NASA NTRS1984

Methods for analysis of cracks in three-dimensional solids

Analytical and numerical methods evaluating the stress-intensity factors for three-dimensional cracks in solids are presented, with reference to fatigue failure in aerospace structures. The exact solutions for embedded elliptical and circular cracks in infinite solids, and the approximate methods, including the finite-element, the boundary-integral equation, the line-spring models, and the mixed methods are discussed. Among the mixed methods, the superposition of analytical and finite element methods, the stress-difference, the discretization-error, the alternating, and the finite element-alternating methods are reviewed. Comparison of the stress-intensity factor solutions for some three-dimensional crack configurations showed good agreement. Thus, the choice of a particular method in evaluating the stress-intensity factor is limited only to the availability of resources and computer programs.

Raju, I. S.↗

NASA NTRS1988

Theoretical studies on flapped delta wings

The effects of leading edge flaps on the aerodynamic characteristics of a low aspect-ratio delta wing are studied theoretically. As an extension of the classical crossflow plane analysis and in order to include separated shear layers, an analogy between three dimensional steady conical and two dimensional unsteady self-similar flows is explored. This analogy provides a simple steady-unsteady relationship. The criteria for the validity of the steady-unsteady analogy are also examined. Two different theoretical techniques are used to represent the separated shear layers based on the steady-unsteady analogy, neglecting the trailing edge effect. In the first approach, each vortex system is represented by a pair of concentrated vortices connected to the separation points by straight feeding sheets. In the second approach, the vortex cloud method is adopted for simulating the flow field in the crossflow plane. The separated shear layers are replaced with a cloud of discrete vortices and the boundary element method is employed to represent the wing trace by a vorticity distribution. A simple merging scheme is used to model the core region of the vortical flow as a single vortex by imposing a restriction on the shear layer rotation angle. The results are compared with experiments and with results from 3-D panel calculations.

Oh, S.↗

NASA NTRS2000

Hypersonic Boundary-Layer Trip Development for Hyper-X

Boundary layer trip devices for the Hper-X forebody have been experimentally examined in several wind tunnels. Five different trip configurations were compared in three hypersonic facilities, the LaRC 20-Inch Mach 6 Air Tunnel, the LaRC 31 -Inch Mach 10 Air Tunnel, and in the HYPULSE Reflected Shock Tunnel at GASL. Heat transfer distributions, utilizing the phosphor thermography and thin-film techniques, shock system details, and surface streamline patterns were measured on a 0.333-scale model of the Hyper-X forebody. Parametric variations include angles-of-attack of 0-deg, 2-deg, and 4-deg; Reynolds numbers based on model length of 1.2 to 15.4 million: and inlet cowl door simulated in both open and closed positions. Comparisons of transition due to discrete roughness elements have led to the selection of a trip configuration for the Hyper-X Mach 7 flight vehicle.

Berry, Scott A.↗

NASA NTRS2001

Visualization of CFD Results in Immersive Virtual Environments

An object-oriented event-driven immersive virtual environment (VE) is described for the visualization of computational fluid dynamics (CFD) results. The VE incorporates the following types of primitive software objects: interface objects, support objects, geometric entities, and finite elements. The fluid domain is discretized using either a multi-block structured grid or an unstructured finite element mesh. The VE allows natural 'fly-through' visualization of the model, the CFD grid, and the model's surroundings. In order to help visualize the flow and its effects on the model, the VE incorporates the following objects: stream objects (lines, surface-restricted lines. ribbons. and volumes); colored surfaces; elevation surfaces; surface arrows; global and local iso-surfaces; vortex cores; and separation/attachment surfaces and lines. Most of these objects can be used for dynamically probing the flow. Particles and arrow animations can be displayed on top of stream objects. Primitive response quantities as well as derived quantities can be used. A recursive tree search algorithm is used for real-time point and value search in the CFD grid.

Wasfy, Tamer M.↗