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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Uniform convergence of an upwind discontinuous Galerkin method for solving scaled discrete-ordinate radiative transfer equations with isotropic scattering

Here we present an error analysis for the discontinuous Galerkin (DG) method applied to the discrete-ordinate discretization of the steady-state radiative transfer equation with isotropic scattering. Under some mild assumptions, we show that the DG method converges uniformly with respect to a scaling parameter which characterizes the strength of scattering in the system. However, the rate is not optimal and can be polluted by the presence of boundary layers. In one-dimensional slab geometries, we demonstrate optimal convergence when boundary layers are not present and analyze a simple strategy for balance interior and boundary layer errors. Some numerical tests are also provided in this reduced setting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗
DOE OSTI2022

Memory-efficient emulation of physical tabular data using quadtree decomposition

Computationally expensive functions are sometimes replaced in simulations with an emulator that approx-imates the true function (e.g., equations of state, wavelength-dependent opacity, or composition-dependent materials properties). For functions that have a constrained domain of interest, this can be done by discretizing the domain and performing a local interpolation on the tabulated function values of each local domain. For these so-called tabular data methods, the method of discretizing the domain and mapping the input space to each subdomain can drastically influence the memory and computational costs of the emulator. This is especially true for functions that vary drastically in different regions. We present a method for domain discretization and mapping that utilizes quadtrees, which results in significant reductions in the size of the emulator with minimal increases to computational costs or loss of global accuracy. We apply our method to the electron-positron Helmholtz free energy equation of state and show over an order of magnitude reduction in memory costs for reasonable levels of numerical accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Vlasov simulations of electron-ion collision effects on damping of electron plasma waves

Collisional effects can play an essential role in the dynamics of plasma waves by setting a minimum damping rate and by interfering with wave-particle resonances. Kinetic simulations of the effects of electron-ion pitch angle scattering on Electron Plasma Waves (EPWs) are presented here. In particular, the effects of such collisions on the frequency and damping of small-amplitude EPWs for a range of collision rates and wave phase velocities are computed and compared with theory. Both the Vlasov simulations and linear kinetic theory find the direct contribution of electron-ion collisions to wave damping significantly reduced from that obtained through linearized fluid theory. To our knowledge, this simple result has not been published before. Simulations have been carried out using a grid-based (Vlasov) approach, based on a high-order conservative finite difference method for discretizing the Fokker-Planck equation describing the evolution of the electron distribution function. Details of the implementation of the collision operator within this framework are presented. Such a grid-based approach, which is not subject to numerical noise, is of particular interest for the accurate measurements of the wave damping rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

A flexible nonlinear diffusion acceleration method for the S N transport equations discretized with discontinuous finite elements

This paper presents a flexible nonlinear diffusion acceleration (NDA) method that discretizes both the S N transport equation and the diffusion equation using the discontinuous finite element method (DFEM). The method is flexible in that the diffusion equation can be discretized on a coarser mesh with the only restriction that it is nested within the transport mesh and the FEM shape function orders of the two equations can be different. The consistency of the transport and diffusion solutions at convergence is defined by using a projection operator mapping the transport into the diffusion FEM space. The diffusion weak form is based on the modified incomplete interior penalty (MIP) diffusion DFEM discretization that is extended by volumetric drift, interior face, and boundary closure terms. In contrast to commonly used coarse mesh finite difference (CMFD) methods, the presented NDA method uses a full FEM discretized diffusion equation for acceleration. Suitable projection and prolongation operators arise naturally from the FEM framework. Via Fourier analysis and numerical experiments for a one-group, fixed source problem the following properties of the NDA method are established for structured quadrilateral meshes: (1) the presented method is unconditionally stable and effective in the presence of mild material heterogeneities if the same mesh and identical shape functions either of the bilinear or biquadratic type are used, (2) the NDA method remains unconditionally stable in the presence of strong heterogeneities, (3) the NDA method with bilinear elements extends the range of effectiveness and stability by a factor of two when compared to CMFD if a coarser diffusion mesh is selected. In addition, the method is tested for solving the C5G7 multigroup, eigenvalue problem using coarse and fine mesh acceleration. Finally, while NDA does not offer an advantage over CMFD for fine mesh acceleration, it reduces the iteration count required for convergence by almost a factor of two in the case of coarse mesh acceleration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

CFD-DEM Modeling of Autothermal Pyrolysis of Corn Stover with a Coupled Particle- and Reactor-Scale Framework

Autothermal operation of fast pyrolysis is an efficient process-intensification technique wherein exothermic oxidation reactions are used to overcome the heat-transfer bottleneck of conventional pyrolysis. The development of accurate, reliable modeling toolsets is imperative to generating a deeper understanding of biomass autothermal pyrolysis systems to support scale-up and industrial deployment. This modeling effort describes the development of single-particle and reactor models which incorporate detailed reaction schemes and simultaneous exothermic oxidation reactions. The particle-scale model was parameterized for corn stover feedstock with particle morphology, density, ash content, and biopolymer composition, all of which impact the emergent conversion characteristics during pyrolysis. Results were then used to parameterize a reactor-scale autothermal pyrolysis model, which was developed using a coarse-grained computational fluid dynamic-discrete element method. The simulation results compared well with experimental results, with the predicted bio-oil, light gas, and biochar yield within 3.0 wt% of the experimental yields. Further analyses were performed to test the influence of equivalence ratio, biomass injection position, and particle size distribution on autothermal pyrolysis. The analysis of the physio-chemical properties of the fluid and solid phase inside the reactor and at the reactor outlet help reveal important process interactions of autothermal pyrolysis.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
DOE OSTI2015

A High-Order Nonlinear Diffusion Acceleration for the SN Equations Discretized with the Discontinuous FEM I: Theory and Numerical Results

Recently, nonlinear diffusion acceleration (NDA) methods have received widespread attention in the radiation transport community because of their inherent efficiency (Willert 2014, Wang 2014, Park 2012). However, NDA is in general only conditionally stable \cite{Larsen2014}, while diffusion synthetic acceleration methods (DSA) are unconditionally stable if suitable discretization methods for the transport and diffusion systems are chosen \cite{Larsen1982}. The presented work introduces a new generalized NDA method for solving the transport equation discretized with $S_N$ (discrete ordinates) and discontinuous Finite Elements. In contrast to common coarse-mesh finite difference methods (CMFD), where the diffusion system is discretized using constant shape functions, the presented NDA method discretizes the diffusion system using the same (higher order) finite element basis functions as those for the transport system. This paper focuses on the solution of the steady-state transport equation with the low-order diffusion equation employing the same mesh and the same spatial shape functions, but the NDA methods currently under development for the Rattle$S_N$ake radiation transport solver\cite{Wang2013} based on the MOOSE\cite{Gaston2009} framework is more general: it allows for the diffusion system to be discretized using an arbitrary group structure and arbitrary mesh nested in their $S_N$ system counterparts, and arbitrary combinations of FEM basis functions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

On the coupling of hyperbolic and parabolic systems: Analytical and numerical approach

The coupling of hyperbolic and parabolic systems is discussed in a domain Omega divided into two distinct subdomains omega(+) and omega(-). The main concern is to find the proper interface conditions to be fulfilled at the surface separating the two domains. Next, they are used in the numerical approximation of the problem. The justification of the interface conditions is based on a singular perturbation analysis, i.e., the hyperbolic system is rendered parabolic by adding a small artifical viscosity. As this goes to zero, the coupled parabolic-parabolic problem degenerates into the original one, yielding some conditions at the interface. These are taken as interface conditions for the hyperbolic-parabolic problem. Actually, two alternative sets of interface conditions are discussed according to whether the regularization procedure is variational or nonvariational. It is shown how these conditions can be used in the frame of a numerical approximation to the given problem. Furthermore, a method of resolution is discussed which alternates the resolution of the hyperbolic problem within omega(-) and of the parabolic one within omega(+). The spectral collocation method is proposed, as an example of space discretization (different methods could be used as well); both explicit and implicit time-advancing schemes are considered. The present study is a preliminary step toward the analysis of the coupling between Euler and Navier-Stokes equations for compressible flows.

Gastaldi, Fabio↗
DOE OSTI2024

A Second Moment Method for k -Eigenvalue Acceleration with Continuous Diffusion and Discontinuous Transport Discretizations

The second moment method is a linear acceleration technique that couples the transport equation to a diffusion equation with transport-dependent additive closures. The resulting low-order diffusion equation can be discretized independent of the transport discretization, unlike diffusion synthetic acceleration, and is symmetric positive definite, unlike quasidiffusion. While this method has been shown to be comparable to quasidiffusion in iterative performance for fixed source and time-dependent problems, it is largely unexplored as an eigenvalue problem acceleration scheme due to the belief that the resulting inhomogeneous source makes the problem ill posed. Recently, a preliminary feasibility study was performed on the second moment method for eigenvalue problems. The results suggested comparable performance to quasidiffusion and more robust performance than diffusion synthetic acceleration. This work extends the initial study to more realistic reactor problems using state-of-the-art discretization techniques. Finally, the results in this paper show that the second moment method is more computationally efficient than its alternatives on complex reactor problems with unstructured meshes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Method to estimate uncertainty associated with parcel size in coarse discrete particle simulation

Coarse grained particle methods significantly reduce the computation cost of large‐scale fluidized bed simulation by lumping many real particles into a computation parcel. This research provides a method to estimate the errors associated with parcel size in large‐scale fluidized bed simulations. This uncertainty is first quantified in small scale domains by comparing results of discrete particle method with that employing coarse parcels of different sizes. Then, this uncertainty is correlated with parcel size and simulation domains consisting of a simple homogeneous cooling system and more complex bubbling and circulating fluidized beds. These correlations allow us to accurately estimate the uncertainty in large‐scale fluidized beds based solely on data obtained in smaller systems. The ability to estimate model‐related uncertainty in larger systems makes this method relevant for industrial applications. © 2018 American Institute of Chemical Engineers AIChE J , 64: 2340–2350, 2018

Lu, Liqiang↗
DOE OSTI2017

An efficient and reliable predictive method for fluidized bed simulation

In past decades, the continuum approach was the only practical technique to simulate large-scale fluidized bed reactors because discrete approaches suffer from the cost of tracking huge numbers of particles and their collisions. This study significantly improved the computation speed of discrete particle methods in two steps: First, the time-driven hard-sphere (TDHS) algorithm with a larger time-step is proposed allowing a speedup of 20-60 times; second, the number of tracked particles is reduced by adopting the coarse-graining technique gaining an additional 2-3 orders of magnitude speedup of the simulations. A new velocity correction term was introduced and validated in TDHS to solve the over-packing issue in dense granular flow. The TDHS was then coupled with the coarse-graining technique to simulate a pilot-scale riser. The simulation results compared well with experiment data and proved that this new approach can be used for efficient and reliable simulations of large-scale fluidized bed systems.

42 ENGINEERING↗
NASA NTRS1980

Newton's method: A link between continuous and discrete solutions of nonlinear problems

Newton's method for nonlinear mechanics problems replaces the governing nonlinear equations by an iterative sequence of linear equations. When the linear equations are linear differential equations, the equations are usually solved by numerical methods. The iterative sequence in Newton's method can exhibit poor convergence properties when the nonlinear problem has multiple solutions for a fixed set of parameters, unless the iterative sequences are aimed at solving for each solution separately. The theory of the linear differential operators is often a better guide for solution strategies in applying Newton's method than the theory of linear algebra associated with the numerical analogs of the differential operators. In fact, the theory for the differential operators can suggest the choice of numerical linear operators. In this paper the method of variation of parameters from the theory of linear ordinary differential equations is examined in detail in the context of Newton's method to demonstrate how it might be used as a guide for numerical solutions.

Thurston, G. A.↗