Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite-Difference Methods”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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 OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Harnessing complexity: Nonlinear optical phenomena in L-shapes, nanocrescents, and split-ring resonators

Here, we conduct systematic studies of the optical characteristics of plasmonic nanoparticles that exhibit C 2v symmetry. In particular, we analyze three distinct geometric configurations: an L-type shape, a crescent, and a split-ring resonator shaped like the Greek letter π. Optical properties are examined using the finite-difference time-domain method. It is demonstrated that all three shapes exhibit two prominent plasmon modes associated with the two axes of symmetry. This is in addition to a wide range of resonances observed at high frequencies corresponding to quadrupole modes and peaks due to sharp corners. Next, to facilitate nonlinear analysis, we employ a semiclassical hydrodynamic model, where the electron pressure term is explicitly accounted for. This model goes beyond the standard Drude description and enables capturing nonlocal and nonlinear effects. Employing this model enables us to rigorously examine the second-order angular resolved nonlinear optical response of these nanoparticles in each of the three configurations. Two pumping regimes are considered, namely, continuous wave (CW) and pulsed excitations. For CW pumping, we explore the properties of the second harmonic generation (SHG). Polarization and angle-resolved SHG spectra are obtained, revealing strong dependence on the nanoparticle geometry and incident wave polarization. The C 2v symmetry is shown to play a key role in determining the polarization states and selection rules of the SHG signal. For pulsed excitations, we discuss the phenomenon of broadband terahertz (THz) generation induced by the difference-frequency generation . It is shown that the THz emission spectra exhibit unique features attributed to the plasmonic resonances and symmetry of the nanoparticles. The polarization of the generated THz waves is also examined, revealing interesting patterns tied to the nanoparticle geometry. To gain deeper insight, we propose an analytical theory that agrees very well with the numerical experiments. The theory shows that the physical origin of the THz radiation is the mixing of various frequency components of the fundamental pulse by the second-order nonlinear susceptibility. An expression for the far-field THz intensity is derived in terms of the incident pulse parameters and the nonlinear response tensor of the nanoparticle. The results presented in this work offer new insights into the linear and nonlinear optical properties of nanoparticles with C 2v symmetry. The demonstrated strong SHG response and efficient broadband THz generation hold great promise for applications in nonlinear spectroscopy, nanophotonics, and optoelectronics. The proposed theoretical framework also provides a valuable tool for understanding and predicting the nonlinear behavior of other related nanostructures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Development of an Optical Library for Coevaporated CdSe x Te 1– x

The conversion efficiency of CdTe solar cells may be improved by bandgap engineering, i.e., changing the bandgap value through the addition of Se in the absorber. The Se alloying enables a short-circuit current density improvement, as it leads to a bandgap energy value decrease. Furthermore, it has been associated with increased minority carrier lifetimes, assuring high open-circuit voltage values. An Se gradient profile control can further optimize the solar cell performance. Thus, an optical model baseline of the CdSe x Te 1–x (CST) compound was developed. Spectroscopic ellipsometry measurements were conducted to accurately extract the optical constants of ten CST layers deposited through coevaporation with x varying from 0 to 1. Using the measured dielectric function spectra from the discrete CST layers with varying x, and considering the composition-induced shift in the critical point energies, an energy-shift model was employed to develop the accurate optical library for the CST compound for any x value to provide data for future modeling and optimization. Furthermore, the library accuracy was validated through optical simulations of the quantum efficiency of a graded CST solar cell using the finite-difference time-domain method by replicating the Se profile in the absorber layer measured through secondary ion mass spectrometry.

14 SOLAR ENERGY↗
DOE OSTI2024

Finite-frequency modeling of regional tropospheric infrasound using realistic atmospheres and terrain

Infrasonic waves have been observed to propagate to regional (greater than 15 km) distances through the troposphere. Infrasound propagation in the geometric acoustics approximation has shown that realistic terrain can scatter acoustic energy from tropospheric ducts; however, ray methods cannot intrinsically capture finite-frequency behavior such as diffraction. A two-dimensional finite-difference time-domain (FDTD) method has been developed to solve linearized equations for infrasound propagation with realistic terrain. Acoustic wave propagation over 100 km with both flat terrain and a Gaussian hill was first simulated in order to compare finite-frequency propagation with ray predictions. The effects of realistic terrain and atmospheres on infrasound signals from a 2012 surface explosion at the Utah Testing and Training Range are then investigated. Propagation through the troposphere is suggested by array processing results, but eigenrays are not predicted due to weak to nonexistent ducting conditions. FDTD modeling suggests that the inclusion of terrain and finite frequency effects helps explain much of the observed signal in a realistic scenario. Furthermore, these results suggest that geometric acoustics may underestimate propagation through the troposphere, and that recorded waveforms at regional distances may be noticeably affected by terrain.

58 GEOSCIENCES↗
DOE OSTI2023

Improved Water Sound Synthesis using Coupled Bubbles

We introduce a practical framework for synthesizing bubble-based water sounds that captures the rich inter-bubble coupling effects responsible for low-frequency acoustic emissions from bubble clouds. We propose coupled-bubble oscillator models with regularized singularities, and techniques to reduce the computational cost of time stepping with dense, time-varying mass matrices. Airborne acoustic emissions are estimated using finite-difference time-domain (FDTD) methods. We propose a simple, analytical surface-acceleration model, and a sample-and-hold GPU wavesolver that is simple and faster than prior CPU wavesolvers. Sound synthesis results are demonstrated using bubbly flows from incompressible, two-phase simulations, as well as procedurally generated examples using single-phase FLIP fluid animations. Our results demonstrate sound simulations with hundreds of thousands of bubbles, and perceptually significant frequency transformations with fuller low-frequency content.

Computer Science↗
DOE OSTI2025

Application of mesh refinement to relativistic magnetic reconnection

During relativistic magnetic reconnection, antiparallel magnetic fields undergo a rapid change in topology, releasing a large amount of energy in the form of non-thermal particle acceleration. This work explores the application of mesh refinement to 2D reconnection simulations to efficiently model the inherent disparity in length-scales. We have systematically investigated the effects of mesh refinement and determined necessary modifications to the algorithm required to mitigate non-physical artifacts at the coarse–fine interface. We have used the ultrahigh-order pseudo-spectral analytical time-domain Maxwell solver to analyze how its use can mitigate the numerical dispersion that occurs with the finite-difference time-domain (or “Yee”) method. Absorbing layers are introduced at the coarse–fine interface to eliminate spurious effects that occur with mesh refinement. We also study how damping the electromagnetic fields and current density in the absorbing layer can help prevent the non-physical accumulation of charge and current density at the coarse–fine interface. Using a mesh refinement ratio of 8 for two-dimensional magnetic reconnection simulations, we obtained good agreement with the high-resolution baseline simulation, using only 36% of the macroparticles and 71% of the node-hours needed for the baseline. The methods presented here are especially applicable to 3D systems where higher memory savings are expected than in 2D, enabling comprehensive, computationally efficient 3D reconnection studies in the future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Mechanism of structural colors in binary mixtures of nanoparticle-based supraballs

Inspired by structural colors in avian species, various synthetic strategies have been developed to produce noniridescent, saturated colors using nanoparticle assemblies. Nanoparticle mixtures varying in particle chemistry and size have additional emergent properties that affect the color produced. For complex multicomponent systems, understanding the assembled structure and a robust optical modeling tool can empower scientists to identify structure-color relationships and fabricate designer materials with tailored color. Here, we demonstrate how we can reconstruct the assembled structure from small-angle scattering measurements using the computational reverse-engineering analysis for scattering experiments method and use the reconstructed structure in finite-difference time-domain calculations to predict color. We successfully, quantitatively predict experimentally observed color in mixtures containing strongly absorbing nanoparticles and demonstrate the influence of a single layer of segregated nanoparticles on color produced. The versatile computational approach that we present is useful for engineering synthetic materials with desired colors without laborious trial-and-error experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Enhancing two-photon spontaneous emission in rare earths using graphene and graphene nanoribbons

The enhancement of two-photon spontaneous emission (2PSE) from trivalent and divalent rare earth ions in proximity to graphene and graphene nanoribbons is calculated for achievable experimental conditions using a combination of finite difference time domain simulations and direct computation of transition rates between energy levels in rare earths. For Er 3+ , we find that the 2PSE rate is initially 8 orders lower than the single-photon spontaneous emission rate but that, with enhancement, 2PSE can reach 2.5% of the overall decay. When graphene nanoribbons are used, we also show that the emission of free-space photon pairs from Er 3+ at 3–3.2 μ⁢m via 2PSE can be increased by ~400. Our calculations show significantly less relative graphene-enhanced 2PSE than previous works, and we attribute this variation to differences in emitter size and assumed graphene mobility. Here, we also show that the internal energy structure of the ion can have an impact on the degree of 2PSE enhancement achievable and find that divalent rare earths are more favorable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Inverse modeling of circular lattices via orbit response measurements in the presence of degeneracy

The number and location of beam position monitors (BPMs) and steerers with respect to the quadrupoles in a circular lattice can lead to degeneracy in the context of fitting linear optics and extracting lattice information from measured closed orbits. Furthermore, the measurement uncertainties due to the imperfection of BPMs and steerers can be propagated by the fitting process in ways that prohibit the successful extraction of discrepancies between lattice elements in the real machine and their description in the corresponding model. We systematically studied the influence of the placement of BPMs and steerers on the reconstruction of linear optics and corresponding lattice information. The derivative of orbit response coefficients with respect to the quadrupole strengths, the Jacobian, is derived as an analytical formula. This analytical version of the Jacobian is used to further derive the theoretical limitations of fitting linear optics from closed orbits in terms of the placement of BPMs and steerers. It is further demonstrated that when evaluating the Jacobian during the fitting procedure, the analytical version can be used in place of the conventional finite-difference computation. This allows for greatly improved efficiency when computing the Jacobian during each iteration of the fitting procedure. The approach is tested with large-scale simulations and the findings are verified by measurement data taken on SIS18 synchrotron at GSI Helmholtz Centre for Heavy Ion Research. The presented methods are of general nature and can be applied to other accelerator lattices as well. The fitting procedure by using the analytical Jacobian is tested in conjunction with various methods for mitigating quasidegeneracy and the results agree with those obtained by using the conventional Jacobian via finite-difference approximation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

An Empirical Quantile Estimation Approach for Chance-Constrained Nonlinear Optimization Problems

We investigate an empirical quantile estimation approach to solve chance-constrained nonlinear optimization problems. Our approach is based on the reformulation of the chance constraint as an equivalent quantile constraint to provide stronger signals on the gradient. In this approach, the value of the quantile function is estimated empirically from samples drawn from the random parameters, and the gradient of the quantile function is estimated via a finite-difference approximation on top of the quantile-function-value estimation. We establish a convergence theory of this approach within the framework of an augmented Lagrangian method for solving general nonlinear constrained optimization problems. The foundation of the convergence analysis is a concentration property of the empirical quantile process, and the analysis is divided based on whether or not the quantile function is differentiable. In contrast to the sampling-and-smoothing approach used in the literature, the method developed in this paper does not involve any smoothing function and hence the quantile-function gradient approximation is easier to implement and there are less accuracy-control parameters to tune. Furthermore, we demonstrate the effectiveness of this approach and compare it with a smoothing method for the quantile-gradient estimation. Numerical investigation shows that the two approaches are competitive for certain problem instances.

Applied Probability↗
DOE OSTI2025

Decoupling Carrier Dynamics and Energy Transport in Ultrafast Near-Field Nanoscopy

Ultrafast near-field optical nanoscopy has emerged as a powerful platform to characterize low-dimensional materials. While analytical and numerical models have been established to account for photoexcited carrier dynamics, quantitative evaluation of the associated pulsed laser heating remains elusive. Here, we decouple the photocarrier density and temperature increase in near-field nanoscopy by integrating the two-temperature model (TTM) with finite-difference time-domain (FDTD) simulations. These results reveal that the electron–phonon coupling in a silicon film after femtosecond laser excitation is most pronounced within approximately 3 ps–substantially shorter than the photocarrier decay time scale at tens of picoseconds. Moreover, the coupled TTM-FDTD method indicates that ultrafast laser heating can cause up to a 14% variation in the near-field signal at a 220 μJ/cm 2 pump pulse fluence. Our numerical results are further validated by transient experiments, highlighting the potential of this method for investigations of carrier and thermal phenomena in emerging nanomaterials and nanodevices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Inverse aqueous transport modeling for emergency response

ALGE is a three-dimensional, finite-difference aqueous transport model that simulates pollutant fate and transport in lakes, rivers, bays, and estuaries by solving the prognostic equations of mass, momentum, and energy. Its current modeling capabilities include transport of dissolved tracer for a series of predefined basins across the continental United States. Recently, an inverse method (also known as backtracking) has been added to ALGE to provide a possible source of a pollutant should one be detected by a sensor in a body of water and a source is not known. This inverse method is a three step process that uses an algorithm to inverse the flow. We demonstrate the new model’s capabilities through simulating the 2021 Piney Point spill in Tampa Bay, Florida (USA). This involves moving tracer backwards from its detection points, encompassing a potential source area, and applying Bayes’ Theorem and $\frac{𝜒}{𝑄}$ to reduce the area within which the true source could be located.

hydrological modeling↗
DOE OSTI2025

HFBTHO-AD: Differentiation of a nuclear energy density functional code

The HFBTHO code implements a nuclear energy density functional solver to model the structure of atomic nuclei. HFBTHO has previously been used to calibrate energy functionals and perform sensitivity analysis by using derivative-free methods. To enable derivative-based optimization and uncertainty quantification approaches, we must compute the derivatives of HFBTHO outputs with respect to the parameters of the energy functional, which are a subset of all input parameters of the code. Here, we use the algorithmic/automatic differentiation (AD) tool Tapenade to differentiate HFBTHO. We compare the derivatives obtained using AD against finite-difference approximation and examine the performance of the derivative computation.

Algorithmic differentiation↗
DOE OSTI2023

Origin of high optical contrast in zinc-zinc oxide electrodeposits for dynamic windows

The control of solar light and heat emission through windows is an important strategy for increasing the energy efficiency of buildings. Reversible Zn electrodeposition has recently emerged as a promising method for constructing electronically tintable robust dynamic windows. In Zn electrodeposits formed from dimethyl sulfoxide (DMSO) electrolytes during device tinting, we observe extraordinary absorption that is in excess of what is predicted by the Beer-Lambert law for a uniform Zn thin film. Here, the charge required to electroplate these films is abnormally low, significantly less than previously reported dynamic windows based on reversible metal electrodeposition, which facilitates the construction of large-area devices that switch uniformly. Finite-difference time-domain (FDTD) simulations are used to investigate the origin of this enhanced absorption, which arises from plasmonic effects among the Zn nanoparticles and ZnO dendrites in the film. The dielectric ZnO dendrites promote absorption via slit-like Zn-dielectric-Zn structures that from hybrid surface plasmon resonance at metal walls. Through these investigations, we provide design principles to construct low-charge and high-contrast metal and metal oxide-based dynamic windows.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
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 OSTI2026

Benchmarking of massively parallel phase-field codes for directional solidification

We present a detailed benchmark comparing two state-of-the-art phase-field implementations for simulating alloy solidification under experimentally relevant conditions. The study investigates the directional solidification of Al-3wt%Cu under high-velocity solidification conditions and SCN-0.46wt% camphor under microgravity conditions from National Aeronautics and Space Administration (NASA) DECLIC-DSI-R experiments. Both codes, one employing finite-difference discretization with uniform mesh and GPU-acceleration (GPU-PF) and the other one employing finite-element discretization with adaptive-mesh and CPU-parallelization (PRISMS-PF), solve the same quantitative phase-field formulation that incorporates an anti-trapping current for the solidification of dilute alloys. We evaluate the predictions of each code for dendritic morphology, primary spacing, and tip dynamics in both 2D and 3D, as well as their numerical convergence and computational performance. While existing benchmark problems have primarily focused on simplified or small-scale simulations, they do not reflect the computational and modeling challenges posed by employing experimentally relevant time and length scales. Our results provide a practical framework for assessing phase-field code performance as well as validating and facilitating their application in integrated computational materials engineering (ICME) workflows that require integration with realistic experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗