Dynamic Programming Approach to Estimate the Capacity Value of Energy Storage
Abstract not provided
SEARCH · Search NASA
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.
Abstract not provided
We present the largest search to date for Y-band dropout galaxies (z ∼ 8 Lyman break galaxies, LBGs) based on 350 arcmin{sup 2} of Hubble Space Telescope observations in the V, Y, J, and H bands from the Brightest of Reionizing Galaxies (BoRG) survey. In addition to previously published data, the BoRG13 data set presented here includes approximately 50 arcmin{sup 2} of new data and deeper observations of two previous BoRG pointings, from which we present 9 new z ∼ 8 LBG candidates, bringing the total number of BoRG Y-band dropouts to 38 with 25.5 ≤ m{sub J} ≤ 27.6 (AB system). We introduce a new Bayesian formalism for estimating the galaxy luminosity function, which does not require binning (and thus smearing) of the data and includes a likelihood based on the formally correct binomial distribution as opposed to the often-used approximate Poisson distribution. We demonstrate the utility of the new method on a sample of 97 Y-band dropouts that combines the bright BoRG galaxies with the fainter sources published in Bouwens et al. from the Hubble Ultra Deep Field and Early Release Science programs. We show that the z ∼ 8 luminosity function is well described by a Schechter function over its full dynamic range with a characteristic magnitude M{sup ⋆}=−20.15{sub −0.38}{sup +0.29}, a faint-end slope of α=−1.87{sub −0.26}{sup +0.26}, and a number density of log{sub 10} ϕ{sup ⋆}[Mpc{sup −3}]=−3.24{sub −0.24}{sup +0.25}. Integrated down to M = –17.7, this luminosity function yields a luminosity density log{sub 10} ϵ[erg s{sup −1} Hz{sup −1} Mpc{sup −3}]=25.52{sub −0.05}{sup +0.05}. Our luminosity function analysis is consistent with previously published determinations within 1σ. The error analysis suggests that uncertainties on the faint-end slope are still too large to draw a firm conclusion about its evolution with redshift. We use our statistical framework to discuss the implication of our study for the physics of reionization. By assuming theoretically motivated priors on the clumping factor and the photon escape fraction we show that the UV luminosity density from galaxy samples down to M = –17.7 can ionize only 10%-50% of the neutral hydrogen at z ∼ 8. Full reionization would require extending the luminosity function down to M = –15. The data are consistent with a substantial fraction of neutral hydrogen at z > 7, in agreement with recent suggestions based on deep spectroscopy of z ∼ 8 LBGs.
This study presents a solution to the optimal control problem of a three degrees-of-freedom (3DOF) wave energy converter (WEC). The three modes are the heave, pitch, and surge. The dynamic model is characterized by a coupling between the pitch and surge modes, while the heave is decoupled. The heave, however, excites the pitch motion through nonlinear parametric excitation in the pitch mode. This paper uses Fourier series (FS) as basis functions to approximate the states and the control. A simplified model is first used where the parametric excitation term is neglected and a closed-form solution for the optimal control is developed. For the parametrically excited case, a sequential quadratic programming approach is implemented to solve for the optimal control numerically. Numerical results show that the harvested energy from three modes is greater than three times the harvested energy from the heave mode alone. Moreover, the harvested energy using a control that accounts for the parametric excitation is significantly higher than the energy harvested when neglecting this nonlinear parametric excitation term.
Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.
The invention provides optimal control of a three-degree-of-freedom wave energy converter using a pseudo-spectral control method. The three modes are the heave, pitch and surge. A dynamic model is characterized by a coupling between the pitch and surge modes, while the heave is decoupled. The heave, however, excites the pitch motion through nonlinear parametric excitation in the pitch mode. The invention can use a Fourier series as basis functions to approximate the states and the control. For the parametric excited case, a sequential quadratic programming approach can be implemented to numerically solve for the optimal control. The numerical results show that the harvested energy from three modes is greater than three times the harvested energy from the heave mode alone. Moreover, the harvested energy using a control that accounts for the parametric excitation is significantly higher than the energy harvested when neglecting this nonlinear parametric excitation term.
The HTV capsule is a High Flux Isotope Reactor (HFIR) target-rod capsule designed to operate at very high temperatures. The graphite containing section of the capsule (in core) is approximately 18 inches (457.2 mm) long and is separated into eight temperature zones. The specimen diameters within each zone are set to achieve the desired gas gap and hence design temperature (900 C, 1200 C or 1500 C). The capsule has five zones containing 0.400 inch (10.16 mm) diameter specimens, two zones containing 0.350 inch (8.89 mm) diameter specimens and one zone containing 0.300 inch (7.62 mm) diameter specimens. The zones have been distributed within the experiment to optimize the gamma heating from the HFIR core as well as minimize the axial heat flow in the capsule. Consequently, there are two 900 C zones, three 1200 C zones, and three 1500 C zones within the HTV capsule. Each zone contains nine specimens 0.210 0.002 inches (5.334 mm) in length. The capsule will be irradiated to a peak dose of 3.17 displacements per atom. The HTV specimens include samples of the following graphite grades: SGL Carbon s NBG-17 and NBG-18, GrafTech s PCEA, Toyo Tanso s IG-110, Mersen s 2114 and the reference grade H-451 (SGL Carbon). As part of the pre-irradiation program the specimens were characterized using ASTM Standards C559 for bulk density, and ASTM C769 for approximate Young s modulus from the sonic velocity. The probe frequency used for the determination of time of flight of the ultrasonic signal was 2.25 MHz. Marked volume (specimen diameter) effects were noted for both bulk density (increased with increasing specimen volume or diameter) and Dynamic Young s modulus (decreased with increasing specimen volume or diameter). These trends are extended by adding the property vs. diameter data for unirradiated AGC-1 creep specimens (nominally 12.5 mm-diameter x 25.4 mm-length). The relatively large reduction in Dynamic Young s Modulus was surprising given the trend for increasing density with increasing volume. The graphite-filler particle size was noted to be influential in the volume dependency data, with finer grained graphites showing the least specimen volume/diameter effect. Here the volume dependency trends are discussed in terms of the graphite s filler-particle size and texture.
We describe the moving wire interface attached to the 1-MV AMS system at LLNL’s Center for Accelerator Mass Spectrometry for the analysis of nonvolatile liquid samples as either discrete drops or from the direct output of biochemical separatory instrumentation, such as high-performance liquid chromatography (HPLC). Discrete samples containing at least a few 10 s of nanograms of carbon and as little as 50 zmol 14 C can be measured with a 3–5% precision in a few minutes. The dynamic range of our system spans approximately 3 orders in magnitude. Sample to sample memory is minimized by the use of fresh targets for each discrete sample or by minimizing the amount of carbon present in a peak generated by an HPLC containing a significant amount of 14 C. As a result, liquid sample AMS provides a new technology to expand our biomedical AMS program by enabling the capability to measure low-level biochemicals in extremely small samples that would otherwise be inaccessible.
Software solves the three-dimensional Poisson equation div(k(grad(u)) = f, by the finite element method for the case when material properties, k, are distributed over hierarchy of edges, facets and tetrahedra in the finite element mesh. Method is described in Weiss, CJ, Finite element analysis for model parameters distributed on a hierarchy of geometric simplices, Geophysics, v82, E155-167, doi:10.1190/GEO2017-0058.1 (2017). A standard finite element method for solving Poisson’s equation is augmented by including in the 3D stiffness matrix additional 2D and 1D stiffness matrices representing the contributions from material properties associated with mesh faces and edges, respectively. The resulting linear system is solved iteratively using the conjugate gradient method with Jacobi preconditioning. To minimize computer storage for program execution, the linear solver computes matrix-vector contractions element-by-element over the mesh, without explicit storage of the global stiffness matrix. Program output vtk compliant for visualization and rendering by 3rd party software. Program uses dynamic memory allocation and as such there are no hard limits on problem size outside of those imposed by the operating system and configuration on which the software is run. Dimension, N, of the finite element solution vector is constrained by the the addressable space in 32-vs-64 bit operating systems. Total storage requirements for the problem. Total working space required for the program is approximately 13*N double precision words.
The Wabash CarbonSAFE project will drill a stratigraphic test well at the Wabash Valley Resources facility near Terre Haute, Indiana. Lithologic, geomechanical, and geophysical data collected from samples testing and an extensive geophysical logging suite will be used to determine the feasibility for the geologic sequestration of 50 million metric tons or more of injected carbon dioxide. The primary target reservoir and caprock to be investigated for the Wabash CarbonSAFE project are the Mt. Simon Sandstone and Eau Claire Shale, respectively. The project team has assessed the planned Wabash #1 drilling site using subsurface geologic data and analyses from both Illinois and Indiana. These data were used to generate geologic structure and isopach thickness maps for target reservoir and caprock layers which were incorporated into a detailed discussion of the nature of the Illinois Basin’s Cambro-Ordovician carbon dioxide storage complex at the drilling site. Formation tops were prognosed for 19 lithologic units to aid drilling and coring procedures. Additional GIS analysis indicates that no known structural features will negatively impact the target storage reservoir or caprock, nor will protected and sensitive areas, groundwater resources, or existing resource development be negatively affected by Wabash #1 drilling activities. The Wabash #1 Geologic/Structure Test well drilling permit was awarded by Indiana Department of Natural Resources (IN-DNR) – Division of Oil and Gas on November 12th, 2019, and the well is expected to be drilled to a total depth of approximately 8,350 feet. Potential impacts on drilling progress include abandoned underground coal mines which mined the Springfield and Survant coals until 1964 and 1954, respectively, and a lost circulation zone in the Potosi Dolomite, which has been demonstrated to have intermittent cavernous porosity from karstification. This document also details the coring and sampling program, proposed logging suite, and well testing program, all of which will be dynamically reviewed as drilling progresses.
Here, in this paper, we propose a tri-level, nested, two-stage price-based demand response (PBDR) framework that considers distribution locational marginal price (DLMP) as DR enabler between load-serving entities (LSE), demand response providers (DRPs), and customers in the day-ahead distribution market. It enables LSE and customer interactions by using multiple DRPs, positioned in-between, and independently optimizes their objectives. The problem is formulated using linear power flow with approximated power losses and its application in DLMP as DR pricing. The tri-level problem is solved using a nested reformulation & decomposition (R&D) method and tested on the real Indian-108 bus distribution system under various dynamic pricings. Further, the temporal–spatial variations in DLMPs are assessed using fairness criteria. Numerical analyses demonstrate that DLMP applications can effectively improve economic efficiency, and transparency in DR programs valuation with a favorable fairness margin. The results show that DLMP as DR pricing signal induces (0-2) % variation in DLMP for DR participation up to 10 %. Further, it gives over 90 % fairness over temporal–spatial variation for all the customers.
This project uses dynamic rupture simulations to investigate high-frequency seismic energy generation. The relevant phenomena (frictional breakdown, shear heating, effective normal-stress fluctuations, material damage, etc.) controlling rupture are strongly interacting and span many orders of magnitude in spatial scale, requiring highresolution simulations that couple disparate physical processes (e.g., elastodynamics, thermal weakening, pore-fluid transport, and heat conduction). Compounding the computational challenge, we know that natural faults are not planar, but instead have roughness that can be approximated by power laws potentially leading to large, multiscale fluctuations in normal stress. The capacity to perform 3D rupture simulations that couple these processes will provide guidance for constructing appropriate source models for high-frequency ground motion simulations. The improved rupture models from our multi-scale dynamic rupture simulations will be used to conduct physicsbased (3D waveform modeling-based) probabilistic seismic hazard analysis (PSHA) for California. These calculation will provide numerous important seismic hazard results, including a state-wide extended earthquake rupture forecast with rupture variations for all significant events, a synthetic seismogram catalog for thousands of scenario events and more than 5000 physics-based seismic hazard curves for California.
This final scientific report documents the Industrial Technology Program (ITP) Stage 2 Concept Development effort on Data Center Energy Reduction and Management Through Real-Time Optimal Control (RTOC). Society is becoming increasingly dependent on information technology systems, driving exponential growth in demand for data center processing and an insatiable appetite for energy. David Raths noted, 'A 50,000-square-foot data center uses approximately 4 megawatts of power, or the equivalent of 57 barrels of oil a day1.' The problem has become so severe that in some cases, users are giving up raw performance for a better balance between performance and energy efficiency. Historically, power systems for data centers were crudely sized to meet maximum demand. Since many servers operate at 60%-90% of maximum power while only utilizing an average of 5% to 15% of their capability, there are huge inefficiencies in the consumption and delivery of power in these data centers. The goal of the 'Recovery Act: Decreasing Data Center Energy Use through Network and Infrastructure Control' is to develop a state of the art approach for autonomously and intelligently reducing and managing data center power through real-time optimal control. Advances in microelectronics and software are enabling the opportunity to realize significant data center power savings through the implementation of autonomous power management control algorithms. The first step to realizing these savings was addressed in this study through the successful creation of a flexible and scalable mathematical model (equation) for data center behavior and the formulation of an acceptable low technical risk market introduction strategy leveraging commercial hardware and software familiar to the data center market. Follow-on Stage 3 Concept Development efforts include predictive modeling and simulation of algorithm performance, prototype demonstrations with representative data center equipment to verify requisite performance and continued commercial partnering agreement formation to ensure uninterrupted development, and deployment of the real-time optimal control algorithm. As a software implementable technique for reducing power consumption, the RTOC has two very desirable traits supporting rapid prototyping and ultimately widespread dissemination. First, very little capital is required for implementation. No major infrastructure modifications are required and there is no need to purchase expensive capital equipment. Second, the RTOC can be rolled out incrementally. Therefore, the effectiveness can be proven without a large scale initial roll out. Through the use of the Impact Projections Model provided by the DOE, monetary savings in excess of $100M in 2020 and billions by 2040 are predicted. In terms of energy savings, the model predicts a primary energy displacement of 260 trillion BTUs (33 trillion kWh), or a 50% reduction in server power consumption. The model also predicts a corresponding reduction of pollutants such as SO2 and NOx in excess of 100,000 metric tonnes assuming the RTOC is fully deployed. While additional development and prototyping is required to validate these predictions, the relative low cost and ease of implementation compared to large capital projects makes it an ideal candidate for further investigation.
The FNAL accelerator complex has been upgrading in increasing beam intensity and beam quality. A new beam halo diagnostic device is required in the beam transport line between booster and Recycler. For this purpose, it was decided to introduce the wide dynamic range monitor technique that was developed in 2012 and has been in operation at the J-PARC beam transport line. The device is a two-dimensional beam profile monitor, and it has a dynamic range of approximately six digits of magnitude by using of Optical Transition Radiation and fluorescence screens. Eliminating harmful beam halos is the most important technique for high-intensity proton accelerators. Therefore, beam halo diagnosis is indispensable and becomes more and more important. New FNAL device has been manufactured in a collaboration between J-PARC and FNAL as a part of U.S.-Japan Science and Technology Cooperation Program in High Energy Physics. The equipment will be manufactured at J-PARC and will be shipped to FNAL in 2025.We designed the device to satisfy FNAL specifications: the beam energy, intensity, and size. Currently, most of the equipments are under construction. The large-aperture optical system has been completed and its optical characteristics are being evaluated at J-PARC. We have been also investigating measurement methods corresponding to FNAL bunch trains. This paper reports on the current status of these developments.
Abstract Objective The Clark Fork–Pend Oreille River basin of northeastern Washington and the Idaho Panhandle historically supported a robust metapopulation of the Westslope Cutthroat Trout (WCT) Oncorhynchus lewisi, a western native salmonid of high cultural and economic value. The construction of impassible hydroelectric dams and smaller instream barriers has prevented the return of migratory WCT to spawning tributaries, leading to the fragmentation of this metapopulation over the past 100 years. One such impassible barrier is Albeni Falls Dam (AFD) near Newport, Washington, which was completed without fish passage capabilities in 1955. We sought to examine large-scale genetic patterns in the study area and determine the most likely spawning tributary of origin for migratory WCT captured below AFD. Methods We created a genetic baseline representative of populations within the Clark Fork–Pend Oreille River basin from upstream and downstream of the dam using 191 biallelic single-nucleotide polymorphism genetic markers. Our data set included 124 collections, which allowed for an examination of population structure and hatchery influence across the study area and provided a robust tool for population assignment. Population assignment tests were conducted using the program RUBIAS. Result Population assignment tests were successful for all pure WCT of unknown origin despite potential influence from hatchery lineages across the study area. Of 83 migratory WCT captured below AFD, approximately 80% were assigned to tributaries upstream of AFD with a posterior assignment probability of at least 90%. Only one fish was assigned to a tributary downstream of AFD. Conclusion Our results indicate that AFD disrupts the natural metapopulation dynamics of WCT populations in the basin. Passage for WCT at this barrier would reestablish metapopulation connectivity within the basin by allowing migratory individuals to make genetic contributions to populations upstream of the dam.
The focus of this research program is on developing density functional methods for the calculation of non-equilibrium phenomena of the transport type (finite wavelength, low frequency) and of the dynamical type (long wavelength, high frequency). One such method is the time-dependent current density functional theory, in which the evolution of the physical system is mimicked in a fictitious non-interacting system subjected to an effective one-particle potential, which we approximate as a function of the equilibrium density and non-equilibrium current density.
This study covers the onset of flow instability (OFI) preliminary results obtained from leveraging correlations, in addition to the preliminary thermal hydraulics results such as pressure, flow velocity, temperature, and oxide layer over the design demonstration experiment (DDE) for the Massachusetts Institute of Technology Reactor (MITR). Current computational fluid dynamics (CFD) models in fluid structure interaction (FSI) have added the capability of assessing margins to onset of nucleate boiling (ONB). This study initiates the capability to model the margin to OFI and ONB presented for the MITR. Such study is supportive of the United States High Performance Research Reactor (USHPRR) program. Previous studies provided preliminary thermal-hydraulic and mechanical analyses of the hydrodynamic effects in the MITR DDE under conservative approximations for plate power distribution. In conclusion, this study focuses on providing insights into the OFI future research direction optimizing the transport of thermal energy, mass-flow rates, flow-channel geometries, and boundary conditions.
Eliminating harmful beam halos is the most important technique for high-intensity proton accelerators. Therefore, beam halo diagnosis is indispensable and becomes more and more important. At J-PARC, a wide dynamic range monitor was installed in the beam transport line in 2012. The device is a two-dimensional beam profile monitor [*, **], and it has a dynamic range of approximately six digits of magnitude by using Optical Transition Radiation and fluorescence screens. The FNAL accelerator complex has been upgrading through increased beam intensity and beam quality. A new beam halo diagnostic device is required in the beam transport line between the booster and recycler. It will be manufactured in a collaboration between J-PARC and FNAL as a part of the U.S.-Japan Science and Technology Cooperation Program in High Energy Physics. We are redesigning the monitor to satisfy FNAL specifications for beam energy, intensity, and size. The equipment will be manufactured at J-PARC and then shipped to FNAL in 2024. In this report, the design of the device will be described.
Purpose: This study aims to evaluate the potential and feasibility of positron emission tomography for dynamic lung tumor tracking during radiation treatment. The authors propose a center of mass (CoM) tumor tracking algorithm using gated-PET images combined with a respiratory monitor and investigate the geometric accuracy of the proposed algorithm. Methods: The proposed PET dynamic lung tumor tracking algorithm estimated the target position information through the CoM of the segmented target volume on gated PET images reconstructed from accumulated coincidence events. The information was continuously updated throughout a scan based on the assumption that real-time processing was supported (actual processing time at each frame ≈10 s). External respiratory motion and list-mode PET data were acquired from a phantom programmed to move with measured respiratory traces (external respiratory motion and internal target motion) from human subjects, for which the ground truth target position was known as a function of time. The phantom was cylindrical with six hollow sphere targets (10, 13, 17, 22, 28, and 37 mm in diameter). The measured respiratory traces consisted of two sets: (1) 1D-measured motion from ten healthy volunteers and (2) 3D-measured motion from four lung cancer patients. The authors evaluated the geometric accuracy of the proposed algorithm by quantifying estimation errors (Euclidean distance) between the actual motion of targets (1D-motion and 3D-motion traces) and CoM trajectories estimated by the proposed algorithm as a function of time. Results: The time-averaged error of 1D-motion traces over all trajectories of all targets was 1.6 mm. The error trajectories decreased with time as coincidence events were accumulated. The overall error trajectory of 1D-motion traces converged to within 2 mm in approximately 90 s. As expected, more accurate results were obtained for larger targets. For example, for the 37 mm target, the average error over all 1D-motion traces was 1.1 mm; and for the 10 mm target, the average error over all 1D-motion traces was 2.8 mm. The overall time-averaged error of 3D-motion traces was 1.6 mm, which was comparable to that of the 1D-motion traces. There were small variations in the errors between the 3D-motion traces, although the motion trajectories were very different. The accuracy of the estimates was consistent for all targets except for the smallest. Conclusions: The authors developed an algorithm for dynamic lung tumor tracking using list-mode PET data and a respiratory motion signal, and demonstrated proof-of-principle for PET-guided lung tumor tracking. The overall tracking error in phantom studies is less than 2 mm.