Search NASA⌕ Search

SEARCH · Search NASA

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
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 OSTI2021

Enhanced Planar Antenna Efficiency Through Magnetic Thin-Films

This work proposes to use magnetic material as the substrate of planar antennas to overcome the platform effect caused by the conducting ground plane. The upper bound of the radiation efficiency of an electric-current-driven low-profile antenna is theoretically derived, which is inversely proportional to the Gilbert damping factor of the magnetic material. Meanwhile, the improvement of radiation due to the use of magnetic material is demonstrated by a three-dimensional (3D) multiphysics and multiscale time-domain model. The simulation results match the theoretical derivation, showing 25% radiation efficiency from a planar antenna backed by a FeGaB thin film with 2.56 m thickness. Furthermore, for conductive ferromagnetic materials, it is shown that the eddy current loss can be well suppressed by laminating the thin film into multiple layers. Here, the radiation efficiency of the modeled antenna with a conductive ferromagnetic substrate is improved from 2.2% to 11.8% by dividing the substrate into 10 layers, with a ferromagnetic material fill factor of 93%.

42 ENGINEERING↗
DOE OSTI2024

Control Algorithms for Multi-Wavefront Sensor Adaptive Optics for Astronomical Exoplanet Imaging

This package consists of simulations of control algorithms being proposed for upcoming adaptive optics projects for exoplanet direct imaging, in which one wavefront affector (deformable mirror) is controlled by two wavefront sensors. Simulations include time-domain behavior under various assumed conditions, parameter optimization routines for the controllers, and stability analysis.

Sengupta, Aditya↗
DOE OSTI2024

GENERATING CORCOS-COHERENT SIGNAL SERIES USING PHASE PERTURBATION​

Phase perturbations are used to simulate the time-domain behavior of pressure loads from turbulent boundary layers while faithfully representing the coherence decay behaviors of turbulence-generated pressure fields. Turbulence "evolves" as it proceeds downstream, and coherence decays as the decaying exponential of the product of a flow-related constant with frequency and distance. Time domain methods for a line, a plane, a closed cyclic surface, and most of a cone were given.

Swift, Stephen Hales↗
DOE OSTI2025

Review of recent WEC-Sim (v6.1) advanced features

WEC-Sim (Wave Energy Converter SIMulator) is an opensource software for modeling the motions, loads and power generation of wave energy converters. WEC-Sim performs simulations in the time domain using hydrodynamic coefficients calculated by boundary element method (BEM) frequency-domain potential flow solvers such as WAMIT, NEMOH, Capytaine, or Ansys AQWA. WEC-Sim development is ongoing, including various features, applications, and example cases used to demonstrate potential use-cases to meet the needs of the growing marine energy industry. Through input from a broad user base and an extensive team of developers and collaborators, new features of WEC-Sim are developed to expand the software’s use cases and improve overall functionality. Three new WEC-Sim features highlighted in this paper include updating WEC-Sim to be compatible with MoorDyn Version 2, incorporation of second order excitation loads (quadratic transfer functions) and allowing for dynamically changing hydrodynamics.

13 HYDRO ENERGY↗
DOE OSTI2020

Influence of Time and Frequency Domain Wave Forcing on the Power Estimation of a Wave Energy Converter Array

Industry-specific tools for analyzing and optimizing the design of wave energy converters (WECs) and associated power systems are essential to advancing marine renewable energy. This study aims to quantify the influence of phase information on the device power output of a virtual WEC array. We run the phase-resolving wave model FUNWAVE-TVD (Total Variation Diminishing) to generate directional waves at the PacWave South site offshore from Newport, Oregon, where future WECs are expected to be installed for testing. The two broad cases presented correspond to mean wave climates during warm months (March–August) and cold months (September–February). FUNWAVE-TVD time series of sea-surface elevation are then used in WEC-Sim, a time domain numerical model, to simulate the hydrodynamic response of each device in the array and estimate their power output. For comparison, WEC-Sim is also run with wave energy spectra calculated from the FUNWAVE-TVD simulations, which do not retain phase information, and with wave spectra computed using the phase-averaged model Simulating WAves Nearshore (SWAN). The use of spectral data in WEC-Sim requires a conversion from frequency to time domain by means of random superposition of wave components, which are not necessarily consistent because of the linear assumption implicit in this method. Thus, power response is characterized by multiple realizations of the wave climates.

13 HYDRO ENERGY↗
DOE OSTI2020

Planck intermediate results

We present the NPIPE processing pipeline, which produces calibrated frequency maps in temperature and polarization from data from the Planck Low Frequency Instrument (LFI) and High Frequency Instrument (HFI) using high-performance computers. NPIPE represents a natural evolution of previous Planck analysis efforts, and combines some of the most powerful features of the separate LFI and HFI analysis pipelines. For example, following the LFI 2018 processing procedure, NPIPE uses foreground polarization priors during the calibration stage in order to break scanning-induced degeneracies. Similarly, NPIPE employs the HFI 2018 time-domain processing methodology to correct for bandpass mismatch at all frequencies. In addition, NPIPE introduces several improvements, including, but not limited to: inclusion of the 8% of data collected during repointing manoeuvres; smoothing of the LFI reference load data streams; in-flight estimation of detector polarization parameters; and construction of maximally independent detector-set split maps. For component-separation purposes, important improvements include: maps that retain the CMB Solar dipole, allowing for high-precision relative calibration in higher-level analyses; well-defined single-detector maps, allowing for robust CO extraction; and HFI temperature maps between 217 and 857 GHz that are binned into 0'.9 pixels (Nside = 4096), ensuring that the full angular information in the data is represented in the maps even at the highest Planck resolutions. The net effect of these improvements is lower levels of noise and systematics in both frequency and component maps at essentially all angular scales, as well as notably improved internal consistency between the various frequency channels. Based on the NPIPE maps, we present the first estimate of the Solar dipole determined through component separation across all nine Planck frequencies. The amplitude is (3366.6 ± 2.7) μK, consistent with, albeit slightly higher than, earlier estimates. From the large-scale polarization data, we derive an updated estimate of the optical depth of reionization of τ = 0.051 ± 0.006, which appears robust with respect to data and sky cuts. There are 600 complete signal, noise and systematics simulations of the full-frequency and detector-set maps. As a Planck first, these simulations include full time-domain processing of the beam-convolved CMB anisotropies. The release of NPIPE maps and simulations is accompanied with a complete suite of raw and processed time-ordered data and the software, scripts, auxiliary data, and parameter files needed to improve further on the analysis and to run matching simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

MOSCATO Development and Integration in Fiscal Year 2023

In FY23, significant developments were made to the MOSCATO code package, enhancing its capabilities in many aspects. The improvements and advancements can be summarized as follows: Code Restructuring: The MOSCATO code underwent a restructuring process, enabling both Nek5000 and NekRS to share the same interface when using MOSCATO. This integration facilitates smoother interactions between the different modules and streamlines the overall simulation process. Automatic Domain Extension: To allow simulations to propagate over longer time frames, an automatic domain extension capability was added to the metal diffusion module. This enhancement enables more accurate and reliable simulations of diffusion processes, especially in scenarios with extended time requirements. Expanded Corrosion Reactions: The updated MOSCATO package now includes additional reactions that contribute to corrosion. This addition enhances the code’s ability to model corrosion effects more comprehensively, providing a more realistic representation of the systems being studied. Example Case Setup: To assist users and showcase MOSCATO’s capabilities, we developed example cases for various applications. Notably, we created example cases for thermal cells and corrosion loops. These pre-configured cases allow users to get started quickly and effectively with their simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Parameterization of generic positive sequence models to represent behavior of inverter based resources in low short circuit scenarios

The connection of large inverter based resources (IBR) in transmission systems is often located geographically and electrically far away from load centers. This, coupled with the displacement of synchronous machine plants, results in a reduction of the network short circuit strength at the point of connection. Under these conditions, state-of-the-art positive sequence simulation platforms and models can have difficulties maintaining numerical stability and/or providing an accurate representation of IBR plant dynamic behavior. As a result, computationally heavy time domain electromagnetic transient (EMT) simulations may be required to evaluate these systems. In this paper, a recently developed improved generic positive sequence model is parameterized to provide a representation of IBR behavior in low short circuit networks. Further, comparisons against generic and detailed EMT models demonstrate the suitability of the improved positive sequence model to study practical stability issues experienced with presently in-service plants. Such a model can provide some of the accuracy of an EMT representation with a much lower computational burden. The performance of the positive sequence model is validated against the behavior shown by both open white box and closed black box EMT domain models from around the world.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Optimal Droop Setting for Congestion Reduction in a 100% Grid-Forming Inverter-based Power System

he high penetration of inverter-based resources (IBRs) introduces new challenges to power systems due to the complex inverter control. However, IBRs can be configured to maximize their benefits to improve system resilience and reliability. This paper proposes a steady-state optimization model that aims to mitigate transmission congestion in a 100% grid- forming (GFM) IBR-based power system. This goal is achieved by determining the optimal droop settings for the GFM IBRs under different congestion conditions due to renewable energy and load variations. The numerical solution is rigorously verified by a high-fidelity model of the IEEE 39-bus test system with detailed GFM IBR control in the time-domain electromagnetic transient (EMT) simulation tool PSCAD. The numerical solution and simulation results show a significant congestion reduction while meeting all other operating requirements. It is also observed that the numerical solving time is substantially less compared to the EMT simulation time.

Nguyen, Quan H.↗
DOE OSTI2022

Charge carrier nonadiabatic dynamics in non-metal doped graphitic carbon nitride

Graphitic carbon nitride (GCN) has attracted significant attention due to its excellent performance in photocatalytic applications. Non-metal doping of GCN has been widely used to improve the efficiency of the material as a photocatalyst. Using a combination of time-domain density functional theory with nonadiabatic molecular dynamics, we study the charge carrier dynamics in oxygen and boron doped GCN systems. The reported simulations provide a detailed time-domain mechanistic description of the charge separation and recombination processes that are of fundamental importance while evaluating the photovoltaic and photocatalytic performance of the material. The appearance of smaller energy gaps due to the presence of dopant states improves the visible light absorption range of the doped systems. At the same time, the nonradiative lifetimes are shortened in the doped systems as compared to the pristine GCN. In the case of boron doped at a carbon (B–C–GCN), the charge recombination time is very long as compared to the other two doped systems owing to the smaller electron–phonon coupling strength between the valence band maximum and the trap state. The results suggest B–C–GCN as the most suitable candidate among three doped systems studied in this work for applications in photocatalysis. Finally, this work sheds light into the influence of dopants on quantum dynamics processes that govern GCN performance and, thus, guides toward building high-performance devices in photocatalysis.

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

Numerical simulation of NSLS-II fast orbit feedback system

We discuss recent work on the construction of a numerical simulation of the fast orbit feedback system at NSLS-II. The simulation operates in the time domain and also includes the spatial domain, i.e. all beam position monitors (BPMs) and correctors. It can accept inputs of the real beam orbit and the measured orbit response matrix. It can also add errors to every stage of the calculation. We present methods to verify the simulation results by comparing the simulation results with experimental data collected at NSLS-II. The results are in very good agreement. The effects of errors from BPMs and correctors, as well as the amplitude of excitation, on the feedback performance, are also explored. This simulation can be used to predict the behavior and performance of the system for a future upgrade.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Efficient derivative computation for unsteady fatigue-constrained nonlinear aero-structural wind turbine blade optimization

Gradient-based optimization offers significant efficiency advantages for wind turbine blade design, but its application has often been limited by the cost and accuracy of finite-difference derivative calculations, especially when fatigue constraints are considered. In this work, we systematically compare and evaluate four differentiation techniques, namely algorithmic differentiation, implicit differentiation, sparsity exploitation, and parallelization, to determine their effectiveness in computing accurate gradients through time-domain aero-structural simulations. By integrating these techniques with unsteady nonlinear aerodynamic and structural models, we develop software designed for accurate gradient computation. We show that combining these techniques addresses memory and runtime challenges associated with long simulations required by design load cases. Specifically, the most effective combination reduces derivative computation wall time by over an order of magnitude compared to finite differencing while maintaining superior accuracy. We demonstrate this approach in a proof-of-concept aero-structural optimization of a wind turbine blade that improves the cost of energy by 12.78 %. This comparative study establishes a viable approach for fatigue-aware blade design that balances computational efficiency with modeling accuracy.

17 WIND ENERGY↗
DOE OSTI2021

Femtosecond damage experiments and modeling of broadband mid-infrared dielectric diffraction gratings

High peak and average power lasers with high wall-plug efficiency, like the Big Aperture Thulium (BAT) laser, have garnered tremendous attention in laser technology. To meet the requirements of the BAT laser, we have developed low-dispersion reflection multilayer dielectric (MLD) gratings suitable for compression of high-energy pulses for operations at 2 micron wavelength. We carried out 10000-on-1 damage tests to investigate the fluence damage thresholds of the designed MLD gratings and mirrors, which were found between 100-230 mJ/cm 2 . An ultrashort pulsed laser (FWHM = 53 fs, λ = 1.9 μ m) operating at 500 Hz was used in the serpentine raster scans. The atomic force microscope images of the damage sites show blister formation of the underlying layers at lower fluences but ablation of the grating pillars at higher fluences. We simulated the dynamic electronic excitation in the MLD optics with a finite-difference in the time domain approach in 2D. The simulation results agree well with the LIDT measurements and the observed blister formation. This model is able to evaluate the absolute LIDT of MLD gratings.

42 ENGINEERING↗
DOE OSTI2020

Tip-Enhanced Raman Nanospectroscopy of Smooth Spherical Gold Nanoparticles

We record nanoscale-resolved chemical images of thiobenzonitrile (TBN)-functionalized smooth gold nanospheres on silicon via tip-enhanced Raman (TER) nano-spectroscopy. The recorded images trace the nascence of the familiar doughnut-shaped scattering profile of nanoparticles on silicon at its origin (the particle surface), which appears as a horseshoe-shaped scattering pattern under our experimental conditions. The local optical field maps are in agreement with their simulated finite-difference time-domain analogues. Analysis of the recorded spectra with the aid of ab initio molecular dynamics-based Raman spectral simulations further suggests that optical rectification and molecular charging take place throughout the course of atomic force microscopy-based TER nanoscale chemical imaging.

Bhattarai, Ashish↗
DOE OSTI2020

Optimum power analysis of a self-reactive wave energy point absorber with mechanically-driven power take-offs

This paper studies a self-reactive ocean wave energy converter (WEC) where energy is absorbed through the relative motion between a floating buoy on the ocean surface and a submerged body. Two types of direct-drive power take off (PTO) systems are examined. One adopts a ball screw system in inverse driving mode with mechanical motion rectifier (MMR) that converts the oscillating bi-directional input heave motion into unidirectional rotation to drive an electromagnetic generator, and another one uses a ball screw system to directly drive the generator (called non-MMR). Dynamic models for both types of PTOs are established and integrated with the overall WEC system model in both the time and frequency domains. Numerical simulation is used to investigate the dynamic performance of both PTOs. The influences of the PTO inerter, which is from the rotational inertia of the flywheel and generator, and the WEC drag damping coefficient are discussed and explored. The analytical closed-form solution for the optimum condition of the non-MMR system is derived based on the criterion of maximum power absorption, and numerical simulation yields the optimum condition for the MMR system. Overall, the results under regular and irregular waves show that the PTO inerter in the MMR system can improve the power absorption for small wave periods and maintain the same performance for large wave periods. The PTO inerter in the non-MMR system shifts the peak frequency. The WEC drag damping negatively influences the power absorption in both PTO systems, especially at large wave periods.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

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

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

Harthcock, Colin [Lawrence Livermore National Labo↗
DOE OSTI2020

Computational modelling and experimental tank testing of the multi float WaveSub under regular wave forcing

A submerged wave device generates energy from the relative motion of floating bodies. In WaveSub, three floats are joined to a reactor; each connected to a spring and generator. Electricity generated damps the orbital movements of the floats. The forces are non-linear and each float interacts with the others. Tuning to the wave climate is achieved by changing the line lengths, so there is a need to understand the performance trade-offs for a large number of configurations. This requires an efficient, large displacement, multidirectional, multi-body numerical scheme. Results from a 1/25 scale wave basin experiment are described. Here, we show that a time domain linear potential flow formulation (Nemoh, WEC-Sim) can match the tank testing provided that suitably tuned drag coefficients are employed. Inviscid linear potential models can match some wave device experiments; however, additional viscous terms generally provide better accuracy. Scale experiments are also prone to mechanical friction, and we estimate friction terms to improve the correlation further. The resulting error in mean power between numerical and physical models is approximately 10%. Predicted device movement shows a good match. Overall, drag terms in time domain wave energy modelling will improve simulation accuracy in wave renewable energy device design. to understand the performance trade-offs for a large number of configurations. This requires an efficient, large displacement, multidirectional, multi-body numerical scheme. Results from a 1/25 scale wave basin experiment are described. Here, we show that a time domain linear potential flow formulation (Nemoh, WEC-Sim) can match the tank testing provided that suitably tuned drag coefficients are employed. Inviscid linear potential models can match some wave device experiments; however, additional viscous terms generally provide better accuracy. Scale experiments are also prone to mechanical friction, and we estimate friction terms to improve the correlation further. The resulting error in mean power between numerical and physical models is approximately 10%. Predicted device movement shows a good match. Overall, drag terms in time domain wave energy modelling will improve simulation accuracy in wave renewable energy device design.

16 TIDAL AND WAVE POWER↗