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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Examining the Feasibility of Modeling the Inner Current Loop of a Grid-Following Inverter in the Phasor-Based Solver

Historically, when modeling a grid-following (GFL) inverter in the phasor domain, the fast inner current loop is typically ignored. To achieve a more accurate modeling, the authors have attempted to model the inner current loop in the phasor domain, but experienced significant numerical stability issues. This paper conducts a detailed analysis and provides insight into the feasibility of modeling the inner current loop of a GFL inverter in the phasor domain. It points out that: a) because of the neglection of the inverter filter dynamic, the bandwidth of the inner current loop becomes infinite with typical parameters in the phasor domain, resulting in an inaccurate representation of the actual inner current loop, whose bandwidth is typically between 600 Hz to 2 kHz, and b) the relatively large simulation time step in the phasor simulation causes numerical stability issues when simulating such an inner current loop with an infinite bandwidth. Detailed electromagnetic (EMT) and phasor modeling, and frequency response analysis are performed to explain this phenomenon. The findings of this paper explain why it is inappropriate to model the fast current loop of a GFL inverter in the phasor domain.

Mohiuddin, Sheik Mohammad↗

DOE OSTI2025

Hard–Soft Acid–Base Theory Explains Photoexcited Carrier Dynamics in Porphyrin/CNT Nanohybrids: Time-Domain Atomistic Analysis

We employ the fundamental chemical concepts of hard− soft acid−base to formulate general principles governing excited-state dynamics in zinc porphyrin (ZnP)/carbon nanotube (CNT) hybrids for energy photoconversion. Atomistic quantum dynamics simulations demonstrate that electron-withdrawing and donating substituents at the ZnP β-pyrrolic position strongly influence the dynamics. ZnP photoexcitation produces subpicosecond electron transfer (ET) from ZnP to CNT, in agreement with the experiment. Substitutions of CN by H and tBu accelerate the ET. The trend is directly related to the hard−soft acid−base concept because the soft−soft interaction between the t Bu- ZnP acid and the mild CNT base enhances the donor−acceptor coupling. Longer coherence and more active vibrational modes facilitate the ET in t Bu-ZnP/CNT. Electron−hole recombination in CN-ZnP/CNT occurs on a hundred picosecond time scale, nicely corroborated by the experiment. The exciton lifetime is extended beyond a nanosecond by the substitutions. The soft−soft interaction in t Bu-ZnP/ CNT increases the splitting between the highest occupied orbitals of the two subsystems, reduces their mixing, and decreases nonadiabatic coupling between the ground and excited states. Rapid decoherence and involvement of low-frequency vibrations favor longer lifetimes. Our investigation reveals that larger p K a of the β-pyrrolic acid gives rapid ET and slow recombination and provides detailed mechanistic information, essential for future optoelectronic applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Localized Exponential Time Differencing Method for Shallow Water Equations: Algorithms and Numerical Study

Here, we explore the performance of the exponential time differencing (ETD) method applied to the rotating shallow water equations. Comparing with explicit time stepping of the same order accuracy in time, the ETD algorithms could reduce the computational time in many cases by allowing the use of large time step sizes while still maintaining numerical stability. To accelerate the ETD simulations, we propose a localized approach that synthesizes the ETD method and overlapping domain decomposition. By dividing the original problem into many subdomain problems of smaller sizes and solving them locally, the proposed approach could speed up the calculation of matrix exponential vector products. Several standard test cases for shallow water equations of one or multiple layers are considered. The results show great potential of the localized ETD method for high-performance computing because each subdomain problem can be naturally solved in parallel at every time step.

58 GEOSCIENCES↗

DOE OSTI2021

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗

DOE OSTI2021

Intermediate scattering functions of a rigid body monoclonal antibody protein in solution studied by dissipative particle dynamic simulation

In the past decade, there was increased research interest in studying internal motions of flexible proteins in solution using Neutron Spin Echo (NSE) as NSE can simultaneously probe the dynamics at the length and time scales comparable to protein domain motions. However, the collective intermediate scattering function (ISF) measured by NSE has the contributions from translational, rotational, and internal motions, which are rather complicated to be separated. Widely used NSE theories to interpret experimental data usually assume that the translational and rotational motions of a rigid particle are decoupled and independent to each other. To evaluate the accuracy of this approximation for monoclonal antibody (mAb) proteins in solution, dissipative particle dynamic computer simulation is used here to simulate a rigid-body mAb for up to about 200 ns. The total ISF together with the ISFs due to only the translational and rotational motions as well as their corresponding effective diffusion coefficients is calculated. The aforementioned approximation introduces appreciable errors to the calculated effective diffusion coefficients and the ISFs. For the effective diffusion coefficient, the error introduced by this approximation can be as large as about 10% even though the overall agreement is considered reasonable. Thus, we need to be cautious when interpreting the data with a small signal change. In addition, the accuracy of the calculated ISFs due to the finite computer simulation time is also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Characterizing patterns and variability of building electric load profiles in time and frequency domains

The rapid development of advanced metering infrastructure provides a new data source—building electrical load profiles with high temporal resolution. Electric load profile characterization can generate useful information to enhance building energy modeling and provide metrics to represent patterns and variability of load profiles. Such characterizations can be used to identify changes to building electricity demand due to operations or faulty equipment and controls. In this study, we proposed a two-path approach to analyze high temporal resolution building electrical load profiles: (1) time-domain analysis and (2) frequency-domain analysis. Furthermore, the commonly adopted time-domain analysis can extract and quantify the distribution of key parameters characterizing load shape such as peak-base load ratio and morning rise time, while a frequency-domain analysis can identify major periodic fluctuations and quantify load variability. We implemented and evaluated both paths using whole-year 15-minute interval smart meter data of 188 commercial office building in Northern California. The results from these two paths are consistent with each other and complementary to represent full dynamics of load profiles. The time- and frequency-domain analyses can be used to enhance building energy modeling by: (1) providing more realistic assumptions about building operation schedules, and (2) validating the simulated electric load profiles using the developed variability metrics against the real building load data.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

PFLOTRAN-SIP: A PFLOTRAN Module for Simulating Spectral-Induced Polarization of Electrical Impedance Data

Spectral induced polarization (SIP) is a non-intrusive geophysical method that collects chargeability information (the ability of a material to retain charge) in the time domain or its phase shift in the frequency domain. Although SIP is a temporal method, it cannot measure the dynamics of flow and solute/species transport in the subsurface over long times (i.e., 10–100 s of years). Data collected with the SIP technique need to be coupled with fluid flow and reactive-transport models in order to capture long-term dynamics. To address this challenge, PFLOTRAN-SIP was built to couple SIP data to fluid flow and solute transport processes. Specifically, this framework couples the subsurface flow and transport simulator PFLOTRAN and geoelectrical simulator E4D without sacrificing computational performance. PFLOTRAN solves the coupled flow and solute-transport process models in order to estimate solute concentrations, which were used in Archie’s model to compute bulk electrical conductivities at near-zero frequency. These bulk electrical conductivities were modified while using the Cole–Cole model to account for frequency dependence. Using the estimated frequency-dependent bulk conductivities, E4D simulated the real and complex electrical potential signals for selected frequencies for SIP. These frequency-dependent bulk conductivities contain information that is relevant to geochemical changes in the system. This study demonstrated that the PFLOTRAN-SIP framework is able to detect the presence of a tracer in the subsurface. SIP offers a significant benefit over ERT in the form of greater information content. It provided multiple datasets at different frequencies that better constrained the tracer distribution in the subsurface. Consequently, this framework allows for practitioners of environmental hydrogeophysics and biogeophysics to monitor the subsurface with improved resolution.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Smoothed boundary method for simulating incompressible flow in complex geometries

Simulating flow through porous media with explicit considerations of complex microstructures is very challenging using conventional sharp-interface methods because of the difficulties in generating meshes conformal to complex geometries. In this work, a diffuse interface embedded boundary method known as the Smoothed Boundary Method (SBM) is utilized to facilitate simulations of fluid dynamics involving complex geometries. In diffuse-interface methods, the geometry is described by a domain parameter. The SBM allows the straightforward reformulation of the time-dependent Navier-Stokes equations in terms of this domain parameter, using only algebraic identities. Thus, enforcing the appropriate boundary conditions at the irregular embedded boundary is greatly simplified. Adaptive mesh refinement is used to increase the accuracy of the diffuse interface method by allowing thinner interfacial thickness to be used in the domain parameter. Here, the SBM-formulated Navier-Stokes equations are solved with the Finite Difference Method on refined mesh systems. Sharp-interface Finite Element Method simulations using the commercial software COMSOL on body-conforming meshes are also provided for comparison. Favorable agreement between the two methods is observed. Since it is no longer necessary for the mesh to conform to the complex geometry, the grid system for the SBM simulations can be generated rapidly and without additional manual interventions, making the entire simulation process more expedient.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Fluorine-Free Precise Polymer Electrolyte for Efficient Proton Transport: Experiments and Simulations

Designing polymers with controlled nanoscale morphologies and scalable synthesis is of great interest in the development of fluorine-free materials for proton-exchange membranes in fuel cells. This study focuses on a precision polyethylene with phenylsulfonic acid branches at every fifth carbon, p5PhSA, with a high ion-exchange capacity (4.2 mmol/g). The polymers self-assemble into hydrophilic and hydrophobic co-continuous nanoscale domains. In the hydrated state, the hydrophilic domain, composed of polar sulfonic acid moieties and water, serves as a pathway for efficient mesoscopic proton conductivity. The morphology and proton transport of p5PhSA are evaluated under hydrated conditions using in situ X-ray scattering and electrochemical impedance spectroscopy techniques. Additionally, at 40 °C and 95% relative humidity, the proton conductivity of p5PhSA is 0.28 S/cm, which is four times greater than Nafion 117 under the same conditions. Atomistic molecular dynamics (MD) simulations are also used to elucidate the interplay between the structure and the water dynamics. The MD simulations show strong nanophase separation between the percolated hydrophilic and hydrophobic domains over a wide range of water contents. The percolated hydrophilic nanoscale domain facilitates the rapid proton transport in p5PhSA and demonstrates the potential of precise hydrocarbon-based polymers as processible and effective proton-exchange membranes.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Variational, stable, and self-consistent coupling of 3D electromagnetics to 1D transmission lines in the time domain

This work presents a new multiscale method for coupling the 3D Maxwell's equations to the 1D telegrapher's equations. While Maxwell's equations are appropriate for modeling complex electromagnetics in arbitrary-geometry domains, simulation cost for many applications (e.g. pulsed power) can be dramatically reduced by representing less complex transmission line regions of the domain with a 1D model. By assuming a transverse electromagnetic (TEM) ansatz for the solution in a transmission line region, we reduce the Maxwell's equations to the telegrapher's equations. Here, we propose a self-consistent finite element formulation of the fully coupled system that uses boundary integrals to couple between the 3D and 1D domains and supports arbitrary unstructured 3D meshes. Additionally, by using a Lagrange multiplier to enforce continuity at the coupling interface, we allow for an absorbing boundary condition to also be applied to non-TEM modes on this boundary. We demonstrate that this feature reduces non-physical reflection and ringing of non-TEM modes off of the coupling boundary. By employing implicit time integration, we ensure a stable coupling, and we introduce an efficient method for solving the resulting linear systems. We demonstrate the accuracy of the new method on two verification problems, a transient O-wave in a rectilinear prism and a steady-state problem in a coaxial geometry, and show the efficiency and weak scalability of our implementation on a cold test of the Z-machine MITL and post-hole convolute.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Conductive Spray Paint Fabricated Cylindrical Coplanar Waveguide for Noncontact Liquid-Level Sensing

A noncontact liquid-level sensor using a cylindrical coplanar waveguide (CCPW) is proposed. The concept is demonstrated using a prototype device made by directly spraying conductive paint on the outer surface of a glass cylinder. The position of the air–liquid interface is measured using the time-domain reflection (TDR) signal from the CCPW. The geometric parameters of the prototype sensor were optimized using a full-wave simulation tool. The experimental results show that the prototype CCPW-based liquid-level sensor has a resolution of 29.8μm over a range of 210 mm at the sweeping frequency bandwidth of 500 MHz–4 GHz. The tradeoff between the resolution and the measurement range was discussed. The proposed sensor provides a low-cost, easy-to-fabricate, and noncontact method for liquid-level measurement, which has potential applications in fuel storage, chemical production, and wastewater treatment.

42 ENGINEERING↗

DOE OSTI2022

GDSA Framework Development and Process Model Integration FY2022

The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). A high priority for SFWST disposal R&D is disposal system modeling (Sassani et al. 2021). The SFWST Geologic Disposal Safety Assessment (GDSA) work package is charged with developing a disposal system modeling and analysis capability for evaluating generic disposal system performance for nuclear waste in geologic media. This report describes fiscal year (FY) 2022 advances of the Geologic Disposal Safety Assessment (GDSA) performance assessment (PA) development groups of the SFWST Campaign. The common mission of these groups is to develop a geologic disposal system modeling capability for nuclear waste that can be used to assess probabilistically the performance of generic disposal options and generic sites. The modeling capability under development is called GDSA Framework (pa.sandia.gov). GDSA Framework is a coordinated set of codes and databases designed for probabilistically simulating the release and transport of disposed radionuclides from a repository to the biosphere for post-closure performance assessment. Primary components of GDSA Framework include PFLOTRAN to simulate the major features, events, and processes (FEPs) over time, Dakota to propagate uncertainty and analyze sensitivities, meshing codes to define the domain, and various other software for rendering properties, processing data, and visualizing results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Electrical Submersible Pumps: A System Modeling Approach for Power Quality Analysis With Variable Frequency Drives

This article proposes a high-level engineering guide to develop an integrated system model for power quality analysis in electrical submersible pump (ESP) applications with variable frequency drives (VFDs). Such analyses are troublesome in the industry due to the complexity of these systems. In this article, simple steps to perform system integration analyses of such arrangements, including torsional analysis, are developed. A simplified VFD-ESP model suitable for coupled electrical and mechanical analysis in steady state is proposed in all their configurations of practical interest. Such a model can be easily implemented in common simulation software, significantly reducing engineering efforts for implementation and analysis. The focus of the proposed model is the prediction of ESP failures that might result from a poor power quality caused by VFDs. Analytical expressions of different types of harmonics in these systems, as well as their accurate locations in the frequency domain, including their interharmonics and common-mode harmonics, are derived for this purpose. Further, the effectiveness of the proposed model is verified through offline and real-time hybrid simulation results. Finally, a comparison between simulation results obtained using the proposed model and measurements collected on a down-scale laboratory prototype is carried out to demonstrate the accuracy of the suggested modeling approach.

42 ENGINEERING↗

DOE OSTI2020

On the character of turbulent-like flows in self-consistent models of core-collapse supernovae

Neutrino-driven convection plays a crucial role in the development of core-collapse supernova (CCSN) explosions. However, the complex mechanism that triggers the shock revival and the subsequent explosion has remained inscrutable for many decades. Multidimensional simulations suggest that the growth of fluid instabilities, especially development of convection, will determine the morphology of the explosion. We have performed 3D simulations using spherical-polar coordinates covering a reduced angular extent (90° computational domain), and with angular resolutions of 2°, 1°, 1/2°, and 1/4°, to study the development of turbulent-like flows in core-collapse supernova explosions on a time scale of order 100 ms. We have employed the multi-physics Chimera code that includes detailed nuclear physics and spectral neutrino transport. Coarse resolution models do not develop an inertial range, presumably due to the bottleneck effect, such that the energy is prevented from cascading down to small scales and tends to accumulate at large scales. High-resolution models instead, start to recover the k −5/3 scaling of Kolmogorov's theory. Stochasticity and few simulation samples limit our ability to predict the development of explosions. Over the simulated time period, our models show no clear trend in improving (or diminishing) conditions for explosion as the angular resolution is increased. However, we find that disordered flow provides an effective pressure, as characterized by the Reynolds stress, behind the shock. This contribution, ~40%–50% of the thermal pressure, aids shock revival and thus the development of the explosion. Finally, we show that the kinetic energy power spectrum of reduced angular extent and full 4π models are consistent, thus indicating that a 90° computational domain is an adequate configuration to study the character of turbulence in CCSNe.

Casanova Bustamante, Jordi↗

DOE OSTI2023

Fracture Intensity Impacts on Reaction Front Propagation and Mineral Weathering in Three‐Dimensional Fractured Media

Abstract Studying reaction front propagation in heterogeneous natural settings is challenging, but numerical simulations can provide insight into the varying spatial and temporal scales of reaction front propagation. Here, the impact of increasing fracture intensity on mineral dissolution rates, and the extent of reaction front propagation is investigated using reactive transport simulations in upscaled discrete fracture network domains with varied fracture intensity. Domain‐averaged dissolution rates vary less than 0.5 log units regardless of the fracture intensity, but the spatial distribution of reactions is controlled by the location and number of dead‐end fractures and the number of connected flowpaths through the domain. Higher fracture intensities lead to more weathering in the domain because of more available mineral for water‐rock interactions. We find that reaction fronts propagate through the primary flowpaths in the first 10,000 years of the simulation for a 10‐m length domain, then propagate into secondary flowpaths and dead‐end fractures between 10,000 and 100,000 years, and finally into the matrix over timescales of hundreds of thousands of years. The domain‐averaged reaction rates decrease through time corresponding to a transition from dissolution in advection‐dominated, fast‐flowing pathways, to dissolution in transport‐limited zones of disconnected fractures and matrix. Matrix dissolution, or dissolution under transport‐limited conditions, is the dominant process at late times in these simulations. The results of these simulations recreate the observed paradox found in nature where highly fractured hillslopes tend to be more weathered but have slower weathering rates, while hillslopes with fewer fractures, are less weathered but have higher dissolution rates.

58 GEOSCIENCES↗

DOE OSTI2021

Cell dynamics simulations of coupled charge and magnetic phase transformation in correlated oxides

In this work, we present a comprehensive numerical study on the kinetics of phase transition that is characterized by two nonconserved scalar order parameters coupled by a special linear-quadratic interaction. This particular Ginzburg-Landau theory has been proposed to describe the coupled charge and magnetic transition in nickelates and the collinear stripe phase in cuprates. The inhomogeneous state of such systems at low temperatures consists of magnetic domains separated by quasimetallic domain walls where the charge order is reduced. By performing large-scale cell dynamics simulations, we find a two-stage phase-ordering process in which a short period of independent evolution of the two order parameters is followed by a correlated coarsening process. The long-time growth and coarsening of magnetic domains is shown to follow the Allen-Cahn power law. We further show that the nucleation-and-growth dynamics during phase transformation to the ordered states is well described by the Kolmogorov-Johnson-Mehl-Avrami theory in two dimensions. On the other hand, the presence of quasimetallic magnetic domain walls in the ordered states gives rise to a very different kinetics for transformation to the high-temperature paramagnetic phase. In this scenario, the phase transformation is initiated by the decay of magnetic domain walls into two insulator-metal boundaries, which subsequently move away from each other. Implications of our findings to recent nano-imaging experiments on nickelates are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Stochastic simulations of temporally and spatially variable Fe cycling within a floodplain aquifer

This data package contains model input and output files for local and global sensitivity analysis (SA) of a reactive transport model simulating redox cycling within a floodplain aquifer. The reaction network is implemented in CrunchFlow and focuses on spatially and temporally heterogeneous Fe cycling, though it contains several other species and redox pathways (37 aqueous species and 7 minerals in total). The package contains data from 600 Monte Carlo simulations (for global SA) and 104 one-at-a-time simulations (for local SA). It includes all input files required to run the CrunchFlow model, as well as limited model output – .csv files containing time series of primary species concentrations at 21 grid cells within the domain. These time series files contain the following species: H+, CO32-, SO42-, Cl-, Ca2+, Mg2+, Na+, K+, Fe(II), Fe(III), HS-, CH3COO-, O2(aq), S(aq), and Br-.These data were used to demonstrate the application of a novel form of global SA, distance-based generalized sensitivity analysis (DGSA) to a reactive transport model. The 2D model simulates the export of reduced species from pockets of fine-grained, organic-rich sediments embedded in a coarse sand aquifer. Stochastic simulations jointly varied 17 key model input parameters across several orders of magnitude, including the spatial variability of each parameter. DGSA on the model results reveal that the amplitude and variability of exported Fe(II) are most sensitive to the interaction between sand permeability and individual reaction rates. By contrast, the propagation of reducing conditions downgradient from the fine-grained lenses depends most heavily on the rate of dissolved organic carbon and sulfate release from the fines.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Terahertz conductivity of two-dimensional materials: a review

Two-dimensional (2D) van der Waals materials are shaping the landscape of next-generation devices, offering significant technological value thanks to their unique, tunable, and layer-dependent electronic and optoelectronic properties. Time-domain spectroscopic techniques at terahertz (THz) frequencies offer noninvasive, contact-free methods for characterizing the dynamics of carriers in 2D materials. They also pave the path toward the applications of 2D materials in detection, imaging, manufacturing, and communication within the increasingly important THz frequency range. In this paper, we overview the synthesis of 2D materials and the prominent THz spectroscopy techniques: THz time-domain spectroscopy, optical-pump THz-probe technique, and optical pump–probe THz spectroscopy. Through a confluence of experimental findings, numerical simulation, and theoretical analysis, we present the current understanding of the rich ultrafast physics of technologically significant 2D materials: graphene, transition metal dichalcogenides, MXenes, perovskites, topological 2D materials, and 2D heterostructures. Finally, we offer a perspective on the role of THz characterization in guiding future research and in the quest for ideal 2D materials for new applications.

2D materials↗