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 1,171 records · Page 65

Control of simulation of a particle accelerator as a spatially distributed system

The work deals with the construction of strategy for finding optimal control of systems for spatially distributed systems on an example of problems in physics beams. The suggested approach is based on use of exact and approximate invariants as a controlling mechanism for control of modelling process of beam lines. The classical problem of optimal control of similar systems as an optimal control problem is formulated in a functional space. This approach is based on algebraic methods Lie and the corresponding matrix formalism that allows you to build not only the formal procedures for finding optimal solutions, but also provide the necessary constructive computational procedures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analysis of PV Advanced Inverter Functions and Setpoints under Time Series Simulation.

Utilities are increasingly concerned about the potential negative impacts distributed PV may have on the operational integrity of their distribution feeders. Some have proposed novel methods for controlling a PV system's grid - tie inverter to mitigate poten tial PV - induced problems. This report investigates the effectiveness of several of these PV advanced inverter controls on improving distribution feeder operational metrics. The controls are simulated on a large PV system interconnected at several locations within two realistic distribution feeder models. Due to the time - domain nature of the advanced inverter controls, quasi - static time series simulations are performed under one week of representative variable irradiance and load data for each feeder. A para metric study is performed on each control type to determine how well certain measurable network metrics improve as a function of the control parameters. This methodology is used to determine appropriate advanced inverter settings for each location on the f eeder and overall for any interconnection location on the feeder.

14 SOLAR ENERGY↗

Comparing event generator predictions and ab initio calculations of ν- 12 C neutral-current quasielastic scattering at 1 GeV

The measurement of neutrino oscillations and exotic physics searches are important parts of the physics program in the near future, with new state-of-the-art experiments planned within the next decade. Future and modern experiments in these fields will make use of nuclear targets. Event generators (EGs) are software used in the analysis of neutrino oscillation experiments. EGs are used to predict kinematic observables for a range of neutrino energies. These simulations make use of simple models of nucleon dynamics. As such, they may fail to capture features of more rigorous theoretical calculations. This work compares EG performance to nuclear theory calculations by comparing observables generated in the two frameworks. We provide a common set of definitions between theory and experiment and assess the adequacy of the implementation of two nuclear physics models in EG simulations. Neutral-current quasielastic (NCQE) scattering events for neutrinos and antineutrinos on a 12 C target are simulated with a specific EG, NEUT, used by the T2K experiment for its analysis. The simulated cross sections are compared to analytic calculations from nuclear theory within the factorization scheme. We compare the NEUT implementation of two different models on nuclear spectral functions: the relativistic Fermi gas (RFG) and the correlated basis spectral function (CBF) to analytic calculations of the same models in the factorization scheme. For both nuclear physics models, we compare the appearance of features in the distributions relevant to experimental analyses. The peak of the cross section dσ/(dΩ dω) is consistent in energy transfer, ω, for RFG and CBF simulations. Qualitatively, the shape of the simulated distribution is similar to the one obtained through theory calculations; however, there are some discrepancies between the theory calculations and the NEUT simulation. While the EG simulations and analytic calculations with the same model of nuclear dynamics show similar overall features, there are still differences between the two. These results demonstrate the importance of benchmarking EGs so their nuclear physics implementations can be improved for the analysis of future experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Determination of Dark Matter Halo Mass from Dynamics of Satellite Galaxies

We propose that the mass of a dark matter halo can be inferred from the dynamical status of its satellite galaxies. Using nine dark matter simulations of halos like the Milky Way (MW), we find that the present-day substructures in each halo follow a characteristic distribution in the phase space of orbital binding energy and angular momentum, and that this distribution is similar from halo to halo, but has an intrinsic dependence on the halo formation history. We build this distribution directly from the simulations for a specific halo and extend the result to halos of similar formation history but different masses by scaling. The mass of an observed halo can then be estimated by maximizing the likelihood in comparing the measured kinematic parameters of its satellite galaxies with these distributions. We test the validity and accuracy of this method with mock samples taken from the simulations. Using the positions, radial velocities, and proper motions of nine tracers and assuming observational uncertainties comparable to those of MW satellite galaxies, we find that the halo mass can be recovered to within $\sim 40 \% $. The accuracy can be improved to within ~25% if 30 tracers are used. Yet, the dependence of the phase-space distribution on the halo formation history sets a minimum uncertainty of $\sim 20 \% $ that cannot be reduced by using more tracers. We believe that this minimum uncertainty also applies to any mass determination for a halo when the phase-space information of other kinematic tracers is used.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamic Building Load Control to Facilitate High Penetration of Solar Photovoltaic Generation (Final Technical Report)

Solar photovoltaic (PV) resources are the most common form of distributed generation in residential and commercial customer premises within electric distribution networks. A higher penetration of PV generation in distribution circuits will impose challenges on maintaining service voltages within the range of industry standards, power quality, and power flow. Buildings consume 74% of the electricity produced in the United States, and a significant portion of the building load is dispatchable, making them responsive to electrical grid needs. Oak Ridge National Laboratory—in collaboration with Southern Company; the University of Tennessee, Knoxville; and the Georgia Institute of Technology—is examining the PV integration issues in distribution-level electrical grids and developing integrated demand-side control and communication systems to enable responsive loads. The proposed responsive loads mechanism performs renewable generation following to increase the penetration of solar PV within each feeder. The specific objectives of this project are to (1) examine distribution-level PV integration scenarios to understand requirements, (2) undertake an end-to-end simulation-based design of a distributed control strategy of loads geographically near the PV generation asset to minimize the effect on the distribution feeder, (3) deploy and demonstrate the control technology developed in partnership with utilities, and (4) perform a scalability analysis at the utility scale. This 3-year integrated project aims to develop, demonstrate, and validate demand-side control technology to enable increased the penetration of renewables while mitigating challenges that arise due to their intermittency. Activities in Budget Period (BP) 1 focused on a literature review and the formal design of a control system for integrating local distribution with generation and loads. The team used modeling and simulation to evaluate the impact of varying buildings loads, variable PV generation, and power flow dynamics on the distribution circuit. The dynamic models developed in BP 1 were used in BP 2 to develop a model-based control design and a test bed. The test bed has enabled the simulation-based testing and comparison of different control designs and formulations applied to different configurations of the distribution grid, PVs, and building loads. The control approaches developed in BP 2 were implemented in BP 3 in the form of hardware deployed at the Central Baptist Church (CBC) in Knoxville, Tennessee, for testing and evaluation. The outcome of this project was the development and demonstration of open-source, low-cost, low-touch sensing and control retrofits to distributed PV generation and building loads that, in a coordinated fashion, provide the load-shaping response needed to integrate high levels of renewable penetration. This research addresses the target metrics by dynamically controlling a load with solar generation variability to minimize the extent of two-way power flow, enhance reliability, facilitate high PV penetration (>100% of peak load in a line segment), and generate scalable software and hardware solutions adaptable to any penetration levels. The research and development activities are focused and designed to be impactful within the relevant 2020 targets time frame.An accurate open-source integration simulation framework for end-to-end control design was developed and deployed at the CBC facility for testing and evaluation. This final report provides a detailed review of the technical results achieved during this 3-year integrated project. A novel spectral analysis of PV data is demonstrated to derive the requirements of the control design. A detailed simulation-based analysis of PV integration at increasing penetration levels is presented using 1 year of PV data to demonstrate the impact on the distribution circuits. Two different control strategies were developed and demonstrated via simulation to track variable PV generation with adaptive load dispatch. The report concludes with a summary of accomplishments and recommendations for a path forward.

14 SOLAR ENERGY↗

Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics

Langevin dynamics has become a popular tool to simulate the Boltzmann equilibrium distribution. When the repartition of the Langevin equation involves the exact realization of the Ornstein-Uhlenbeck noise, in addition to the conventional density evolution, there exists another type of discrete evolution that may not correspond to a continuous, real dynamical counterpart. This virtual dynamics case is also able to produce the desired stationary distribution. Different types of repartition lead to different numerical schemes, of which the accuracy and efficiency are investigated through studying the harmonic oscillator potential, an analytical solvable model. By analyzing the asymptotic distribution and characteristic correlation time that are derived by either directly solving the discrete equations of motion or using the related phase space propagators, it is shown that the optimal friction coefficient resulting in the minimum characteristic correlation time depends on the time interval chosen in the numerical implementation. In conclusion, when the recommended “middle” scheme is employed, both analytical and numerical results demonstrate that, for good numerical performance in efficiency as well as accuracy, one may choose a friction coefficient in a wide range from around the optimal value to the high friction limit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparative Analysis of the Powers of the Two-Sample Kolmogorov–Smirnov and Anderson–Darling Tests Under Various Alternatives

In this paper we conduct a comparative analysis of the powers of the two-sample Kolmogorov–Smirnov and Anderson–Darling tests under various alternatives using simulation. We consider two examples. In the first example the alternatives to the standard normal distribution are the distributions of the so-called contaminated normal model. We study the influence of a small contamination with a positive shift on the powers of the test. In the second example the alternatives are the logistic and the Laplace distributions, which are symmetric and differ in shape from the normal distribution having a larger kurtosis coefficient and heavier tails.

97 MATHEMATICS AND COMPUTING↗

A Mixed integer linear programming‐based distributed energy management for networked microgrids considering network operational objectives and constraints

Abstract Mixed integer linear programming (MILP)–based distributed energy management for networked microgrids embedded modern distribution systems is proposed. Considering the diverse ownership of microgrids, distributed energy resources (DERs) that interface directly with utilities and responsive loads, an alternating direction method of multipliers–based distributed framework was formulated for the scheduling of networked microgrids embedded modern distribution systems by adjusting nodal price signals iteratively. In addition, to make the formulated optimization problems resolvable through more accessible and popular MILP solvers, different linearisation techniques were employed to transform the nonlinear terms into linear or mixed integer linear formats. The proposed MILP‐based distributed method preserves all participants' autonomy (e.g., microgrids, DERs that interface directly with utilities and responsive loads), while incentivising them to actively participate in the distribution system operation with price signals. The proposed method is validated with results of numerical simulation using a modern distribution system consisting of multiple networked microgrids, DERs that interface directly with utilities, as well as responsive loads.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Effects of Magnetic and Kinetic Helicities on the Growth of Magnetic Fields in Laminar and Turbulent Flows by Helical Fourier Decomposition

We present a numerical and analytical study of incompressible homogeneous conducting fluids using a helical Fourier representation. We analytically study both small- and large-scale dynamo properties, as well as the inverse cascade of magnetic helicity, in the most general minimal subset of interacting velocity and magnetic fields on a closed Fourier triad. We mainly focus on the dependency of magnetic field growth as a function of the distribution of kinetic and magnetic helicities among the three interacting wavenumbers. By combining direct numerical simulations of the full magnetohydrodynamics equations with the helical Fourier decomposition, we numerically confirm that in the kinematic dynamo regime the system develops a large-scale magnetic helicity with opposite sign compared to the small-scale kinetic helicity, a sort of triad-by-triad α -effect in Fourier space. Concerning the small-scale perturbations, we predict theoretically and confirm numerically that the largest instability is achived for the magnetic component with the same helicity of the flow, in agreement with the Stretch–Twist–Fold mechanism. Vice versa, in the presence of Lorentz feedback on the velocity, we find that the inverse cascade of magnetic helicity is mostly local if magnetic and kinetic helicities have opposite signs, while it is more nonlocal and more intense if they have the same sign, as predicted by the analytical approach. Our analytical and numerical results further demonstrate the potential of the helical Fourier decomposition to elucidate the entangled dynamics of magnetic and kinetic helicities both in fully developed turbulence and in laminar flows.

79 ASTRONOMY AND ASTROPHYSICS↗

Topology of Large-scale Structures of Galaxies in Two Dimensions—Systematic Effects

We study the two-dimensional topology of galactic distribution when projected onto two-dimensional spherical shells. Using the latest Horizon Run 4 simulation data, we construct the genus of the two-dimensional field and consider how this statistic is affected by late-time nonlinear effects—principally gravitational collapse and redshift space distortion (RSD). We also consider systematic and numerical artifacts, such as shot noise, galaxy bias, and finite pixel effects. We model the systematics using a Hermite polynomial expansion and perform a comprehensive analysis of known effects on the two-dimensional genus, with a view toward using the statistic for cosmological parameter estimation. We find that the finite pixel effect is dominated by an amplitude drop and can be made less than 1% by adopting pixels smaller than 1/3 of the angular smoothing length. Nonlinear gravitational evolution introduces time-dependent coefficients of the zeroth, first, and second Hermite polynomials, but the genus amplitude changes by less than 1% between z = 1 and z = 0 for smoothing scales R{sub G}>9 Mpc/h. Non-zero terms are measured up to third order in the Hermite polynomial expansion when studying RSD. Differences in the shapes of the genus curves in real and redshift space are small when we adopt thick redshift shells, but the amplitude change remains a significant ∼O(10%) effect. The combined effects of galaxy biasing and shot noise produce systematic effects up to the second Hermite polynomial. It is shown that, when sampling, the use of galaxy mass cuts significantly reduces the effect of shot noise relative to random sampling.

79 ASTRONOMY AND ASTROPHYSICS↗

Heat flux modeling using ion drift effects in DIII-D H-mode plasmas with resonant magnetic perturbations

The heat flux patterns measured in low-collisionality DIII-D H-mode plasmas strongly deviate from simultaneously measured CII emission patterns, used as indicator of particle flux, during applied resonant magnetic perturbations. While the CII emission clearly shows typical striations, which are similar to magnetic footprint patterns obtained from vacuum field line tracing, the heat flux is usually dominated by one large peak at the strike point position. The vacuum approximation, which only considers applied magnetic fields and neglects plasma response and plasma effects, cannot explain the shape of the observed heat flux pattern. One possible explanation is the effect of particle drifts. This is included in the field line equations and the results are discussed with reference to the measurement. Electrons and ions show different drift motions at thermal energy levels in a guiding center approximation. While electrons hardly deviate from the field lines, ions can drift several centimetres away from field line flux surfaces. A model is presented in which an ion heat flux, based on the ion drift motion from various kinetic energies as they contribute to a thermal Maxwellian distribution, is calculated. The simulated heat flux is directly compared to measurements with a varying edge safety factor q{sub 95}. This analysis provides evidence for the dominate effect of high-energy ions in carrying heat from the plasma inside the separatrix to the target. High-energy ions are deposited close to the unperturbed strike line, while low-energy ions can travel into the striated magnetic topology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Metal–Carbodithioate-Based 3D Semiconducting Metal–Organic Framework: Porous Optoelectronic Material for Energy Conversion

Solar energy conversion requires the working compositions to generate photoinduced charges with high potential and the ability to deliver charges to the catalytic sites and/or external electrode. These two properties are typically at odds with each other and call for new molecular materials with sufficient conjugation to improve charge conductivity but not as much conjugation as to overly compromise the optical band gap. In this work, we developed a semiconducting metal-organic framework (MOF) prepared explicitly through metal-carbodithioate"(-CS 2 ) n M" linkage chemistry, entailing augmented metal-linker electronic communication. The stronger ligand field and higher covalent character of metal-carbodithioate linkages-when combined with spirofluorene-derived organic struts and nickel-(II) ion-based nodes-provided a stable, semiconducting 3D-porous MOF, Spiro-CS 2 Ni. This MOF lacks long-range ordering and is defined by a flexible structure with non-aggregated building units, as suggested by reverse Monte Carlo simulations of the pair distribution function obtained from total scattering experiments. The solvent-removed "closed pore"material recorded a Brunauer-Emmett-Teller area of similar to 400m 2 /g, where the "open pore" form possesses90 wt % solvent-accessible porosity. Electrochemical measurements suggest that Spiro-CS 2 Ni possesses a band gap of 1.57 eV(sigma = 10 -7 S/cm at -1.3 V bias potential),which can be further improved by manipulating the d-electron configuration through an axial coordination (ligand/substrate), the latter of which indicates usefulness as an electrocatalyst and/or a photoelectrocatalyst(upon substrate binding). Transient-absorption spectroscopy reveals a long-lived photo-generated charge-transfer state (τ CR = 6.5 mu s) capable of chemical transformation under a biased voltage. Spiro-CS 2 Ni can endure a compelling range of pH(1-12 for weeks) and hours of electrochemical and photoelectrochemical conditions in the presence of water and organic acids. We believe this work provides crucial design principles for low-density, porous,light-energy-conversion materials.

36 MATERIALS SCIENCE↗

The BayREN Integrated Commercial Retrofits (BRICR) Project: An Introduction and Preliminary Results

BayREN Integrated Commercial Retrofits (BRICR) is a DOE-funded project which aims to enhance the capacity of energy efficiency programs to recruit participants, develop retrofits, and measure outcomes in small and medium-sized commercial buildings, a sector notoriously hard to reach and expensive to serve, that accounts for 2/3 of US commercial floor space. To address these barriers, BRICR leverages existing incentives, financing, data, and open source software to facilitate two paths for comprehensive improvements: a deep energy retrofit, or serial upgrades integrated into capital improvement and maintenance cycles. BRICR is developing an integrated workflow for iterative energy modeling of commercial buildings for city energy program managers and auditors - starting with mass building-scale simulation based on public records and proceeding through audit, retrofit, and measurement and verification stages. BRICR builds on existing tools including LBNL’s CBES, NREL’s OpenStudio, PNNL’s Audit Template, and DOE’s BuildingSync and SEED PlatformTM. At each stage, BRICR uses available information to inform simulations (starting with public records but iteratively augmented with observations from energy program staff) to improve the quality of the models that inform decision making. This paper presents initial results from energy models of 1699 office and retail buildings in San Francisco. Building stock data from public records were translated to create a BuildingSync file for each building which was stored in SEED. Each BuildingSync file was then translated to multiple OpenStudio Workflow files for EnergyPlus simulation to estimate energy savings of energy conservation measures (ECM). Energy savings predictions for each ECM were written back to an updated BuildingSync file for each building and re-uploaded to SEED. The distribution of baseline energy simulations was calibrated against the publicly disclosed distribution of energy benchmarking data to increase confidence in results.

Piette, M↗

An Islanding Detection Test Platform for Multi-Inverter Islands Using Power HIL

When an unintentional island is formed on the electric power system, distributed energy resources (DERs) are typically required to detect and de-energize the island. This requirement may become more challenging as the number of DERs in an island rises. Thus, it is of interest to experimentally verify whether DERs can successfully detect and de-energize islands containing many DERs connected at different points. This paper presents a power hardware-in-the-loop (PHIL) platform for testing the duration of islands containing multiple inverters connected at multiple points in a network. The PHIL platform uses real-time simulation to represent islanded distribution circuits, with DER inverters connected in hardware. This allows efficient testing of a large number of island configurations simply by changing the distribution circuit model in the real-time simulator. A method for calculating the quality factor of an arbitrary distribution circuit in real time, designed for anti-islanding tests, is also presented. Experimental results are included demonstrating the use of the PHIL method to test a variety of three-inverter islands.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Visualisation of structural inhomogeneities in strongly scattering media using the method of spatially-resolved reflectometry: Monte Carlo simulation

Two-dimensional spatial intensity distributions of diffuse scattering of near-infrared laser radiation from a strongly scattering medium, whose optical properties are close to those of skin, are obtained using Monte Carlo simulation. The medium contains a cylindrical inhomogeneity with the optical properties, close to those of blood. It is shown that stronger absorption and scattering of light by blood compared to the surrounding medium leads to the fact that the intensity of radiation diffusely reflected from the surface of the medium under study and registered at its surface has a local minimum directly above the cylindrical inhomogeneity. This specific feature makes the method of spatially-resolved reflectometry potentially applicable for imaging blood vessels and determining their sizes. It is also shown that blurring of the vessel image increases almost linearly with increasing vessel embedment depth. This relation may be used to determine the depth of embedment provided that the optical properties of the scattering media are known. The optimal position of the sources and detectors of radiation, providing the best imaging of the vessel under study, is determined. (biophotonics)

60 APPLIED LIFE SCIENCES↗

Evaluating Transactive Controls of Integrated Transmission and Distribution Systems using the Framework for Network Co-Simulation

With an ever-evolving power grid, concerns regarding how to maintain system stability, efficiency, and reliability remain constant because of increasing uncertainties and decreasing rotating inertia. To alleviate some of these concerns, demand response represents a viable solution and is virtually an untapped resource in the current power grid. This work describes a hierarchical control framework that allows coordination between distributed energy resources and demand response. This control framework is composed of two control layers: a coordination layer that ensures aggregations of resources are coordinated to achieve system objectives and a device layer that controls individual resources to assure the predetermined power profile is tracked in real time. Large-scale simulations are executed to study the hierarchical control, requiring advancements in simulation capabilities. Technical advancements necessary to investigate and answer control interaction questions, including the Framework for Network Co-Simulation platform and Arion modeling capability, are detailed. Insights into the interdependencies of controls across a complex system and how they must be tuned, as well as validation of the effectiveness of the proposed control framework, are yielded using a large-scale integrated transmission system model coupled with multiple distribution systems.

demand response↗

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗

Coherent control of the vibrational state population in a nonpolar molecule

A coherent control scheme for the population distribution in the vibrational states of nonpolar molecules is proposed. Our theoretical analysis and results of numerical simulations for the interaction of the hydrogen molecular ion in its electronic ground state with an infrared laser pulse reveal a selective two-photon transition between the vibrational states via a coupling with the first excited dissociative state. We demonstrate that for a given temporal intensity profile the population transfer between vibrational states, or a superposition of vibrational states, can be made complete for a single chirped pulse or a train of chirped pulses, which accounts for the accumulated phase difference due to the ac Stark effect. Effects of a spatial intensity (or focal) averaging are discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗