Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
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 OSTI2020

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

Direct E&M Field Solver for a Plane Wave without Sources

A three dimensional direct field solver for the electric and magnetic fields using the Finite Difference Time Domain (FDTD) method is described. Example simulations using periodic and absorbing boundary conditions in a vacuum medium for a plane wave propagating in three dimensions are explored and are compared to theory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Numerical and Analytical Modeling of the Effect of Cracks on the Self-Inductance of a COTS YJ-41003-TC Toroid

COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μ r is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μ core >> μ air ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.

36 MATERIALS SCIENCE↗
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 OSTI2023

Plasmonic hybrid core-shell (HyCoS) AgPt NP template hybridized with GQDs for SERS enhancement of 4-MBA and BT

Surface-enhanced Raman spectroscopy (SERS) is an attractive vibrational spectroscopic technique that can enable a non-destructive and ultra-sensitive detection down to the single-molecule level. Herein, a novel hybrid SERS platform is developed based on hybrid core-shell (HyCoS) AgPt nanoparticles (NPs) and graphene quantum dots (GQDs) for the enhancement of Raman vibration of 4-mercaptobenzoic acid (4-MBA) and benzenethiol (BT). The unique design of HyCoS AgPt NPs induces strong electromagnetic mechanism (EM) enhancement through the amplification of electromagnetic fields by the excitation of high-density surface plasmons and hot spots. Superior localized surface plasmon resonance (LSPR) is generated by the AgPt core-shell and background Ag NP coupling, which is systematically investigated by the optical properties and FDTD simulations. The background Ag NPs can further increase the coverage of metallic NPs, leading to higher-density hot spots and enhanced SERS response. At the same time, GQDs can provide plentiful accessible edges for the charge transfer to the HOMO and LUMO of 4-MBA and BT based on the chemical mechanism (CM) enhancement. The mixing approach of GQDs and target molecules on the HyCoS AgPt NPs can significantly amplify the Raman signals via the strong adsorption of probe molecules by the π–π interaction. The enhancement factors of proposed SERS platform can reach ~10 7 and ~10 5 for the 4-MBA and BT respectively.

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

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2022

Photonic Luminescent Solar Concentrator Design for High Efficiency, Low Cost Multijunction Photovoltaics

Despite the extraordinary advances in solar cell efficiency in laboratory settings, the deployment of solar cells continues to be limited to low efficiency (<25%) silicon cells because of cost. In this work, we take advantage of the extraordinary optical properties afforded by nanophotonic structures to create a photonic luminescent solar concentrator for an InGaP-Si multijunction concentrator cell. Finite difference time domain (FDTD) simulations demonstrated a concentrator that could effectively capture, downconvert, and guide concentrated light to an InGaP subcell while still transmitting longer wavelengths to a Si subcell. We fabricated the photonic luminescent solar concentrator, which was comprised of CdSe/CdS quantum dots embedded within alternating layers of Si 3 N 4 and SiO 2 , and experimentally verified the optical performance, showing a 40% increase in light guiding and a significant reduction in reabsorption losses in the plane of the luminescent concentrator as compared to traditional designs. Finally, we utilized modified detailed balance calculations that accounted for cell and optical losses and showed >30% efficiencies are possible with this design, demonstrating the potential to meet the demands for high efficiency, inexpensive solar modules.

Eisler, Carissa N.↗
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↗
NASA NTRS1991

A real scale simulator for high frequency LEMP

The real scale simulator is described which was designed by the Centre d'Etudes de Gramat (CEG) to study the coupling of fast rise time Lightning Electromagnetic pulse in a fighter aircraft. The system capability of generating the right electromagnetic environment was studied using a Finite Difference Time Domain (FDTD) computer program. First, data of inside stresses are shown. Then, a time domain and a frequency domain approach is exposed and compared.

Gauthier, D.↗
NASA NTRS1991

A theoretical analysis of the electromagnetic environment of the AS330 super Puma helicopter external and internal coupling

Numerical techniques such as Finite Difference Time Domain (FDTD) computer programs, which were first developed to analyze the external electromagnetic environment of an aircraft during a wave illumination, a lightning event, or any kind of current injection, are now very powerful investigative tools. The program called GORFF-VE, was extended to compute the inner electromagnetic fields that are generated by the penetration of the outer fields through large apertures made in the all metallic body. Then, the internal fields can drive the electrical response of a cable network. The coupling between the inside and the outside of the helicopter is implemented using Huygen's principle. Moreover, the spectacular increase of computer resources, as calculations speed and memory capacity, allows the modellization structures as complex as these of helicopters with accuracy. This numerical model was exploited, first, to analyze the electromagnetic environment of an in-flight helicopter for several injection configurations, and second, to design a coaxial return path to simulate the lightning aircraft interaction with a strong current injection. The E field and current mappings are the result of these calculations.

Flourens, F.↗
NASA NTRS1991

A Fourier collocation time domain method for numerically solving Maxwell's equations

A new method for solving Maxwell's equations in the time domain for arbitrary values of permittivity, conductivity, and permeability is presented. Spatial derivatives are found by a Fourier transform method and time integration is performed using a second order, semi-implicit procedure. Electric and magnetic fields are collocated on the same grid points, rather than on interleaved points, as in the Finite Difference Time Domain (FDTD) method. Numerical results are presented for the propagation of a 2-D Transverse Electromagnetic (TEM) mode out of a parallel plate waveguide and into a dielectric and conducting medium.

Shebalin, John V.↗
NASA NTRS1991

A two-dimensional time domain near zone to far zone transformation

In a previous paper, a time domain transformation useful for extrapolating 3-D near zone finite difference time domain (FDTD) results to the far zone was presented. In this paper, the corresponding 2-D transform is outlined. While the 3-D transformation produced a physically observable far zone time domain field, this is not convenient to do directly in 2-D, since a convolution would be required. However, a representative 2-D far zone time domain result can be obtained directly. This result can then be transformed to the frequency domain using a Fast Fourier Transform, corrected with a simple multiplicative factor, and used, for example, to calculate the complex wideband scattering width of a target. If an actual time domain far zone result is required it can be obtained by inverse Fourier transform of the final frequency domain result.

Luebbers, Raymond J.↗
NASA NTRS1991

Wideband finite difference time domain implementation of surface impedance boundary conditions for good conductors

Surface impedance boundary conditions are employed to reduce the solution volume during the analysis of scattering from lossy dielectric objects. In a finite difference solution, they also can be utilized to avoid using small cells, made necessary by shorter wavelengths in conducting media, throughout the solution volume. A 1-D implementation for a surface impedance boundary condition for good conductors in the Finite Difference Time Domain (FDTD) technique.

Beggs, John H.↗
NASA NTRS1991

Time domain scattering and radar cross section calculations for a thin, coated perfectly conducting plate

Radar cross section (RCS) calculations for flat, perfectly conducting plates are readily available through the use of conventional frequency domain techniques such as the Method of Moments (MOM). However, if the plate is covered with a dielectric material that is relatively thick in comparison with the wavelength in the material, these frequency domain techniques become increasingly difficult to apply. We present the application of the Finite Difference Time Domain (FDTD) Technique to the problem of electromagnetic scattering and RCS calculations from a thin, perfectly conducting plate that is coated with a thick layer of lossless dielectric material. Both time domain and RCS calculations are presented and disclosed.

Luebbers, Raymond J.↗
NASA NTRS1991

A two-dimensional time domain near zone to far zone transformation

A time domain transformation useful for extrapolating three dimensional near zone finite difference time domain (FDTD) results to the far zone was presented. Here, the corresponding two dimensional transform is outlined. While the three dimensional transformation produced a physically observable far zone time domain field, this is not convenient to do directly in two dimensions, since a convolution would be required. However, a representative two dimensional far zone time domain result can be obtained directly. This result can then be transformed to the frequency domain using a Fast Fourier Transform, corrected with a simple multiplicative factor, and used, for example, to calculate the complex wideband scattering width of a target. If an actual time domain far zone result is required, it can be obtained by inverse Fourier transform of the final frequency domain result.

Luebbers, Raymond J.↗