Data-Driven Modeling of Grid-Forming Inverter Dynamics Using Power Hardware-in-the-Loop Experimentation
Not Available
24 POWER TRANSMISSION AND DISTRIBUTION↗
SEARCH · Engineering Papers
Results for “analytical model”
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.
Combined Normal and Disruption, Electromagnetic Transient, Thermal, and Structural Analysis of COMPASS Upgrade
We report that COMPASS Upgrade is a replacement for the COMPASS device at the Institute of Plasma Physics (IPP) of the Czech Academy of Sciences. It is a cryogenic copper machine with a major radius of 0.89 m, with 5 T at the plasma centerline and 2-MA plasma current. It is intended to develop ITER relevant plasma densities and high power fluxes in the divertor region. The entire vacuum vessel is planned to reach 500 °C. The COMPASS-U TF coil is a bitter plate-like design with radial conductor plates that extend the width of the TF inner leg. Currents redistribute radially based on resistive and inductive effects. This necessitates an electromagnetic (EM) transient analysis of the coil current distribution. For copper coils like COMPASS Upgrade, the solution must be coupled with a thermal solution to properly include the temperature-dependent resistive effects. The EM transient solution includes the TF magnetic loads, and a stress pass on the results produces in-plane loads. The EM solution chosen requires inclusion of all the conducting and nonconducting regions surrounding the TF coil. With the TF current distribution solved, the addition of the poloidal coils and resulting background fields allows the determination of the out-of-plane loads. Integration of the Joule heating produced from the EM solution produces temperature distributions throughout the shot, which can be read in to a structural pass to include thermal stresses in the evaluation. This is true of both TF and PF temperature prediction with appropriate packing fractions. Currents in the PF coils, with computed background fields, produce Lorentz loads. Time transients of the PF currents produce startup eddy currents in the structures included in the model; inclusion of the vessel and other passive structures allows these to be included in a stress pass. The EM transient solution includes all the elements needed for a disruption analysis with prescribed motions and quench of the plasma. From a model originally intended to simulate current diffusion in the TF coil, the analysis can be refined to be used to predict the performance of most of the components of the tokamak or provide boundary conditions for more detailed submodels. In this article, the use of the current diffusion model for bounding assessments of the COMPASS Upgrade TF, PF, vessel, external structure, and disruption analyses will be presented.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A Multi-Model, Multi-Scale Research Program in Stressors, Responses, and Coupled Systems Dynamics at the Energy-Water-Land Nexus and for Concentrated, Interdependent Infrastructures: Toward Next Generation Capabilities in Integrated Impacts, Adaptation, and Vulnerability (I-IAV) Modeling and a Community of Practice
The goal of this research program was to build a next generation integrated suite of science-driven modeling and analytic capabilities, and a more expanded and connected community of practice, for analyses of the stressors, impacts, adaptations and vulnerabilities of global and regional change. The emphasis was on understanding energy-water-land interactions and feedbacks and interdependent infrastructures at appropriate regional and temporal scales. Although the scope spans many complex facets of data, modeling, and analysis, as well as scales appropriate for integrated impacts and adaptation research, the focus of this effort was the development of multi-model, multi-scale capabilities spanning the domains of Multi-Sector Dynamics (MSD) models; Impact, Adaptation, and Vulnerability (IAV) models; and Earth System Models (ESMs).
54 ENVIRONMENTAL SCIENCES↗
Semi-Analytical Hierarchical Bayesian Inference of Nonlinear Model Structure in Stochastic Dynamics: Applied to Compartmental Models of Infectious Diseases
A Bayesian computational framework for parsimonious inference in stochastic nonlinear dynamical systems is presented. This framework enables the concurrent estimation of system states, time-varying parameters, time-invariant parameters, and the optimal sparsity structure of the model parameters. Because differential equation-based models are often simplified mechanistic or phenomenological representations, robust inference from noisy measurement data requires explicit treatment of model error and uncertainty. Model error and time-varying parameters can be represented as random processes, enabling inference while making minimal assumptions about the underlying sources of discrepancy and variability. Adopting stochastic differential equation representations affords the model significant flexibility, but can also render it susceptible to overfitting during statistical inversion, where the inferred model may track noise rather than the underlying signal. To alleviate the effects of overfitting and to enable the discovery of the optimal sparse representation of the time-invariant parameters, a Bayesian sparse learning algorithm is embedded within the framework. This sparse learning framework adopts an approximate hierarchical Bayesian setting defined by a series of semi-analytical expressions. The model structure inference framework is validated using a stochastic compartmental model for tracking and forecasting active cases of an infectious disease. Compartmental models describe population-level infectious disease dynamics through interactions among population fractions grouped by disease state. Mathematically, such models consist of a system of coupled ordinary differential equations. This example adopts an expressive compartmental model that includes multiple possible interactions between disease states, motivated by early uncertainty surrounding COVID-19 reinfection dynamics and their implications for long-term epidemic forecasting. The sparse learning exercise permits the inference of a priori unknown epidemiological dynamics from simulated public health data, discovering the nested compartmental model that optimizes the trade-off between average data-fit and model complexity. It is shown that inducing sparsity among the model parameters eliminates redundant interactions between compartments, equivalently revealing the optimal coupling structure between differential equations.
97 MATHEMATICS AND COMPUTING↗
gcamfaostat: An R package to prepare, process, and synthesize FAOSTAT data for global agroeconomic and multisector dynamic modeling
The gcamfaostat R package is designed for the preparation, processing, and synthesis of the Food and Agriculture Organization (FAO) Statistics (FAOSTAT) agroeconomic data. The primary purpose is to facilitate FAOSTAT data use in global economic and multisector dynamic models while ensuring transparency, traceability, and reproducibility. Here, we provide an overview of the development of gcamfaostat (v1.0) and demonstrate its capabilities in generating and maintaining agroeconomic data required for the Global Change Analysis Model (GCAM). Our initiative seeks to enhance the quality and accessibility of data for the global agroeconomic modeling community, with the aim of fostering more robust and harmonized outcomes in a collaborative, efficient, and open-source framework. One of the important features of the package is the possibility to construct the FAO Food Balance Sheets at the disaggregated commodity level (with over 500 commodities), which provides a comprehensive and detailed data input for a variety of analytical and modeling applications. The processed data and visualizations offered by gcamfaostat can also be valuable to a broader audience interested in gaining insights into the intricacies of global agriculture.
97 MATHEMATICS AND COMPUTING↗
Kinetic, metabolic, and statistical analytics: addressing metabolic transport limitations among organelles and microbial communities
Microbial organisms engage in a variety of metabolic interactions. A crucial part of these interactions is the exchange of molecules between different organelles, cells, and the environment. The main forces mediating this metabolic exchange are transporters. This transport can be difficult to measure experimentally because several transport mechanisms remain opaque. However, theoretical calculations about the inputs and outputs of cells via metabolic exchanges have enabled the successful inference of the workings of intra-organismal and inter-organismal systems. Kinetic, metabolic, and statistical modeling approaches in combination with omics data are enhancing our knowledge and understanding about metabolic exchange and mass resource allocation. Furthermore, this model-driven analytics approach can guide effective experimental design and yield new insights into biological function and control.
59 BASIC BIOLOGICAL SCIENCES↗
Comparison of Polarized Radiative Transfer Codes Used by the EHT Collaboration
Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT Collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I , Q , U , and V , respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties.
79 ASTRONOMY AND ASTROPHYSICS↗
Strong Dark Matter Self-interactions Diversify Halo Populations within and surrounding the Milky Way
Abstract We perform a high-resolution cosmological zoom-in simulation of a Milky Way (MW)–like system, which includes a realistic Large Magellanic Cloud analog, using a large differential elastic dark matter self-interaction cross section that reaches ≈100 cm 2 g −1 at relative velocities of ≈10 km s −1 , motivated by the diverse and orbitally dependent central densities of dwarf galaxies within and surrounding the MW. We explore the effects of dark matter self-interactions on satellite, splashback, and isolated halos through their abundance, central densities, maximum circular velocities, orbital parameters, and correlations between these variables. We use an effective constant cross section model to analytically predict the stages of our simulated halos’ gravothermal evolution, demonstrating that deviations from the collisionless R max – V max relation can be used to select deeply core-collapsed halos, where V max is a halo’s maximum circular velocity, and R max is the radius at which it occurs. We predict that a sizable fraction (≈20%) of subhalos with masses down to ≈10 8 M ⊙ is deeply core collapsed in our SIDM model. Core-collapsed systems form ≈10% of the isolated halo population down to the same mass; these isolated, core-collapsed halos would host faint dwarf field galaxies with extremely steep central density profiles. Finally, most halos with masses above ≈10 9 M ⊙ are core-forming in our simulation. Our study thus demonstrates how self-interactions diversify halo populations in an environmentally dependent fashion within and surrounding MW-mass hosts, providing a compelling avenue to address the diverse dark matter distributions of observed dwarf galaxies.
79 ASTRONOMY AND ASTROPHYSICS↗
Applications of Multibody for Everybody (M4E) in Marine Energy
Multibody for Everybody (M4E) is an open-source symbolic dynamics modeling framework designed to automate the derivation of equations of motion and simulation of constrained multibody systems using the joint coordinate formulation.
16 TIDAL AND WAVE POWER↗
Amorphization and siliconization of silicon carbide as a first wall material
The understanding and prediction of silicon carbide (SiC) material evolution exposed to SOL plasma conditions is of prime interest because SiC represents a promising main chamber wall plasma-facing material for next-step fusion devices (low hydrogenic diffusion, good mechanical and thermal properties under neutron irradiation). Gross and net Si erosion rates from SiC surfaces in contact with a well-diagnosed L-mode plasma in the DIII-D tokamak have been simulated and the surface concentrations of impurities have been tracked as a function of time. Coupled simulation of surface model and impurity transport demonstrates amorphization of crystalline SiC exposed to L-mode plasma due to the accumulation of displacement damages under ion irradiation. This affects the lifetime of SiC plasma facing components. Surface evolution is tightly coupled to impurity transport in the plasma and therefore needs to be integrated with impurity transport simulations to effectively predict Si erosion rates and sub-surface concentrations as a function of time. The simulation workflow couples a semi-analytical surface model to the impurity transport code GITR. The surface model is a homogeneous mixed-material model that tracks physical & chemical sputtering and reflection of impurities. Gross erosion is primarily influenced by the background plasma parameters and redeposition patterns are mainly influenced by the prompt redeposition due to the gyro-orbits of impurity ions. Although crystalline form of SiC is preferable for fusion wall applications because of resistance to neutron irradiation, this work indicates that crystalline SiC will undergo amorphization under D plasma contact with implications of higher sputtering and fuel retention. These results direct us to explore the effects of amorphization on crystalline SiC and further the physics basis of SiC usage as first wall material for fusion environments.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Steady-State Analysis of Microgrid Distributed Control Under Denial of Service Attacks
Distributed control is an effective method to coordinate multiple components in a microgrid. It usually has mandatory requirements for communication graph topology; changing the topology unexpectedly could disrupt the operation of a distributed controller. A denial of service (DoS) attack blocks communication channels to cause variations in communication graph. Unforeseen topology changes render a controller to receive incomplete system information. Control decisions made based on incorrect information can lead to significantly reduced grid supporting capability. As a result, the vital power balance between generation and load may be challenged. Determining power balance requires steady-state power flow analysis; however, the cyber-induced steady-state analysis is not well studied in the literature. Here, in this article, we show that a cyber–physical microgrid with unknown DoS attacks has augmented uncertain power flow equations, whose special structure creates new computational issue to determining power balance. To deal with the issue, we develop a novel cyber-induced microgrid model: the cyber–physical system is transformed into a pure physical model, where the disruption effects of DoS in the cyber layer are shown to be equivalent to the variations of the physical model parameters. With the transformed system, we develop a sufficient condition to ascertain the system power balance under unknown DoS attacks. In addition, practical insights for resilient distributed control design are shown as well.
24 POWER TRANSMISSION AND DISTRIBUTION↗
Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling
Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.
58 GEOSCIENCES↗
Single-photon absorption and emission from a natural photosynthetic complex
Photosynthesis is generally assumed to be initiated by a single photon from the Sun, which, as a weak light source, delivers at most a few tens of photons per nanometre squared per second within a chlorophyll absorption band1. Yet much experimental and theoretical work over the past 40 years has explored the events during photosynthesis subsequent to absorption of light from intense, ultrashort laser pulses. Here, we use single photons to excite under ambient conditions the light-harvesting 2 (LH2) complex of the purple bacterium Rhodobacter sphaeroides, comprising B800 and B850 rings that contain 9 and 18 bacteriochlorophyll molecules, respectively. Excitation of the B800 ring leads to electronic energy transfer to the B850 ring in approximately 0.7 ps, followed by rapid B850-to-B850 energy transfer on an approximately 100-fs timescale and light emission at 850–875 nm. Using a heralded single photon source20,21 along with coincidence counting, we establish time correlation functions for B800 excitation and B850 fluorescence emission and demonstrate that both events involve single photons. We also find that the probability distribution of the number of heralds per detected fluorescence photon supports the view that a single photon can upon absorption drive the subsequent energy transfer and fluorescence emission and hence, by extension, the primary charge separation of photosynthesis. An analytical stochastic model and a Monte Carlo numerical model capture the data, further confirming that absorption of single photons is correlated with emission of single photons in a natural light-harvesting complex.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Universal function for heavy-ion fusion cross sections
A universal function for heavy-ion fusion cross sections, Y = $\sqrt{π}$ Xerfc (- X ) + exp(- X 2 ) is proposed. By scaling both the cross section σ(E) and the energy E, heavy-ion fusion cross section data are found to follow closely a universal function over the whole energy range. The scaling is developed from either a simple, empirical single-Gaussian barrier distribution model for the representation d 2 (σE)/dE 2 , or the modified Siwek-Wilczynski model. The cross section expressions of these models are analytical, and can be easily used for all heavy-ion fusion excitation functions. Thus a bench marking of heavy-ion fusion excitation functions has been achieved. Finally, a general discussion regarding the universal function is given.
low & intermediate energy heavy-ion reactions↗
WaterTAP3 (The Water Technoeconomic Assessment Pipe-Parity Platform)
The Water Technoeconomic Assessment Pipe-Parity Platform (WaterTAP3) was developed under the National Alliance for Water Innovation (NAWI) to facilitate consistent technoeconomic assessments of desalination treatment trains. The WaterTAP3 is an analytically robust modeling tool that can be used to evaluate water technology cost, energy, environmental, and resiliency tradeoffs across different water sources, sectors, and scales. The model simulates steady-state water treatment train performance and costs including flow and constituent mass balance across unit processes, based on source water conditions, configurations of treatment technologies, and system-level techno-economic assumptions. Users can build a new treatment train by connecting any number of unit processes, specific for their context and system, or selecting a train from the treatment train library. The model contains various technical and cost parameter options for a range of treatment processes and a library of influent water quality characteristics for a variety of source waters and case studies. Users can customize water quality parameters to evaluate the technology performance in their context. The model can be set up for different assessment needs including simulation, optimization, and uncertainty and sensitivity analyses. The results from WaterTAP3 can help identify trade-offs among the different system performance metrics, with insight on how particular technologies or systems promote pipe-parity. The flexibility and comprehensive scope of the tool makes it a promising solution to industry-wide water technoeconomic evaluations, leading to more informed water investment decisions and technologies. As a user-friendly, open-source platform, WaterTAP3 can be used by industry, academia, policymakers, planners, and those with or without extensive analytical experience.
Miara, Ariel↗
Water Technoeconomic Assessment Pipe-Parity Platform (WaterTAP3)
The Water Technoeconomic Assessment Pipe-Parity Platform (WaterTAP3) was developed under the National Alliance for Water Innovation (NAWI) to facilitate consistent technoeconomic assessments of desalination treatment trains. The WaterTAP3 is an analytically robust modeling tool that can be used to evaluate water technology cost, energy, environmental, and resiliency tradeoffs across different water sources, sectors, and scales. The model simulates steady-state water treatment train performance and costs including flow and constituent mass balance across unit processes, based on source water conditions, configurations of treatment technologies, and system-level techno-economic assumptions. Users can build a new treatment train by connecting any number of unit processes, specific for their context and system, or selecting a train from the treatment train library. The model contains various technical and cost parameter options for a range of treatment processes and a library of influent water quality characteristics for a variety of source waters and case studies. Users can customize water quality parameters to evaluate the technology performance in their context. The model can be set up for different assessment needs including simulation, optimization, and uncertainty and sensitivity analyses. The results from WaterTAP3 can help identify trade-offs among the different system performance metrics, with insight on how particular technologies or systems promote pipe-parity. The flexibility and comprehensive scope of the tool makes it a promising solution to industry-wide water technoeconomic evaluations, leading to more informed water investment decisions and technologies. As a user-friendly, open-source platform, WaterTAP3 can be used by industry, academia, policymakers, planners, and those with or without extensive analytical experience. A publicly available graphical user interface is currently under development.
Miara, Ariel↗
Matter power spectra in modified gravity: a comparative study of approximations and N -body simulations
ABSTRACT Testing gravity and the concordance model of cosmology, $\Lambda$CDM, at large scales is a key goal of this decade’s largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon, and K-mouflage. In particular, we compare the predictions from N-body simulations solving the full scalar field equation, two N-body codes with approximate time integration schemes, a parametrized modified N-body implementation, and the analytic halo model reaction approach. We find the modification to the $\Lambda$CDM spectrum is in 2 per cent agreement at $z\le 1$ and $k\le 1~h\,{\rm Mpc}^{-1}$ over all gravitational models and codes, in accordance with many previous studies, indicating these modelling approaches are robust enough to be used in forthcoming survey analyses under appropriate scale cuts. We further make public the new code implementations presented, specifically the halo model reaction K-mouflage implementation and the relativistic Cubic Galileon implementation.
Bose, B. (ORCID:0000000319658614)↗
Massive islands
We comment on the role of the graviton mass in recent calculations of the Page curve using holographic ideas. All reliable calculations of the Page curve in more than 2+1 spacetime dimensions have been performed in systems with massive gravitons. A crucial ingredient in these calculations is the formation of islands, regions that contribute to the entropy of degrees of freedom located elsewhere. While most often simply ignored, it is indeed true that mass of the graviton does not appear to significantly affect the calculations that appeared in the literature. We use the freedom to change the graviton mass to give an extremely simple model of analytically tractable island formation in general dimensions. We do however note that if one attempts to take the limit of zero graviton mass, any contribution from the islands disappears. This raises the question to what extent entanglement islands can play a role in standard massless gravity.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗