Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite-difference time-domain (FDTD)”

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 19 records

Field–potential finite-difference time-domain (FiPo FDTD) technique for computational electromagnetics

Modeling light–matter interactions at the nanoscale requires accurate handling of coupled quantum and electromagnetic systems. This coupling requires information about the electric scalar potential Φ and the magnetic vector potential A, which are not typically calculated in standard computational electromagnetics implementations. To that end, we have developed a field–potential finite-difference time-domain (FiPo FDTD) algorithm, which solves a set of first-order equations for Φ and A alongside equations for the electric and magnetic fields E and H. The FiPo Basic code is essentially conventional FDTD, but with an added module that calculates the potentials. The FiPo Hybrid code self-consistently calculates both fields and potentials and is particularly suitable for coupling with quantum electronic transport solvers because it can be sourced by the potentials themselves. To terminate the domain and mimic infinite space, we have derived and implemented a convolutional perfectly matched layer (CPML) absorbing boundary condition for FiPo FDTD whose performance is on par with state-of-the-art CPMLs for standard FDTD. We present FiPo simulation results on several example systems.

Avazpour, L. [University of Wisconsin-Madison, WI ↗

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↗

Investigating relationship between surface topography and emissivity of metallic additively manufactured parts

Due to the direct relationship between thermal history and mechanical behavior, in situ thermal monitoring is key in gauging quality of parts produced with additive manufacturing (AM). Accurate monitoring of temperatures in an AM process requires knowledge of environment and object parameters including object emissivity. The emissivity is dependent on several variables, including: wavelength, material composition, temperature, and surface topography. Researchers have been concerned with the thermal emissivity dependence on temperature since large ranges are seen in metal powder bed processes, but there is also an extensive range of surfaces produced by AM. This work focused on discovering what roughness characteristics control thermal emissivity through investigation of prototypic 316 stainless steel AM samples produced with a range of build conditions on a laser powder bed fusion machine. Through experimental measurements of emissivity using hemispherical directional reflectance (HDR), guided by simulations using a finite-difference time-domain (FDTD) Maxwell solver, it was found that combinations of existing roughness parameters describing both height and slope of the surface correlate well with emissivity changes. Furthermore, these parameters work well due to their apt description of surface features encouraging internal reflection, which is the phenomenon that increases emissivity when a surface falls under the geometric optical region conditions.

36 MATERIALS SCIENCE↗

Hybrid SERS platform by adapting both chemical mechanism and electromagnetic mechanism enhancements: SERS of 4-ATP and CV by the mixture with GQDs on hybrid PdAg NPs

Surface-enhanced Raman spectroscopy (SERS) is a versatile analytical technique widely adapted for the identification of target analytes even at an extremely low concentration. Herein, a unique hybrid SERS platform is demonstrated by utilizing the graphene quantum dots (GQDs) and hybrid PdAg NPs for the detection of 4-amino-thiophenol (4-ATP) and crystal violet (CV). GQDs offer the rich charge transfer to the LUMO and HOMO of 4-ATP and CV by the CM enhancement. A significant EM enhancement is achieved by the strong localized surface plasmon resonance (LSPR) of uniquely designed hybrid core-shell PdAg NP, possessing the narrow nano-gaps with the highly dense small Ag NPs and Pd core Ag shell configuration. As a result, orders of enhancement of originally weak 4-ATP and CV Raman signal is witnessed by the combination of electromagnetic mechanism (EM) and chemical mechanism (CM) enhancements on a single SERS substrate. The hybrid GQD/HNP SERS substrate demonstrates significantly improved hot spots and strongly localized electromagnetic fields as confirmed by the finite-difference time-domain (FDTD) simulations. There for a proper mixture ratio of GQDs and probe molecules presents largely enhanced SERS signals by a strong adsorption of probe molecules through π - π interaction and the enhancement factors vary by the mixture ratio depending on the adsorption kinetics and Raman cross-section of molecules.

36 MATERIALS SCIENCE↗

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↗

Scalable freeform optimization of wide-aperture 3D metalenses by zoned discrete axisymmetry

We introduce a novel framework for design and optimization of 3D freeform metalenses that attains nearly linear scaling of computational cost with diameter, by breaking the lens into a sequence of radial “zones” with 𝑛-fold discrete axisymmetry, where 𝑛 increases with radius. This allows vastly more design freedom than imposing continuous axisymmetry, while avoiding the compromises of the locally periodic approximation (LPA) or scalar diffraction theory. Using a GPU-accelerated finite-difference time-domain (FDTD) solver in cylindrical coordinates, we perform full-wave simulation and topology optimization within each supra-wavelength zone. We validate our approach by designing millimeter and centimeter-scale, poly-achromatic, 3D freeform metalenses which outperform the state of the art. By demonstrating the scalability and resulting optical performance enabled by our “zoned discrete axisymmetry” (ZDA) and supra-wavelength domain decomposition, we highlight the potential of our framework to advance large-scale meta-optics and next-generation photonic technologies.

Sun, Mengdi [Wesleyan University]↗

Finite-difference time-domain methods

The finite-difference time-domain (FDTD) method is a widespread numerical tool for full-wave analysis of electromagnetic fields in complex media and for detailed geometries. Applications of the FDTD method cover a range of time and spatial scales, extending from subatomic to galactic lengths and from classical to quantum physics. Technology areas that benefit from the FDTD method include biomedicine — bioimaging, biophotonics, bioelectronics and biosensors; geophysics — remote sensing, communications, space weather hazards and geolocation; metamaterials — sub-wavelength focusing lenses, electromagnetic cloaks and continuously scanning leaky-wave antennas; optics — diffractive optical elements, photonic bandgap structures, photonic crystal waveguides and ring-resonator devices; plasmonics — plasmonic waveguides and antennas; and quantum applications — quantum devices and quantum radar. This Primer summarizes the main features of the FDTD method, along with key extensions that enable accurate solutions to be obtained for different research questions. Additionally, hardware considerations are discussed, plus examples of how to extract magnitude and phase data, Brillouin diagrams and scattering parameters from the output of an FDTD model. Furthermore, the Primer ends with a discussion of ongoing challenges and opportunities to further enhance the FDTD method for current and future applications.

42 ENGINEERING↗

Formation and optical properties of indium nanoparticle arrays for deep-UV plasmonics

We utilize a combined computational-experimental approach to examine the influence of indium nanoparticle (NP) array distributions on deep-ultraviolet (UV) plasmon resonances. For photon energies < 5.7 eV, analysis of ellipsometric spectra reveals an increase in silicon reflectance induced by indium NP arrays on silicon. For various energies in the range 5.7–7.0 eV, a decrease in reflectance is induced by the NP arrays. Similar trends in reflectance are predicted from finite-difference time-domain (FDTD) simulations using NP size distributions extracted from atomic-force micrographs as input. In addition, in the energy range of 7.4–9.2 eV, the FDTD simulations reveal reflectance minima, characteristic of localized surface plasmon resonances. Here, electron energy-loss spectroscopy collected from individual indium NPs reveals the presence of LSPR at ≈ 8 eV, further supporting the promise of indium NP arrays on silicon for deep-UV plasmonics.

36 MATERIALS SCIENCE↗

A transient near to far field transformation method and verification benchmarking procedure

The numerical calculation of electromagnetic far fields in the time-domain requires a near to far field transformation (NTFF) method. While time-domain NTFF methods for popular finite-difference time-domain (FDTD) approaches are well established, there is little discourse on NTFF methods for finite-element time-domain (FETD) codes. Here, this work is concerned with the development of an NTFF method for the Empire FETD code, which utilizes curl and divergence conforming elements. This discretization presents a difficulty in obtaining the equivalent electric current for the NTFF. Straightforward finite element interpolation of the fields is shown to give poor accuracy. Alternative interpolation methods are recommended. An expanding magnetic quadrupole pulse benchmark problem, which is fully developed in the appendices, provides the basis for quantitative comparison.

FETD↗

Photonic ring resonator notch filters for astronomical OH line suppression

Photonic ring resonator arrays used as notch filters are a promising novel solution to improve the signal-to-noise ratio of ground-based astronomical observations by suppressing OH emission lines in the near-infrared (NIR) wavelength range (0.9-2.5 μm). We aim to fabricate a series of ring resonators connected by a waveguide, each with its resonance wavelength and full-width-half-maximum (FWHM) matched with one of the OH emission lines. The fundamental structure of ring resonator-based notch filters is a ring resonator evanescently coupled to a waveguide. Due to the optical loss of the ring resonator, the transmission of light through the waveguide will be suppressed at the resonance wavelength of the ring resonator. The degree of suppression can be controlled by varying the coupling strength between the waveguide and the ring resonator, while the resonance wavelength can be controlled by the geometry and size of the ring. However, since a single ring resonator has multiple modes each with a different wavelength separated by its free spectral range (FSR), the additional suppression at other wavelengths may undermine the signal. Therefore, the ideal FSR off the ring resonators needs to be larger than the measurement window. For this application, we aim the FSR to be 1 larger than 30 nm . Our initial measurement results revealed two major challenges for this application. First, the FSR of the ring resonators is limited by the bending loss of ring resonators. In a standard 220 nm silicon-on-insulator material platform, the FSR of ring resonators cannot exceed 30 nm; Second, due to the wafer thickness variation and fabrication inaccuracy, the achievable accuracy in the resonance wavelength of a fabricated single ring resonator is only ~1.5 nm. In this work, we report our novel solutions to the above two challenges. For the first one, we have developed a novel structure based on the Vernier effect to greatly raise the FSR of ring resonators. Serially coupled Vernier ring structures is a typical approach to increase the FSR without encountering bending loss, but it requires restrictive matching of the coupling and optical loss of the rings and any imperfect matching can result in even noisier signals. On the contrary, our structure, as shown in shown in Fig. 1, does not require the restrictive matching. The signals of two output ports in our structure marked in Fig. 1 will be added together to form the final output signal, and Fig. 2 shows the finite-difference time-domain (FDTD) simulation results of the final output signal. As can be seen, the FSR of the entire structure is well-beyond 60 nm with little noise at other wavelengths. For the second one, we have developed an easy-to-implement method to permanently tune the resonance wavelengths of ring resonators within 0.1 nm accuracy. The tuning mechanism is based on the property of hydrogen silsesquioxane (HSQ) resist. As a negative tone electron beam resist, HSQ upon exposure will cross-link and eventually form SiO2. The refractive 2 index of HSQ resist and formed SiO2 are slightly different , enabling the tuning. By having a thin layer of HSQ resist on top of ring resonator structures followed by SiO2 cladding, the resonance wavelengths of ring resonators can be post-tuned Advances in Optical Astronomical Instrumentation 2019, edited by Simon Ellis, Céline d'Orgeville, Downloaded From: https://www.spiedigitallibrary.org/conference-proceedings-of-spie on 23 Sep 2020 Terms of Use: https://www.spiedigitallibrary.org/terms-of-use through electron beam exposure. We have found that, within 5 nm tuning range, the resonance wavelength shift and dosage applied obeys a linear relationship. So a simple calibration beforehand can result in resonance wavelength tuning within 0.1 nm accuracy. The transmission spectra of a 3-ring device before and after tuning is shown in Fig. 3.

Astronomical optics↗

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↗

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↗

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.↗

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)↗

Parametric study of optical transmission through plasmonic hole arrays modulated by the phase transition of vanadium dioxide

We have performed comprehensive electromagnetic simulations and preliminary experiments to explore the effects of geometrical and material parameters on the extraordinary optical transmission (EOT) through periodic arrays of subwavelength holes in a bilayer stack consisting of a gold or silver film atop a vanadium dioxide film (Au/Ag + VO 2 ), where the latter undergoes a semiconductor-to-metal phase transition. Using the finite-difference time-domain (FDTD) and finite-element methods (FEM), we vary iteratively the array periodicity, VO 2 film thickness and hole diameters, as well as the refractive index inside the VO 2 -layer holes and the VO 2 optical constants. For each variation, we compare the metallic-to-semiconducting ratios of the zero-order transmission ( T 00 ) peaks and find sharp maxima in these ratios within narrow parameter ranges. The maxima arise from Fabry-Perot and Fano-type resonances that minimize T 00 in the semiconducting phase of the perforated bilayers. At a fixed array period, the primary factors controlling the VO 2 -enabled EOT modulation are the VO 2 thickness, diameter of the VO 2 -layer holes, and absorption in the two VO 2 phases. Besides uncovering the origins of the higher metallic-phase T 00 , this study provides a protocol for optimizing the performance of the bilayer hole arrays for potential uses as dynamically tunable nano-optical devices.

Donev, Eugenii U. (ORCID:000000030555219X)↗

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↗

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↗

Plasmonic Core–Shell–Satellites with Abundant Electromagnetic Hotspots for Highly Sensitive and Reproducible SERS Detection

In this work, we develop a Ag@Al2O3@Ag plasmonic core–shell–satellite (PCSS) to achieve highly sensitive and reproducible surface-enhanced Raman spectroscopy (SERS) detection of probe molecules. To fabricate PCSS nanostructures, we employ a simple hierarchical dewetting process of Ag films coupled with an atomic layer deposition (ALD) method for the Al2O3 shell. Compared to bare Ag nanoparticles, several advantages of fabricating PCSS nanostructures are discovered, including high surface roughness, high density of nanogaps between Ag core and Ag satellites, and nanogaps between adjacent Ag satellites. Finite-difference time-domain (FDTD) simulations of the PCSS nanostructure confirm an enhancement in the electromagnetic field intensity (hotspots) in the nanogap between the Ag core and the satellite generated by the Al2O3 shell, due to the strong core–satellite plasmonic coupling. The as-prepared PCSS-based SERS substrate demonstrates an enhancement factor (EF) of 1.7 × 107 and relative standard deviation (RSD) of ~7%, endowing our SERS platform with highly sensitive and reproducible detection of R6G molecules. We think that this method provides a simple approach for the fabrication of PCSS by a solid-state technique and a basis for developing a highly SERS-active substrate for practical applications.

47 OTHER INSTRUMENTATION↗