SEARCH · Search NASA

Results for “difference in differences”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A GPU Accelerated Mixed‐Precision Finite Difference Informed Random Walker (FDiRW) Solver for Strongly Inhomogeneous Diffusion Problems

In nature, many complex multi‐physics coupling problems exhibit significant diffusivity inhomogeneity, where one process occurs several orders of magnitude faster than others temporally. Simulating rapid diffusion alongside slower processes demands intensive computational resources due to the necessity for small time steps. To address these computational challenges, we have developed an efficient numerical solver named Finite Difference informed Random Walker (FDiRW). In this study, we propose a GPU‐accelerated, mixed‐precision configuration for the FDiRW solver to maximize efficiency through GPU multi‐threaded parallel computation and lower precision computation. Numerical evaluation results reveal that the proposed GPU‐accelerated mixed‐precision FDiRW solver can achieve a 117× speedup over the CPU baseline, while an additional 1.75× speedup is achieved by employing lower precision GPU computation. Notably, for large model sizes, the GPU‐accelerated mixed‐precision FDiRW solver demonstrates strong scaling with the number of nodes used in simulation. When simulating radionuclide absorption processes by porous wasteform particles with a medium‐sized model of 192 × 192 × 192, this approach reduces the total computational time to 10 min, enabling the simulation of larger systems with strongly inhomogeneous diffusivity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

A Minkowski difference-based advancing front packing technique for generating convex noncircular particles in complex domains

In this work, a Minkowski difference-based advancing front approach is proposed to generate convex and non-circular particles in a predefined computational domain. Two specific algorithms are developed to handle the contact conformity of generated particles with the boundaries of the computational domain. The first, called the open form, is used to handle the smooth contact of generated particles with (external) boundaries, while the other, called the closed form, is proposed to handle the internal boundaries of a computational domain with a complex cavity. The Gilbert-Johnson-Keerthi (GJK) method is used to efficiently solve the contact detection between the newly generated particle at the front and existing particles. Furthermore, the problem of one-sided particle lifting, which can cause some defects in the packing structure in existing advancing front methods during packing generation, is highlighted and an effective solution is developed. Several examples of increasing complexity are used to demonstrate the efficiency and applicability of the proposed packing generation approach. The numerical results show that the generated packing is not only more uniform, but also achieves a higher packing density than existing advancing front methods.

42 ENGINEERING↗

DOE OSTI2022

High-order accurate finite difference discretisations on fully unstructured dual quadrilateral meshes

Here, we present a novel approach for high-order accurate numerical differentiation on unstructured meshes of quadrilateral elements. To differentiate a given function, an auxiliary function with greater smoothness properties is defined which when differentiated provides the derivatives of the original function. The method generalises traditional finite difference methods to meshes of arbitrary topology in any number of dimensions for any order of derivative and accuracy. We demonstrate the accuracy of the numerical scheme using dual quadrilateral meshes and a refinement method based on subdivision surfaces. The scheme is applied to the solution of a range of partial differential equations, including both linear and nonlinear, second and fourth order equations, and a time-dependent first order equation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Flow and transport in three-dimensional discrete fracture matrix models using mimetic finite difference on a conforming multi-dimensional mesh

Here, we present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

A variational mimetic finite difference method for elliptic interface problems on non-matching polytopal meshes with geometric interface inconsistencies

A new variational mimetic finite difference method for elliptic interface problems with perfect and imperfect thermal contacts on non-matching polytopal meshes with geometric interface inconsistencies is developed and analyzed theoretically and numerically. The method is defined on multiple non-matching submeshes with gaps and overlaps along their interfaces. The discrete equations are derived from a minimization problem for the augmented Dirichlet functional. For a perfect thermal contact, the functional uses a modified mimetic gradient with extended stencil which couples unknowns from both sides of an interface, as well as penalty terms to enforce weak continuity of temperature across the interface. The method leads to a symmetric positive definite matrix for any scaling of the penalty terms. For an imperfect thermal contact, the Dirichlet functional is supplemented with a quadratic jump term along the interface related to the interface thermal resistance. We prove that the method conserves the total heat flux across each interface. In conclusion, the obtained results are verified with numerical experiments showing convergence in the discrete L 2 and L ∞ norms.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Spatially-resolved dynamic sampling of different phasic magnetic resonances of nanoparticle ensembles in a magnetotactic bacterium Magnetospirillum magnetotacticum

Abstract Nanoscaled magnetic particle ensembles are promising building blocks for realizing magnon based binary logic. Element-specific real-space monitoring of magnetic resonance modes with sampling rates in the GHz regime is imperative for the experimental verification of future complex magnonic devices. Here we present the observation of different phasic magnetic resonance modes using the element-specific technique of time-resolved scanning transmission x-ray microscopy within a chain of dipolarly coupled Fe 3 O 4 nanoparticles (40–50 nm particle size) inside a single cell of a magnetotactic bacterium Magnetospirillum magnetotacticum . The particles are probed with 25 nm resolution at the Fe L 3 x-ray absorption edge in response to a microwave excitation of 4.07 GHz. A plethora of resonance modes is observed within multiple particle segments oscillating in- and out-of-phase, well resembled by micromagnetic simulations.

time-resolved scanning transmission x-ray microsco↗

DOE OSTI2025

Solving disorder in (3D) real space: a comparative study of the three-dimensional difference pair distribution function and atomic resolution holography reconstructions

The quantitative analysis of local ordering principles in disordered crystalline systems has gained much attention over the past few years, as it is often considered crucial for optimizing material functionality. This development has been driven by significant advancements in computational and experimental methods, which have led to the establishment and widespread use of various analytical techniques. In this study, we perform model calculations to compare the effectiveness of atomic resolution holography and three-dimensional difference pair distribution function analysis (3D-ΔPDF). Using Cu 3 Au as a model system, we demonstrate an approach to derive local order parameters quantitatively and show that both techniques are well suited to quantifying chemical short-range order correlations and local bond-distance variations. By evaluating the strengths and limitations of both techniques, we advocate for their combined use to solve complex short-range order problems accurately.

3D-ΔPDF↗

DOE OSTI2026

Rank-Limiting Strategies for Optimizing Tensor-Train Finite-Difference Time-Domain Simulations

We introduce rank-limiting strategies to optimize tensor-train decompositions for three-dimensional finite-difference time-domain simulations using the relationship between the tensors and their specific dimensionality. These include the use of hard caps on the inner ranks of the tensor train decomposition and the use of a group rounding algorithm taking into account all field components simultaneously. Here, several numerical examples are considered to verify the efficacy of the proposed optimization strategies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Bichromatic Wave Selection for Validation of the Difference-Frequency Transfer Function for the OC6 Validation Campaign: Preprint

The focus of the Offshore Code Comparison Collaboration, Continuation, with Correlation and unCertainity (OC6) project, which operates under the International Energy Agency (IEA) Wind Task 30, is to refine the accuracy of engineering tools used to design offshore wind turbines. In support of this work, a new validation campaign is being developed that seeks to better understand the nonlinear wave loading that excites floating wind systems at their low-frequency, rigid-body modes in surge and pitch. The validation data will be used in a three-way validation between simplified engineering tools and higher-fidelity tools, such as computational fluid dynamics (CFD). Irregular wave spectrums, which are traditionally used to examine the nonlinear wave interaction with offshore structures, are too computationally expensive to be simulated in CFD tools, and so bichromatic wave cases will be used instead. This paper reviews the process used to choose the bichromatic wave pairs to be used in the campaign to validate the second-order difference-frequency quadratic and potential loads at the surge and pitch natural frequencies of a floating semisubmersible.

17 WIND ENERGY↗

DOE OSTI2024

Lattice Green’s Functions for High-Order Finite Difference Stencils

Lattice Green's Functions (LGFs) are fundamental solutions to discretized linear operators, and as such they are a useful tool for solving discretized elliptic PDEs on domains that are unbounded in one or more directions. The majority of existing numerical solvers that make use of LGFs rely on a second-order discretization and operate on domains with free-space boundary conditions in all directions. Under these conditions, fast expansion methods are available that enable precomputation of 2D or 3D LGFs in linear time, avoiding the need for brute-force multi-dimensional quadrature of numerically unstable integrals. Here we focus on higher-order discretizations of the Laplace operator on domains with more general boundary conditions, by (1) providing an algorithm for fast and accurate evaluation of the LGFs associated with high-order dimension-split centered finite differences on unbounded domains, and (2) deriving closed-form expressions for the LGFs associated with both dimension-split and Mehrstellen discretizations on domains with one unbounded dimension. Through numerical experiments we demonstrate that these techniques provide LGF evaluations with near machine-precision accuracy, and that the resulting LGFs allow for numerically consistent solutions to high-order discretizations of the Poisson's equation on fully or partially unbounded 3D domains.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Coarse Mesh Finite Difference Acceleration for Pebble Tracking Transport in Griffin

We implemented a coarse mesh finite difference (CMFD) for accelerating transport calculations with PTT (pebble tracking transport) in the Griffin code. More specifically, extensions for transport update with the consideration of scattering operator and CMFD projection were implemented for PTT. The implementation was verified with a simplified PBR (pebble bed reactor) benchmark problem and significant performance improvements in CPU time was observed.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Effect of Substituent Location on the Relationship between the Transition Dipole Moments, Difference Static Dipole, and Hydrophobicity in Squaraine Dyes for Quantum Information Devices

Aggregates of organic dyes that exhibit excitonic coupling have a wide array of applications, including medical imaging, organic photovoltaics, and quantum information devices. The optical properties of a dye monomer, as a basis of dye aggregate, can be modified to strengthen excitonic coupling. Squaraine (SQ) dyes are attractive for those applications due to their strong absorbance peak in the visible range. While the effects of substituent types on the optical properties of SQ dyes have been previously examined, the effects of various substituent locations have not yet been investigated. In this study, density functional theory (DFT) and time-dependent density functional theory (TD-DFT) were used to investigate the relationships between SQ substituent location and several key properties of the performance of dye aggregate systems, namely, difference static dipole (Δd), transition dipole moment (μ), hydrophobicity, and the angle (θ) between Δd and μ. We found that attaching substituents along the long axis of the dye could increase μ while placement off the long axis was shown to increase Δd and reduce θ. The reduction in θ is largely due to a change in the direction of Δd as the direction of μ is not significantly affected by substituent position. Hydrophobicity decreases when electron-donating substituents are located close to the nitrogen of the indolenine ring. These results provide insight into the structure–property relationships of SQ dyes and guide the design of dye monomers for aggregate systems with desired properties and performance.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2022

Application of linear prolongation to coarse mesh finite difference acceleration in CASMO5

The Coarse-Mesh Finite Difference (CMFD) method has been used for over a decade to accelerate the convergence of the Method of Characteristics (MOC) solution to the two- dimensional particle transport equation in CASMO5. Numerical testing, along with widespread use in production-level calculations, have shown that the current CMFD implementation provides stability and robustness for a wide range of realistic reactor physics problems. However, the recent development of linear prolongation has attracted attention from the community as a way to further improve the performance and stability of CMFD. Two interpolation methods for linear prolongation are presented in this work and implemented into a test version of CASMO5. The performance of the proposed interpolations, referred to as the bilinear and linear directional schemes, is evaluated in terms of runtime relative to the default constant or uniform scaling update. Numerical results indicate that the use of linear prolongation can reduce the transport solver runtime on average by approximately 10% when tested with two hundred randomly selected cases. The new directional linear interpolation, combined with default constant boundary updates, is found to provide the highest reduction in runtime for the cases analyzed. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Facile Synthesis of Palladium-Based Nanocrystals with Different Crystal Phases and a Comparison of Their Catalytic Properties

A relatively unexplored aspect of noble-metal nanomaterials is polymorphism, or their ability to crystallize in different crystal phases. Here, in this work, a method is reported for the facile synthesis of Ru@Pd core–shell nanocrystals featuring polymorphism, with the core made of hexagonally close-packed (hcp)-Ru while the Pd shell takes either an hcp or face-centered cubic (fcc) phase. The polymorphism shows a dependence on the shell thickness, with shells thinner than ≈1.4 nm taking the hcp phase whereas the thicker ones revert to fcc. The injection rate provides an experimental knob for controlling the phase, with one-shot and drop-wise injection of the Pd precursor corresponding to fcc-Pd and hcp-Pd shells, respectively. When these nanocrystals are tested as catalysts toward formic acid oxidation, the Ru@Pd hcp nanocrystals outperform Ru@Pd fcc in terms of both specific activity and peak potential. Density functional theory calculations are also performed to elucidate the origin of this performance enhancement.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Mechanistic Studies on the Insertion of Carbonyl Substrates into Cu‐H: Different Rate‐Limiting Steps as a Function of Electrophilicity

Abstract We report mechanistic studies on the insertion reactions of [(NHC)Cu(μ‐H)] 2 complexes with carbonyl substrates by UV‐vis and 1 H NMR spectroscopic kinetic studies, H/D isotopic labelling, and X‐ray crystallography. The results of these comprehensive studies show that the insertion of Cu‐H with an aldehyde, ketone, activated ester/amide, and unactivated amide consist of two different rate limiting steps: the formation of Cu‐H monomer from Cu‐H dimer for more electrophilic substrates, and hydride transfer from a transient Cu‐H monomer for less electrophilic substrates. We also report spectroscopic and crystallographic characterization of rare Cu‐hemiacetalate and Cu‐hemiaminalate moieties from the insertion of an ester or amide into the Cu−H bond.

Tran, Ba L.↗

DOE OSTI2021

Bimetallic Copper/Ruthenium/Osmium Complexes: Observation of Conformational Differences Between the Solution Phase and Solid State by Atomic Pair Distribution Function Analysis

Abstract High‐energy X‐ray scattering and pair distribution function analysis (HEXS/PDF) is a powerful method to reveal the structure of materials lacking long‐range order, but is underutilized for molecular complexes in solution. We demonstrate the application of HEXS/PDF with 0.26 Å resolution to uncover the solution structure of five bimetallic Cu I /Ru II /Os II complexes. HEXS/PDF of each complex in acetonitrile solution confirms the pairwise distances in the local coordination sphere of each metal center as well as the metal⋅⋅⋅metal distances separated by over 12 Å. The metal⋅⋅⋅metal distance detected in solution is compared with that from the crystal structure and molecular models to confirm that distortions to the metal bridging ligand are unique to the solid state. This work presents the first example of observing sub‐Ångström conformational differences by direct comparison of solution phase and solid‐state structures and shows the potential for HEXS/PDF in the determination of solution structure of single molecules.

Xie, Zhu‐Lin↗

DOE OSTI2020

Mechanistic Studies on the Insertion of Carbonyl Substrates into Cu-H: Different Rate-Limiting Steps as a Function of Electrophilicity

In this paper, we report mechanistic studies on the insertion reactions of [(NHC)Cu(H)] 2 complexes with carbonyl substrates by UV-vis and 1 H NMR spectroscopic kinetic studies, H/D isotopic labelling, and X-ray crystallography. The results of these comprehensive studies show that the insertion of Cu-H with an aldehyde, ketone, activated ester/amide, and unactivated amide consist of two different rate limiting steps: the formation of Cu-H monomer from Cu-H dimer for more electrophilic substrates and hydride transfer from a transient Cu-H monomer for less electrophilic substrates. We also report spectroscopic and crystallographic characterization of rare Cu-hemiacetalate and Cu-hemiaminalate moieties from the insertion of an ester and amide into the Cu-H bond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Bimetallic Copper/Ruthenium/Osmium Complexes: Observation of Conformational Differences Between the Solution Phase and Solid State by Atomic Pair Distribution Function Analysis

High-energy X-ray scattering and pair distribution function analysis (HEXS/PDF) is a powerful method to reveal the structure of materials lacking long-range order, but is underutilized for molecular complexes in solution. Here we demonstrate the application of HEXS/PDF with 0.26 Å resolution to uncover the solution structure of five bimetallic Cu I /Ru II /Os II complexes. HEXS/PDF of each complex in acetonitrile solution confirms the pairwise distances in the local coordination sphere of each metal center as well as the metal...metal distances separated by over 12 Å. The metal…metal distance detected in solution is compared with that from the crystal structure and molecular models to confirm that distortions to the metal bridging ligand are unique to the solid state. This work presents the first example of observing sub-Ångström conformational differences by direct comparison of solution phase and solid-state structures and shows the potential for HEXS/PDF in the determination of solution structure of single molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗