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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A family of independent Variable Eddington Factor methods with efficient preconditioned iterative solvers

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2015

2nd-Order CESE Results For C1.4: Vortex Transport by Uniform Flow

The Conservation Element and Solution Element (CESE) method was used as implemented in the NASA research code ez4d. The CESE method is a time accurate formulation with flux-conservation in both space and time. The method treats the discretized derivatives of space and time identically and while the 2nd-order accurate version was used, high-order versions exist, the 2nd-order accurate version was used. In regards to the ez4d code, it is an unstructured Navier-Stokes solver coded in C++ with serial and parallel versions available. As part of its architecture, ez4d has the capability to utilize multi-thread and Messaging Passage Interface (MPI) for parallel runs.

Space-time CE/SE method↗
NASA NTRS2010

Modeling and Characterization of Near-Crack-Tip Plasticity from Micro- to Nano-Scales

Methodologies for understanding the plastic deformation mechanisms related to crack propagation at the nano-, meso- and micro-length scales are being developed. These efforts include the development and application of several computational methods including atomistic simulation, discrete dislocation plasticity, strain gradient plasticity and crystal plasticity; and experimental methods including electron backscattered diffraction and video image correlation. Additionally, methodologies for multi-scale modeling and characterization that can be used to bridge the relevant length scales from nanometers to millimeters are being developed. The paper focuses on the discussion of newly developed methodologies in these areas and their application to understanding damage processes in aluminum and its alloys.

Glaessgen, Edward H.↗
NASA NTRS2011

Modeling and Characterization of Near-Crack-Tip Plasticity from Micro- to Nano-Scales

Methodologies for understanding the plastic deformation mechanisms related 10 crack propagation at the nano, meso- and micro-length scales are being developed. These efforts include the development and application of several computational methods including atomistic simulation, discrete dislocation plasticity, strain gradient plasticity and crystal plasticity; and experimental methods including electron backscattered diffraction and video image correlation. Additionally, methodologies for multi-scale modeling and characterization that can be used to bridge the relevant length scales from nanometers to millimeters are being developed. The paper focuses on the discussion of newly developed methodologies in these areas and their application to understanding damage processes in aluminum and its alloys.

Glaessgen, Edward H.↗
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

SAM Theory Manual

The System Analysis Module (SAM) is an advanced and modern system analysis tool being developed at Argonne National Laboratory under the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. SAM development aims for advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability for reactor transient analyses. To facilitate the code development, SAM utilizes an object-oriented application framework (MOOSE), and its underlying meshing and finite-element library (libMesh) and linear and non-linear solvers (PETSc), to leverage modern advanced software environments and numerical methods. SAM focuses on modeling advanced reactor concepts such as SFRs (sodium fast reactors), LFRs (lead-cooled fast reactors), and FHRs (fluoride-salt-cooled high temperature reactors) or MSRs (molten salt reactors). These advanced concepts are distinguished from light-water reactors in their use of single-phase, low-pressure, high-temperature, and low Prandtl number (sodium and lead) coolants. As a new code development, the initial effort has been focused on modeling and simulation capabilities of heat transfer and single-phase fluid dynamics responses in Sodium-cooled Fast Reactor (SFR) systems. The system-level simulation capabilities of fluid flow and heat transfer in general engineering systems and typical SFRs have been verified and validated. This document provides the theoretical and technical basis of the code to help users understand the underlying physical models (such as governing equations, closure models, and component models), system modeling approaches, numerical discretization and solution methods, and the overall capabilities in SAM. As the code is still under ongoing development, this SAM Theory Manual will be updated periodically to keep it consistent with the state of the development.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

A CFD model for biomass combustion in a packed bed furnace

Climate change has now become an important issue which is affecting environment and people around the world. Global warming is the main reason of climate change which is increasing day by day due to the growing demand of energy in developed countries. Use of renewable energy is now an established technique to decrease the adverse effect of global warming. Biomass is a widely accessible renewable energy source which reduces CO{sub 2} emissions for producing thermal energy or electricity. But the combustion of biomass is complex due its large variations and physical structures. Packed bed or fixed bed combustion is the most common method for the energy conversion of biomass. Experimental investigation of packed bed biomass combustion is difficult as the data collection inside the bed is challenging. CFD simulation of these combustion systems can be helpful to investigate different operational conditions and to evaluate the local values inside the investigation area. Available CFD codes can model the gas phase combustion but it can’t model the solid phase of biomass conversion. In this work, a complete three-dimensional CFD model is presented for numerical investigation of packed bed biomass combustion. The model describes the solid phase along with the interface between solid and gas phase. It also includes the bed shrinkage due to the continuous movement of the bed during solid fuel combustion. Several variables are employed to represent different parameters of solid mass. Packed bed is considered as a porous bed and User Defined Functions (UDFs) platform is used to introduce solid phase user defined variables in the CFD. Modified standard discrete transfer radiation method (DTRM) is applied to model the radiation heat transfer. Preliminary results of gas phase velocity and pressure drop over packed bed have been shown. The model can be useful for investigation of movement of the packed bed during solid fuel combustion.

09 BIOMASS FUELS↗
NASA NTRS1989

On equivalence of discrete-discrete and continuum-discrete design sensitivity analysis

Developments in design sensitivity analysis (DSA) method have been made using two fundamentally different approaches as shown. In the first approach, a discretized structural finite element model is used to carry out DSA. There are three different methods in the discrete DSA approach: finite difference, semi-analytical, and analytical methods. The finite difference method is a popular one due to its simplicity, but a serious shortcoming of the method is the uncertainty in the choice of a perturbation step size of design variables. In the semi-analytical method, the derivatives of stiffness matrix is computed by finite differences, whereas in the analytical method, the derivatives are obtained analytically. For the shape design variable, computation of analytical derivative of stiffness matrix is quite costly. Because of this, the semi-analytical method is a popular choice in discrete shape DSA approach. However, recently, Barthelemy and Haftka presented that the semi-analytical method can have serious accuracy problems for shape design variables in structures modeled by beam, plate, truss, frame, and solid elements. They found that accuracy problems occur even for a simple cantilever beam. In the second approach, a continuum model of the structure is used to carry out DSA.

Choi, Kyung K.↗
DOE OSTI2013

Thick diffusion limit boundary layer test problems

We develop two simple test problems that quantify the behavior of computational transport solutions in the presence of boundary layers that are not resolved by the spatial grid. In particular we study the quantitative effects of 'contamination' terms that, according to previous asymptotic analyses, may have a detrimental effect on the solutions obtained by both discontinuous finite element (DFEM) and characteristic-method (CM) spatial discretizations, at least for boundary layers caused by azimuthally asymmetric incident intensities. Few numerical results have illustrated the effects of this contamination, and none have quantified it to our knowledge. Our test problems use leading-order analytic solutions that should be equal to zero in the problem interior, which means the observed interior solution is the error introduced by the contamination terms. Results from DFEM solutions demonstrate that the contamination terms can cause error propagation into the problem interior for both orthogonal and non-orthogonal grids, and that this error is much worse for non-orthogonal grids. This behavior is consistent with the predictions of previous analyses. We conclude that these boundary layer test problems and their variants are useful tools for the study of errors that are introduced by unresolved boundary layers in diffusive transport problems. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Extracting Hydrocarbon from Shale: An Investigation of the Factors That Influence the Decline and the Tail of the Production Curve

Understanding physical processes that control the long-term production of hydrocarbon from shale formations is important for both predicting the yield and increasing it. In this work, we explore the processes that could control the tail of the production curve by using a discrete fracture network method to calculate the total travel time from the rock matrix to small-scale fractures to the primary hydraulic fracture network. The factors investigated include matrix diffusion, extent of the small-scale fracture zone (or tributary fracture zone/TFZ) consisting of natural, reactivated and induced fractures, and the percentage of free hydrocarbon in the primary fracture network. Individual and combined parameter spaces are explored for each of these to understand the limits of these parameters as well as any systematic correlations between pairs of parameters. Although recent studies have shown that the matrix diffusion in virgin shale influences the production tail only after nearly 20 years, we demonstrate that matrix diffusion in the region of the TFZ significantly impacts production within the first year itself. Additionally, we found that the depth of TFZ fracturing region had no effect on the shape of the production curves although the total mass of the hydrocarbon produced increases with the depth. We also show that one can fit the production data using a site-specific set of parameters representing the diffusion in the TFZ, depth of the TFZ and the free hydrocarbon in the large-scale fractures.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The role of rotation-vibration coupling in symmetric and asymmetric isotopomers of ozone

A theoretical framework and a computer code (SpectrumSDT) are developed for accurate calculations of coupled rotational–vibrational states in triatomic molecules using hyper-spherical coordinates and taking into account the Coriolis coupling effect. Concise final formulas are derived for the construction of the Hamiltonian matrix using an efficient combination of the variational basis representation and discrete variable representation methods with locally optimized basis sets and grids. First, the new code is tested by comparing its results with those of the APH3D program of Kendrick et al. [Kendrick, Pack, Walker, and Hayes, J. Chem. Phys. 110, 6673 (1999)]. Then, accurate calculations of the rovibrational spectra are carried out for doubly substituted symmetric ( 18 O 16 O 18 O) and asymmetric ( 18 O 18 O 16 O) ozone isotopomers for the total angular momentum up to J = 5. Together with similar data recently reported for the singly substituted symmetric ( 16 O 18 O 16 O) and asymmetric ( 16 O 16 O 18 O) ozone isotopomers, these calculations quantify the role of the Coriolis coupling effect in the large mass-independent isotopic enrichment of ozone, observed in both laboratory experiments and the atmosphere of the Earth. It is found that the Coriolis effect in ozone is relatively small, as evidenced by deviations of its rotational constants from the symmetric-top-rotor behavior, magnitudes of parity splittings (Λ-doubling), and ratios of rovibrational partition functions for asymmetric vs symmetric ozone molecules. It is concluded that all of these characteristics are influenced by the isotopic masses as much as they are influenced by the overall symmetry of the molecule. It is therefore unlikely that the Coriolis coupling effect could be responsible for symmetry-driven mass-independent fractionation of oxygen isotopes in ozone.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Can Electrical Resistance Tomography be used for imaging unsaturated moisture flow in cement-based materials with discrete cracks?

Previously, it has been shown that Electrical Resistance Tomography (ERT) can be used for monitoring moisture flow in undamaged cement-based materials. In this work, we investigate whether ERT could be used for imaging three-dimensional (3D) unsaturated moisture flow in cement-based materials that contain discrete cracks. Novel computational methods based on the so-called absolute imaging framework are developed and used in ERT image reconstructions, aiming at a better tolerance of the reconstructed images with respect to the complexity of the conductivity distribution in cracked material. ERT is first tested using specimens with physically simulated cracks of known geometries, and corroborated with numerical simulations of unsaturated moisture flow. Next, specimens with loading-induced cracks are imaged; here, ERT reconstructions are evaluated qualitatively based on visual observations and known properties of unsaturated moisture flow. Results indicate that ERT is a viable method of visualizing 3D unsaturated moisture flow in cement-based materials with discrete cracks. - Highlights: • 3D EIT is developed to visualize water ingress in cracked mortar. • Mortar with different size discrete cracks are used. • The EIT results are corroborated with numerical simulations. • EIT results accurately show the temporal and spatial variation of water content. • EIT is shown to be a viable method to monitor flow in cracks and matrix.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized demoralization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin demoralization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS, and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobean linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Discrete maximum principle theory will be presented for general finite volume approximations on unstructured meshes. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc, will. be addressed as needed.

Barth, Timothy J.↗
DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part I: the Linear Setting

Fully implicit Runge--Kutta (IRK) methods have many desirable properties as time integration schemes in terms of accuracy and stability, but high-order IRK methods are not commonly used in practice with numerical PDEs due to the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic preconditioning framework for solving the systems of equations that arise from IRK methods applied to linear numerical PDEs (without algebraic constraints). Additionally, this framework also naturally applies to discontinuous Galerkin discretizations in time. Under quite general assumptions on the spatial discretization that yield stable time integration, the preconditioned operator is proven to have condition number bounded by a small, order-one constant, independent of the spatial mesh and time-step size, and with only weak dependence on number of stages/polynomial order; for example, the preconditioned operator for 10th-order Gauss IRK has condition number less than two, independent of the spatial discretization and time step. The new method can be used with arbitrary existing preconditioners for backward Euler-type time-stepping schemes and is amenable to the use of three-term recursion Krylov methods when the underlying spatial discretization is symmetric. The new method is demonstrated to be effective on various high-order finite-difference and finite element discretizations of linear parabolic and hyperbolic problems, demonstrating fast, scalable solution of up to 10th-order accuracy. The new method consistently outperforms existing block preconditioning approaches, and in several cases, the new method can achieve 4th-order accuracy using Gauss integration with roughly half the number of preconditioner applications and wallclock time as required using standard diagonally IRK methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A mass conservative well-balanced reconstruction at wet/dry interfaces for the Godunov-type shallow water model

This paper presents a novel mass conservative, positivity preserving wetting and drying treatment for Godunov-type shallow water models with second-order bed elevation discretization. The novel method allows to compute water depths equal to machine accuracy without any restrictions on the time step or any threshold that defines whether the finite volume cell is considered to be wet or dry. The resulting scheme is second-order accurate in space and keeps the C-property condition at fully flooded area and also at the wet/dry interface. For the time integration, a second-order accurate Runge–Kutta method is used. The method is tested in two well-known computational benchmarks for which an analytical solution can be derived, a C-property benchmark and in an additional example where the experimental results are reproduced. Overall, the presented scheme shows very good agreement with the reference solutions. The method can also be used in the discontinuous Galerkin method.

47 OTHER INSTRUMENTATION↗