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 793 records · Page 44

Electron Beam Welding of Nimonic 80A Superalloy: Microstructure Evolution and EBSD Study After Aging Heat Treatment

The purpose of this research was to evaluate integrity and microstructure changes of the electron beam-welded Nimonic 80A and the base metal after aging heat treatment process, i.e., 700 °C/72 h/air cooling. Here, microstructure, microhardness and tensile were investigated by various tools, i.e., scanning electron microscope with energy-dispersive spectroscopy and electron backscatter diffraction, x-ray diffraction analysis and Vickers hardness. Evaluation of microstructure shows that the weld zone mainly consists of columnar dendrites. The effect of aging heat treatment on the microstructure morphology and phase structure is very limited. Also, it was observed that there are obvious segregation and precipitates in intracrystalline and along the grain boundaries after aging heat treatment. Electron backscatter diffraction results show that the weld zone and base metal have a clear difference in the grain boundary characterization and kernel average misorientation at room temperature and after aging heat treatment, which leads to a difference in microhardness of the weld zone. In addition, the relationship between microhardness and tensile strength is discussed.

36 MATERIALS SCIENCE↗

Mass transport analysis for tritium removal in FHRs

Tritium control and mitigation is one of the most significant issues in Fluoride Salt-cooled Hightemperature Reactors (FHRs). To address this issue, a cross-flow tritium removal facility has been proposed to remove molecular tritium, T 2 , from the primary coolant. To model the performance of the design, a method is developed based on the logarithmic mean difference of the square root of the tritium partial pressures on the upstream side and the purging gas side. A computer code has been developed based on this global analysis method to study the effectiveness of tritium removal in such a cross-flow tritium removal facility. The code is benchmarked against experimental data from a hydrogen separation experiment in the literature and compared with the results from a code based on traditional finite volume method. Evaluation of the cross-flow tritium removal facility designed for FHR systems has also been performed. The results show how the key factors, i.e., the fluid flow rates and tritium inlet concentration affect the performance of the facility.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Half-closed discontinuous Galerkin discretisations

Here we introduce the concept of half-closed nodes for nodal discontinuous Galerkin (DG) discretisations. Unlike more commonly used closed nodes in DG, where on every element nodes are placed on all of its boundaries, half-closed nodes only require nodes to be placed on a subset of the element's boundaries. The effect of using different nodes on DG operator sparsity is studied and we find in particular for there to be no difference in the sparsity pattern of the Laplace operator whether closed or half-closed nodes are used. On quadrilateral/hexahedral elements we use the Gauss-Radau points as the half-closed nodes of choice, which we demonstrate is able to speed up DG operator assembly in addition to leverage previously known superconvergence results. We also discuss in this work some linear solver techniques commonly used for Finite Element or discontinuous Galerkin methods such as static condensation and block-based methods, and how they can be applied to half-closed DG discretisations.

97 MATHEMATICS AND COMPUTING↗

Impact of bulk-edge coupling on observation of anyonic braiding statistics in quantum Hall interferometers

Abstract Quantum Hall interferometers have been used to probe fractional charge and statistics of quasiparticles. We present measurements of a small Fabry–Perot interferometer in which the electrostatic coupling constants which affect interferometer behavior can be determined experimentally. Near the center of the ν = 1/3 state this device exhibits Aharonov–Bohm interference interrupted by a few discrete phase jumps, and Φ 0 oscillations at higher and lower magnetic fields, consistent with theoretical predictions for detection of anyonic statistics. We estimate the electrostatic parameters K I and K I L by two methods: using the ratio of oscillation periods in compressible versus incompressible regions, and from finite-bias conductance measurements. We find that the extracted K I and K I L can account for the deviation of the phase jumps from the theoretical anyonic phase θ a = 2 π /3. At integer states, we find that K I and K I L can account for the Aharonov–Bohm and Coulomb-dominated behavior of different edge states.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Microstructure-Sensitive Uncertainty Quantification for Crystal Plasticity Finite Element Constitutive Models Using Stochastic Collocation Methods

Uncertainty quantification (UQ) plays a major role in verification and validation for computational engineering models and simulations, and establishes trust in the predictive capability of computational models. In the materials science and engineering context, where the process-structure-property-performance linkage is well known to be the only road mapping from manufacturing to engineering performance, numerous integrated computational materials engineering (ICME) models have been developed across a wide spectrum of length-scales and time-scales to relieve the burden of resource-intensive experiments. Within the structure-property linkage, crystal plasticity finite element method (CPFEM) models have been widely used since they are one of a few ICME toolboxes that allows numerical predictions, providing the bridge from microstructure to materials properties and performances. Several constitutive models have been proposed in the last few decades to capture the mechanics and plasticity behavior of materials. While some UQ studies have been performed, the robustness and uncertainty of these constitutive models have not been rigorously established. In this work, we apply a stochastic collocation (SC) method, which is mathematically rigorous and has been widely used in the field of UQ, to quantify the uncertainty of three most commonly used constitutive models in CPFEM, namely phenomenological models (with and without twinning), and dislocation-density-based constitutive models, for three different types of crystal structures, namely face-centered cubic (fcc) copper (Cu), body-centered cubic (bcc) tungsten (W), and hexagonal close packing (hcp) magnesium (Mg). Our numerical results not only quantify the uncertainty of these constitutive models in stress-strain curve, but also analyze the global sensitivity of the underlying constitutive parameters with respect to the initial yield behavior, which may be helpful for robust constitutive model calibration works in the future.

36 MATERIALS SCIENCE↗

A System-Level Model for Estimating Residual Strain and Life of Nuclear Reactor Coolant System Components Under Connected-System-Thermal–Mechanical Boundary Conditions

Background: Environmental-assisted fatigue (EAF) is a major issue for the long-term survival of nuclear power plant fleets in the U.S. and worldwide. Multi-material welded regions (e.g., nozzles) and other high-stress regions of reactor coolant system (RCS) components are prone to EAF-related damage. Objective: The discussed work describes a system-level finite element (FE) model of RCS components of a pressurized water reactor (PWR). This is with the goal of predicting the stress hotspots, strain residuals, strain amplitudes and the resulting fatigue lives. Methods: Herein, the FE model was developed considering system-level loading conditions (under connected system thermal-mechanical boundary conditions). Thermal-mechanical stress analysis was performed considering thermal stratification and a design-basis reactor loading cycle. Based on the FE model results, the strain residuals, strain amplitudes and resulting fatigue lives of RCS components were predicted. Results: The results show that some of the RCS components can have significantly different strain amplitudes, residual strain, and fatigue lives, despite having similar geometry and material. Higher residual strain can lead to accelerated cyclic hardening of material and the associated effect of EAF. The simulated component-level strain profile (under realistic multi-axial-multi-physics loading cycle) can guide the selection of appropriate test inputs for conducting laboratory-scale EAF tests, which is a focus of future works. Conclusions: Despite similar geometry and material the RCS component can have significantly different strain profiles and resulting fatigue lives. The discussed approach can help to identify and prioritize the RCS components for conducting expensive nondestructive evaluation (NDE) inspections.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Bridging the time scale in exascale computing of chemical systems (Final Technical Report)

This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0019441. The theme of this project was to develop and apply methods that allowed for the acceleration of atomistic calculations, particularly in challenging areas such as multiphase systems, electrified interfaces, uncertainty estimation, and applications requiring chemical accuracy, which tend to be applications where simulation time is severely bottlenecked by the computational time requirements. Much of the focus was on the application of emerging machine-learning methodologies, although a wide range of methodologies were employed. This report has two major sections. The first focuses on the methodological advances themselves. Within this part, we report a number of major advances, a few examples of which are described here. We report the first machine-learning scheme for the acceleration of electronically grand-canonical calculations (that is, those applicable to electrochemistry). We report new methods of performing transfer learning, in which physics-based priors can be used to provide predictions, often with uncertainty estimates, of images well outside of training sets; we also offer ways to fine-tune these transfer-learning models. We provide a new systematic means to generate and apply minimal training data sets to very large (10,000’s of atoms) systems, with only small training sets appropriate for electronic structure. We developed new methodologies to integrate surface vibrations into surface adsorption calculations. We made advances to the applicability of diffusion Monte Carlo methods to allow (learned) force prediction, finite-size error correction, and force-free means of searching for transition states. We integrated machine-learned atomistic predictions into mechanism generation codes. Additionally, we released new software including AmpTorch, a modernized version of our original atomistic machine-learning code Amp. The second part of this report focuses on the scientific applications that accompanied, and were often enabled by, the methodological advances described earlier. A few examples follow, but full details are in the individual chapters of the report. For example, we developed a general theory of phonon-induced friction on molecular adsorbates. We showed fundamentally how solvent influences the adsorption and desorption process and how it differs from the processes typically involved at the solid–gas interface, making aqueous-phase and electrocatalysis different from traditional thermocatalysis. We examined how metal–insulator and magnetic transitions can be probed, and accelerated exciton dynamics via Frenkel Hamiltonian parameters. We showed that the nearsighted force-training approach, developed within this project, can predict both the stability and reactivity of large nanoparticles, and can also lead to insights on catalyst coverage on binding energies and entropies. These applied studies, which generally integrated with our method development, allowed us to push forward the theoretical understanding of several reaction classes.

08 HYDROGEN↗

A concurrent irradiation-mechanics multiscale coupling model

The vast majority of our current knowledge regarding the basic mechanisms controlling irradiation effect on mechanical properties is based almost entirely on the results of post-irradiation experiments or theoretical models. However, the concurrent effects of irradiation, mechanical stress, and thermal damage on the failure phenomena of materials and components remain largely unexplored due to its internal multiscale-multiphysics coupling nature. Here, we present here a concurrent irradiation-mechanics multiscale coupling model. The concurrent evolutions of nanoscale irradiation defect clusters, microscale dislocation configurations, and mechanical responses are well captured based on coupling cluster dynamics, discrete dislocation dynamics, and the finite element methods using an effective time marching scheme. Model predictions of defect densities and size are in general agreement with experimental observations. Irradiation hardening is shown to take place also in samples undergoing concurrent irradiation-mechanical loading, similar to samples tested post-irradiation. However, the occurrence of plastic flow localization and dislocation channel formation is not accompanied with apparent yield drop (softening) under concurrent irradiation-mechanical loading conditions, which is different from the post-irradiation case.

42 ENGINEERING↗

Quantitative correlation between the void morphology of niobium-tin wires and their irreversible critical current degradation upon mechanical loading

Understanding the critical current performance variation of Nb 3 Sn superconducting wires under mechanical loading is a crucial issue for the design of next generation accelerator and fusion magnets. In these applications, the mechanical properties of the conductors may become a limiting factor due to the strong electro-magnetic forces resulting from the combination of large magnets and intense magnetic fields. In particular, the presence of voids in the superconducting filament structure, which are formed during the fabrication and the reaction heat treatment, determines localized stress concentrations and possibly the formation of cracks. In this work, we demonstrate a quantitative correlation between the void morphology and the electro-mechanical limits measured on different Bronze route Nb 3 Sn wires. Hot Isostatic Pressing (HIP) prior to the reaction heat treatment is utilized to partially eliminate the voids. The wires’ void distributions - with and without HIP treatment - are detected and statistically analyzed using high energy X-ray micro tomography. The stress concentration due to the shape and distribution of the voids as well as their impact on the electro-mechanical properties are determined through finite element method modeling. Lastly, the results are quantitatively correlated with the experimentally determined limits of the irreversible critical current degradation upon mechanical loading.

36 MATERIALS SCIENCE↗

Analytical Modeling of a Double-Sided Flux Concentrating E-Core Transverse Flux Machine with Pole Windings

In this paper, a nonlinear analytical model based on the Magnetic Equivalent Circuit (MEC) method is developed for a double-sided E-Core Transverse Flux Machine (TFM). The proposed TFM has a cylindrical rotor, sandwiched between E-core stators on both sides. Ferrite magnets are used in the rotor with flux concentrating design to attain high airgap flux density, better magnet utilization, and higher torque density. The MEC model was developed using a series-parallel combination of flux tubes to estimate the reluctance network for different parts of the machine including air gaps, permanent magnets, and the stator and rotor ferromagnetic materials, in a two-dimensional (2-D) frame. An iterative Gauss-Siedel method is integrated with the MEC model to capture the effects of magnetic saturation. A single phase, 1 kW, 400 rpm E-Core TFM is analytically modeled and its results for flux linkage, no-load EMF, and generated torque, are verified with Finite Element Analysis (FEA). The analytical model significantly reduces the computation time while estimating results with less than 10 percent error.

30 DIRECT ENERGY CONVERSION↗

A Performance-Portable MultiGPU Implementation of 3D Euler Equations using ProtoX and IRIS

Computational scientists often face challenges when developing and optimizing code for high-performance computing (HPC), especially when trying to leverage GPUs. Given the heterogeneity of the nodes that comprise many modern HPC facilities, considerable demand exists for performance portable solutions for the core computational kernels used in many scientific computing libraries. In this work, we demonstrate a fourth-order finite volume method–based implementation of the Euler equations, which are an integral part of computational fluid dynamics. Our performance-portable multiGPU implementation for Euler equations uses ProtoX to generate kernels and IRIS for portability. ProtoX is a domain-specific language that uses a structured-grid partial differential equation library called Proto as its front end and the SPIRAL code generation system as its back end to generate optimized kernels for different architectures. Optimized kernels generated by ProtoX are orchestrated through the IRIS intelligent runtime system to provide portability. Two levels of optimizations within the IRIS runtime— directed acyclic graph fusion and task fusion—are explored to efficiently utilize computing resources in a multiGPU environment. Performance improvement through these optimizations is showcased by comparing the base ProtoX-IRIS implementation on AMD GPUs (Frontier node) and on NVIDIA GPUs (NVIDIA DGX-1).

Mankad, Het↗

Semiclassical shell-structure micro-macroscopic approach for the level density

Level density ρ(E,A) is derived for a one-component nucleon system with a given energy E and particle number A within the mean-field semiclassical periodic-orbit theory beyond the saddle-point method of the Fermi gas model. We obtain ρ∝I ν (S)/S ν , with I ν (S) being the modified Bessel function of the entropy S. Within the micro-macro-canonical approximation (MMA), for a small thermal excitation energy U, with respect to rotational excitations E rot , one obtains ν = 3/2 for ρ(E,A). In the case of excitation energy U larger than E rot but smaller than the neutron separation energy, one finds a larger value of ν = 5/2. A role of the fixed spin variables for rotating nuclei is discussed. The MMA level density ρ reaches the well-known grand-canonical ensemble limit (Fermi gas asymptote) for large S related to large excitation energies, and also reaches the finite micro-canonical limit for small combinatorial entropy S at low excitation energies (the constant “temperature” model). Fitting the ρ(E,A) of the MMA to the experimental data for low excitation energies, taking into account shell and, qualitatively, pairing effects, one obtains for the inverse level density parameter K a value which differs essentially from that parameter derived from data on neutron resonances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hygro-thermo-mechanical modeling of partially saturated air-entrained concrete containing dissolved salt and exposed to freeze-thaw cycles

Highlights: • A model for cyclic freezing of air-entrained concrete containing salt is proposed. • The freeze-thaw hysteresis, also including the effect of dissolved salt, is modeled. • The diffusive and convective flux of salt is coupled to the freeze-thaw processes. • The freeze-thaw response of air-entrained concrete containing salt is reproduced. • Liquid absorption from reservoirs increases with an increased concentration of salt. In cold regions, understanding the freeze-thaw behavior of air-entrained concrete is important for designing durable structures and assessing the remaining service life of existing structures. This study presents a hygro-thermo-mechanical multiphase model that describes the cyclic freeze-thaw behavior of partially saturated air-entrained concrete containing dissolved salt. An equilibrium and a non-equilibrium approach are adopted to model the ice formation, including the freeze-thaw hysteresis, inside the porous network. The model also considers the diffusive and convective transport of the dissolved salt coupled to the freeze-thaw processes. Two examples are presented to verify and highlight the capabilities of the model. The first example shows that the model is capable of reproducing the experimentally observed mechanical response of specimens containing NaC1-solutions of different concentrations. In the second example, a larger absorption of liquid from an external reservoir is obtained with an increasing salt concentration in the reservoir, which is consistent with experimental observations.

36 MATERIALS SCIENCE↗

Measuring the Error in Approximating the Sub-Level Set Topology of Sampled Scalar Data

This paper studies here the influence of the definition of neighborhoods and methods used for creating point connectivity on topological analysis of scalar functions. It is assumed that a scalar function is known only at a finite set of points with associated function values. In order to utilize topological approaches to analyze the scalar-valued point set, it is necessary to choose point neighborhoods and, usually, point connectivity to meaningfully determine critical-point behavior for the point set. Two distances are used to measure the difference in topology when different point neighborhoods and means to define connectivity are used: (i) the bottleneck distance for persistence diagrams and (ii) the distance between merge trees. Usually, these distances define how different scalar functions are with respect to their topology. These measures, when properly adapted to point sets coupled with a definition of neighborhood and connectivity, make it possible to understand how topological characteristics depend on connectivity. Noise is another aspect considered. Five types of neighborhoods and connectivity are discussed: (i) the Delaunay triangulation; (ii) the relative neighborhood graph; (iii) the Gabriel graph; (iv) the k-nearest-neighbor (KNN) neighborhood; and (v) the Vietoris–Rips complex. It is discussed in detail how topological characterizations depend on the chosen connectivity.

97 MATHEMATICS AND COMPUTING↗

A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries

Herein, we present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an hp adaptive mesh refinement (hp-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The hp-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The hp-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.

97 MATHEMATICS AND COMPUTING↗

A discrete-continuum coupled finite element modelling approach for fibre reinforced concrete

Highlights: • An easy-to-implement but effective meso-scale modelling approach is developed for FRC. • It is capable of simulating and interpreting all main failure mechanisms in FRC. • It is critically validated by experiments with single or multiple random fibres under various loading conditions. • It can be used to optimise key design factors of FRC at meso/material scale by parametric studies. - Abstract: Fibre reinforced concrete (FRC) exhibits complicated failure modes such as fibre breakage, mortar cracking, crushing and spalling, fibre-mortar interfacial debonding, depending on many material properties, geometric dimensions, boundary and loading conditions. Most existing numerical models are unable to reproduce these failure modes that may occur simultaneously or sequentially in a specimen, mainly due to difficulties in generating finite element meshes with a large number of randomly-oriented fibres. Herein we develop a discrete-continuum coupled finite element modelling approach for FRC materials capable of effectively simulating all the major failure modes. The continuum damaged plasticity model is used to simulate damage and fracture behaviour of the mortar, while debonding of fibre-mortar interfaces is modelled by nonlinear cohesive interfacial elements. Unique techniques are devised to generate conforming meshes between fibres and the surrounding mortar so that the randomly-oriented fibres are easily modelled. The modelling approach is validated by simulating single fibre pullout tests with different inclination angles, notched and non-notched direct tensile tests and three-point bending beam tests with randomly-distributed multiple fibres.

36 MATERIALS SCIENCE↗

Fatigue life prediction of powder bed fused–laser beam AlSi10Mg: Incorporating critical defects via crystal plasticity modelling

The current study provides a microstructurally-based computational framework to predict the fatigue life of additive manufactured (AM), i.e., powder bed fused–laser beam (PBF-LB), AlSi10Mg specimens using the crystal plasticity finite element method (CPFEM). The fractography analysis, electron backscatter diffraction (EBSD), uniaxial and cyclic responses, and fatigue life of specimens were used to inform the computational framework. CPFE simulation was used to compute fatigue indicator parameters (FIPs) as fatigue driving forces. A new fatigue criterion is introduced based on FIPs, which was calibrated using experimental fatigue data. The proposed fatigue measure was evaluated versus the specimens with critical defects of various sizes and locations subjected to different stress amplitudes. The results show that the developed framework can capture the fatigue life of samples with different critical defect locations and sizes along with different stress amplitudes for both high-cycle fatigue (HCF) and very high-cycle fatigue (VHCF) regimes.

Additive manufacturing↗

Massively parallel axisymmetric fluid model for streamer discharges

A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗