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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Elevated temperature mechanical properties of Inconel 617 surface oxide using nanoindentation

Inconel 617 is a principal candidate material for helium gas cooled very-high-temperature reactors with outlet temperatures of 700–950 °C. Recent findings showed that this alloy develops unique surface oxide layers especially at high temperature (HT) helium environment with distinctive wear, friction and contact properties. This study investigates the elevated temperature mechanical properties of Inconel 617 top surface layers aged in HT helium environment. Nanoindentation technique is used to obtain load-displacement graphs of the alloy top surface oxide in temperatures ranging from 25 °C up to 600 °C. In addition, using finite element analysis along with an iterative regression technique, a semi-numerical method is developed to further measure and quantify the material parameters and, in particular, time-independent and creep characteristics of the oxide. Although Young's modulus of the oxide is found to be relatively close to the bulk for the tested temperatures, the yield strength and hardness, in comparison to the bulk material, increase significantly as the material is oxidized after aging. The oxide exhibits significant softening as the temperature increases to 600 °C. Unlike the bulk material, diffusion through the grains is found to be the dominant creep mechanism for the oxide. Considerable difference between the mechanical properties of the oxide and the bulk material shows the need for accurate measurements of near surface mechanical properties, if reliable predictive contact and tribological models are sought at elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Anderson acceleration stability in NDA-accelerated k-eigenvalue problems

Anderson acceleration (AA) has been used to improve the stability and convergence rate of multiphysics iterative methods for reactor analysis. Most applications studied assume a tightly converged solution for the different physics problems, and AA is usually applied to state variables like temperature, density, and heat generation rate. In this paper, we study the theoretical performance of AA in NDA-accelerated k-eigenvalue problems. The problems and algorithms studied are simplified from the coupled iteration scheme adopted by MPACT and many other high-fidelity whole-core reactor codes. Compared to previous analyses of AA for these iteration schemes, we study the case with a partially converged neutronics solution and possibly partially converged nonlinear diffusion acceleration (NDA)/coarse mesh finite difference (CMFD) solutions. We observe that the performance of the iteration scheme with AA is very sensitive to the initial guess and is affected by the partially converged CMFD solutions. When the NDA solution is fully converged, using AA cannot achieve the optimal convergence rate in large-sized problems. Conversely, if the NDA solution is partially converged, the iteration scheme with AA can diverge or converge extremely slowly. It is found that the loss of robustness for AA is due to the fact that it is applied to the iterative subspace of state variables rather than the fundamental unknowns of the governing equations. To improve the robustness, the scalar flux should also be considered in the implementation of AA. After considering the residuals of flux, we observe that the stability is regardless of the partial convergence of NDA solutions. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Validation of Hole-Drilling Residual Stress Measurements in Workpieces of Various Thickness

A recent revision to the ASTM E837 standard for near-surface residual stress measurement by the hole-drilling method describes a new thickness-dependent stress calculation procedure applicable to “thin” and “intermediate” workpieces for which strain versus depth response depends on workpiece thickness. This new calculation procedure differs from that of the prior standard, which applies only to thick workpieces with strain versus depth response independent of thickness. Herein we assess the new calculation procedures by performing hole-drilling residual stress measurements in samples with a range of thickness. Near-surface residual stress is measured in a thick aluminum plate containing near-surface residual stress from a uniform shot peening treatment, and in samples of different thickness removed from the plate at the peened surface. A finite element (FE) model is used to assess consistency between measured residual stress across the range of sample thickness. Measured residual stress varies with sample thickness, with thinner samples exhibiting smaller near-surface compressive stress and a larger gradient of subsurface stress. These trends are consistent with both observed bending (curvature) of the removed samples and the trend in FE-calculated expected residual stress. The measured and expected residual stresses are in good agreement for samples of intermediate thickness, but the agreement decreases with sample thickness. Measured residual stress is invariant with gage circle diameter. The new thickness-dependent stress calculation procedure for hole-drilling provides meaningful improvement compared to thick-workpiece calculations.

42 ENGINEERING↗
DOE OSTI2025

Identification of $^{3}$He–$^{3}$H clusters in the $^{6}$Li+$^{89}$Y experiment using particle-$\gamma$ coincidence measurement

The 6 Li+ 89 Y experiment was performed to explore the reaction mechanism induced by a weakly bound nucleus 6 Li and its cluster configuration. Here, the particle-$\gamma$ coincidence method was used to identify the different reaction channels. The $\gamma$-rays coincident with 3 He/ 3 H indicate that the 3 H/ 3 He stripping reaction plays a significant role in the formation of Zr/Nb isotopes. The obtained results support the existence of a 3 He- 3 H cluster in 6 Li. Direct and sequential transfer reactions are adequately discussed, and the FRESCO code is used to perform precise finite-range cyclic redundancy check calculations. In the microscopic calculation, direct cluster transfer is more predominant than sequential transfer in 3 H transfer. However, the direct cluster transfer is of comparable magnitude to the sequential transfer in the 3 He transfer.

CRC calculations↗
DOE OSTI2021

Fidelity and entanglement entropy for infinite-order phase transitions

Here, we study the fidelity and the entanglement entropy for the ground states of quantum systems that have infinite-order quantum phase transitions. In particular, we consider the quantum O(2) model with a spin-S truncation, where there is an infinite-order Gaussian (IOG) transition for S = 1 and there are Berezinskii-Kosterlitz-Thouless (BKT) transitions for S ≥ 2. We show that the height of the peak in the fidelity susceptibility (χ F ) converges to a finite thermodynamic value as a power law of 1 / L for the IOG transition and as 1 / ln (L) for BKT transitions. The peak position of χ F resides inside the gapped phase for both the IOG transition and BKT transitions. On the other hand, the derivative of the block entanglement entropy with respect to the coupling constant $S^{'}_{vN}$ has a peak height that diverges as ln 2 (L) for S = 1 and ln 3 (L) for S ≥ 2 and can be used to locate both kinds of transitions accurately. We include higher-order corrections for finite-size scalings and obtain the value of the central charge consistent with c = 1 predicted by conformal field theory. The crossing point of χ F between different system sizes is at the IOG point for S = 1 but is inside the gapped phase for S ≥ 2, while those of $S^{'}_{vN}$ are at the phase-transition points for all S truncations. Our work elaborates on how to use the finite-size scaling of χ F or $S^{'}_{vN}$ to detect infinite-order quantum phase transitions and discusses the efficiency and accuracy of the two methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Bipolar Membranes with an Electrospun 3D Junction

Freestanding bipolar membranes (BPMs) with an entirely new and transformative morphology were fabricated and characterized. The key disruptive design element was a water-splitting/water-generation junction layer of finite thickness composed of interwoven and interlocking electrospun anion-exchange polymer (AEP) fibers and cation-exchange polymer (CEP) fibers. Methods were developed to sandwich the 3D junction layer between two pre-formed dense AEP and CEP films., with the closure of all interfiber voids, where the processing steps are easily incorporated into a commercial roll-to-roll membrane manufacturing scheme. A series of membranes were made using different polymers and junction layers to identify the optimum morphology and composition for water-splitting and water-generation applications. For water splitting, the 3D junction BPMs worked remarkably well with operating current densities at/above 1.0 A/cm2 (10X greater than commercial BPMs) at a transmembrane voltage drop of only 1.1 V. In water generation mode, the 3D junction BPMs operated at 0.5 A/cm2, a world-record current density. The membranes were found to exhibit outstanding durability and can be manufactured at scale for low cost. The high operating current densities and stable long-term operation are due to the 3D junction layer design, where there is a high interfacial area for water splitting or water generation reactions and where the interlocking fibers prevent delamination of the outer films. These new BPMs are ideal candidate materials for both existing water-splitting electrodialysis separations and new electrochemical processes, such as large-scale direct air capture of CO 2 , reactors for CO 2 reduction, self-hydrating fuel cells, and redox flow batteries.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2019

Finite Element Modeling of Orthogonal Machining of Brittle Materials Using an Embedded Cohesive Element Mesh

Machining of brittle materials is common in the manufacturing industry, but few modeling techniques are available to predict materials’ behavior in response to the cutting tool. The paper presents a fracture-based finite element model, named embedded cohesive zone–finite element method (ECZ–FEM). In ECZ–FEM, a network of cohesive zone (CZ) elements are embedded in the material body with regular elements to capture multiple randomized cracks during a cutting process. The CZ element is defined by the fracture energy and a scaling factor to control material ductility and chip behavior. The model is validated by an experimental study in terms of chip formation and cutting force with two different brittle materials and depths of cut. The results show that ECZ–FEM can capture various chip forms, such as dusty debris, irregular chips, and unstable crack propagation seen in the experimental cases. For the cutting force, the model can predict the relative difference among the experimental cases, but the force value is higher by 30–50%. The ECZ–FEM has demonstrated the feasibility of brittle cutting simulation with some limitations applied.

Takabi, Behrouz↗
DOE OSTI2018

High-Temperature Sapphire Pressure Sensors for Harsh Environments

The project aimed to develop and characterize an optical pressure sensor capable of withstanding high temperatures, up to 1000 °C. Survivability in harsh environments was enabled using sapphire for both the substrate and the optical fibers, therefore mitigating thermal expansion mismatch. Such material cannot be machined through conventional microfabrication techniques, and unique additive and subtractive methods had to be studied and implemented. Subtractive methods included laser-machining the sapphire substrate using a picosecond laser machining tool. At first, the parameters of importance were theoretically and experimentally evaluated such that the fluence, the pulse spacing, and the number of passes were found to strongly influence the laser ablation process. A simulation tool was designed to predict the final pattern when inputting a set of parameters and a profile. This model was complemented with laser ablation material physics models and Bayesian uncertainty analysis for one dimensional and three dimensional ablation processes. An additive process called thermocompression bonding was studied in order to seal the diaphragm substrate to the back cavity. At first, bonding strength using this technique was examined using the chevron-notch method that allow to determine the fracture toughness of this bond. This parameter is independent of the geometry used and therefore can be extended to other geometries used in actual sensor bonding. A theoretical and experimental study of internal damage caused by the laser micromachining of the sapphire material was carried out. The research identified enhanced strength and fracture toughness along laser ablated surfaces. This was determined using bend bar specimens and Vicker’s microindentations. Softening of the plasticity behavior (as indicated by nanoindentation experiments) was also investigated using finite deformation plasticity modeling via a contact mechanics finite element simulation. Good correlations between experiments and the model were achieved to determine material parameters that can be used in future MEMS based sensor design and modeling to minimize sensor failure loads. The enhanced strength and fracture toughness was determined to be due to compressive residual stress near laser ablated surfaces. This was confirmed through x-ray spectra taken along different crystal structure planes in combination and Bayesian uncertainty analysis. The Bayesian analysis inferred the entire residual strain tensor (six components) and the residual stress tensor. The results illustrated compressive in-plane stress was the most likely cause of enhanced strength and toughness after laser machining. Due to major equipment failures, a parallel route for fabricating the sensor was taken, using commercially available sapphire removal method (subtractive) and ceramic adhesive bonding (additive). The fully-packaged sensor utilized laser micromachining for the sapphire substrate fabrication and was comprised only of high-temperature-compatible materials. Calibration tests of the sensor in both room-temperature and high-temperature plane-wave tubes showed good linearity and consistent frequency response, in the tested frequency range. The sensor also showed good thermal stability in tests up to 300 °C, proving the fiber-optic lever and its packaging process to be a promising approach for high-temperature measurements. The calibration from 23 °C to 300 °C gave a raw sensitivity of 3.9 × 10 -7 V/Pa and a normalized sensitivity of 9.8 × 10 -8 Pa -1 which were consistent within the temperature range. The sensor had a resonance at 16.3 kHz and a flat band (±3 dB) up to 10.3 kHz. The noise floor of the sensor at 1 kHz was 316 nV/Hz 1/2 , resulting in a MDP of 0.8 Pa. Previously, diaphragm delamination due to CTE mismatch was found at 500 °C, causing packaging failure. Unfortunately high temperature validation in hot jet environment could not be made in time. In the future, it is suggested that the two sapphire substrates be bonded by thermal compression bonding with an intermediate platinum layer instead of using an alumina-based adhesive, and therefore achieve higher temperature capacity up to 1000 °C.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2017

Enhanced Performance of a Self-Powered Organic/Inorganic Photodetector by Pyro-Phototronic and Piezo-Phototronic Effects

Self-powered photodetectors (PDs) have long been realized by utilizing photovoltaic effect and their performances can be effectively enhanced by introducing the piezo-phototronic effect. Recently, a novel pyro-phototronic effect is invented as an alternative approach for performance enhancement of self-powered PDs. Here, a self-powered organic/inorganic PD is demonstrated and the influences of externally applied strain on the pyro-phototronic and the photovoltaic effects are thoroughly investigated. Under 325 nm 2.30 mW cm -2 UV illumination and at a -0.45% compressive strain, the PD's photocurrent is dramatically enhanced from ≈14.5 to ≈103 nA by combining the pyro-phototronic and piezo-phototronic effects together, showing a significant improvement of over 600%. Theoretical simulations have been carried out via the finite element method to propose the underlying working mechanism. Moreover, the pyro-phototronic effect can be introduced by applying a -0.45% compressive strain to greatly enhance the PD's response to 442 nm illumination, including photocurrent, rise time, and fall time. This work provides in-depth understandings about the pyro-phototronic and the piezo-phototronic effects on the performances of self-powered PD to light sources with different wavelengths and indicates huge potential of these two effects in optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Three-body resonances in the φ 4 theory

We study the properties of three-body resonances using a lattice complex scalar φ4 theory with two scalars, with parameters chosen such that one heavy particle can decay into three light ones. We determine the two- and three-body spectra for several lattice volumes using variational techniques, and then analyze them with two versions of the three-particle finite-volume formalism: the Relativistic Field Theory approach and the Finite-Volume Unitarity approach. We find that both methods provide an equivalent description of the energy levels, and we are able to fit the spectra using simple parametrizations of the scattering quantities. By solving the integral equations of the corresponding three-particle formalisms, we determine the pole position of the resonance in the complex energy plane and thereby its mass and width. We find very good agreement between the two methods at different values of the coupling of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Formulation and computation of dynamic, interface-compatible Whitney complexes in three dimensions

A discrete De Rham complex enables compatible, structure-preserving discretizations for a broad range of partial differential equations problems. Such discretizations can correctly reproduce the physics of interface problems, provided the grid conforms to the interface. However, large deformations, complex geometries, and evolving interfaces makes generation of such grids difficult. We develop and demonstrate two formally equivalent approaches that, for a given background mesh, dynamically construct an interface-conforming discrete De Rham complex. Both approaches start by dividing cut elements into interface-conforming subelements but differ in how they build the finite element basis on these subelements. The first approach discards the existing non-conforming basis of the parent element and replaces it by a dynamic set of degrees of freedom of the same kind. The second approach defines the interface-conforming degrees of freedom on the subelements as superpositions of the basis functions of the parent element. These approaches generalize the Conformal Decomposition Finite Element Method (CDFEM) and the extended finite element method with algebraic constraints (XFEM-AC), respectively, across the De Rham complex.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

First-Principles Simulation of the (Li–Ni–Vacancy)O Phase Diagram and Its Relevance for the Surface Phases in Ni-Rich Li-Ion Cathode Materials

Despite several reports on the surface phase transformations from a layered to a disordered spinel and a rock-salt structure at the surface of the Ni-rich cathodes, the precise structures and compositions of these surface phases are unknown. The phenomenon, in itself, is complex and involves the participation of several contributing factors. Of these factors, transition metal (TM) ion migration toward the interior of the particle and hence formation of TM-densified surface layers, triggered by oxygen loss, is thermodynamically probable. Here in this study, we simulate the thermodynamic phase equilibria as a function of TM ion content in the cathode material in the context of lithium nickel oxides, using a combined approach of first-principles density functional calculations, the cluster expansion method, and grand canonical Monte Carlo simulations. We developed a unified lattice Hamiltonian that accommodates not only rock-salt like structures but also topologically different spinel-like structures. Also, our model provides a foundation to investigate metastable cation compositions and kinetics of the phase transformations. Our investigations predict the existence of several Ni-rich phases that were, to date, unknown in the scientific literature. Our simulated phase diagrams at finite temperature show a very low solubility range of the prototype spinel phase. Finally, we find a partially disordered spinel-like phase with far greater solubility that is expected to show very different Li diffusivity compared to that of the prototype spinel structure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A theoretical approach for transient shock strengthening in high-energy-density laser compression experiments

In high-energy-density shock compression experiments, the desired state of compression is typically achieved by passing shock waves through a sequence of different materials. Here, a theoretical approach for transiently strengthening such shocks passing through interfaces of experimentally relevant impedance ratios is examined. A semi-analytical method based on characteristics analysis is used to solve the problem of a shock passing from one material to another through an intermediate region of non-uniform impedance; by appropriately designing this region, a greater shock strength can be achieved in the second material for a finite duration than in the absence of this region. When a shock passes into a material of higher impedance, the shock strengthening increases to a maximum before decreasing to an asymptotic value as the strength of the incident shock is increased. For shocks passing into materials of decreasing impedance, the shock strengthening increases monotonically as the strength of the incident shock is increased when the impedance ratio is above a critical threshold. Incorporating multiple intermediate materials can further increase the strength of the transmitted shock, with an exponential discretization of intermediate material impedances being the most effective distribution for strengthening strong shocks. The results suggest that up to 25% and 9% increases in pressure behind the leading shock can be achieved for materials of increasing and decreasing impedance, respectively. The technique is applied to the design of laser-driven dynamic compression experiments, and the results of the analysis are verified via comparison to simulations performed with the HYADES hydrodynamics code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Different measures of speckle and coherence at the output of a multimode optical fiber

Excitation of a multimode fiber with a focused spatially coherent light of finite bandwidth results in a partially coherent light at the output of the fiber. In this report, we observe the properties of speckle and classical coherence of such light with analytical theory, numerical modeling, and experimentally. Of particular interest is the relationship between measures of coherence and speckle and their dependence on input source bandwidth and fiber length. Speckle contrast is easy to measure experimentally and there exist at least two different methods to generate ensembles of random speckles. We demonstrate that speckle contrast evaluated over the ensemble of external diffusers is related to the number of effective modes—one of the characteristics of beam global coherence. The other speckle contrast measure evaluated over the ensemble of random bends and twists of the fiber is related to residual coherence, which is the pedestal on the average modulus of the complex degree of the coherence function on the output endface of the fiber.

47 OTHER INSTRUMENTATION↗
DOE OSTI2019

Attosecond Streaking Time Delays: Finite-Range Interpretation and Applications

We review theoretical studies of the attosecond streaking time delay concept in photoionization via the investigation of the electron dynamics in the streaking field after the transition of the photoelectron into the continuum upon absorption of an extreme ultraviolet photon. Based on the results, a so-called finite range interpretation was introduced, which highlighted that the delay is accumulated until the streaking pulse ends and, hence, over a finite range of the potential of the parent ion. Following a discussion of the analysis leading to this interpretation, we summarize a few applications which provide insights into different aspects of the streaking time delay concept in photoionization. Besides a review of previously presented results, we give an analysis of the relevance of the first half-cycle of the streaking field and an outlook regarding the perspective of using the streaking method to resolve dynamical changes in the potential that the photoelectron explores during its propagation in the continuum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗
DOE OSTI2019

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↗