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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Degradation of ILW Drums Due to Internal Metallic Corrosion - 18036

Within the British nuclear inventory, metallic Intermediate Level Waste (ILW) is encapsulated in grout and stored in stainless steel drums. Based on recent inspections held in the relevant storage facilities (Sellafield, Cumbria), a proportion of these drums have exhibited considerable distortion, raising concerns about the long-term safety of the packages. Distortion may be caused by the generation of voluminous and potentially flammable substances resulting from metallic corrosion occurring within the drums' core. Potential oxygen influx in the event of steel fracture and consequently drum's failure could potentially cause release or even ignition of the encapsulated radioactive material. The mechanical degradation of the encapsulants (stainless steel and grout) induced from the progressive corrosion of the encapsulated metallic components (and the associated volume expansion) was investigated using both Finite Element Modelling (FEM) techniques as well as experimental trials. A parametric study using simplified 3D FE models was performed to determine their structural integrity under different scenarios. Synchrotron X-Ray tomography (XRT) was used to monitor the behaviour of miniaturized ILW drum-simulants, whilst the latter were subjected to controlled accelerated corrosion in situ. (authors)

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

Declamped Piezoelectric Coefficients in Patterned 70/30 Lead Magnesium Niobate–Lead Titanate Thin Films

Lateral subdivision of blanket piezoelectric thin films increases the functional properties through both increased domain wall mobility and declamping of the intrinsic response. This work presents the local effects of substrate declamping on the piezoelectric coefficient d 33,f of 300 nm thick, rhombohedral, {001}-oriented lead magnesium niobate–lead titanate thin films at the 70/30 composition (70PMN–30PT). Films grown by chemical solution deposition on platinized Si substrates are patterned into strip structures ranging from 0.75 to 9 µm in width. The longitudinal piezoelectric coefficient, d 33,f , is interrogated as a function of position across the patterned structures by three approaches: finite element modeling, piezoresponse force microscopy, and nanoprobe synchrotron X-ray diffraction. It is found that d 33,f increases from the clamped value of 40–50 to ≈160 pm V -1 at the free sidewall under 200 kV cm -1 excitation. The sidewalls partially declamp the piezoelectric response 500–600 nm into the patterned structure, raising the piezoelectric response at the center of features with lateral dimensions less than 1 µm (3:1 width to thickness aspect ratio). The normalized data from all three methods are in excellent agreement, with quantitative differences providing insight to the field dependence of the piezoelectric coefficient and its declamping behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulated effect of defect volume and location on very high cycle fatigue of laser beam powder bed fused AlSi10Mg

This study quantifies the interaction between volumetric defect location and size on the very high cycle fatigue (VHCF) of laser beam powder bed fused (LB-PBF) AlSi10Mg. Crystal plasticity finite element method (CPFEM) simulations were used to investigate the effects of defect location and size on the driving force for crack initiation. The CPFEM model was calibrated against uniaxial and cyclic experimental data of LB-PBF AlSi10Mg. Defect characteristics were informed by experimental data from the specimens produced in various geometries to create realistic representative volume elements (RVEs) with equivalent volume fractions of defects. By embedding defects of varying sizes and locations within the RVEs, fatigue indicator parameters (FIPs) were calculated to analyze the impact of defects’ characteristics on fatigue performance. Different combinations of defect volume and locations were generated for various microstructure instantiations, providing insight into extreme value fatigue responses. Larger defect volumes located on free surfaces consistently generated the highest FIPs, suggesting defect size and boundary proximity intensify stress concentration effects. RVEs with multiple smaller defects produced lower FIPs than those with single large critical defects. These findings underscore the critical role of defect characteristics on fatigue life, providing a foundation for future predictive modeling in fatigue-sensitive AM applications.

AlSi10Mg↗
DOE OSTI2020

Superfluid condensate fraction and pairing wave function of the unitary Fermi gas

The unitary Fermi gas is a many-body system of two-component fermions with zero-range interactions tuned to infinite scattering length. Despite much activity and interest in unitary Fermi gas and its universal properties, there have been great difficulties in performing accurate calculations of the superfluid condensate fraction and pairing wave function. In this paper, we present auxiliary-field lattice Monte Carlo simulations using a lattice interaction which accelerates the approach to the continuum limit, thereby allowing for robust calculations of these difficult observables. As a benchmark test, we compute the ground-state energy of 33 spin-up and 33 spin-down particles. As a fraction of the free Fermi gas energy $E_{\text{FG}}$, we find $E_0/E_{\text{FG}}$ = 0.369(2), 0.372(2), using two different definitions of the finite-system energy ratio, in agreement with the latest theoretical and experimental results. We then determine the condensate fraction by measuring off-diagonal long-range order in the two-body density matrix. We find that the fraction of condensed pairs is α = 0.43(2). Further, we also extract the pairing wave function and find the pair correlation length to be $ζ_pk_F$ = 1.8(3)ℏ, where $k_F$ is the Fermi momentum. Provided that the simulations can be performed without severe sign oscillations, the methods we present here can be applied to superfluid neutron matter as well as more exotic $\textit{P}$-wave and $\textit{D}$-wave superfluids.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Efficient real space formalism for hybrid density functionals

We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn–Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly competitive with the fast Fourier transform scheme while not imposing any restrictions on the boundary conditions. We implement this formalism for both the unscreened and range-separated variants of hybrid functionals. We verify its accuracy and efficiency through comparisons with established planewave codes for isolated as well as bulk systems. In particular, we demonstrate up to an order-of-magnitude speedup in time to solution for the real space method. We also apply the framework to study the structure of liquid water using ab initio molecular dynamics, where we find good agreement with the literature. Overall, the current formalism provides an avenue for efficient real-space DFT calculations with hybrid density functionals.

Chemistry↗
DOE OSTI2021

Damage development in neutron-irradiated concrete in a test reactor: Hygro-thermal and mechanical simulations

This paper reports the development of a 3D mesoscale hygro-thermal-mechanical simulation approach to predict damage in concrete irradiated in a test reactor. This framework, developed in MOOSE, considers the effects of elevated temperature, moisture content, and high neutron fluence (energy threshold, E > 0.1 MeV) on the mortar and aggregates separately. The first-stage simulation implements hygro-thermal analysis to determine the temperature and RH inside the specimen as a function of imposed radiation energy. These are used as inputs to the second stage, which considers radiation-induced volumetric expansion (RIVE) of aggregates, and creep, shrinkage, and stress-strain response of mortar to predict the expansion, stresses, and damage in specimens made using different coarse aggregates and subjected to different irradiation times. The irradiation time-dependent damage in the mortar is expressed using an isotropic damage parameter. This multi-physics model serves as a predictive tool for damage quantification in concrete due to neutron irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Numerical Simulation of Critical Factors Controlling Heat Extraction from Geothermal Systems Using a Closed-Loop Heat Exchange Method

Closed-loop heat exchange for geothermal energy production involves injecting working fluid down a well that extends through the geothermal resource over a significant length to absorb heat by conduction through the well pipe. The well then needs to return to the surface for energy recovery and fluid re-injection to complete the cycle. We have carried out mixed convective-conductive fluid-flow modeling using a wellbore flow model for TOUGH2 called T2Well to investigate the critical factors that control closed-loop geothermal energy recovery. T2Well solves a mixed explicit-implicit set of momentum equations for flow in the pipe with full coupling to the implicit three-dimensional integral finite difference equations for Darcy flow in the porous medium. T2Well has the option of modeling conductive heat flow from the porous medium to the pipe by means of a semi-analytical solution which makes the computation very efficient because the porous medium does not have to be discretized. When the fully three-dimensional option is chosen, the porous medium is discretized and heat flow to the pipe is by conduction and convection, depending on reservoir permeability and other factors. Simulations of the closed-loop system for a variety of parameter values have been carried out to elucidate the heat recovery process. To the extent that convection may occur to aid in heat delivery to the pipe, the permeability of the geothermal reservoir, whether natural or stimulated, is an important property in heat extraction. The injection temperature and flow rate of the working fluid strongly control the ultimate energy recovery. Pipe diameter also plays a strong role in heat extraction, but is correlated with flow rate. Similarly, the choice of working fluid plays an important role, with water showing better heat extraction than CO2 for certain flow rates, while the CO2 has higher pressure at the production wellhead which can aid in surface energy recovery. In general, we find complex interactions between the critical factors that will require advanced computational approaches to fully optimize.

Oldenburg, Curtis↗
DOE OSTI2018

Modeling the diffusion of heat energy within composites of homogeneous materials using the uncertainty principle

The goal of this paper is to develop a highly accurate and efficient numerical method for the solution of a time-dependent partial differential equation with a piecewise constant coefficient, on a finite interval with periodic boundary conditions. The resulting algorithm can be used, for example, to model the diffusion of heat energy in one space dimension, in the case where the spatial domain represents a medium consisting of two homogeneous materials. The resulting model has, to our knowledge, not yet been solved in closed form through analytical methods, and is difficult to solve using existing numerical methods, thus suggesting an alternative approach. The approach presented in this paper is to represent the solution as a linear combination of wave functions that change frequencies at the interfaces between different materials. It is demonstrated through numerical experiments that using the Uncertainty Principle to construct a basis of such functions, in conjunction with a spectral method, a mathematical model for heat diffusion through different materials can be solved much more efficiently than with conventional time-stepping methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Determining Stability Margins in Adiabatic Superconducting Magnets with 3-D Finite Element Analysis

Superconducting magnets play a key role in the development of experiments at Fermilab; understanding the operating stability of these can allow us to utilize more potent magnets for future experiments (like the proposed Muon Collider), optimize the design of magnets in more immediate experiments (like Mu2e), and research the future use of more exotic materials (like high-temperature superconductors). This summer, I developed a 3-D parametric FEA program in ANSYS Mechanical APDL that simulates quench in superconducting magnets, and I also developed a parametric MATLAB program that predicts thermal behavior in magnet quench using the MIITS method. These programs can provide useful quenching parameters (like minimum quench energy and normal zone propagation velocity) for different cases of quench, leading to the previously mentioned objective of magnet design optimization. To test the programs, preliminary cases were run and the data produced was compared and analyzed. The results of these analyses, as well as the program operating methods, are discussed in this project.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Incorporating point defect generation due to jog formation into the vector density-based continuum dislocation dynamics approach

During plastic deformation of crystalline materials, point defects such as vacancies and interstitials are generated by jogs on moving dislocations. A detailed model for jog formation and transport during plastic deformation was developed within the vector density-based continuum dislocation dynamics framework. Here as a part of this model, point defect generation associated with jog transport was formulated in terms of the volume change due to the non-conservative motion of jogs. Balance equations for the vacancies and interstitials including their rate of generation due to jog transport were also formulated. A two-way coupling between point defects and dislocation dynamics was then completed by including the stress contributed by the eigen-strain of point defects. A jog drag stress was further introduced into the mobility law of dislocations to account for the energy dissipation during point defects generation. A number of test problems and a fully coupled simulation of dislocation dynamics and point defects generation and diffusion were performed. The results show that there is an asymmetry of vacancy and interstitial generation due to the different formation energies of the two types of defects. The results also show that a higher hardening rate and a higher dislocation density are obtained when the point defect generation mechanism is coupled to dislocation dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Characterization of latent image of electron beam resist via critical-dimension resonant soft X-ray scattering

Metrology plays a crucial role in semiconductor manufacturing by providing accurate and precise measurement and characterization of critical parameters. With the development of high-resolution extreme ultraviolet lithography (EUVL) processes, critical dimensions are shrinking to sub-10 nm. Resist materials encounter the challenge of providing heightened sensitivity and a handle on exacerbating stochastic variations. A comprehensive understanding of the chemical profile of the latent image is pivotal for mitigating stochastic effects and optimizing pattern quality. However, the subtle differences in chemistry between the exposed and unexposed regions of the resists make it extremely challenging to characterize the latent images with sub-nanometer precision. Here, we develop the metrology with critical-dimension resonant soft X-ray scattering (CD-RSoXS) to probe the chemical profiles of latent images stored in resist after exposure. The combination of absorption spectroscopy and enhanced scattering contrast makes it possible to characterize the subtle structural and chemical variations in the latent image. Moreover, the results of the measurements are compared with the simulations with a finite element method–based Maxwell solver to extract a detailed profile of the latent and developed images. We demonstrate that the CD-RSoXS technique can provide valuable insights into the high spatial resolution and local chemical sensitivity simultaneously, which is crucial to understanding the resolution limits and stochastic effects in EUVL processes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring. We have deployed strainmeters and then conducted injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters (30 to 40m depth) was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 n/d over a few days. The signal at two strainmeters at shallow depth were consistent, although the magnitude of the horizontal strains were different reflecting the different radial directions from the well. The signal at a deep strainmeter deployed at reservoir depth was much different, with tensile vertical strains and compressive horizontal strains. These data can be interpreted by inverting poroelastic forward models developed using numerical and analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. This model is fast and can be inverted to estimate reservoir stiffness and geometry. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. The proxy model is periodically updated and refined using the finite element model to ensure accuracy. This approach significantly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. We have shown that the strain tensor in the caprock is sensitive to pressure in the reservoir, boundaries in the reservoir, and pressure in the caprock caused by leaks. These results indicate that coupling strain tensor data with inversion has the potential to help evaluate reservoirs during initial characterization, and to monitor them during the CO2 injection and storage process.

Murdoch, Larry↗
DOE OSTI2018

Polymer Nanosheet Containing Star‐Like Copolymers: A Novel Scalable Controlled Release System

Abstract Poly(ε‐caprolactone) (PCL)‐based nanomaterials, such as nanoparticles and liposomes, have exhibited great potential as controlled release systems, but the difficulties in large‐scale fabrication limit their practical applications. Among the various methods being developed to fabricate polymer nanosheets (PNSs) for different applications, such as Langmuir–Blodgett technique and layer‐by‐layer assembly, are very effort consuming, and only a few PNSs can be obtained. In this paper, poly(ε‐caprolactone)‐based PNSs with adjustable thickness are obtained in large quantity by simple water exposure of multilayer polymer films, which are fabricated via a layer multiplying coextrusion method. The PNS is also demonstrated as a novel controlled guest release system, in which release kinetics are adjustable by the nanosheet thickness, pH values of the media, and the presence of protecting layers. Theoretical simulations, including Korsmeyer–Peppas model and Finite‐element analysis, are also employed to discern the observed guest‐release mechanisms.

Cao, Peng‐Fei↗
DOE OSTI2016

Unconditionally stable, second-order accurate schemes for solid state phase transformations driven by mechano-chemical spinodal decomposition

Here, we consider solid state phase transformations that are caused by free energy densities with domains of non-convexity in strain-composition space; we refer to the non-convex domains as mechano-chemical spinodals. The non-convexity with respect to composition and strain causes segregation into phases with different crystal structures. We work on an existing model that couples the classical Cahn-Hilliard model with Toupin’s theory of gradient elasticity at finite strains. Both systems are represented by fourth-order, nonlinear, partial differential equations. The goal of this work is to develop unconditionally stable, second-order accurate time-integration schemes, motivated by the need to carry out large scale computations of dynamically evolving microstructures in three dimensions. We also introduce reduced formulations naturally derived from these proposed schemes for faster computations that are still second-order accurate. Although our method is developed and analyzed here for a specific class of mechano-chemical problems, one can readily apply the same method to develop unconditionally stable, second-order accurate schemes for any problems for which free energy density functions are multivariate polynomials of solution components and component gradients. Apart from an analysis and construction of methods, we present a suite of numerical results that demonstrate the schemes in action.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A new numerically stable sequential algorithm for coupled finite-strain elastoplastic geomechanics and flow

In this work, we propose a new numerically stable sequential implicit method for coupled flow and finite-strain multiplicative elastoplastic geomechanics. We find from stability analysis that the sequential method that solves the flow problem first by fixing the first Piola-Kirchhoff total stress, solving the geomechanics (solid deformation) problem at the next step, is unconditionally stable (contractive and B-stable). In this sense, this method named the fixed first Piola-Kirchhoff stress method is an extension of the fixed stress method in coupled flow and infinitesimal geomechanics. We also study the fixed second Piola-Kirchhoff stress method, comparing it with the fixed first Piola-Kirchhoff stress method, because the constitutive relations are formulated by the second Piola-Kirchhoff total stress, although fixing the second Piola-Kirchhoff total stress field does not provide theoretical unconditional stability. In space discretization, we use the finite element method for the geomechanics problem with the total Lagrangian approach, while employing the finite volume method for the flow problem. Geometrical nonlinearity from the total Lagrangian approach results in full-tensor permeability even if the initial permeability is isotropic. To consider the full-tensor permeability accurately, we employ the multipoint flux approximation in solving the flow problem. In time discretization, the backward Euler method is used. Then, we perform numerical experiments of the two fixed-stress sequential methods with various scenarios, and find superiority of their numerical stability. Specifically, the two fixed-stress sequential methods provide almost identical results for elasticity. For J2 plasticity, the two sequential methods yield slightly different numerical results although the trends of the results are similar. For the Drucker-Prager plasticity, we find that the results from the two methods are almost identical. For all the test cases, the two sequential methods are numerically stable.

58 GEOSCIENCES↗
DOE OSTI2020

In situ indentation and high cycle tapping deformation responses in a nanolaminate crystalline/amorphous metal composite

The incorporation of nanostructured and amorphous metals into modern applications is reliant on the understanding of deformation and failure modes in constrained conditions. To study this, a 105 nm crystalline Cu/160 nm amorphous Cu 45 Zr 55 (at.%) multilayer structure was fabricated with the two crystalline layers sputter deposited between the top-middle-bottom amorphous layers and prepared to electron transparency. The multilayer was then in situ indented either under a single load to a depth of ~ 100 nm (max load of ~ 100 μN) or held at 20 μN and then repeatedly indented with an additional 5 μN up to 20,000 cycles in a transmission electron microscope to compare the deformation responses in the nanolaminate. For the single indentation test, the multilayer showed serrated load-displacement behavior upon initial indentation inductive of shear banding. Additionally, at an indentation depth of ~ 32 nm, the multilayer exhibited perfect plastic behavior and no strain hardening. Both indented and fatigue-indented films revealed diffraction contrast changes with deformation. Subsequent Automated Crystal Orientation Mapping (ACOM) measurements confirmed and quantified global texture changes in the crystalline layers with specifically identified grains revealing rotation. Using a finite element model, the in-plane displacement vectors under the indent mapped conditions where ACOM determined grain rotation was observed, indicating the stress flow induced grain rotation. Furthermore, the single indented Cu layers also exhibited evidence of deformation induced grain growth, which was not evident in the fatigue-indented Cu based multilayer. Finally, the single indented multilayer retained a significant plastic crater in the upper most amorphous layer that directly contacted the indenter; a negligible crater impression in the same region was observed in the fatigued tested multilayer. These differences are explained by the different loading methods, applied load, and deformation mechanisms experienced in the multilayers.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Increased Defect Resistance and Ordering in MnBi 2 (Se 1– x Te x ) 4 via Accurate Diffusion Monte Carlo

Stabilizing materials and controlling defect formation remain key challenges in materials science, particularly for theory, where small energy differences must be resolved for accurate predictions. Here, we applied state-of-the-art theoretical methods to topological materials, focusing on MnBi 2 Te 4 (MBT), which is a promising intrinsic magnetic topological insulator. Antisite defects in MBT alter its electronic structure and magnetism, degrading topological properties and causing experimental inconsistencies. Using diffusion Monte Carlo and density functional theory, we investigated the thermodynamic stability and defect formation in MBT, MnBi 2 Se 4 (MBS), and MnBi 2 (Se 1–x Te x ) 4 . We found that MnBi 2 Se 2 Te 2 can be stable at finite temperatures, with higher defect formation energies due to stronger Mn–Se bonding and reduced internal strain. Se preferentially substitutes Te near Mn instead of Te in the outer layer, encouraging long-range ordering when incorporated. For MnBi 2 (Se 1–x Te x ) 4 , cluster expansion phase diagrams revealed solid solution behavior when x <0.5 and phase separation for larger x. MBT and MBS are topological insulators; therefore, the MnBi 2 (Se 1–x Te x ) 4 family could offer tunable topological behavior and improved stability.

MnBi2Te4↗
DOE OSTI2020

Near, far, wherever you are: simulations on the dose efficiency of holographic and ptychographic coherent imaging

Different studies in X-ray microscopy have arrived at conflicting conclusions about the dose efficiency of imaging modes involving the recording of intensity distributions in the near (Fresnel regime) or far (Fraunhofer regime) field downstream of a specimen. A numerical study is presented on the dose efficiency of near-field holography, near-field ptychography and far-field ptychography, where ptychography involves multiple overlapping finite-sized illumination positions. Unlike what has been reported for coherent diffraction imaging, which involves recording a single far-field diffraction pattern, it is found that all three methods offer similar image quality when using the same fluence on the specimen, with far-field ptychography offering slightly better spatial resolution and a lower mean error. These results support the concept that (if the experiment and image reconstruction are done properly) the sample can be near or far; wherever you are, photon fluence on the specimen sets one limit to spatial resolution.

42 ENGINEERING↗