Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “deterministic method”

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 325 records · Page 18

Offline Maximizing Minimally Invasive Proper Orthogonal Decomposition for Reduced-Order Modeling of S n Radiation Transport

Deterministic solutions to the Sn radiation transport equation can be computationally expensive to calculate. Reduced-order modeling enables efficient approximation of the full-order model (FOM) solution. We propose a novel method for constructing reduced-order models (ROMs) of the S n radiation transport equation, offline maximizing minimally invasive (OMMI) proper orthogonal decomposition (POD). POD uses the method of snapshots to create a reduced-order basis for constructing an ROM. Minimally invasive POD leverages the sweep infrastructure existing in deterministic transport codes to create a POD-based ROM, even when infeasible by traditional methods. Offline maximizing minimally invasive proper orthogonal decomposition (OMMI-POD) extends minimally invasive POD by performing sweeps offline, therefore maximizing the potential speedup. OMMI-POD does so by creating a library of reduced systems from a training set. This library of reduced systems is then interpolated to provide a rapid approximate solution of the S n radiation transport equation. The model is evaluated on a set of test problems, achieving a low error with a 466 times speedup over the FOM. Also presented is a study of the effect of sampling method on the performance of OMMI-POD, specifically comparing naive uniform sampling to the more accurate and computationally expensive greedy sampling.

97 MATHEMATICS AND COMPUTING↗

Comparison of ADVANTG Simulations to Simple Shielding Measurements - Paper 36

Monte Carlo (MC) radiation transport simulations, like those performed using Monte Carlo N-Particle (MCNP), can be highly accurate but require time to reduce the statistical uncertainties in the calculated detector response(s). Many different 'variance reduction' techniques devised over the years increase the efficiency of MC calculations, but these methods can be more art than science. Many of these methods are problem specific or require careful tuning of parameters for good performance. One of the best variance reduction methods, in terms of both reliable performance increases and ease of use, is the use of weight windows. Weight windows are derived from an importance map in space and energy that allows the transport code to focus more simulation time on particles that have a higher likelihood of being important to the final detector response. Less time is spent simulating particles that have little chance of making a contribution to the detector response. MCNP already includes the abilities to generate and use weight windows. Unfortunately, generating weight windows can take a very long time, perhaps longer than the original problem. Oak Ridge National Laboratory has led the development of hybrid methods using a discrete ordinates adjoint calculation to quickly develop an importance map (weight windows) for use in a Monte Carlo code. The Consistent Adjoint Driven Importance Sampling (CADIS) method uses an adjoint calculation to develop both an importance map and a biased source distribution that work together. This method is reasonably automatic, is easy to use, and has been implemented in the Automated Variance reduction Generator (ADVANTG) package that works with MCNP5. ADVANTG reads an MCNP input file and a small amount of extra user information and then runs an adjoint deterministic calculation using Denovo. The importance map and biased source distribution derived from the CADIS method can then be incorporated back into the original MCNP input file to enhance the figure-of-merit (FOM) of the calculation. The FOM is a measure of the efficiency of the calculation using the amount of run time and corresponding level of statistical uncertainty in the detector result. ADVANTG-enhanced MCNP input files have been shown to increase the FOM of many types of problems, obtaining the same level of statistical uncertainty in a shorter time compared to standard MCNP. ADVANTG can also generate variance reduction parameters (importance map and biased source) using the Forward-Weighted CADIS (FW-CADIS) method for optimizing mesh tallies. The CADIS method creates variance reduction parameters that optimize a tally's total quantity (flux, dose rate, interaction rate, etc.) and will not generally increase performance for all parts of a segmented tally. To optimize a segmented tally, the FW-CADIS method uses two discrete ordinates calculations to make variance reduction parameters that help compute a segmented tally with uniform relative uncertainties in each bin. This report will demonstrate that using ADVANTG for neutron and photon radiation shielding problems will (1) preserve the standard MCNP results and (2) obtain those results either with less statistical uncertainty or in shorter times than the standard MCNP calculation. The four demonstration problems in this report are based on simple neutron or photon transmission measurements through different materials. All of the problems here used the CADIS method and still saw improvement across the flux tallies segmented by particle energy. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Phase 2: Updated Geothermal Play Fairway Analysis of the Tularosa Basin, New Mexico

Play fairway analysis (PFA) methodologies, both deterministic and stochastic, are being applied for medium to low enthalpy geothermal resources in the Tularosa Basin, in south-central New Mexico. These PFA methods use datasets representative of heat, groundwater and fault-fracture related permeability to locate areas with the lowest risk of geothermal energy development. Our Phase 1 PFAs predicted several plays throughout the study area, with a range of certainty and risk, based upon legacy data. In Phase 2, new data were collected for high priority plays and integrated into the PFAs to enhance certainty and further reduce risk. Field geology, gravity surveys and alluvial fan geomorphology supplemented the understanding of geologic structure and potential permeability, while additional well sampling, geothermometry and temperature logging provided valuable information about heat in the subsurface. In order to statistically compare and contrast the layers of evidence in the Tularosa Basin with the broader Great Basin, a probabilistic Weights of Evidence analysis was conducted in parallel with the inclusive and deterministic petroleum industry logic PFA. Both methods yielded 3 similar plays, while Weights of Evidence identified 2 additional plays. A sixth play was inferred around the NASA-Aerospace Data Facility Southwest site, north-east of Las Cruces, based on a promising shallow temperature survey, geothermometers, and structural interpretations of geophysical data. Here, the White Sands Missile Range main cantonment (headquarters) play is now given high priority, along with McGregor Range, Fort Bliss, due to the discovery of relatively high geothermometer values.

15 GEOTHERMAL ENERGY↗

Convergence Analysis of Fixed Point Chance Constrained Optimal Power Flow Problems

For optimal power flow problems with chance constraints, a particularly effective method is based on a fixed point iteration applied to a sequence of deterministic power flow problems. However, a priori, the convergence of such an approach is not necessarily guaranteed. Here this article analyses the convergence conditions for this fixed point approach, and reports numerical experiments including for large IEEE networks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Allocation and Operation of A Hydropneumatic Energy Storage with Building Microgrid

Behind-meter application for a recently invented hydropneumatic Ground-Level Integrated Diverse Energy Storage (GLIDES) is focused in this research. Considering possible uncertainties from solar radiation and electricity load in a building microgrid, two-stage stochastic programming is adopted for 15-min operation and a mixed-integer nonlinear stochastic model is built with integrated energy systems of GLIDES, solar panel and power generating unit. In order to make the model computationally tractable, we use sample average approximate method to decrease associated complexity with sample size and improve solving efficiency. The experimental results in deterministic and stochastic operation indicate that best candidate solution from SAA are superior to the solution from expected value problem in terms of optimality gap and variance.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Some Results of Multicriterial and Robust Optimization of Nuclear Power Systems Structure

The possibility of applying multicriterial and robust optimization methods as well as stochastic methods of determining the effect of uncertainties in the parameters of a system to structural optimization of nuclear power is discussed for the example of a hypothetical nuclear-power system consisting of thermal and fast reactors. It is shown that the optimal structures of nuclear power systems obtained by such methods are balanced in terms of conflicting criteria and technologically more diversified than those determined by single-criterion deterministic optimization. Cost-based innovative technologies, which can be adopted only by imposing additional external limitations in a single-criterion problem, appear in such structures. The examined methods make it possible to determine on a systematic basis the structures of nuclear power that are balanced in terms of a set of conflicting criteria, since factors having both positive and negative effects on the development of nuclear-power systems are included in the calculations, and to find economically efficient directions of structural changes in the system in order to improve stability and efficiency.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Operational Probabilistic Tools for Solar Uncertainty (OPTSUN) (Final Project Report for DOE Solar Forecasting II Project)

Increasing levels of solar PV can challenge system operations and may require novel methods to operate the power system reliably and efficiently. Power system operating plans generally use deterministic forecasts, in which the variable energy resources are represented by the expected value for each interval of the decision horizon. Probabilistic forecasts are relatively new but have the potential to address the shortfalls of deterministic forecasts. However, understanding how best to use such forecasts is still a key gap in industry and was the focus of this project. The project had three workstreams. In a forecasting workstream, improvements were made to baseline probabilistic forecasts using a number of new approaches such as machine learning methods and improved input data. In a design workstream, advanced simulation tools used these forecasts to investigate newly proposed reserve determination methods. Lastly, in a demonstration workstream a scheduling management platform (SMP) was developed to leverage probabilistic forecasts in a modular and customizable manner. In order to study the benefits that could be accrued, the project team collaborated with three utility partners (Duke Energy, Southern Company and Hawaiian Electric) to deliver improved probabilistic forecasts for each region and to model each region in case studies using advanced production cost modeling tools. Different methods to determine operating reserve requirements from probabilistic forecasts were developed, simulated, and tested across each region. The benefits of using these newly proposed methods varied by utility, but, in general, using probabilistic forecasts as well as historical data to set the reserve requirements seems to improve reliability related results, with less risk of reserve or supply shortfalls. The cost implications were not always straightforward; in some cases the new methods could show a reduction in expected operating costs, but often the increase in reserves associated with better risk mitigation using probabilistic forecasts could result in an increase in operating costs in the simulations. The SMP tool was developed to process probabilistic forecasts from their initial receipt through to scheduling decisions. This open-source tool consists of several modules for scenario development, reserve requirements calculation, and visualization. The SMP tool was demonstrated to a wide range of operators and stakeholders at all three utilities and further improved based on their feedback. The tool will be available on www.epri.com/optsun. The proposed probabilistic information-based reserve determination approaches have the potential to be implemented by different regions to ensure an economic and reliable power system operation on power systems integrating increasing levels of variable renewable resources. The innovative yet practical methods developed in this project demonstrated tangible benefits from using probabilistic forecasts beyond just study-based assessments to include three unique balancing areas. The demonstrated benefits across the multiple utility environments, are expected to provide system operators in all regions the confidence required and a platform to adopt the new forecasting and operating methods.

14 SOLAR ENERGY↗

Tensor decompositions for count data that leverage stochastic and deterministic optimization

There is growing interest to extend low-rank matrix decompositions to multi-way arrays, or tensors. One fundamental low-rank tensor decomposition is the canonical polyadic decomposition (CPD). The challenge of fitting a low-rank, nonnegative CPD model to Poisson-distributed count data is of particular interest. Several popular algorithms use local search methods to approximate the maximum likelihood estimator (MLE) of the Poisson CPD model. Here, this work presents two new algorithms that extend state-of-the-art local methods for Poisson CPD. Hybrid GCP-CPAPR combines Generalized Canonical Decomposition (GCP) with stochastic optimization and CP Alternating Poisson Regression (CPAPR), a deterministic algorithm, to increase the probability of converging to the MLE over either method used alone. Restarted CPAPR with SVDrop uses a heuristic based on the singular values of the CPD model unfoldings to identify convergence toward optimizers that are not the MLE and restarts within the feasible domain of the optimization problem, thus reducing overall computational cost when using a multi-start strategy. We provide empirical evidence that indicates our approaches outperform existing methods with respect to converging to the Poisson CPD MLE.

CPAPR↗

Integral Kernel Methods for Nonlinear Parabolic-Elliptic Systems

Nonlinear parabolic-elliptic systems arise in many physical, biological, and chemical phenomena such as chemotaxis, ion transport, self-gravitating particles, and Brownian vortices. Existing methods struggle with the strong coupling and high nonlinearity and nonlocality of some of these systems, especially the ill-conditioned, convection-dominated problems. To overcome numerical difficulties, current approaches rely on initial guesses, preconditioning, or iterative techniques with no convergence guarantees. They might suffer from poor scalability, large memory usage, and difficulty to parallelize. Inspired by the connection of parabolic-elliptic systems to stochastic processes, we introduce a novel meshless, monolithic, and fully explicit method that naturally encapsulates the elliptic and parabolic operators into a single step which updates each node deterministically with global information. By being fully quadrature-based, it avoids solving systems of discretized equations and does not utilize initial guesses or preconditioning, while requiring little memory and being easy to parallelize. We first derive the method in an integral kernel formulation with quadratic complexity in the number of integration nodes and then leverage kernel-independent fast multipole methods (FMM) to present a scalable algorithm with linear complexity. We provide numerical examples for the Poisson-Nernst-Planck equations in one, two, and three dimensions, together with the derivation of the integral kernel for each case. Furthermore, the examples demonstrate the fast convergence and scalability of the FMM-accelerated algorithm, as well as its suitability for convection-dominated problems, making it competitive against traditional PDE solvers.

PDE systems↗

Deterministic fabrication of highly reproducible monochromatic quantum emitters in hexagonal boron nitride

Quantum emitters in hexagonal boron nitride are important room temperature single-photon sources. However, conventional fabrication methods yield quantum emitters with dispersed and inconsistent spectral profiles, limiting their potential for practical quantum applications, which demand reproducible high quality single-photon sources. Here, we report the deterministic creation of highly reproducible monochromatic quantum emitters by applying carbon-ion implantation on freestanding hexagonal boron nitride flakes, while a carbon mask with suitable thickness was adapted to optimize the implantation results. Quantum emitters fabricated using this approach exhibited thermally limited monochromaticity, with an emission center wavelength of 590.7 ± 2.7 nm, a narrow full width at half maximum of 7.1 ± 1.7 nm, an emission rate of 1 MHz without optical engineering, and exceptional stability under ambient conditions. Density functional theory calculations and scanning transmission electron microscopy suggest that these emitters are comprised of boron centered carbon tetramers. This method provides a reliable single-photon source for optical quantum computing and potential future industry-scale applications.

Hua, Muchuan [Argonne National Laboratory (ANL), A↗

Effects of surface runoff and infiltration partition methods on hydrological modeling: A comparison of four schemes in two watersheds in the Northeastern US

The ratio according to which the precipitation/snowmelt is divided into surface and subsurface water is of primary importance to watershed scale hydrological modeling. The present study compared four different surface runoff and infiltration partition methods in the Richards-equation-based SWAT model (RSWAT), i.e., daily curve number method (DC-RSWAT), hourly curve number method (HC-RSWAT), Green-Ampt method (GA-RSWAT), and effective infiltration capacity method (EIC-RSWAT). These four versions of RSWAT and a daily-curve-number-based SWAT (DC-SWAT) were used to simulate daily flow rate (from 2001 to 2015) at the outlets of the Tuckahoe Creek Watershed (TCW) and Greensboro Watershed (GW) of the Chesapeake Bay in Maryland. Parameter sensitivity analysis, calibration, and prediction uncertainty analysis were conducted. The results indicate that surface runoff and infiltration partition methods have pronounced impacts on model performance. Comparative analysis showed that DC-RSWAT outperformed other versions of RSWAT for both calibration (from 2001 to2007) and validation (from 2008 to 2015) periods in the two test watersheds. In general, DC-RSWAT achieved more accurate deterministic flow simulation and reliable uncertainty estimation than other versions. As such, daily-curve-number method is recommended for surface runoff and infiltration partition in RSWAT for simulating stream flow.

hydrology, modeling, streamflow↗

Deterministic stabilization of eight-way 2D diffractive beam combining using pattern recognition

We demonstrate a new method for controlling diffractive, high-power beam combination, sensing phase errors by analyzing the intensity pattern of uncombined side beams at the output. A square array of eight beams is combined with <0.3% stability and 84.6% efficiency. As the channel count is increased, so does the usable information, enabling scaling to large channel counts without significant slowing of control loop response time, an advantage over single-input algorithms.

47 OTHER INSTRUMENTATION↗

Economic evaluation on CO₂-EOR of onshore oil fields in China

Carbon dioxide enhanced oil recovery (CO₂-EOR) and sequestration in depleted oil reservoirs is a plausible option for utilizing anthropogenic CO₂ to increase oil production while storing CO₂ underground. Evaluation of the storage resources and cost of potential CO₂-EOR projects is an essential step before the commencement of large-scale deployment of such activities. In this paper, a hybrid techno-economic evaluation method, including a performance model and cost model for onshore CO₂-EOR projects, has been developed based on previous studies. Total 296 onshore oil fields, accounting for about 70% of total mature onshore oil fields in China, were evaluated by the techno-economic method. The key findings of this study are summarized as follows: (1) deterministic analysis shows there are approximately 1.1 billion tons (7.7 billion barrels) of incremental crude oil and 2.2 billion tons CO₂ storage resource for onshore CO₂-EOR at net positive revenue within the Chinese oil fields reviewed under the given operating strategy and economic assumptions. (2) Sensitivity study highlights that the cumulative oil production and cumulative CO₂ storage resource are very sensitive to crude oil price, CO₂ cost, project lifetime, discount rate and tax policy. High oil price, short project lifetime, low discount rate, low CO₂ cost, and low tax policy can greatly increase the net income of the oil enterprise, incremental oil recovery and CO₂ storage resource. (3) From this techno-economic evaluation, the major barriers to large-scale deployment of CO₂-EOR include complex geological conditions, low API of crude oil, high tax policy, and lack of incentives for the CO₂-EOR project.

02 PETROLEUM↗

Hamming distance and the onset of quantum criticality

Simulating models for quantum correlated matter unveils the inherent limitations of deterministic classical computations. In particular, in the case of quantum Monte Carlo methods, this is manifested by the emergence of negative weight configurations in the sampling, that is, the sign problem (SP). There have been several recent calculations which exploit the SP to locate underlying critical behavior. Here, utilizing a metric that quantifies phase-space ergodicity in such sampling, the Hamming distance, we suggest a significant advance on these ideas to extract the location of quantum critical points in various fermionic models, in spite of the presence of a severe SP. Combined with other methods, exact diagonalization in our case, it elucidates both the nature of the different phases as well as their location, as we demonstrate explicitly for the honeycomb and triangular Hubbard models, in both their U(1) and SU(2) forms. Furthermore, our approach exemplifies a possible path allowing the exploration of the phase diagram of a variety of fermionic quantum models hitherto considered to be impractical via quantum Monte Carlo simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nuclear Software Validation for 7% Enriched UO 2 Fuel Lattices

This analysis performs validation for the Virtual Environment for Reaction Application’s (VERA’s) neutronics codes MPACT and Shift against measured data from the Seven Percent Critical Experiment (7uPCX) at Sandia National Laboratory. The benchmarking supports the future application of VERA for analysis of high-assay low-enriched uranium (HALEU) fuel in commercial pressurized water reactors (PWRs). VERA demonstrates very good agreement with measured data, and the deterministic code MPACT also agrees very well with higher-fidelity stochastic methods KENO-VI and Shift in both total neutron flux and fission rate distribution comparisons. Despite some significant differences between the 7uPCX and existing commercial PWRs that challenge the neutron cross section data and transport solver methods used by MPACT, all of the results are excellent. The agreement increases the confidence in the applicability and accuracy of VERA for reactor analysis of HALEU fuels in commercial reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Estimating the epidemic threshold on networks by deterministic connections

For many epidemic networks some connections between nodes are treated as deterministic, while the remainder are random and have different connection probabilities. By applying spectral analysis to several constructed models, we find that one can estimate the epidemic thresholds of these networks by investigating information from only the deterministic connections. Nonetheless, in these models, generic nonuniform stochastic connections and heterogeneous community structure are also considered. The estimation of epidemic thresholds is achieved via inequalities with upper and lower bounds, which are found to be in very good agreement with numerical simulations. Since these deterministic connections are easier to detect than those stochastic connections, this work provides a feasible and effective method to estimate the epidemic thresholds in real epidemic networks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Robust Operation of Soft Open Points in Active Distribution Networks with High Penetration of Photovoltaic Integration

Distributed generators (DGs) including photovoltaic panels (PVs) have been integrated dramatically in active distribution networks (ADNs). Due to the strong volatility and uncertainty, the high penetration of PV generation immensely exacerbates the conditions of voltage violation in ADNs. However, the emerging flexible interconnection technology based on soft open points (SOPs) provides increased controllability and flexibility to the system operation. For fully exploiting the regulation ability of SOPs to address the problems caused by PV, this paper proposes a robust optimization method to achieve the robust optimal operation of SOPs in ADNs. A two-stage adjustable robust optimization model is built to tackle the uncertainties of PV outputs, in which robust operation strategies of SOPs are generated to eliminate the voltage violations and reduce the power losses of ADNs. A column-and-constraint generation (C&CG) algorithm is developed to solve the proposed robust optimization model, which are formulated as second-order cone program (SOCP) to facilitate the accuracy and computation efficiency. Here, case studies on the modified IEEE 33-node system and comparisons with the deterministic optimization approach are conducted to verify the effectiveness and robustness of the proposed method.

14 SOLAR ENERGY↗