Search NASA⌕ Search

SEARCH · Search NASA

Results for “discretization methods”

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

Investigation of particle–wall interaction in a pseudo-2D fluidized bed using CFD-DEM simulations

Here, we report on discrete element method simulations of a pseudo-two-dimensional (pseudo-2D) fluidized bed to investigate particle–wall interactions. Detailed information on macroscopic flow field variables, including solids pressure, granular temperature, and normal and tangential wall stresses are analyzed. The normal wall stress differs from the solids pressure because of the strong anisotropic flow behavior in the pseudo-2D system. A simple linear relationship exists between normal wall stress and solids pressure. In addition, an effective friction coefficient can be derived to characterize particle–wall flow interaction after evaluating the normal and tangential wall stresses. The effects of inter-particle and particle–wall friction coefficients are evaluated. Strong anisotropic flow behavior in the pseudo-2D system needs to be considered to validate the two-fluid model where the boundary condition is usually developed based on an isotropic assumption. The conclusion has been confirmed by simulation with different particle stiffnesses. Assumptions in the newly developed model for 2D simulation are further examined against the discrete element method simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗
NASA NTRS1977

Monte Carlo turbulence simulation

The paper describes turbulence simulation experiments based on the principles of control system theory, that is, the construction of a system characterized by a system function such that upon exciting the system with prescribed noise processes the output of the system is a realization of a random processing the desired statistical attributes of turbulence. An experimental autocorrelation of Jimsphere measurements of wind velocity was approximated to simulate turbulent wind. From the approximate autocorrelation function, the required system function is obtained, and a discrete time system is designed. Another method of simulation is to solve the convolution integral by filter techniques. Other methods include discrete Fourier simulation and self-similar simulation.

Fichtl, G. H.↗
NASA NTRS1990

What can formal methods offer to digital flight control systems design

Formal methods research begins to produce methods which will enable mathematic modeling of the physical behavior of digital hardware and software systems. The development of these methods directly supports the NASA mission of increasing the scope and effectiveness of flight system modeling capabilities. The conventional, continuous mathematics that is used extensively in modeling flight systems is not adequate for accurate modeling of digital systems. Therefore, the current practice of digital flight control system design has not had the benefits of extensive mathematical modeling which are common in other parts of flight system engineering. Formal methods research shows that by using discrete mathematics, very accurate modeling of digital systems is possible. These discrete modeling methods will bring the traditional benefits of modeling to digital hardware and hardware design. Sound reasoning about accurate mathematical models of flight control systems can be an important part of reducing risk of unsafe flight control.

Good, Donald I.↗
DOE OSTI2018

Modeling of rock scour using coupled 3-D discrete element and Lattice Boltzmann Methods

Copyright © 2018 ARMA, American Rock Mechanics Association. Scour of rock in unlined rock spillway channels is a critical issue facing many of the world’s dams. From a modeling point of view this poses a challenging and interesting problem that combines rock mechanics and hydraulics of turbulent flow. We analyze this interaction between the blocky rock mass and water by directly modeling the solid and fluid phases—the individual polyhedral blocks are modeled using the Discrete Element Method (DEM) while the water is modeled using the Lattice Boltzmann Method (LBM). The LBM mesh is entirely independent of the DEM discretization, making it possible to refine the LBM mesh such that transient and varied fluid pressures acting of the rock surface are directly modeled. This provides the capability to investigate the effect of water pressure inside the fractured rock mass, along potential sliding planes, and can be extended to rock falls and slides into standing bodies of water such as lakes and reservoirs. Herein we present preliminary results to demonstrate the capabilities of the methodology.

Gardner, MH↗
NASA NTRS1986

Sensitivity analysis of discrete structural systems

Methods for calculating sensitivity derivatives for discrete structural systems are surveyed, primarily covering literature published during the past two decades. Methods are described for calculating derivatives of static displacements and stresses, eigenvalues and eigenvectors, transient structural response, and derivatives of optimum structural designs with respect to problem parameters. The survey is focused on publications developed in nonstructural fields such as electronics, controls, and physical chemistry which are directly applicable to structural problems. Most notable among the nonstructural-based methods are the adjoint variable technique from control theory, and the Green's function and FAST methods from physical chemistry.

Adelman, H. M.↗
DOE OSTI2021

The divergence of nearby trajectories in soft-sphere DEM

The n-body instability is investigated with the soft-sphere discrete element method. The divergence of nearby trajectories is quantified by the dynamical memory time. Using the inverse proportionality between the dynamical memory time and the largest Lyapunov exponent, the soft-sphere discrete element method results are compared to previous hard-sphere molecular dynamics data for the first time. Good agreement is observed at low concentrations and the degree of instability is shown to increase asymptotically with increasing spring stiffness. At particle concentrations above 30%, the soft-sphere Lyapunov exponents increase faster than the corresponding hard-sphere data. Finally, this paper concludes with a demonstration of how this case study may be used in conjunction with regression testing and code verification activities.

42 ENGINEERING↗
DOE OSTI2020

An error-controlled adaptive time-stepping method for particle advancement in coupled CFD-DEM simulations

Coupled Computational-Fluid-Dynamics (CFD) and Discrete-Element-Method (DEM) models provide an accurate description of multiphase physical systems where a solid granular particle phase exists in an underlying gaseous continuous medium. The time integration of the granular phase in these simulations is typically handled using an explicit scheme with a constant time-step among all particles that is invariant in time to resolve inter-particle collisions. A locally third-order accurate adaptive time integration technique for particles that employs an embedded locally second-order scheme for error determination is presented in this work. The particle time-step size is dynamically adapted based on solution error, thus leading to significant savings in computational time. The efficacy of our scheme is quantified using four test cases of varying complexity (binary collision, homogeneous cooling system, fluidized bed and hopper discharge). The adaptive time-stepping method exhibits improved performance (~ 2–3 times in most of the cases studied) compared to three commonly used non-adaptive time-step methods (first-order Euler-explicit, second-order Adams-Bashforth and third-order Runge-Kutta schemes), while maintaining the same level of accuracy and parallel scalability.

42 ENGINEERING↗
DOE OSTI2015

Year End Progress Report on Rattlesnake Improvements

Rattlesnake is a MOOSE-based radiation transport application developed at INL to support modern multi-physics simulations. At the beginning of the last year, Rattlesnake was able to perform steady-state, transient and eigenvalue calculations for the multigroup radiation transport equations. Various discretization schemes, including continuous finite element method (FEM) with discrete ordinates method (SN) and spherical harmonics expansion method (PN) for the self-adjoint angular flux (SAAF) formulation, continuous FEM (CFEM) with SN for the least square (LS) formulation, diffusion approximation with CFEM and discontinuous FEM (DFEM), have been implemented. A separate toolkit, YAKXS, for multigroup cross section management was developed to support Rattlesnake calculations with feedback both from changes in the field variables, such as fuel temperature, coolant density, and etc., and in isotope inventory. The framework for doing nonlinear diffusion acceleration (NDA) within Rattlesnake has been set up, and both NDA calculations with SAAF-SN-CFEM scheme and Monte Carlo with OpenMC have been performed. It was also used for coupling BISON and RELAP-7 for the full-core multiphysics simulations. Within the last fiscal year, significant improvements have been made in Rattlesnake. Rattlesnake development was migrated into our internal GITLAB development environment at the end of year 2014. Since then total 369 merge requests has been accepted into Rattlesnake. It is noted that the MOOSE framework that Rattlesnake is based on is under continuous developments. Improvements made in MOOSE can improve the Rattlesnake. It is acknowledged that MOOSE developers spent efforts on patching Rattlesnake for the improvements made on the framework side. This report will not cover the code restructuring for better readability and modularity and documentation improvements, which we have spent tremendous effort on. It only details some of improvements in the following sections.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1980

Hybrid vortex method for lifting surfaces with free-vortex flow

A Nonlinear Hybrid Vortex method (NHV-method) has been developed for predicting the aerodynamic characteristics of wings exhibiting leading- and side-edge separations. This method alleviates the drawbacks of the Nonlinear Discrete Vortex method (NDV-method, also known as the multiple line vortex method.) The NHV-method combines continuous-vorticity and vortex-line representations of the wing and its separated free shear layers. Continuous vorticity is used in the near-field calculations, while discrete vortex-lines are used in the far-field calculations. The wing and its free shear layers are divided into quadrilateral vortex panels having second-order vorticity distributions. The aerodynamic boundary conditions and continuity of the vorticity distributions are satisfied at certain nodal points on the vortex panels. An iterative technique is used to satisfy these conditions in order to obtain the vorticity distribution and the wake shape. Distributed and total aerodynamic loads are then calculated.

Kandil, O. A.↗
NASA NTRS1978

Three dimensional steady and unsteady asymmetric flow past wings of arbitrary planforms

The nonlinear discrete vortex method is extended to treat the problem of asymmetric flows past a wing with leading edge separation, including steady and unsteady flows. The problem is formulated in terms of a body fixed frame of reference and the nonlinear-discrete vortex method is modified accordingly. Although the method is general, only examples of flows past delta wings are presented due to the availability of experimental data as well as approximate theories. Comparison of results with experimental results for a delta wing undergoing a steady rolling motion at zero angle of attack demonstrate the superiority of the present method over existing approximate theories in obtaining highly accurate loads. Numerical results for yawed wings at large angles of attack are also presented. In all cases, total load coefficients, pressure distributions, and shapes of the free vortex sheets are shown.

Kandil, O. A.↗
NASA NTRS1977

Three dimensional steady and unsteady asymmetric flow past wings of arbitrary planforms

The nonlinear discrete vortex method was extended to treat the problem of asymmetric flows past a wing with leading-edge separation, including steady and unsteady flows. The problem was formulated in terms of a body-fixed frame of reference, and the nonlinear discrete vortex method was modified accordingly. Only examples of flows past delta wings are presented. Comparison of these results with experimental results for a delta wing undergoing a steady rolling motion at zero angle of attack demonstrates the superiority of the present method in obtaining highly accurate loads. Numerical results for yawed wings at large angles of attack are also presented. In all cases, total load coefficients, pressure distributions and shapes of the free-vortex sheets are shown.

Kandil, O. A.↗
DOE OSTI2025

Gaussian integral method for void fraction

Here, a novel method, the Gaussian Integral Method (GIM), is presented for calculating void fractions in Computational Fluid Dynamics–Discrete Element Method (CFD-DEM) simulations. GIM is versatile and applicable to various grid types, including structured and unstructured polyhedral meshes, without requiring special boundary treatments. An optimization technique is introduced to make GIM independent of grid resolution and type. The method is validated against experimental data from a fluidized bed, demonstrating that GIM produces realistic simulations closely resembling experimental observations. Additionally, unstructured polyhedral grids using GIM outperform structured grids of equivalent resolution, yielding results more aligned with experimental data. The gradient of the void fraction is computed in the CFD solver and utilized in the DEM solver for precise estimation at particle locations. Overall, GIM provides an effective solution for void fraction calculations in particulate media simulations with complex geometries, enhancing the accuracy and applicability of CFD-DEM simulations for industrial processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2011

Analysis of the Radiative Transfer Equation with Highly Asymmetric Phase Function

This paper considers a scalar radiative transfer problem with high scattering anisotropy, Two computational methods are presented based on decomposition of the diffuse light field into a regular and anisotropic part. The first algorithm (DOMAS) singles out the anisotropic radiance in the forward scattering peak using the Small-Angle Modification of RTE. The second algorithm (DOM2+) separates the single scattering radiance as an anisotropic part, which largely defines the fine detail of the total radiance in the backscattering directions. In both cases, the anisotropic part is represented analytically. With anisotropy subtraction, the regular part of the signal. which requires a numerical solution, is essentially smoothed as a function of angles. Further, the transport equation is obtained for the regular part that contains an additional source function from the anisotropic part of the signal. This equation is solved with the discrete ordinates method. A conducted numerical analysis of this work showed that algorithm DOMAS has a strong advantage as compared to the standard discrete ordinates method for simulation of the radiance transmission, and DOM2 + is the best of the three for the reflection computations. Both algorithms offer at least a factor of three acceleration of convergence of the azimuthal series for highly anisotropic phase functions.

Korkin, Sergey V.↗
NASA NTRS1986

Ideal shrinking and expansion of discrete sequences

Ideal methods are described for shrinking or expanding a discrete sequence, image, or image sequence. The methods are ideal in the sense that they preserve the frequency spectrum of the input up to the Nyquist limit of the input or output, whichever is smaller. Fast implementations that make use of the discrete Fourier transform or the discrete Hartley transform are described. The techniques lead to a new multiresolution image pyramid.

Watson, Andrew B.↗
NASA NTRS1985

A new approach to flow problems past a screen

A new method to simulate flows past a two-dimensional porous plate is proposed. The method is a simple application of a discrete vortex method which has been well established for flows past an inclined flat plate. In this method, a flow past a porous plate is described as superposition of the velocity potential of a uniform flow on that of a flow past a flat plate. The effect of porosity is replaced by that of a mass flow rate of the approaching freestream passing through the plate. Calculated flow features are compared with experiments. It is shown that the present method is effective in the simulation of the flow past a porous plate.

Inoue, O.↗
NASA NTRS2013

Inverse Modeling of Texas NOx Emissions Using Space-Based and Ground-Based NO2 Observations

Inverse modeling of nitrogen oxide (NOx) emissions using satellite-based NO2 observations has become more prevalent in recent years, but has rarely been applied to regulatory modeling at regional scales. In this study, OMI satellite observations of NO2 column densities are used to conduct inverse modeling of NOx emission inventories for two Texas State Implementation Plan (SIP) modeling episodes. Addition of lightning, aircraft, and soil NOx emissions to the regulatory inventory narrowed but did not close the gap between modeled and satellite observed NO2 over rural regions. Satellitebased top-down emission inventories are created with the regional Comprehensive Air Quality Model with extensions (CAMx) using two techniques: the direct scaling method and discrete Kalman filter (DKF) with Decoupled Direct Method (DDM) sensitivity analysis. The simulations with satellite-inverted inventories are compared to the modeling results using the a priori inventory as well as an inventory created by a ground-level NO2 based DKF inversion. The DKF inversions yield conflicting results: the satellite based inversion scales up the a priori NOx emissions in most regions by factors of 1.02 to 1.84, leading to 3-55% increase in modeled NO2 column densities and 1-7 ppb increase in ground 8 h ozone concentrations, while the ground-based inversion indicates the a priori NOx emissions should be scaled by factors of 0.34 to 0.57 in each region. However, none of the inversions improve the model performance in simulating aircraft-observed NO2 or ground-level ozone (O3) concentrations.

Nitrogen dioxide↗
DOE OSTI2018

A three-dimensional meso-scale approach to concrete fracture based on combined DEM with X-ray μCT images

Highlights: • We investigate fracture in concrete beams subjected to 3-point bending. • We use a discrete element method (DEM) to capture fracture. • Concrete is described as 4-phase material. • Evolution of micro-structural phenomena is investigated. • The 3D results are carefully compared with the 2D ones. - Abstract: The paper shows three-dimensional (3D) numerical mesoscopic results on progressive concrete fracture at aggregate level in a notched concrete beam under bending. A classical particle discrete element method was used to capture fracture in concrete. Heterogeneity was taken into account with a 4-phase concrete description: aggregate, cement matrix, macro-voids and interfacial transitional zones. Both the real shape and location of aggregate particles in concrete were assumed in DEM analyses based on X-ray μCT scans. A satisfactory agreement in terms of the vertical force versus crack mouth opening displacement evolution and crack geometry was achieved between DEM analyses and own laboratory tests. The evolution of broken contacts, contact normal forces, particle rotations, internal energies was studied at the aggregate level. The 3D results were also carefully compared with the 2D ones. The 3D model proved to be a promising simulation tool in predicting concrete fracture at the mesoscopic and macroscopic level.

36 MATERIALS SCIENCE↗