Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed 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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Electric Power Distribution System Model Simplification Using Segment Substitution

Quasi-static time-series (QSTS) simulation is used to simulate the behavior of distribution systems over long periods of time (typically hours to years). The technique involves repeatedly solving the load-flow problem for a distribution system model and is useful for distributed energy resource (DER) planning. When a QSTS simulation has a small time step and a long duration, the computational burden of the simulation can be a barrier to integration into utility workflows. One way to relieve the computational burden is to simplify the system model. The segment substitution method of simplifying distribution system models introduced in this paper offers model bus reduction of up to 98% with a simplification error as low as 0.2% (0.002 pu voltage). In contrast to existing methods of distribution system model simplification, which rely on topological inspection and linearization, the segment substitution method uses black-box segment data and an assumed simplified topology.

Reiman, Andrew P. (ORCID:0000000302042960)↗
DOE OSTI2021

Local structure analysis of low-temperature neutron pair distribution function coupled with molecular dynamics simulations of CH 4 and CO 2 hydrates from 2 to 210 K

CH 4 hydrates occur naturally and are an abundant potential fuel source with a corresponding risk of potent greenhouse gas release, due to their stability conditions at low temperature and high pressure. Byproduct CO 2 can be exchanged with CH 4 in these natural deposits and can potentially stabilize the hydrates at higher temperatures. CH 4 , CO 2 , and mixed CH 4 -CO 2 hydrates are studied with in situ neutron pair distribution function experiments from 2 to 210 K to investigate the impact of varying the CH 4 -CO 2 guest composition in the gas hydrate structure. These experiments combined with Reverse Monte Carlo analysis allow for the characterization of intermolecular CH 4 and CO 2 interactions with the water molecules which form the hydrate lattice and how they impact the local structure of the lattice itself. Results indicate that when CH 4 and CO 2 co-occupy the hydrate, the host is more strongly distorted than in the pure CH 4 and CO 2 hydrates, but this becomes less defined with increasing temperature. The presence of CO 2 in mixed hydrate increases the stability range and creates a barrier for CH 4 to completely leave the structure.

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

Electromagnetic Transient (EMT) Simulation Algorithms for Evaluation of Large-Scale Extreme Fast Charging Systems (T&D Models)

Simulation of high-fidelity models of extreme fast charging (XFC) systems and large-area power grids with many XFCs can be time consuming in traditional simulators. Traditional simulators use a single method of discretization for all the components that results in imposing a large computational burden of inverting a large matrix as well as increased computations related to single method of discretization (that is typically a trapezoidal method). To overcome the problem of simulating large-area power grids with many XFCs, in this paper, advanced numerical simulation algorithms are applied for the first time together to reduce the dimension of matrix inversion. Here, the algorithms include numerical stiffness-based segregation, time constant-based segregation, clustering and aggregation on differential algebraic equations (DAEs), and multi-order integration approaches. These algorithms apply multiple discretization algorithms rather than a single discretization algorithm that further reduces the computational burden. The approaches mentioned here have resulted in speed-up of up to 18x in the simulation of a single distribution system with 15 XFCs and of up to 271x in the simulation of a transmission-distribution system with 300 XFCs in multiple distribution feeders with respect to conventional simulators (like power systems computer aided design [PSCAD]).

42 ENGINEERING↗
DOE OSTI2015

Multidimensional dosimetry of {sup 106}Ru eye plaques using EBT3 films and its impact on treatment planning

Purpose: The purpose of this study was to establish a method to perform multidimensional radiochromic film measurements of {sup 106}Ru plaques and to benchmark the resulting dose distributions against Monte Carlo simulations (MC), microdiamond, and diode measurements. Methods: Absolute dose rates and relative dose distributions in multiple planes were determined for three different plaque models (CCB, CCA, and COB), and three different plaques per model, using EBT3 films in an in-house developed polystyrene phantom and the MCNP6 MC code. Dose difference maps were generated to analyze interplaque variations for a specific type, and for comparing measurements against MC simulations. Furthermore, dose distributions were validated against values specified by the manufacturer (BEBIG) and microdiamond and diode measurements in a water scanning phantom. Radial profiles were assessed and used to estimate dosimetric margins for a given combination of representative tumor geometry and plaque size. Results: Absolute dose rates at a reference depth of 2 mm on the central axis of the plaque show an agreement better than 5% (10%) when comparing film measurements (MCNP6) to the manufacturer’s data. The reproducibility of depth-dose profile measurements was <7% (2 SD) for all investigated detectors and plaque types. Dose difference maps revealed minor interplaque deviations for a specific plaque type due to inhomogeneities of the active layer. The evaluation of dosimetric margins showed that for a majority of the investigated cases, the tumor was not completely covered by the 100% isodose prescribed to the tumor apex if the difference between geometrical plaque size and tumor base ≤4 mm. Conclusions: EBT3 film dosimetry in an in-house developed phantom was successfully used to characterize the dosimetric properties of different {sup 106}Ru plaque models. The film measurements were validated against MC calculations and other experimental methods and showed a good agreement with data from BEBIG well within published tolerances. The dosimetric information as well as interplaque comparison can be used for comprehensive quality assurance and for considerations in the treatment planning of ophthalmic brachytherapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Comprehensive study of the interactions between the critical dimensionless numbers associated with multiphase flow in 3D porous media

Multiphase flow in porous media is of signifcant interest in many engineering applications, such as geologic carbon sequestration, enhanced oil recovery, and groundwater contamination and remediation. In order to advance the fundamental understanding of multiphase flow in complex three-dimensional (3D) porous media, the interactions between the critical dimensionless numbers, including the contact angle, viscosity ratio, and capillary (Ca) number, were explored using X-ray micro-computed tomography (micro-CT) scanning and lattice Boltzmann (LB) modeling. In this study, the 3D pore structure information was extracted from micro-CT images and then used as interior boundary conditions of flow modeling in a pore-scale LB simulator to simulate multiphase flow within the pore space. A Berea sandstone sample was scanned and then two-phase flow LB simulations were performed based on the micro-CT images. The LB-simulated water/CO 2 distributions agreed well with the micro-CT scanned images. Simulation findings showed that a decreasing contact angle causes a decrease in wetting-fluid relative permeability and an increase in non-wetting fluid relative permeability. A rising Ca number increases both wetting and non-wetting fluid relative permeabilities. An increasing viscosity ratio (the ratio of non-wetting fluid viscosity to wetting fluid viscosity) facilitates the increase of non-wetting fluid relative permeability and mitigates the reduction of wetting fluid relative permeability, when the contact angle decreases continuously. The primary novel finding of this study is that the viscosity ratio affects the rate of change of the relative permeability curves for both phases when the contact angle changes continuously. To date, we believe it is the first time that comprehensive interactions between these dimensionless numbers are demonstrated in a sandstone sample based on real 3D structures. We also investigated the role of the changes of flow direction and sample location on relative permeability curves. Simulation results showed that the change in non-wetting fluid relative permeability was larger when the flow direction was switched from vertical to horizontal, which indicated that there was stronger anisotropy in larger pore networks that were primarily occupied by the non-wetting fluid.

42 ENGINEERING↗
DOE OSTI2016

ERRATIC FLARING OF BL LAC IN 2012–2013: MULTIWAVELENGTH OBSERVATIONS

BL Lac, the eponymous blazar, flared to historically high levels at millimeter, infrared, X-ray, and gamma-ray wavelengths in 2012. We present observations made with Herschel, Swift, NuSTAR, Fermi, the Submillimeter Array, CARMA, and the VLBA in 2012–2013, including three months with nearly daily sampling at several wavebands. We have also conducted an intensive campaign of 30 hr with every-orbit observations by Swift and NuSTAR, accompanied by Herschel, and Fermi observations. The source was highly variable at all bands. Time lags, correlations between bands, and the changing shapes of the spectral energy distributions can be explained by synchrotron radiation and inverse Compton emission from nonthermal seed photons originating from within the jet. The passage of four new superluminal very long baseline interferometry knots through the core and two stationary knots about 4 pc downstream accompanied the high flaring in 2012–2013. The seed photons for inverse Compton scattering may arise from the stationary knots and from a Mach disk near the core where relatively slow-moving plasma generates intense nonthermal radiation. The 95 spectral energy distributions obtained on consecutive days form the most densely sampled, broad wavelength coverage for any blazar. The observed spectral energy distributions and multi-waveband light curves are similar to simulated spectral energy distributions and light curves generated with a model in which turbulent plasma crosses a conical shock with a Mach disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

FACET Emittance Growth

FACET, the Facility for Advanced Accelerator and Experimental Tests, is a new facility being constructed in sector 20 of the SLAC linac primarily to study beam driven plasma wakefield acceleration. The FACET beamline consists of a chicane and final focus system to compress the 23 GeV, 3.2 nC electron bunches to {approx}20 {micro}m long and {approx}10 {micro}m wide. Simulations of the FACET beamline indicate the short-duration and large, 1.5% rms energy spread beams may suffer a factor of four emittance growth from a combination of chromaticity, incoherent synchrotron radiation (ISR), and coherent synchrotron radiation (CSR). Emittance growth is directly correlated to head erosion in plasma wakefield acceleration and is a limiting factor in single stage performance. Studies of the geometric, CSR, and ISR components are presented. Numerical calculation of the rms emittance can be overwhelmed by long tails in the simulated phase space distributions; more useful definitions of emittance are given. A complete simulation of the beamline is presented as well, which agrees with design specifications.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Cosmic Reionization on Computers: Properties of the Post-reionization IGM

We present a comparison between several observational tests of the post-reionization intergalactic medium and the numerical simulations of reionization completed under the Cosmic Reionization On Computers (CROC) project. The CROC simulations match the gap distribution reasonably well, and also provide a good match for the distribution of peak heights, but there is a notable lack of wide peaks in the simulated spectra and the flux-probability distribution functions are poorly matched in the narrow redshift interval 5.5<z<5.7, with the match at other redshifts being significantly better, albeit not exact. Both discrepancies are related: simulations show more opacity than the data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Evaluating seasonal and regional distribution of snowfall in regional climate model simulations in the Arctic

In this study, we investigate how the regional climate model HIRHAM5 reproduces the spatial and temporal distribution of Arctic snowfall when compared to CloudSat satellite observations during the examined period of 2007–2010. For this purpose, both approaches, i.e., the assessments of the surface snowfall rate (observation-to-model) and the radar reflectivity factor profiles (model-to-observation), are carried out considering spatial and temporal sampling differences. The HIRHAM5 model, which is constrained in its synoptic representation by nudging to ERA-Interim, represents the snowfall in the Arctic region well in comparison to CloudSat products. The spatial distribution of the snowfall patterns is similar in both identifying the southeastern coast of Greenland and the North Atlantic corridor as regions gaining more than twice as much snowfall as the Arctic average, defined here for latitudes between 66 and 81°N. Excellent agreement (difference less than 1%) in the Arctic-averaged annual snowfall rate between HIRHAM5 and CloudSat is found, whereas ERA-Interim reanalysis shows an underestimation of 45% and significant deficits in the representation of the snowfall rate distribution. From the spatial analysis, it can be seen that the largest differences in the mean annual snowfall rates are an overestimation near the coastlines of Greenland and other regions with large orographic variations as well as an underestimation in the northern North Atlantic Ocean. To a large extent, the differences can be explained by clutter contamination, blind zone or higher resolution of CloudSat measurements, but clearly HIRHAM5 overestimates the orographic-driven precipitation. The underestimation of HIRHAM5 within the North Atlantic corridor south of Svalbard is likely connected to a poor description of the marine cold air outbreaks which could be identified by separating snowfall into different circulation weather type regimes. By simulating the radar reflectivity factor profiles from HIRHAM5 utilizing the Passive and Active Microwave TRAnsfer (PAMTRA) forward-modeling operator, the contribution of individual hydrometeor types can be assessed. Looking at a latitude band at 72–73°N, snow can be identified as the hydrometeor type dominating radar reflectivity factor values across all seasons. The largest differences between the observed and simulated reflectivity factor values are related to the contribution of cloud ice particles, which is underestimated in the model, most likely due to the small sizes of the particles. The model-to-observation approach offers a promising diagnostic when improving cloud schemes, as illustrated by comparison of different schemes available for HIRHAM5.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Network Reduction Algorithm for Developing Distribution Feeders for Real-Time Simulators

As advanced grid-support functions (AGF) become more widely used in grid-connected photovoltaic (PV) inverters, utilities are increasingly interested in their impacts when implemented in the field. These effects can be understood by modeling feeders in real-time simulators and test PV inverters using power hardware-in-the-loop (PHIL) techniques. This paper presents a novel feeder model reduction algorithm using a ruin & reconstruct methodology that enables large feeders to be solved and operated on real-time computing platforms. Two Hawaiian Electric feeder models in Synergi Electric's load flow software were converted to reduced order models in OpenDSS, and subsequently implemented in the OPAL-RT real-time digital testing platform. Smart PV inverters were added to the realtime model with AGF responses modeled after characterizing commercially available hardware inverters. Finally, hardware inverters were tested in conjunction with the real-time model using PHIL techniques so that the effects of AGFs on the feeders could be analyzed.

14 SOLAR ENERGY↗
DOE OSTI2017

Network Reduction Algorithm for Developing Distribution Feeders for Real-Time Simulators: Preprint

As advanced grid-support functions (AGF) become more widely used in grid-connected photovoltaic (PV) inverters, utilities are increasingly interested in their impacts when implemented in the field. These effects can be understood by modeling feeders in real-time systems and testing PV inverters using power hardware-in-the-loop (PHIL) techniques. This paper presents a novel feeder model reduction algorithm using a Monte Carlo method that enables large feeders to be solved and operated on real-time computing platforms. Two Hawaiian Electric feeder models in Synergi Electric's load flow software were converted to reduced order models in OpenDSS, and subsequently implemented in the OPAL-RT real-time digital testing platform. Smart PV inverters were added to the real-time model with AGF responses modeled after characterizing commercially available hardware inverters. Finally, hardware inverters were tested in conjunction with the real-time model using PHIL techniques so that the effects of AGFs on the choice feeders could be analyzed.

14 SOLAR ENERGY↗
DOE OSTI2011

SIMULATING THE EFFECTS OF INITIAL PITCH-ANGLE DISTRIBUTIONS ON SOLAR FLARES

In this work, we model both the thermal and non-thermal components of solar flares. The model we use, HYLOOP, combines a hydrodynamic equation solver with a non-thermal particle tracking code to simulate the thermal and non-thermal dynamics and emission of solar flares. In order to test the effects of pitch-angle distribution on flare dynamics and emission, a series of flares is simulated with non-thermal electron beams injected at the loop apex. The pitch-angle distribution of each beam is described by a single parameter and allowed to vary from flare to flare. We use the results of these simulations to generate synthetic hard and soft X-ray emissions (HXR and SXR). The light curves of the flares in Hinode's X-ray Telescope passbands show a distinct signal that is highly dependent on pitch-angle distribution. The simulated HXR emission in the 3-6 keV bandpass shows the formation and evolution of emission sources that correspond well to the observations of pre-impulsive flares. This ability to test theoretical models of thermal and non-thermal flare dynamics directly with observations allows for the investigation of a wide range of physical processes governing the evolution of solar flares. We find that the initial pitch-angle distribution of non-thermal particle populations has a profound effect on loop top HXR and SXR emission and that apparent motion of HXR is a natural consequence of non-thermal particle evolution in a magnetic trap.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Pore-scale multiphase flow modeling and imaging of CO 2 exsolution in Sandstone

This study utilizes synchrotron X-ray micro-tomography and pore scale modeling to investigate the process of gas exsolution and how it affects non-wetting phase relative permeability. Exsolved gas distributions are measured on Domengine and Boise sandstone samples using synchrotron X-ray micro-tomography. Observed gas phase distributions are compared to a new model that simulates the growth and distribution of exsolved gas phase at the pore-scale. Water relative permeability curves are calculated using a Stokes flow simulator with modeled and observed gas distributions, under various conditions, such as rock geometry, and pressure depletion rates. By comparing the actual bubble distributions with modeled distributions, we conclude that exsolved gas is more likely to form and accumulate at locations with higher water velocities. This suggests that convective delivery of CO2 to the gas bubble is a primary mechanism for bubble growth, as compared to diffusive transport through the aqueous phase. For carbonated brine flowing up a fault at half a meter per day, with 5% exsolved gas, the water relative permeability is estimated to be 0.6~0.8 for various sandstones. The reduction of water mobility reduces upward brine migration when even a small amount of exsolution occurs.

58 GEOSCIENCES↗
DOE OSTI2012

Structure and Transformation of Amorphous Calcium Carbonate: A Solid-State 43Ca NMR and Computational Molecular Dynamics Investigation

Amorphous calcium carbonate (ACC) is a metastable precursor to crystalline CaCO{sub 3} phases that precipitates by aggregation of ion pairs and prenucleation clusters. We use {sup 43}Ca solid-state NMR spectroscopy to probe the local structure and transformation of ACC synthesized from seawater-like solutions with and without Mg{sup 2+} and computational molecular dynamics (MD) simulations to provide more detailed molecular-scale understanding of the ACC structure. The {sup 43}Ca NMR spectra of ACC collected immediately after synthesis consist of broad, featureless resonances with Gaussian line shapes (FWHH = 27.6 {+-} 1 ppm) that do not depend on Mg{sup 2+} or H{sub 2}O content. A correlation between {sup 43}Ca isotropic chemical shifts and mean Ca-O bond distances for crystalline hydrous and anhydrous calcium carbonate phases indicates indistinguishable maximum mean Ca-O bond lengths of {approx}2.45 {angstrom} for all our samples. This value is near the upper end of the published Ca-O bond distance range for biogenic and synthetic ACCs obtained by Ca-X-ray absorption spectroscopy. It is slightly smaller than the values from the structural model of Mgfree ACC by Goodwin et al. obtained from reverse Monte Carlo (RMC) modeling of X-ray scattering data and our own computational molecular dynamics (MD) simulation based on this model. An MD simulation starting with the atomic positions of the Goodwin et al. RMC model using the force field of Raiteri and Gale shows significant structural reorganization during the simulation and that the interconnected carbonate/water-rich channels in the Goodwin et al. model shrink in size over the 2 ns simulation time. The distribution of polyhedrally averaged Ca-O bond distances from the MD simulation is in good agreement with the {sup 43}Ca NMR peak shape, suggesting that local structural disorder dominates the experimental line width of ACC.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Data-Driven Distribution System Load Modeling for Quasi-Static Time-Series Simulation

This paper presents a data-driven distribution system load modeling methodology targeting quasi-static time-series (QSTS) simulation. The proposed methodology is appropriate for modeling down to the level of the customer transformer, and it has three main features: 1) both load pattern diversity and intra-second load variability are considered, 2) the load profiles can be populated for multiple nodes on a circuit in such a way that the diversity factor of the feeder can be defined and tuned, and 3) the load aggregation method can be used to populate the profiles for different nodes at various load aggregation levels. Furthermore, as the foundation of the modeling methodology, variability and diversity libraries have been established based on high-resolution load data collected on customer transformers from real utility feeders. The proposed modeling methodology has been used to build load data sets for both the IEEE 123-bus feeder model and a realistic utility feeder model. The QSTS simulation results on the two evaluation feeders have demonstrated that the load data sets established from the proposed modeling methodology can effectively capture the load impact on feeder operations. For the realistic utility feeder, the effectiveness of the proposed methodology has also been validated by comparing the voltage characteristics of the feeder with modeled loads and the voltage characteristics of realistic voltage data from the same feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Bulk cloud microphysical properties as seen from numerical simulation and remote sensing products: case study of a hailstorm event over the La Plata Basin

Hailstorms develop over the La Plata Basin, in south-eastern South America, more often during later winter and early austral spring, between September and October. These systems have significant socioeconomic impacts over the region. Thus, a better understanding of how atmospheric drivers modulate the formation of hailstorms is important to improve the forecast of such phenomena. In this study, we selected a hailstorm event observed over the eastern La Plata Basin during 14–15 July 2016 to evaluate the performance of the Brazilian developments on the Regional Atmospheric Modelling System (BRAMS) model. The ability of the model in simulating cloud microphysical properties was evaluated by comparing simulations driven by different global forcings against in situ and remote sensing observations. The model results showed good skill in capturing the basic characteristics of the thunderstorm, particularly in terms of the spatial distribution of hydrometeors. The simulated spatial distribution of hail covers locations where hail fall was reported. The BRAMS simulations suggest that, despite relatively low values of the convective available potential energy (CAPE) (700–1000 J kg -1 ), environments with strong 0–8-km bulk shear (60–70 kt, ~30.9–36.0 m s –1 ) can promote the formation of ice clouds and hail fall over the eastern La Plata Basin. To be more conclusive, however, further research is needed to understand how different combinations of CAPE and shear affect hail formation over the region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2010

Understanding Spatial and Spectral Morphologies of Ultracompact H II Regions

The spatial morphology, spectral characteristics, and time variability of ultracompact H II regions provide strong constraints on the process of massive star formation. We have performed simulations of the gravitational collapse of rotating molecular cloud cores, including treatments of the propagation of ionizing and non-ionizing radiation. We here present synthetic radio continuum observations of H II regions from our collapse simulations, to investigate how well they agree with observation, and what we can learn about how massive star formation proceeds. We find that intermittent shielding by dense filaments in the gravitationally unstable accretion flow around the massive star leads to highly variable H II regions that do not grow monotonically, but rather flicker, growing and shrinking repeatedly. This behavior appears able to resolve the well-known lifetime problem. We find that multiple ionizing sources generally form, resulting in groups of ultracompact H II regions, consistent with observations. We confirm that our model reproduces the qualitative H II region morphologies found in surveys, with generally consistent relative frequencies. We also find that simulated spectral energy distributions (SEDs) from our model are consistent with the range of observed H II region SEDs, including both regions showing a normal transition from optically thick to optically thin emission, and those with intermediate spectral slopes. In our models, anomalous slopes are solely produced by inhomogeneities in the H II region, with no contribution from dust emission at millimeter or submillimeter wavelengths. We conclude that many observed characteristics of ultracompact H II regions appear consistent with massive star formation in fast, gravitationally unstable, accretion flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Elucidation of Iron Gettering Mechanisms in Boron-Implanted Silicon Solar Cells

To facilitate cost-effective manufacturing of boron-implanted silicon solar cells as an alternative to BBr 3 diffusion, we performed a quantitative test of the gettering induced by solar-typical boron-implants with the potential for low saturation current density emitters (< 50 fA/cm 2 ). We show that depending on the contamination level and the gettering anneal chosen, such boron-implanted emitters can induce more than a 99.9% reduction in bulk iron point defect concentration. The iron point defect results as well as synchrotron-based Nano-X-ray-fluorescence investigations of iron precipitates formed in the implanted layer imply that, with the chosen experimental parameters, iron precipitation is the dominant gettering mechanism, with segregation-based gettering playing a smaller role. We reproduce the measured iron point defect and precipitate distributions via kinetics modeling. First, we simulate the structural defect distribution created by the implantation process, and then we model these structural defects as heterogeneous precipitation sites for iron. Unlike previous theoretical work on gettering via boron- or phosphorus-implantation, our model is free of adjustable simulation parameters. The close agreement between the model and experimental results indicates that the model successfully captures the necessary physics to describe the iron gettering mechanisms operating in boron-implanted silicon. Furthermore, this modeling capability allows high-performance, cost-effective implanted silicon solar cells to be designed.

14 SOLAR ENERGY↗