Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite-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.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Seismic H2: Version 1.0

Seismic-H2 is an integrated software package for geological hydrogen reservoir simulation, optimization, and leakage monitoring. The package includes multiple components: (1) code used for modeling seismic wave propagation in 3D heterogeneous elastic media based on finite-difference method to support detection of geological hydrogen storage reservoir leakage; (2) 3D reservoir simulations of leaks from an underground reservoir and 3D simulations of saline aquifers and depleted gas reservoirs; (3) seismic monitoring costs of passive and active seismic monitoring required for UHS; (4) rock physics calculations and interpolations for converting the reservoir simulations from part (2) into the elastic media models in part (1); (5) pre-processing seismic data; and lastly (6), a GUI interface that combines these different components.

Creasy, Neala↗
DOE OSTI2025

Photochemically Induced Acousto-optics Fluid Simulations

PIAFS is a finite-difference code to solve the compressible Navier-Stokes equations with chemical heating on Cartesian grids. It models chemical reactions of air (oxygen and carbon dioxide) with ozone subject to radiation. It uses a high-order WENO spatial discretization and explicit Runge-Kutta time integration. It is capable of parallel simulations using MPI. The code is written in C/C++.

Oudin, AlbertineN [Lawrence Livermore National Lab↗
DOE OSTI2022

A Statistical Interpolation Code for Ocean Analysis and Forecasting

Abstract We present a data assimilation package for use with ocean circulation models in analysis, forecasting, and system evaluation applications. The basic functionality of the package is centered on a multivariate linear statistical estimation for a given predicted/background ocean state, observations, and error statistics. Novel features of the package include support for multiple covariance models, and the solution of the least squares normal equations either using the covariance matrix or its inverse—the information matrix. The main focus of this paper, however, is on the solution of the analysis equations using the information matrix, which offers several advantages for solving large problems efficiently. Details of the parameterization of the inverse covariance using Markov random fields are provided and its relationship to finite-difference discretizations of diffusion equations are pointed out. The package can assimilate a variety of observation types from both remote sensing and in situ platforms. The performance of the data assimilation methodology implemented in the package is demonstrated with a yearlong global ocean hindcast with a 1/4° ocean model. The code is implemented in modern Fortran, supports distributed memory, shared memory, multicore architectures, and uses climate and forecasts compliant Network Common Data Form for input/output. The package is freely available with an open source license from www.tendral.com/tsis/ .

Srinivasan, Ashwanth↗
DOE OSTI2024

ZFP: A compressed array representation for numerical computations

HPC trends favor algorithms and implementations that reduce data motion relative to FLOPS. We investigate the use of lossy compressed data arrays in place of traditional IEEE floating point arrays to store the primary data of calculations. Simulation is fundamentally an exercise in controlled approximation, and error introduced by finite-precision arithmetic (or lossy compression) is just one of several sources of error that need to be managed to ensure sufficient accuracy in a computed result. We describe ZFP, a compressed numerical format designed for in-memory storage of multidimensional arrays, and summarize theoretical results that demonstrate that the error of repeated lossy compression can be bounded and controlled. Furthermore, we establish a relationship between grid resolution and compression-induced errors and show that, contrary to conventional floating point, ZFP reduces finite-difference errors with finer grids. We present example calculations that demonstrate data reduction by 4x or more with negligible impact on solution accuracy. Our results further demonstrate several orders-of-magnitude increase in accuracy using ZFP over IEEE floating point and Posits for the same storage budget.

Lindstrom, Peter↗
DOE OSTI2025

Laser damage of crazed electron-beam high-reflectors following infrared and ultraviolet irradiation in the nanosecond pulse regime

Laser damage of optical components can be a limiting factor in scaling the energetics of high-peak and average power laser systems. Specifically for optical coatings, damage under nanosecond pulsed irradiation is initiated by pre-existing defects in the coating layers, including those that cause discontinuities in the structure, like craze lines. Crazing or cracking in a multilayer dielectric optical coating is induced when the overall coating stress is sufficiently tensile, and is an occasionally observed issue when employing more porous deposition techniques like electron-beam evaporation. Here, in this study, electron-beam high-reflectors were fabricated utilizing process parameters that are known to induce crazing based on prior processing history to systematically evaluate the impact of crazing on reflector damage performance for 1064 and 355 nm lasers. The crazing that was observed was apparently nucleated at nodular defects. When the cross-section of these nodules was investigated, it was observed that there were cracks into the fused silica substrate of approximately 5 µm in depth. The craze lines were irradiated with 1064 and 355 nm light at fluences slightly above the onset of damage initiation fluence of the coating. The 1064 nm irradiated sub-apertures exhibit laser damage but with no spatial correlation with the craze line, whereas the 355 nm irradiated area exhibited many damage sites along the craze line. Finite-difference time-domain electric-field simulations were conducted, and ∼2× field amplification in hafnia was observed for the 355 nm wavelength case. The laser damage can be attributed to a slight electric-field intensification coincidental with an area where UV damage-prone precursors are known to occur. The 355 nm laser damage in uncoated fused silica substrates has been previously correlated to initiate through localized UV absorption at the broken silica bonds in the tips of fractures.

Harthcock, Colin [Lawrence Livermore National Labo↗
DOE OSTI2024

Simple strategy for the simulation of axially symmetric large-area metasurfaces

Metalenses are composed of nanostructures for focusing light and have been widely explored in many exciting applications. However, their expanding dimensions pose simulation challenges. We propose a method to simulate metalenses in a timely manner using vectorial wave and ray tracing models. We sample the metalens’s radial phase gradient and locally approximate the phase profile by a linear phase response. Each sampling point is modeled as a binary blazed grating, employing the chosen nanostructure, to build a transfer function set. The metalens transmission or reflection is then obtained by applying the corresponding transfer function to the incoming field on the regions surrounding each sampling point. Fourier optics is used to calculate the scattered fields under arbitrary illumination for the vectorial wave method, and a Monte Carlo algorithm is used in the ray tracing formalism. We validated our method against finite-difference time domain simulations at 632 nm, and we were able to simulate metalenses larger than 3000 wavelengths in diameter on a personal computer.

Martins, Augusto (ORCID:0000000295546481)↗
DOE OSTI2024

Interplay between gain and loss in arrays of nonlinear plasmonic nanoparticles: toward parametric downconversion and amplification

With the help of a theoretical model and finite-difference time-domain (FDTD) simulations based on the hydrodynamic–Maxwell model, we examine the effect of difference-frequency generation (DFG) in an array of L-shaped metal nanoparticles (MNPs) characterized by intrinsic plasmonic nonlinearity. The outcomes of the calculations reveal the spectral interplay between gain and loss in the vicinity of the fundamental frequency of the localized surface plasmon resonances. Subsequently, we identify different array thicknesses and pumping regimes facilitating parametric amplification and spontaneous parametric downconversion. Our results suggest that the parametric amplification regime becomes feasible on a scale of hundreds of nanometers and spontaneous parametric downconversion on the scale of tens of nanometers, opening up new exciting opportunities for developing building blocks of photonic metasurfaces.

Shah, Syed A.↗
DOE OSTI2025

Pulse duration dependent effects of ultrafast laser induced damage on a 1030 nm multi-layer dielectric mirror for high repetition rate, high average power laser systems

High repetition rate, high peak, and average power laser systems are crucial for next-generation particle accelerators, inertial confinement fusion, and secondary particle sources. These applications demand durable laser optics, particularly interference coatings on optics lasting millions of shots at high fluence. This study focuses on designing, testing, and simulating multi-layer dielectric (MLD) mirrors for pulse durations of 260 fs, 77 fs, and 25 fs at 1030 nm wavelength and 45-degree incidence angle with p -polarization. S-on-1 laser-induced damage thresholds (LIDT) for varying pulse numbers were determined, with single-shot LIDT values of 0.98 Jcm -2 , 1.63 Jcm -2 , and 2.3 Jcm -2 for 25 fs, 77 fs, and 260 fs respectively. A strong correlation between blister shape and local fluence was observed, implying that the layer expansion in a blister depends on local fluence. We have also examined mechanisms responsible for laser-induced stress generation and energy release rates in blister formation. Damage mechanisms are further explored by finite-difference time-domain (FDTD) simulations, incorporating Keldysh strong field ionization, whose predictions were in excellent agreement with the onset of damage determined experimentally. These findings offer insights for enhancing MLD coating technology, promising more efficient and resilient laser systems for diverse scientific and industrial applications.

Noor, Mohamed Yaseen (ORCID:000000021036644X)↗
DOE OSTI2022

Angle- and polarization-resolved luminescence from suspended and hexagonal boron nitride encapsulated MoSe 2 monolayers

The polarized photoluminescence from atomically thin transition metal dichalcogenides is a frequently applied tool to scrutinize optical selection rules and valley physics, yet it is known to sensibly depend on a variety of internal and external material and sample properties. In this work, we apply combined angle- and polarization-resolved spectroscopy to explore the interplay of excitonic physics and phenomena arising from the commonly utilized encapsulation procedure on the optical properties of atomically thin M o S e 2 . We probe monolayers prepared in both suspended and encapsulated manners. We show that the hBN encapsulation significantly enhances the linear polarization of exciton photoluminescence emission at large emission angles. This degree of linear polarization of excitons can increase up to ∼ <#comment/> 17 % <#comment/> in the hBN encapsulated samples. As we confirm by finite-difference time-domain simulations, it can be directly connected to the optical anisotropy of the hBN layers. In comparison, the linear polarization at finite exciton momenta is significantly reduced in a suspended M o S e 2 monolayer, and becomes notable only in cryogenic conditions. This phenomenon strongly suggests that the effect is rooted in the k-dependent anisotropic exchange coupling inherent in 2D excitons. Our results have strong implications on further studies on valley contrasting selection rules and valley coherence phenomena using standard suspended and encapsulated samples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Design of a Nanophotonic Hybrid Coupler

This project aims to design and simulate a nanophotonic hybrid coupler for detecting vacuum ultraviolet (VUV) light with higher efficiency. VUV light is of significant interest because of its applications ranging from high energy physics to space science and electronic industry. However, current VUV photodetectors have very low efficiencies because most materials absorb in the VUV region (100-200 nm). To address this challenge, I have designed a plasmonic-dielectric coupler with wavelength shifting properties, enabling high Purcell enhancement and far-field confinement. The coupler consists of a cylindrical cSi resonator on top of a thin ZnO layer, and Ag and SiO2 are used as substrates. Simulations are done using the finite-difference time domain method to solve Maxwell’s equations, and specifically with Meep and Lumerical. The effectiveness of the coupler is characterized by the Purcell factor and far-field radiation. The Purcell factor reaches a maximum near 390 nm and is in the o rder of 10^4. Furthermore, the coupler also allows partial control of the directionality of the far-field radiation. With the presence of the resonator, the radiation is directed towards the coupled metalens. Overall, our design combines the advantages of plasmonic and dielectric couplers, along with the intrinsic wavelength shifting properties of materials to enhance VUV detection significantly, and such a device can be used in, for example, the DUNE project.

Zhou, Jiayang↗
DOE OSTI2022

In-situ Characterizations of Microstructural Degradation of Perovskite Solar Cells

Rapid progress has been achieved in perovskite solar cells (PSCs), and their efficiencies have improved from 3.8 % to 24.2 % in less than a decade. With low-cost processing, PSCs have shown exciting photovoltaic (PV) properties, such as effective optical absorption, a long carrier lifetime, and unique defect tolerance. While recent studies demonstrated improved stability up to 100 days, PSC technology is still challenged to meet the stringent industry requirements for commercialization. Despite considerable efforts, the underlying physical mechanisms for the inferior stability of PSCs are not well understood. One reason for this divergence is that many established measurement techniques (e.g., quantum efficiency, photoluminescence) probe the properties on length scales far greater than that of electronic and/or structural inhomogeneity (i.e., < 1 μm near grain boundaries) and therefore characterize convoluted and/or averaged properties. Ion/electron beam-based techniques have been extensively used to access the microstructures of PSCs, enabling atomic/nanoscale structural, chemical, optical, and electrical characterizations. For example, focused ion beam (FIB) milling produces an atomically smooth surface that minimizes the artifacts attributed to the surface roughness. FIB techniques can also create a well-defined cross-section of PSCs without mechanical damage in a physical cleaving sample preparation. While powerful, there are some concerns about possible beam damage of inorganic-organic perovskites via chemical-bond breakage and local heating. This project aims to comprehensively understand how the microstructural/interfacial properties of PSCs (e.g., Methylammonium Lead Iodide [MAPbI 3 ]) are modified under the irradiating ion beams. Specifically, we investigate the sub-surface properties of PSCs before and after Ar-ion beam injections. Kelvin probe force microscopy (KPFM) measures the contact potential differences (CPDs). Photoluminescence (PL) microscopy in conjunction with Finite-Difference Time-Domain (FDTD) simulations infers the formation of a “dead layer” (< 15 nm) on the subsurface of MAPbI 3 during Ar+ milling processes while preserving the initial bulk properties. The x-ray photoemission spectroscopy (XPS) confirms this modified surface is a lead-rich and iodine-deficient surface. We initiate customizing in-situ measurement setup while measuring the local optical and electrical properties of PSC under thermal (cooling, heating) and light stressors. Our results provide in-depth knowledge of the ion-beam impact on metal-halide perovskites and how this modified sub-surface impacts their properties under accelerated stressors of light and heat. Intensive Monte Carlo simulations of an electron beam interacting with PSCs provide the beam energy distribution in PSCs, proposing possible measurement conditions of using e-beam with minimizing beam damage. Our in-situ measurement platform can accommodate the diverse architecture of PSC devices for studying deterioration mechanisms under mixed environmental stressors.

14 SOLAR ENERGY↗
DOE OSTI2024

Pore-scale visualization of natural hydrate-bearing sediments

Accurate modeling of gas hydrate reservoir productivity and geomechanical risks associated with subsurface dissociation of natural gas hydrates (NGH) requires the determination of model parameters through physical testing on natural hydrate-bearing sediments (HBS). This involves investigating the hydro-mechanical behavior of undisturbed hydrate samples from nature under in situ conditions using pressure core characterization and analysis, which provides a unique opportunity for research. By employing state-of-the-art micro computed tomography imagery on cryogenically preserved, hydrate-bearing sediment samples, we can determine hydrate saturation as well as permeability with and without the presence of hydrates in the sediment. Furthermore, utilizing a machine learning based image segmentation technique, it is possible to extract pore space and grain information. Subsections of the entire image volume were used to determine anisotropic permeabilities using a finite-difference method Stokes solver (FDMSS). Additionally, permeability measurements on whole pressure and temperature preserved hydrate-bearing core were analyzed by utilizing the National Energy Technology Laboratory’s (NETL) Pressure Core Characterization and X-ray CT Visualization Tool (PCXT) to manipulate, cut, and analyze pressure preserved sediment. Permeabilities were measured under a broad range of vertical stress states to simulate expected pressure changes during production scenarios, and the results show that permeabilities derived from images are in agreement with those from traditional core derived experiments. The collected stress-dependent permeability, permeability anisotropy, and corresponding gas hydrate saturations provide valuable input into numerical simulations of reservoir productivity. These properties have been proven to be key parameters determining a long-term reservoir response under depressurization.

Liu, Mengwei [Oak Ridge Institute for Science and ↗
DOE OSTI2026

Determination of Oscillator Strength of CdSe Nanoplatelets from Strong Light-Matter Coupling

We report exciton-polaritons formed from core-only, core-crown, and core-shell NPLs, which, alongside finite-difference time-domain (FDTD) simulations of Fabry-Pérot cavity mode volume, allow for the verification of core-only and determination of core-crown and core-shell oscillator strengths. Additionally, we show that polaritons can be used not only as a means of altering the properties of coupled emitters, but as an analytical tool to better understand the properties of coupled emitters.

Hancock, Benjamin↗
DOE OSTI2023

Asymmetric Dual-Grating Dielectric Laser Accelerator Optimization

Although hundreds of keV in energy gain have already been demonstrated in dielectric laser accelerators (DLAs), the challenge of creating structures that can confine electrons for multiple millimeters remains. We focus here on dual gratings with single-sided drive, which have experimentally demonstrated energy modulation numerous times. Using a Finite-Difference Time-Domain simulation to find the fields within various DLA structures and correlating these results with particle tracking simulation, we look at the impact of teeth height and width, as well as gap and offset, on the performance of these structures. We find a tradeoff between electron throughput and acceleration; however, we also find that for any given grating geometry, there is a gap and offset that will allow some charge acceleration. For our 780 nm laser wavelength, this results in a 1200 nm optimal gap size for most gratings.

Crisp, Sophie (ORCID:0000000327503727)↗
DOE OSTI2022

Ultrafast Laser Material Damage Simulation—A New Look at an Old Problem

The chirped pulse amplification technique has enabled the generation of pulses of a few femtosecond duration with peak powers multi-Tera and Peta–Watt in the near infrared. Its implementation to realize even shorter pulse duration, higher energy, and higher repetition rate laser systems relies on overcoming the limitations imposed by laser damage of critical components. In particular, the laser damage of coatings in the amplifiers and in post-compression optics have become a bottleneck. The robustness of optical coatings is typically evaluated numerically through steady-state simulations of electric field enhancement in multilayer stacks. However, this approach cannot capture crucial characteristics of femtosecond laser induced damage (LID), as it only considers the geometry of the multilayer stack and the optical properties of the materials composing the stack. This approach neglects that in the interaction of an ultrashort pulse and the materials there is plasma generation and associated material modifications. Here, we present a numerical approach to estimate the LID threshold of dielectric multilayer coatings based on strong field electronic dynamics. In this dynamic scheme, the electric field propagation, photoionization, impact ionization, and electron heating are incorporated through a finite-difference time-domain algorithm. We applied our method to simulate the LID threshold of bulk fused silica, and of multilayer dielectric mirrors and gratings. The results are then compared with experimental measurements. The salient aspects of our model, such as the implementation of the Keldysh photoionization model, the impact ionization model, the electron collision model for ‘low’-temperature, dense plasma, and the LID threshold criterion for few-cycle pulses are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hot Electron-Driven Photocatalysis Using Sub-5 nm Gap Plasmonic Nanofinger Arrays

Semiconductor photocatalysis has received increasing attention because of its potential to address problems related to the energy crisis and environmental issues. However, conventional semiconductor photocatalysts, such as TiO2 and ZnO, can only be activated by ultraviolet light due to their wide band gap. To extend the light absorption into the visible range, the localized surface plasmon resonance (LSPR) effect of noble metal nanoparticles (NPs) has been widely used. Noble metal NPs can couple incident visible light energy to strong LSPR, and the nonradiative decay of LSPR generates nonthermal hot carriers that can be injected into adjacent semiconductor material to enhance its photocatalytic activity. Here we demonstrate that nanoimprint-defined gap plasmonic nanofinger arrays can function as visible light-driven plasmonic photocatalysts. The sub-5 nm gaps between pairs of collapsed nanofingers can support ultra-strong plasmon resonance and thus boost the population of hot carriers. The semiconductor material is exactly placed at the hot spots, providing an efficient pathway for hot carrier injection from plasmonic metal to catalytic materials. This nanostructure thus exhibits high plasmon-enhanced photocatalytic activity under visible light. The hot carrier injection mechanism of this platform was systematically investigated. The plasmonic enhancement factor was calculated using the finite-difference time-domain (FDTD) method and was consistent with the measured improvement of the photocatalytic activity. This platform, benefiting from the precise controllable geometry, provides a deeper understanding of the mechanism of plasmonic photocatalysis.

42 ENGINEERING↗
DOE OSTI2025

Jipole: A Differentiable ipole-based Code for Radiative Transfer in Curved Spacetimes

Recent imaging of supermassive black holes by the Event Horizon Telescope has relied on exhaustive parameter-space searches, matching observations to large, precomputed libraries of theoretical models. As observational data become increasingly precise, the limitations of this computationally expensive approach grow more acute, creating a pressing need for more efficient methods. In this work, we present Jipole, an automatically differentiable (AD), ipole-based code for radiative transfer in curved spacetimes, designed to compute image gradients with respect to underlying model parameters. These gradients quantify how parameter changes—such as the black hole’s spin or the observer’s inclination—affect the image, enabling more efficient parameter estimation and reducing the number of required images. We validate Jipole against ipole in two analytical tests and then compare pixelwise intensity derivatives from AD with those from finite-difference methods. We then demonstrate the utility of these gradients by performing parameter recovery for an analytical model in three increasingly complex cases for the injected image: ideal, blurred, and blurred with added noise. In most cases, high-accuracy fits are obtained in only a few optimization steps, failing only in cases with extremely low signal-to-noise ratios. These results highlight the potential of AD-based methods to accelerate robust, high-fidelity model-data comparisons in current and future black hole imaging efforts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Tree hydrodynamic modelling of the soil–plant–atmosphere continuum using FETCH3

Modelling the water transport along the soil–plant–atmosphere continuum is fundamental to estimating and predicting transpiration fluxes. A Finite-difference Ecosystem-scale Tree Crown Hydrodynamics model (FETCH3) for the water fluxes across the soil–plant–atmosphere continuum is presented here. The model combines the water transport pathways into one vertical dimension, and assumes that the water flow through the soil, roots, and above-ground xylem can be approximated as flow in porous media. This results in a system of three partial differential equations, resembling the Richardson–Richards equation, describing the transport of water through the plant system and with additional terms representing sinks and sources for the transfer of water from the soil to the roots and from the leaves to the atmosphere. The numerical scheme, developed in Python 3, was tested against exact analytical solutions for steady state and transient conditions using simplified but realistic model parameterizations. The model was also used to simulate a previously published case study, where observed transpiration rates were available, to evaluate model performance. With the same model setup as the published case study, FETCH3 results were in agreement with observations. Through a rigorous coupling of soil, root xylem, and stem xylem, FETCH3 can account for variable water capacitance, while conserving mass and the continuity of the water potential between these three layers. FETCH3 provides a ready-to-use open access numerical model for the simulation of water fluxes across the soil–plant–atmosphere continuum.

54 ENVIRONMENTAL SCIENCES↗