Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time-domain simulation”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Bayesian inference for plasmonic nanometrology

Here, we introduce a Bayesian method for the characterization of plasmonic nanoparticles, which is applicable to both near- and far-field problems. Designed to combine data generated from any photon-plasmon interaction experiment with physically motivated theoretical models, our approach leverages state-of-the-art Markov chain Monte Carlo sampling techniques and returns parameter estimates on nanometric scales. Simulated spectral data sets, describing resonant scattering of photons from ellipsoidal and toroidal nanoparticles, are explored as concrete examples of our approach, with the resulting Bayesian estimates showing excellent agreement with the ground truth, even under conditions of high statistical noise. By incorporating Bayes factors into the method as well, we reveal how model selection can determine which one of competing geometric shapes better explains the observed data. Our comprehensive nanometrology procedure can be tailored to a variety of light-particle interaction models, and its reliance on Bayesian inference furnishes automatic uncertainty quantification. In addition to applicability to a host of plasmonic configurations such as nanoparticle dimers, trimers, and array studies, it is proposed that the presented analysis can be extended to the quantum regime, where nonclassical photon statistics may provide additional insight for inference of scatterer properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2025

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

Tuning charge carrier dynamics through spacer cation functionalization in layered halide perovskites: an ab initio quantum dynamics study

Dion–Jacobson (DJ) phase two-dimensional (2D) hybrid halide perovskites are promising for environmentally stable optoelectronic device applications due to their attractive photophysical properties, including long charge carrier lifetime and partially suppressed non-radiative losses. However, the atomistic details of the structure–property relationship are severely limited, which substantially restricts the strategic materials designed for these layered halide perovskites (LHPs). Here, we combine nonadiabatic molecular dynamics and time-domain density functional theory to understand the effects of spacer cation functionalization on the ground and excited-state charge carrier dynamics in DJ phase LHPs. Our in-depth study reveals that the fluorination of spacer cations considerably restricts the thermal motions of LHPs at ambient conditions. The compact structure weakens the overall electron–phonon interactions and reduces the thermal fluctuations of the band edges in real-time. These dynamic modifications partially mitigate the non-radiative recombination, prolonging the lifetime of photogenerated charge carriers. The suppressed carrier loss mechanism strongly suggests that the fluorinated spacer cation-based LHPs would exhibit enhanced performance as optoelectronic materials. These systematic simulations of excited state carrier dynamics elaborate that chemically viable functionalization of the spacer cations is a robust approach to improve the photophysical properties of 2D halide perovskites strategically. These insights also guide us to propose a few potential strategies to design highly beneficial spacer cations of LHPs that can be introduced in next-generation optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Elucidating Norrish type I reactive pathways by ultrafast X-ray absorption spectroscopy

Norrish type I reactions selectively cleave carbon–carbon bonds directly adjacent to carbonyl groups. Despite their broad use in combination with aromatic carbonyls for additive manufacturing and dental UV curing applications, the nature of the photochemically active state and its population mechanism remain insufficiently understood. Detailed mechanistic insight requires mapping of the photoexcited population flow involving internal conversion and intersystem crossing. We present a time-domain study of gas phase acetophenone as a prototypical aromatic carbonyl combining soft X-ray time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS) spectroscopy at the oxygen K-edge with ab initio multiple spawning (AIMS) simulations. Exploiting the specific sensitivity of TR-NEXAFS spectroscopy to states with nπ* character, we observe population transfer from the initially excited 1 ππ* state to the 1 nπ* state with a time constant of (0.13 ± 0.02) ps after an initial induction period of (0.12 ± 0.02) ps without population transfer, in quantitative agreement with the AIMS simulations. The population in the 1 nπ* state subsequently decays via intersystem crossing, likely mediated by a 3 ππ* state, within (3.17 ± 0.66) ps to a long-lived 3 nπ* state, which is presumed to be active towards Norrish type I chemistry.

Graßl, Martin [SLAC National Accelerator Laborator↗
DOE OSTI2024

ATAT: Astronomical Transformer for time series and Tabular data

Context. The advent of next-generation survey instruments, such as theVera C. RubinObservatory and its Legacy Survey of Space and Time (LSST), is opening a window for new research in time-domain astronomy. The Extended LSST Astronomical Time-Series Classification Challenge (ELAsTiCC) was created to test the capacity of brokers to deal with a simulated LSST stream. Aims. Our aim is to develop a next-generation model for the classification of variable astronomical objects. We describe ATAT, the Astronomical Transformer for time series And Tabular data, a classification model conceived by the ALeRCE alert broker to classify light curves from next-generation alert streams. ATAT was tested in production during the first round of the ELAsTiCC campaigns. Methods. ATAT consists of two transformer models that encode light curves and features using novel time modulation and quantile feature tokenizer mechanisms, respectively. ATAT was trained on different combinations of light curves, metadata, and features calculated over the light curves. We compare ATAT against the current ALeRCE classifier, a balanced hierarchical random forest (BHRF) trained on human-engineered features derived from light curves and metadata. Results. When trained on light curves and metadata, ATAT achieves a macro F1 score of 82.9 ± 0.4 in 20 classes, outperforming the BHRF model trained on 429 features, which achieves a macro F1 score of 79.4 ± 0.1. Conclusions. The use of transformer multimodal architectures, combining light curves and tabular data, opens new possibilities for classifying alerts from a new generation of large etendue telescopes, such as theVera C. RubinObservatory, in real-world brokering scenarios.

Astronomy & Astrophysics↗
DOE OSTI2021

Spatially resolved x-ray detection with photonic crystal scintillators

We study the self-collimation phenomenon in photonic crystals (PhC) of wide bandgap materials for ultra-fast and high spatial resolution x-ray detection. We work on various heavy inorganic scintillators: BaF 2 , GaN, ZnO, CsI:Tl, NaI:Tl, LYSO, WO 4 compounds, and plastic scintillators. Conventional scintillator detectors do not rely on a direct detection mechanism; hence, they require intricate design and fabrication processes. We offer a PhC design to observe self-collimation phenomena and overcome the ongoing spatial resolution challenges with these types of materials. We investigate the photonic band diagrams and iso-frequency contours. Fourier transforms based on finite-difference time-domain and frequency domain simulations are done for verifying and analyzing the self-collimation with the selected material. Light extraction efficiency at the PhC–air interface, depending on the truncation distance from the excitation point, is measured. Beam divergence values are calculated at 1 mm propagation distance. The vertical field profiles are obtained to observe the confinement. For the spatial resolution analysis, cross-sectional beam profiles have been examined. Furthermore, Gaussian envelopes are fitted to beam profiles for a consistent data analysis, and full-width-at-half-maximum values are considered. As a result, we theoretically prove and demonstrate the spatially resolved x-ray detection at the sub-micrometer level for a wide range of scintillator materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Manipulating ultrafast even-order nonlinear chiral responses of L-tryptophan by polarization pulse shaping

Molecular chirality has long been monitored in the frequency domain in the ultraviolet, visible, and infrared regimes. Recently developed time-domain approaches can detect time-dependent chiral dynamics by enhancing intrinsically weak chiral signals. Even-order nonlinear signals in chiral molecules have gained attention thanks to their existence in the electric dipole approximation, without relying on the weaker higher-order multipole interactions. We illustrate the optimization of temporal polarization pulse-shaping in various frequency ranges (infrared/optical and optical/X ray) to enhance chiral nonlinear signals. These signals can be recast as an overlap integral of matter and field pseudoscalars which contain the relevant chiral information. Simulations are carried out for second- and fourth-order nonlinear spectroscopies in L-tryptophan.

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

Constraining Type Ia supernova explosions and early flux excesses with the Zwicky Transient Factory

In the new era of time-domain surveys, Type Ia supernovae are being caught sooner after explosion, which has exposed significant variation in their early light curves. Two driving factors for early-time evolution are the distribution of 56 Ni in the ejecta and the presence of flux excesses of various causes. We perform an analysis of the largest young SN Ia sample to date. We compare 115 SN Ia light curves from the Zwicky Transient Facility to the turtls model grid containing light curves of Chandrasekhar mass explosions with a range of 56 Ni masses, 56 Ni distributions, and explosion energies. We find that the majority of our observed light curves are well reproduced by Chandrasekhar mass explosion models with a preference for highly extended 56 Ni distributions. We identify six SNe Ia with an early-time flux excess in our gr-band data (four ‘blue’ and two ‘red’ flux excesses). We find an intrinsic rate of 18 ± 11 percent of early flux excesses in SNe Ia at z < 0.07, based on three detected flux excesses out of 30 (10 percent) observed SNe Ia with a simulated efficiency of 57 percent. This is comparable to rates of flux excesses in the literature but also accounts for detection efficiencies. Two of these events are mostly consistent with circumstellar material interaction, while the other four have longer lifetimes in agreement with companion interaction and 56 Ni-clump models. We find a higher frequency of flux excesses in 91T/99aa-like events (44 ± 13 percent).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Detecting strongly lensed type Ia supernovae with LSST

ABSTRACT Strongly lensed supernovae are rare and valuable probes of cosmology and astrophysics. Upcoming wide-field time-domain surveys, such as the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST), are expected to discover an order-of-magnitude more lensed supernovae than have previously been observed. In this work, we investigate the cosmological prospects of lensed type Ia supernovae (SNIa) in LSST by quantifying the expected annual number of detections, the impact of stellar microlensing, follow-up feasibility, and how to best separate lensed and unlensed SNIa. We simulate SNIa lensed by galaxies, using the current LSST baseline v3.0 cadence, and find an expected number of 44 lensed SNIa detections per year. Microlensing effects by stars in the lensing galaxy are predicted to lower the lensed SNIa detections by ∼8 per cent. The lensed events can be separated from the unlensed ones by jointly considering their colours and peak magnitudes. We define a ‘gold sample’ of ∼10 lensed SNIa per year with time delay >10 d, >5 detections before light curve peak, and sufficiently bright (mi < 22.5 mag) for follow-up observations. In 3 yr of LSST operations, such a sample is expected to yield a 1.5 per cent measurement of the Hubble constant.

Astronomy & Astrophysics↗
DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Characterizing particle-like charge-migration dynamics with high-order harmonic sideband spectroscopy

Here, we introduce high-order harmonic sideband spectroscopy and show that it can be a robust probe of charge migration (CM) in a halogenated carbon-chain molecule. We simulate both the CM and high-order harmonic generation (HHG) dynamics using ab initio time-dependent density-functional theory. We find that CM dynamics initiated along the molecular backbone induces sidebands in the HHG spectrum driven by a delayed laser pulse that is polarized perpendicular to the molecular axis. Monitoring the spectrum as either the HHG laser frequency or the relative delay is scanned allows for the extraction of detailed information about the time-domain characteristics of the CM process.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Pseudospectral particle-in-cell formulation with arbitrary charge and current-density time dependencies for the modeling of relativistic plasmas

This paper introduces a formulation of the particle-in-cell (PIC) method for the modeling of relativistic plasmas, that leverages the ability of the pseudospectral analytical time-domain solver (PSATD) to handle arbitrary time dependencies of the charge and current densities during one PIC cycle (applied to second-order polynomial dependencies here). Here, the formulation is applied to a modified set of Maxwell's equations that was proposed earlier in the context of divergence cleaning, and to recently proposed extensions of the PSATD-PIC algorithm. Detailed analysis and testings revealed that, under some condition, the formulation can expand the range of numerical parameters under which PIC simulations are stable and accurate when modeling relativistic plasmas such as, e.g., plasma-based particle accelerators.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

DQ Admittance Extraction for Inverter-Based Resources

The power grid industry is pushing for electromagnetic transient (EMT)-based studies for generation interconnection and planning process due to high penetrations of inverter-based resource (IBRs). Vendor-specific and site-specific black-box IBR models are preferred in those simulation studies. For small-signal analysis, measurement-based admittance models are necessary. In this paper, we demonstrate the extraction of frequency-domain dq frame IBR models. These linear models are specific to operating conditions. We demonstrate two extraction methods: frequency scan and a step response-based method. The latter relies on converting time-domain responses to Laplace-domain expressions via eigensystem realization analysis (ERA). Both can lead to dq admittance representation and the latter is time saving.

admittance model↗
DOE OSTI2024

A Time-Domain Protection Approach for AC Transmission Systems With Grid-Forming Resources

Ac transmission protection must reliably detect, classify, and locate short-circuit faults from voltage and current measurements. At present, these functionalities, which have been classically engineered using phasors approaches, are being challenged by the dynamic behavior and fault-current limits of converter-based generation. This paper tackles these challenges by engineering a time-domain protection approach that leverages the classical Bergeron model in a new manner. Low- and high-impedance faults are detected and classified by ascertaining how well line voltage and current measurements match the Bergeron equations. Faults are located by posing a novel one-variable optimization problem, whereas voltage and current waveforms at the fault location are estimated by unveiling rigorous relationships. The proposed elements are secure against external faults, measurement errors, and variation of line parameters and sampling time. Furthermore, these advances are tested via electromagnetic transient simulations and are significant to satisfy IEEE and North American Electric Reliability Corporation requirements.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Data-Driven Cyber-Attack Detection for PV Farms via Time-Frequency Domain Features

The internetworking of grid-connected power electronics converters (PECs) in photovoltaic (PV) farms has inevitably expanded the cyber-attack surfaces. Here this paper presents a comprehensive study on cyber-attack detection and diagnosis for PEC-enabled PV farms via single waveform sensor to distinguish between normal conditions, open-circuit faults, short-circuit faults, and cyber-attacks. To our knowledge, this has not been attempted before. Firstly, we propose frequency-domain magnitude-based residuals to identify short-circuit faults and a time-domain mean current vector-based feature to distinguish open-circuit faults from other threats. These features can fully reflect the specific physical characteristics of PV farms during threat duration. Secondly, unlike micro phasor measurement units (µPMU) and raw electric waveform-based methods, the proposed innovative features can address novel cyber-attacks that are excluded from the training process. Thirdly, an online hardware-in-the-loop (HIL) testbed using the OPAL-RT real-time digital simulator has verified the effectiveness. The monitoring system runs in real-time while using HIL as an operational solar farm and a National Instruments (NI) data acquisition card as the electric waveform sensor at the point of coupling.

42 ENGINEERING↗
DOE OSTI2026

Enabling Grid-Forming Control with Fault Ride-Through in Unbalanced Distribution Networks

Distribution networks are often unbalanced, causing oscillatory responses in inverter control designed for balanced conditions. Here, to address this problem, this paper proposes a novel time-domain transformation appropriate for inverter control and enables the decomposition of three-phase unbalanced signals into constant positive and negative components. Relations useful for calculating unbalanced active and reactive power are derived from first principle, providing insight into vector products of unbalanced three-phase signals. Furthermore, a grid-forming control effective under unbalanced conditions is developed, which delivers superior performance while meeting UNIFI1 specifications for grid-forming control under unbalanced conditions. specifications applicable to category 4 inverter-based resource, like setting and regulating frequency/voltage, providing voltage support, sharing active power, injecting negative sequence current, and riding through faults. A current limiter is proposed for safe fault ride-through and integrates with the grid-forming control featuring frequency/voltage droop controllers and current and voltage control loops. The transformation of interconnected inverters is formulated and stability of the proposed control analyzed to support robust parameter selections. The effectiveness of the proposed transformation and grid-forming control is demonstrated through analytical results and real-time simulation of a IEEE 123 distribution network on the Real-Time Digital Simulator. Comparison with existing methods shows that the proposed strategy satisfies the UNIFI specifications with a much better performance.

24 - POWER TRANSMISSION AND DISTRIBUTION↗