Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite difference method”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Mechanics of the Ring Tension Test (RTT): A Finite Element-based Investigation

The ring tension test (RTT) is an experimental method for determining mechanical behavior in a material’s circumferential or hoop direction. It is a crucial test for testing anisotropic materials with tube geometry, such as nuclear fuel cladding or irradiated pipes. Several RTT configurations exist, each with their own advantages and disadvantages. However, this test is significantly more complex than traditional tensile testing and can be especially sensitive to small differences and inconsistencies in the test setup and geometry, ultimately affecting the derived mechanical properties. Previous research has focused on method development, and little work has been done on understanding the subtle differences between an ideal test and experiments, specifically when the tests are performed on highly irradiated materials in hot cells. In this work, a finite element-based investigation of the RTT is conducted. Two promising test configurations are investigated, comparing their ability to determine accurate material strengths through plastic deformation. Several non-ideal conditions and uncontrollable effects which are likely to occur during experimental testing such as machining tolerances, variations of specimen geometry from nominal dimensions, rotation of specimens and fixturing, and other test setup discrepancies are studied. The sensitivity of measured strengths to these conditions is presented. A mechanics-based approach to describing and correcting raw data to determine actual strengths is also included for one of the configurations, resulting in a robust correction method with highly accurate material strength measurements. Based on these analyses, the hemicylindrical mandrel configuration is recommended with a gauge region oriented at a 45° angle.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Computational modeling of continuum scale constitutive equations to improve biomass feedstock material handling and conveying systems

Existing design practices for common feedstock handling equipment such as hoppers and screw feeders have not been successful in many situations. The lack of an appropriate toolset to characterize and model biomass flow behavior is deemed the root cause of plugging, stalls, and flow-related issues in biomass handling systems design and operations. New standardized characterization protocols in conjunction with appropriate constitutive models will lead to successful biomass handling system design and operations. This study investigates constitutive models for bulk solid flow based on the first principles to overcome issues in biomass handling. The focus is on two major issues in biomass handling, i.e., initiating and sustaining of flow in the hopper that is widely used in biomass conveying operation. Accordingly, this article focuses on analytical models of biomass hopper flow. Mohr-Coulomb, Drucker-Prager, and modified Cam-clay models are investigated as constitutive bulk solid flow models. Parameters of those models are determined using a set of triaxial tests. Triaxial tests are conducted using the cubical triaxial tester (CTT) in the Food and Biomaterials Lab of the Pennsylvania State University. Milled 1mm corn stover and 2mm wood chip biomass materials having different flow characteristics are examined, and their respective flow behavior out of a hopper are simulated using the Finite Element method implementing identified constitutive models and determined fundamental mechanical properties. It is found that, with the same hopper geometry, ground corn stover 1mm tends to exhibit funnel flow whereas ground Douglas fir 2 mm tends to exhibit mass flow. These simulation results are validated with hopper flow experiments. This study demonstrates how biomass flow problems can be predicted and mitigated by numerically modeled fundamental constitutive equations.

09 BIOMASS FUELS↗
DOE OSTI2019

Thermal behavior characterization of a kilowatt-power-level cryogenically cooled Yb:YAG active mirror laser amplifier

We present an analysis of the thermal behavior of a high-energy kilowatt-average-power diode-pumped cryogenically cooled Yb:YAG active mirror laser amplifier based on measurements and simulations. Maps of the temperature distribution of the laser material at pump powers up to 1 kW were obtained for the first time by spatially and spectrally resolving the fluorescence induced by a scanning probe beam. The wavefront distortion resulting from the front surface deformation and the overall deformation of the gain medium assembly were measured using a Mach–Zehnder interferometer. The measured deformations agree well with the results of thermomechanical modeling using finite element method simulations, and with the results of focal length shift measurements. The relative contributions to the optical path difference (OPD) of the mechanical deformations, refractive index changes, and electronic contribution are discussed. We propose that the Cr 4+ :YAG cladding plays a significant role in both the temperature distribution and the overall OPD changes. Moreover, the pump-induced mechanical deformations of the assembly dominate the OPD changes in this kilowatt-average-pump-power cryogenically cooled Yb:YAG active mirror laser.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Novel methods for the time-dependent Maxwell's equations and their applications

This dissertation investigates three different mathematical models based on the time-domain Maxwell's equations using three different numerical methods: a Yee scheme using a non-uniform grid, a nodal discontinuous Galerkin (nDG) method, and a newly developed discontinuous Galerkin method named the weak Galerkin (WG) method. The non-uniform Yee scheme is first applied to an electromagnetic metamaterial model. Stability and superconvergence error results are proved for the method, which are then confirmed through numerical results. Additionally, a numerical simulation of backwards wave propagation through a negative-index metamaterial is given using the presented method. Next, the nDG method is used to simulate signal propagation through a corrugated coaxial cable through the use of axisymmetric Maxwell's equations. Stability and error analysis are performed for the semi-discrete method, and are verified through numerical results. The nDG method is then used to simulate signal propagation through coaxial cables with a number of different corrugations. Finally, the WG method is developed for the standard time-domain Maxwell's equations. Similar to the other methods, stability and error analysis are performed on the method and are verified through a number of numerical experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Optimal Polynomial Smoothers and One‐Sided V‐Cycles for Poisson Problems

The solution to the Poisson equation arising from the spectral element discretization of the incompressible Navier‐Stokes equations needs robust preconditioning strategies. One such strategy is multigrid. To realize the potential of multigrid methods, effective smoothing strategies are needed. Chebyshev polynomial smoothers, in conjunction with pointwise Jacobi or additive Schwarz methods (ASMs), prove to be an effective smoother. Other polynomial smoothers, however, may provide superior convergence to the multigrid preconditioner. The authors compare the standard Chebyshev polynomial smoothers to both the novel fourth‐kind Chebyshev polynomial smoothers proposed by Lottes as well as smoothers based on the polynomial of best uniform approximation to as proposed by Kraus, Vassilevski, and Zikatanov. At the cost of symmetry, further improvements may be made. For example, a order polynomial smoother on both sides of the V‐cycle may be substituted with an order polynomial smoother on one side at no additional cost. The choice of omitting the postsmoother in favor of higher‐order polynomial presmoothing is advantageous in cases where the multigrid approximation property constant is large. The authors consider a 2D model problem based on finite differences to motivate the choice of polynomial smoother, order, and whether to apply postsmoothing for the target application of high‐order ‐geometric multigrid methods for GPU architectures. Results from both domains demonstrate the substantial improvement of these approaches over the standard Chebyshev polynomial smoother with a symmetric V‐cycle.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Customizable two-species kinetic equilibria for nonuniform low-beta plasmas

Here, two methods are developed for constructing self-consistent two-species kinetic equilibria for low-beta electrostatic plasmas, in which the magnetic field is uniform and fixed in time. The first method uses Taylor series approximations to construct distribution functions that can be specified analytically. The second method relies on numerically solving a nonlinear ordinary differential equation and produces exact—to numerical precision—equilibria. In both methods, the equilibrium distribution functions for ions and electrons are expressed in terms of constants of motion and satisfy the steady-state Vlasov-Poisson equation system. Provided that the ion drift speed does not exceed the ion thermal speed, the equilibria can be specified with customizable density and electrostatic potential profiles. The methods can thereby be tailored to different applications and are successfully applied to construct kinetic equilibria for cross-field plasmas with sheared flows, large density variations, and different levels of magnetization. The equilibria are used to initialize fourth-order finite-volume Vlasov-Poisson simulations in (x, v x , v y ) coordinates and the associated temporal evolution is used to assess the accuracy of each method. Lastly, the low-amplitude deviations observed in these simulations demonstrate that the kinetic equilibria are robust and that they provide a valuable means of studying the dynamics of nonuniform magnetized plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Heat and Mass Transfer in Cesium Evaporation from Graphite Surface in Argon

The aim of the present work is mathematical modeling of the heat and mass transfer of {sup 137}Cs under the action of low-temperature equilibrium plasma flow on sections of graphite contaminated by radionuclides. The methods of investigation are calculation of the temperature and mass rate of evaporation of radiocesium by solving a system of quasi-homogeneous nonstationary heat conduction equations with boundary conditions of the first and third kind, which take account of mass transfer, by means of finite-difference elements using the ANSYS package.

42 ENGINEERING↗
DOE OSTI2024

A Characteristics Approach to the Finite Element Method

Herein, we present a new method for solving the linear Boltzmann transport equation. Two commonly used and well-understood methods for solving partial differential equations are the method of characteristics (MOC) and the finite element method (FEM). We propose a new method that combines the fundamental concept of the FEM with the analytic solution from the MOC to obtain coefficients for the FEM basis function expansion. Traditionally, coefficients for the FEM basis function expansion are obtained via matrix inversion. Instead, we solve for the coefficients with the MOC and represent the underlying fields with the basis function expansion using these coefficients. We provide a convergence study for our method with results from two sets of FEM basis functions: Gauss-Legendre and Gauss-Lobatto sets. We also compare two different variations of our method categorized as short characteristics and intermediate characteristics.

42 ENGINEERING↗
DOE OSTI2018

FELIX-2.0: New version of the finite element solver for the time dependent generator coordinate method with the Gaussian overlap approximation

The time-dependent generator coordinate method (TDGCM) is a powerful method to study the largeamplitude collective motion of quantum many-body systems such as atomic nuclei. Under the GaussianOverlap Approximation (GOA), the TDGCM leads to a local, time-dependent Schrödinger equation in amulti-dimensional collective space. In this paper, we present the version 2.0 of the code FELIX that solvesthe collective Schrödinger equation in a finite element basis. This new version features: (i) the ability tosolve a generalized TDGCM+GOA equation with a metric term in the collective Hamiltonian, (ii) supportfor new kinds of finite elements and different types of quadrature to compute the discretized Hamiltonianand overlap matrices, (iii) the possibility to leverage the spectral element scheme, (iv) an explicit Krylovapproximation of the time propagator for time integration instead of the implicit Crank–Nicolson methodimplemented in the first version, (v) an entirely redesigned workflow. We benchmark this release on ananalytic problem as well as on realistic two-dimensional calculations of the low-energy fission of240Puand256Fm. Low to moderate numerical precision calculations are most efficiently performed with simplexelements with a degree 2 polynomial basis. Higher precision calculations should instead use the spectralelement method with a degree 4 polynomial basis. We emphasize that in a realistic calculation of fissionmass distributions of240Pu, FELIX-2.0 is about 20 times faster than its previous release (within a numericalprecision of a few percents).

Regnier, D.↗
DOE OSTI2018

Optimal Control of a Degenerate PDE for Surface Shape

Controlling the shapes of surfaces provides a novel way to direct self-assembly of colloidal particles on those surfaces and may be useful for material design. This motivates the investigation of an optimal control problem for surface shape in this paper. Specifically, we consider an objective (tracking) functional for surface shape with the prescribed mean curvature equation in graph form as a state constraint. The control variable is the prescribed curvature. We prove existence of an optimal control, and using improved regularity estimates, we show sufficient differentiability to make sense of the first order optimality conditions. This allows us to rigorously compute the gradient of the objective functional for both the continuous and discrete (finite element) formulations of the problem. Numerical results are shown to illustrate the minimizers and optimal controls on different domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

DPC Disposal Thermal Scoping Analysis

This is a progress report on thermal modeling for dual-purpose canister (DPCs) direct disposal that covers several available calculation methods and addresses creep and temperature-dependent properties in a salt repository. Three modeling approaches are demonstrated: A semi-analytical calculation method that uses linear solutions with superposition and imaging, to represent a central waste package in a larger array; A finite difference model of coupled thermal creep, implemented in FLAC2D; and An integrated finite difference thermal-hydrologic modeling approach for repositories in different generic host media, implemented in PFLOTRAN. These approaches are at different levels of maturity, and future work is expected to add refinements and establish the best applications for each.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

Matrix Methods for Solving Hartree-Fock Equations in Atomic Structure Calculations and Line Broadening

Atomic structure of N-electron atoms is often determined by solving the Hartree-Fock equations, which are a set of integro-differential equations. The integral part of the Hartree-Fock equations treats electron exchange, but the Hartree-Fock equations are not often treated as an integro-differential equation. The exchange term is often approximated as an inhomogeneous or an effective potential so that the Hartree-Fock equations become a set of ordinary differential equations (which can be solved using the usual shooting methods). Because the Hartree-Fock equations are an iterative-refinement method, the inhomogeneous term relies on the previous guess of the wavefunction. In addition, there are numerical complications associated with solving inhomogeneous differential equations. This work uses matrix methods to solve the Hartree-Fock equations as an integro-differential equation. It is well known that a derivative operator can be expressed as a matrix made of finite-difference coefficients; energy eigenvalues and eigenvectors can be obtained by using linear-algebra packages. The integral (exchange) part of the Hartree-Fock equation can be approximated as a sum and written as a matrix. The Hartree-Fock equations can be solved as a matrix that is the sum of the differential and integral matrices. We compare calculations using this method against experiment and standard atomic structure calculations. This matrix method can also be used to solve for free-electron wavefunctions, thus improving how the atoms and free electrons interact. This technique is important for spectral line broadening in two ways: it improves the atomic structure calculations, and it improves the motion of the plasma electrons that collide with the atom.

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

Effective sensor location for detection of change in structural dynamic response

This paper investigates the sensitivity of structural system response to the sensor location by investigating consequences of small changes in the location to the structural system response. The paper discusses how maximum observability (based on mode shape and the participation of that mode in the input provided) drives optimal location. The structural responses were investigated in terms of the g-rms response for various low-frequency inputs (pure sinusoids and real-life inputs such as an earthquake and trains). Results were then analyzed in the context of Modal Contributions Factors (MCF) and changes to the Force-to-response Transfer Functions (TRFs). A modal-matching process is first presented using a Matlab TM -based Finite Element Method (FEM) model of a cantilever beam and instrumentation to determine the location of a small mass based on three different criteria. Subsequently, the structural response is investigated using experiments and the FEM model. The accelerometer of small mass (at 1/3 height) was moved up or down to obtain changes in the structural response (TRF) to various realistic low-frequency inputs. Modal Contribution Factor (MCF) and derivative (slope) of the associated mode-shapes were correlated to the observed changes in TRFs. Results show how optimal sensor locations for detecting change in structural response can be based on the MCFs and the associated mode-shapes.

Woodall, J.↗
DOE OSTI2017

Flexible Spatial Partitions in MPACT Through Module-Based Data Passing

As part of the Consortium for Advanced Simulation of Light Water Reactors (CASL), the Virtual Environment for Reactor Applications (VERA) is being developed to provide high-fidelity multiphysics simulations of nuclear reactor cores. The MPACT code being developed collaboratively by Oak Ridge National Laboratory and the University of Michigan is the primary deterministic neutron transport solver available within VERA. MPACT employs the 2D/1D method to solve 3D problems using 2D method of characteristics (MOC) radial transport solvers with 1D nodal expansion method-based P3 (NEM-P3) axial transport solvers, which are coupled using axial and radial transverse leakages. A 3D coarse mesh finite difference (CMFD) solver is used to provide accelerated convergence and stability. Quarter-core problems in MPACT are typically run on moderately sized clusters ({approx}500-10,000 cores). As such, it relies heavily on spatial decomposition, both radially and axially, using MPI data communication along spatial domain boundaries to pass angular flux boundary condition and neighboring data. The initial spatial decomposition strategy implemented into MPACT overlaid a Cartesian grid over the core to form the domain boundaries. With this approach, a single buffer spanning across the entire boundary between two spatial domains is passed to a domain's neighbor. Additionally, only one neighboring domain index was allowed in each direction. For example, the entire eastern boundary of a domain must be the entire western boundary of the corresponding neighbour's domain. This implementation helped to execute cases on several thousand cores, where substantial radial decomposition was used. The default for several years was to use assembly decomposition, which yielded 73 radial partitions in the Watts Bar Unit 1 (WBN1) models being analyzed. However, as CASL formed stronger connections with industry partners, there was a strong push to allow MPACT to run on fewer cores, preferably in the 500-1,000 core range. Such cases require much less radial decomposition, typically with 8-16 radial domains on each plane. Unfortunately, the previous decomposition strategy is very limited in its ability to yield well-balanced domains, and it does not permit partitions with an arbitrary number of domains, particularly with small numbers of radial domains (3, 5, 7, 11, 13, etc.). With these limitations, a more flexible approach that provided more even-load balancing was needed. The details of this new strategy are outlined here, along with insights to its advantages. The examples provided are for the WBN1 core, which has been the focus of several initial CASL analyses. It is worth noting that spatial decomposition strategies have been employed in many codes with focus on a number of different transport methodologies. The references cited are not exhaustive as this has been very important research topic across many areas.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Residually Stressed Bimaterial Beam Specimen for Measuring Environmentally Assisted Crack Growth

Background: Subcritical crack growth can occur in a brittle material when the stress intensity factor is smaller than the fracture toughness if an oxidizing agent (such as water) is present at the crack tip. Objective: Here we present a novel bi-material beam specimen which can measure environmentally assisted crack growth rates. The specimen is “self-loaded” by residual stress and requires no external loading. Methods: Two materials with different coefficient of thermal expansion are diffusion bonded at high temperature. After cooling to room temperature a subcritical crack is driven by thermal residual stresses. A finite element model is used to design the specimen geometry in terms of material properties in order to achieve the desired crack tip driving force. Results: The specimen is designed so that the crack driving force decreases as the crack extends, thus enabling the measurement of the crack velocity versus driving force relationship with a single test. The method is demonstrated by measuring slow crack growth data in soda lime silicate glass and validated by comparison to previously published data. Conclusions: The self-loaded nature of the specimen makes it ideal for measuring the very low crack velocities needed to predict brittle failure at long lifetimes.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Generalization of NESTLE into a Multi-Energy N-Group Formulation

The NESTLE code is a few-group, diffusion-based core simulation code that employs the Nodal Expansion Method (NEM). NESTLE has the capability to solve both steady-state and transient problems, using either Cartesian or hexagonal geometry, and to apply thermal-hydraulic feedback including single and two-phase flow. A modern version of the NESTLE nodal-diffusion simulator is currently maintained and developed at the University of Tennessee. The original version of the code was limited to only two or four neutron energy groups, therefore, this article describes recent work carried out to expand NESTLE into a general multi-group solver. Generalization of a fixed group finite-difference diffusion theory into arbitrary N groups is a well-known problem [3], however, this paper describes the intricacies associated with the generalization of the Nodal Expansion Method, and the methodical and computational challenges encountered and managed along the way.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A level set approach for the computational study of a yield stress fluid filling a thin mold

Many important engineering and scientific applications such as cement slurries, foams, crude oil, and granular avalanches involve the concept of yield stress. Therefore, modeling yield stress fluids in different flow configurations, including the accurate prediction of the yield surface, is important. In this paper, we present a computational model based on the finite element method to study the flow of yield stress fluids in a thin mold and compare the results with data from flow visualization experiments. Here, we use the level set method to describe the interface between the filling fluid and air. We use polypropylene glycol as a model Newtonian fluid and Carbopol for the model yield stress fluid, as the Carbopol solution demonstrates yielding without thixotropy. To describe the yielding and shear-thinning behavior, we use a generalized Newtonian constitutive equation with a Bingham–Carreau–Yasuda form. We compare the results obtained from the mold filling experiments with the results from the three-dimensional (3D) model and from a reduced-order Hele-Shaw (HS) model that is two-dimensional, including the effect of shear-thinning along the thin direction only approximately. We show that both the 3D and the HS model can capture the experimental meniscus shape reasonably well for all the fluids considered at three different flow rates. This indicates that the shape evolution is insensitive to the dimensionality of the model. However, the viscosity and yield surfaces predicted by the 3D and HS models are different. The HS model underestimates the high viscosity and unyielded regions compared to the estimation by the 3D model.

42 ENGINEERING↗
DOE OSTI2020

CFD-DEM and PR-DNS studies of low-temperature densely packed beds

Over the past few decades, granular media is gaining attention as a viable option for heat transfer fluids (HTFs). Several research efforts are studying the use of particle-based heat transfer fluids in a wide variety of applications. With this motivation, the current work focusses on analyzing the different heat transfer mechanisms in low-temperature mono-sized densely packed granular media. To study the heat transfer behavior of granular media at different scales, the current work employs a two-way coupled computational strategy. The motion of particles is solved using the Discrete Element Method (DEM) and the interstitial air is solved using a Finite-Volume (CFD) approach. The Open-Source library CFDEM Coupling® is used in the current study to join the Finite Volume PISO solver of OpenFOAM® and the DEM solver of LIGGGHTS®. Typically, particle-particle contact conduction and particle-air convection are the most popular closure models. But recent research identified a different heat transfer phenomenon in packed beds that cannot be identified by conduction or convection models. While closure models were developed to implement this on a CFD-DEM framework, they did not capture the effect of intra-particulate thermal gradients on this phenomenon. Hence the current work also employs Particle-Resolved Direct Numerical Simulations (PR-DNS) to gain valuable insights allowing for the modification of existing models. A new closure model is then proposed here and is implemented in the CFD-DEM framework. This model provides key insights into the different heat transfer mechanism of packed beds.

42 ENGINEERING↗