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.
AI-Based Faster-Than-Real-Time Stability Assessment of Large Power Systems with Applications on WECC System
Achieving clean energy goals will require significant advances in regard to addressing the computational needs for next-generation renewable-dominated power grids. One critical obstacle that lies in the way of transitioning today’s power grid to a renewable-dominated power grid is the lack of a faster-than-real-time stability assessment technology for operating a fast-changing power grid. This paper proposes an artificial intelligence (AI) -based method that predicts the system’s stability margin information (e.g., the frequency nadir in the frequency stability assessment and the critical clearing time (CCT) value in the transient stability assessment) directly from the system operating conditions without performing the conventional time-consuming time-domain simulations over detailed dynamic models. Since the AI method shifts the majority of the computational burden to offline training, the online evaluation is extremely fast. This paper has tested the AI-based stability assessment method using multiple dispatch cases that are converted and tuned from actual dispatch cases of the Western Electricity Coordinating Council (WECC) system model with more than 20,000 buses. The results show that the AI-based method could accurately predict the stability margin of such a large power system in less than 0.2 milliseconds using the offline-trained AI agent. Therefore, the proposed method has great potential to achieve faster-than-real-time stability assessment for practical large power systems while preserving sufficient accuracy.
Numerical Modeling and Optimization of the iProTech Pitching Inertial Pump (PIP) Wave Energy Converter (WEC) (Cooperative Research and Development Final Report, CRADA Number: CRD-22-22968)
This work generated a first-of-its-kind automated workflow to couple time-domain simulations of wave energy converters written in one software language with a set of design generation and evaluation scripts written in another software language. This automated workflow used an existing optimization package to analyze the sensitivity of different design parameters on the power output of a specific WEC, iProTech’s Pitching Inertial Pump (PIP). Geometric, inertial, and power take-off variables were all varied and optimized to find values that produced the highest amount of power generated over varying wave conditions. The findings on these parameter sensitivity studies are used to inform future design iterations of the PIP WEC. Including more design variables in the optimizations will only increase computational run time and further software development is needed to analyze a larger optimization.
Simple strategy for the simulation of axially symmetric large-area metasurfaces
Metalenses are composed of nanostructures for focusing light and have been widely explored in many exciting applications. However, their expanding dimensions pose simulation challenges. We propose a method to simulate metalenses in a timely manner using vectorial wave and ray tracing models. We sample the metalens’s radial phase gradient and locally approximate the phase profile by a linear phase response. Each sampling point is modeled as a binary blazed grating, employing the chosen nanostructure, to build a transfer function set. The metalens transmission or reflection is then obtained by applying the corresponding transfer function to the incoming field on the regions surrounding each sampling point. Fourier optics is used to calculate the scattered fields under arbitrary illumination for the vectorial wave method, and a Monte Carlo algorithm is used in the ray tracing formalism. We validated our method against finite-difference time domain simulations at 632 nm, and we were able to simulate metalenses larger than 3000 wavelengths in diameter on a personal computer.
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.
kynema-fmb [SWR-23-07]
Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.
Imaging Photonic Resonances within an All‐Dielectric Metasurface via Photoelectron Emission Microscopy
Dielectric metasurfaces, through volume‐type photonic resonances, enable precise control of light‐matter interactions for applications including imaging, holography, and sensing. The application space of dielectric metasurfaces has extended from infrared to visible wavelengths by incorporating high refractive index materials, such as titanium dioxide (TiO 2 ). Understanding the fundamental and fabrication limits for these applications requires metrology with nanoscale resolution, sensitivity to electromagnetic fields within the meta‐atom volume, and far‐field excitation. In this work, photoelectron emission microscopy (PEEM) is used to image field distributions of photonic resonances in a TiO 2 metasurface excited with far‐field, visible‐wavelength illumination. The local volumetric field variations within the meta‐atoms are analyzed as a function of illumination angle and polarization by comparing photoelectron images to finite‐difference time‐domain simulations. This study determines the inelastic mean free path of very low‐energy (<1 eV) photoelectrons to be 35 ± 10 nm, which is comparable to the meta‐atom height thereby highlighting PEEM sensitivity to resonances within the volume. Additionally, the simulations reveal high sensitivity of PEEM images to an in‐plane component of the illumination k ‐vector. These results demonstrate that photoelectron imaging with subwavelength resolution offers unique advantages for examining light‐matter interactions in volume‐type (as opposed to surface) photonic modes within dielectric nanophotonic structures.
Non-close-packed plasmonic Bravais lattices through a fluid interface-assisted colloidal assembly and transfer process
Abstract The assembly of colloids at fluid interfaces followed by their transfer to solid substrates represents a robust bottom-up strategy for creating colloidal monolayers over large, macroscopic areas. In this study, we showcase how subtle adjustments in the transfer process, such as varying the contact angle of the substrate and controlling deposition speed and direction, enable the realization of all five two-dimensional Bravais lattices. Leveraging plasmonic core–shell microgels as the building blocks, we successfully engineered non-close-packed plasmonic lattices exhibiting hexagonal, square, rectangular, centered rectangular, and oblique symmetries. Beyond characterizing the monolayer structures and their long-range order, we employed extinction spectroscopy alongside finite difference time domain simulations to comprehensively investigate and interpret the plasmonic response of these monolayers. Additionally, we probed the influence of the refractive index environment on the plasmonic properties by two methods: first, by plasma treatment to remove the microgel shells, and second, by overcoating the resulting gold nanoparticle lattices with a homogeneous refractive index polymer film. Graphical Abstract
Performance analysis and tank test validation of a hybrid ocean wave-current energy converter with a single power takeoff
This paper introduces a hybrid ocean wave-current energy converter (HWCEC) that harvests energy from ocean waves and current simultaneously with a single power takeoff (PTO). Specifically, the wave energy is extracted through relative heaving motion between a floating buoy and a submerged second body, while the current energy is extracted using a marine current turbine (MCT). Energy from both sources are integrated by a hybrid PTO whose concept is based on a mechanical motion rectifier (MMR). The hybrid PTO with two one-way clutches converts bidirectional, up-and-down motion from the waves into unidirectional rotation of the generator. Meanwhile a third one-way clutch couples the MCT with the same generator. The wave and current can simultaneously or separately drive the same generator through different engagement and disengagement statuses of the one-way clutches. Time-domain simulation is conducted with hydrodynamic coefficients obtained from computational fluid dynamics analysis and boundary element method. Tank tests were conducted for a HWCEC under co-existing wave and current inputs. For comparison, separate baseline tests of a turbine and a two-body point absorber, each acting in isolation, are conducted. Experimental results validate the dynamic modeling and show that a HWCEC can increase the output power with a range between 29-87% over either current turbine and wave energy converter acting individually, and it can reduce by up to 70% the peak-to-average power ratio compared with the wave energy converter on the tested conditions.
Joint scheduling of energy, fast and primary frequency response reserves in integrated transmission–distribution networks
Inverter-based distributed energy resources (DERs) connected to distribution networks (DNs) can provide fast frequency support, but their reserve deliverability depends on feeder constraints and differs from synchronous primary frequency response (PFR). Existing transmission–distribution coordination studies usually treat reserve generically or neglect feeder-level feasibility, while frequency-security scheduling studies rarely represent distribution feeders explicitly. This paper develops a bi-level day-ahead scheduling framework for integrated transmission–distribution networks that jointly clears energy, transmission-side PFR, and distribution-side fast frequency response (FFR) under exogenous hourly inertia and largest-loss inputs from an external unit commitment (UC) schedule. The transmission problem is modeled with DC-optimal power flow (OPF) and closed-form second-order cone (SOC) frequency-security constraints, whereas each DN is represented by a reserve-aware branch-flow AC-OPF so that scheduled fast reserves remain deliverable during activation. The bi-level problem is reformulated through Karush–Kuhn–Tucker (KKT) conditions into a mixed-integer SOC program, and a penalty term is used to tighten the distribution-network relaxation. In the reduced test system, lower exogenous inertia increased the required primary response from 179.64 MW to 191.08 MW, distribution-side fast response reduced total frequency-response procurement by up to 4.9%, and neglecting distribution constraints overstated the combined distribution-side energy and reserve award by up to 18%. In the expanded study, the largest case was solved in 2.02 s with a 0.00% optimality gap. Time-domain simulations kept the frequency nadir above 59.0 Hz in all tested hours. These results demonstrate the value of fast-response modeling and distribution-feasible reserve delivery in coordinated market clearing.
Process window estimation in manufacturing through Entropy-Sigma active learning
In manufacturing, there exist boundary identification problems for defining parameter spaces that meet desired thresholds on outcomes. This paper presents an Entropy-Sigma acquisition function for active learning of the process window/map in manufacturing using a Gaussian Process surrogate. Here, the method is applied to identify the stability boundary for the stability process map in machining using time-domain simulations with a periodic sampling stability metric. Results show that the proposed Entropy-Sigma method significantly outperforms Latin hypercube sampling or grid-based methods. The described method can be applied to identify the process window/map for any manufacturing application using a quantitative process outcome metric.
Cutting force and stability for inserted cutters using structured light metrology
This paper describes force and stability prediction for inserted end and face mills using a reverse engineering approach. Structured light scanning is used to identify the spatial coordinates of the points that define the multiple insert cutting edges. These points are used to extract the cutting edge radius and angle and are then incorporated in a time domain simulation that predicts cutting force and tool/workpiece deflection for user-selected operating parameters. Good agreement between predicted and measured cutting forces is demonstrated. The process stability is then predicted using the validated force model and measured structural dynamics.
Cutting force and stability prediction for inserted cutters
Here, this paper presents an approach for inserted end/face mill force and stability predictions using a reverse engineering approach. Structured light scanning is implemented to identify the spatial coordinates of the points that define the multiple insert cutting edges. These points are analyzed to extract the cutting edge radius and angle values, which are then incorporated in a time domain simulation that predicts cutting force and tool/workpiece deflection for user-selected operating parameters. Good agreement between predicted and measured cutting forces is first demonstrated. The process stability is then predicted using the validated force model and measured structural dynamics. Observations about the stability dependence on insert geometry are finally presented for a selected end mill.
Digital force prediction for milling
This paper presents a fully digital, integrated approach for milling force prediction with arbitrary end mill-work material combinations. The approach includes: 1) structured light scanning to identify the end mill’s cutting edge macro-geometry along the tool axis; 2) structured light scanning to measure the cutting edge cross-sectional rake and relief profiles; 3) finite element analysis of orthogonal cutting to determine the force model coefficients that relate the force to chip area using the work material’s constitutive model and tool’s rake and relief profiles; and 4) time domain simulation with inputs that include the measured cutting edge macro-geometry, finite element-based force model, and measured structural dynamics. Milling force predictions are compared to in-process measurements to validate the method.
Cost-optimal wave-powered persistent oceanographic observation
Not Available
Comparable Enhancement of TERS Signals from WSe 2 on Chromium and Gold
Plasmonic tip-sample junctions, at which the incident and scattered optical fields are localized and optimally enhanced, are often exploited to achieve ultrasensitive and highly spatially localized tip-enhanced Raman scattering (TERS). Recent work has demonstrated that the sensitivity and spatial resolution that are required to probe single molecules are attainable in such platforms. In this work, we observe and rationalize comparable TERS from few-layer WSe 2 single crystals exfoliated onto Au- and Cr-coated Si substrates, using a plasmonic TERS probe excited with a 638-nm laser. Our experimental observations are supported by finite-difference time-domain simulations that illustrate that the attainable field enhancement factors at the Au-Au and Au-Cr tip-sample nanojunctions are comparable in magnitude. Furthermore, through a combined experimental and theoretical analysis, we propose that besides Au/Ag, several metallic substrates may be used to record bright TERS spectral images.
Interference-Limited Absorption in Dense Molecular Nanolayers Near Reflecting Surfaces
We investigate linear resonant absorption by a dense ensemble of molecules confined to a sub-wavelength layer in two geometries: (i) a free-standing film in a homogeneous space and (ii) the same film placed at a controlled distance from a reflecting surface. In both cases, increasing the effective light–matter coupling (via molecular density/oscillator strength) produces a nonmonotonic response: absorption rises to an optimum and then decreases as the film becomes increasingly radiatively bright and reflective. Finite-difference time-domain simulations and analytical transfer-matrix calculations agree quantitatively and yield compact ridge conditions for the optimum. We interpret the trends using a scattering/port picture: the isolated film is a symmetric two-port system (reflection and transmission), which bounds single-sided resonant absorption to ≤50% in the ultrathin limit (reflecting transition saturation), whereas adding a mirror suppresses transmission and converts the structure into an effectively one-port absorber. In the mirror-backed geometry, interference can cancel reflection, and unity absorption is obtained at critical coupling, when radiative leakage is balanced by intrinsic molecular loss. These results clarify fundamental limits and design rules for collective absorption in dense molecular layers near dielectric or metallic boundaries.
Tuning Hybrid exciton–Photon Fano Resonances in Two-Dimensional Organic–Inorganic Perovskite Thin Films
As easy-to-grow quantum wells with narrow excitonic features at room temperature, two-dimensional (2D) Ruddleson–Popper perovskites are promising for realizing novel nanophotonic devices based on exciton–photon interactions. Here, we demonstrate a distinct hybrid exciton–photon Fano resonance in (C 4 H 9 NH 3 ) 2 PbI 4 thin films prepared via spin coating. Using a classical coupled-oscillator model and finite-difference time-domain simulations, we link the Fano interference to the coupling of the exciton with the Rayleigh-like scattering of the film microstructure. Combining colloidal plasmonic cavities with the 2D perovskite films, we demonstrate tuning of the Fano resonance. In combination with silver nanoparticles, the exciton–photon Fano interference couples to the in-plane plasmonic modes with indications of Rabi splitting. By creating a nanoparticle on mirror geometry, we address the out-of-plane excitonic component, reaching an intermediate coupling regime. Furthermore, these structures suggest possible photonic targets for biomolecular self-assembly applications.