Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stochastic processes”

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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Chance-Constrained System of Systems Based Operation of Power Systems

In this paper, a chance-constrained system of systems (SoS) based decision-making approach is presented for stochastic scheduling of power systems encompassing active distribution grids. Based on the concept of SoS, the independent system operator (ISO) and distribution companies (DISCOs) are modeled as self-governing systems. These systems collaborate with each other to run the entire power system in a secure and economic manner. Each self-governing system accounts for its local reserve requirements and line flow constraints with respect to the uncertainties of load and renewable energy resources. A set of chance constraints are formulated to model the interactions between the ISO and DISCOs. The proposed model is solved by using analytical target cascading (ATC) method, a distributed optimization algorithm in which only a limited amount of information is exchanged between collaborative ISO and DISCOs. In this paper, a 6-bus and a modified IEEE 118-bus power systems are studied to show the effectiveness of the proposed algorithm.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

Stochastic Multi-Timescale Power System Operations With Variable Wind Generation

This paper describes a novel set of stochastic unit commitment and economic dispatch models that consider stochastic loads and variable generation at multiple operational timescales. The stochastic model includes four distinct stages: stochastic day-ahead security-constrained unit commitment (SCUC), stochastic real-time SCUC, stochastic real-time security-constrained economic dispatch (SCED), and deterministic automatic generation control (AGC). These sub-models are integrated together such that they are continually updated with decisions passed from one to another. The progressive hedging algorithm (PHA) is applied to solve the stochastic models to maintain the computational tractability of the proposed models. Comparative case studies with deterministic approaches are conducted in low wind and high wind penetration scenarios to highlight the advantages of the proposed methodology, one with perfect forecasts and the other with current state-of-the-art but imperfect deterministic forecasts. The effectiveness of the proposed method is evaluated with sensitivity tests using both economic and reliability metrics to provide a broader view of its impact.

17 WIND ENERGY↗
DOE OSTI2021

Distribution System State Estimation With High Penetration of Demand Response Enabled Loads

Demand-side operations incentivize utility customers to take part in various grid services. A demand response enabled load (DREL) is a flexible grid asset that schedules electricity consumption in response to a time-of-use (TOU) energy price. Consequently, its energy profile differs from that of a conventional load that is insensitive to price. This difference may cause new challenges for distribution system state estimation (DSSE). It is well known that DSSE often needs to use pseudo-measurements based on historic load profiles to increase system observability. However, historic profiles of conventional loads are not representative of DREL behaviors. Additionally, the inaccuracy impacts DSSE results and other DSSE-dependent operations. In this paper, we propose an online pseudo-measurement generation approach for DSSE with DRELs. We formulate an optimization model to represent DRELs self-adjusting actions. Sampling-based stochastic optimization techniques are proposed to account for uncertainties in DRELs. A set of representative DREL behavior data corresponding to the samples are used to characterize DREL pseudo-measurements. Case studies with modified IEEE 123-bus test system verify the validity of the proposed work.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Observability Analysis of a Power System Stochastic Dynamic Model Using a Derivative-Free Approach

Serving as a prerequisite to power system dynamic state estimation, the observability analysis of a power system dynamic model has recently attracted the attention of many power engineers. However, because this model is typically nonlinear and large-scale, the analysis of its observability is a challenge to the traditional derivative-based methods. Indeed, the linear-approximation-based approach may provide unreliable results while the nonlinear-technique-based approach inevitably faces extremely complicated derivations. Furthermore, because power systems are intrinsically stochastic, the traditional deterministic approaches may lead to inaccurate observability analyses. In this work, facing these challenges, we propose a novel polynomial-chaos-based derivative-free observability analysis approach that not only is free of any linear approximations, but also accounts for the stochasticity of the dynamic model while bringing a low implementation complexity. Furthermore, this approach enables us to quantify the degree of observability of a stochastic model, what conventional deterministic methods cannot do. The excellent performance of the proposed method has been demonstrated by performing extensive simulations using a synchronous generator model with IEEE-DC1A exciter and the TGOV1 turbine governor.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Chance-Constrained Energy Management System for Power Grids With High Proliferation of Renewables and Electric Vehicles

This paper proposes a two-stage energy management system (EMS) for power grids with massive integration of electric vehicles (EVs) and renewable energy resources. The first stage economic dispatch determines the optimal operating points of charging stations and battery swapping stations (BSS) for EVs under plug-in and battery swapping modes, respectively. The proposed stochastic model predictive control (SMPC) problem in this stage is characterized through a chance-constrained optimization formulation that can effectively capture the system and the forecast uncertainties. A distributed algorithm, the alternating direction method of multipliers (ADMM), is applied to accelerate the optimization computation through parallel computing. The second stage is aimed in coordinating the EV charging mechanisms to continuously follow the first-stage solutions, i.e., the target operating points, and meeting the EV customers' charging demands captured via the Advanced Metering Infrastructure (AMI). In this work, the proposed solution offers a holistic control strategy for large-scale centralized power grids in which the aggregated individual parameters are predictable and the system dynamics do not vary sharply within a short time-interval.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Bi-Level Adaptive Storage Expansion Strategy for Microgrids Using Deep Reinforcement Learning

Battery energy storage (BES) is a versatile resource for the secure and economic operation of microgrids (MGs). Prevailing stochastic optimization-based approaches for BES expansion planning for MGs are computationally complicated. This work proposes a data-driven bi-level multi-period BES expansion planning framework to determine the siting, sizing, and timing of BES installations. The proposed planning framework unifies deep reinforcement learning (DRL) and linear programming, thereby decoupling the determinations for the integer and continuous decision variables in two time scales, respectively. In the upper level, a rainbow DRL agent with quantile regression is trained to provide dynamic planning policies to accommodate stochastic renewable energy resources (RESs), load, and battery price changes efficiently. Further, the lower level computes the optimal operation of MGs with frequency constraints to hedge the islanding contingency. The two levels communicate with one another by exchanging storage configuration and operating expenses in order to accomplish the shared goal of minimizing investment and operation costs. Comparative case studies on an MG are carried out to demonstrate the superiority of the proposed DRL-based solution to the mixed-integer linear programming counterpart on efficiency, scalability, and adaptability.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2012

The effects of additive outliers on time Series Components and robust estimation: a case study on the oymapinar dam, turkey

In this paper, the behaviour of the body of the Oymapinar Dam (Antalya, Turkey) resulting from reservoir water level changes was investigated by time series analysis. First, the possible best fitting methods for determining the existing additive outliers in series were investigated. For this purpose, a time series with known parameters consisting of trend, periodical and stochastic components was obtained. Then, a new time series with additive outlier(s) was created by adding outlier(s) to series with known parameters. The parameters of the components were determined by the least squares method and it was seen that outlier(s) affect the components of the series. It was found that the correct data were contaminated by an outlier to a similar degree as the autoregressive model. The outliers affecting the components were determined by Andrews, Bisquare, Cauchy, Huber and Welsch estimators in trend analysis. The results show that the parameters of the trend were not affected by outliers and the results produced by the estimators are numerically close to each other. Thus, series of the outlier(s) of the Oymapinar Dam were determined by Bisquare estimator in a trend analysis. As a result of series analysis it was determined that the linear trend resulted from a long-term periodical of changes and the dam body responded periodically depending on the water level changes.

42 ENGINEERING↗
DOE OSTI2011

Reference air kerma and kerma-area product as estimators of peak skin dose for fluoroscopically guided interventions

Purpose: To determine more accurate regression formulas for estimating peak skin dose (PSD) from reference air kerma (RAK) or kerma-area product (KAP). Methods: After grouping of the data from 21 procedures into 13 clinically similar groups, assessments were made of optimal clustering using the Bayesian information criterion to obtain the optimal linear regressions of (log-transformed) PSD vs RAK, PSD vs KAP, and PSD vs RAK and KAP. Results: Three clusters of clinical groups were optimal in regression of PSD vs RAK, seven clusters of clinical groups were optimal in regression of PSD vs KAP, and six clusters of clinical groups were optimal in regression of PSD vs RAK and KAP. Prediction of PSD using both RAK and KAP is significantly better than prediction of PSD with either RAK or KAP alone. The regression of PSD vs RAK provided better predictions of PSD than the regression of PSD vs KAP. The partial-pooling (clustered) method yields smaller mean squared errors compared with the complete-pooling method.Conclusion: PSD distributions for interventional radiology procedures are log-normal. Estimates of PSD derived from RAK and KAP jointly are most accurate, followed closely by estimates derived from RAK alone. Estimates of PSD derived from KAP alone are the least accurate. Using a stochastic search approach, it is possible to cluster together certain dissimilar types of procedures to minimize the total error sum of squares.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2011

An approach to assessing stochastic radiogenic risk in medical imaging

Purpose: This letter suggests a formalism, the medical effective dose (MED), that is suitable for assessing stochastic radiogenic risks in diagnostic medical procedures. Methods: The MED is derived from radiobiological and probabilistic first principals, including: (1) The independence of radiation-induced biological effects in neighboring voxels at low doses; (2) the linear no-threshold assumption for stochastic radiation injury (although other dose-response relationships could be incorporated, instead); (3) the best human radiation dose-response data currently available; and (4) the built-in possibility that the carcinogenic risk to an irradiated organ may depend on its volume. The MED involves a dose-risk summation over irradiated voxels at high spatial resolution; it reduces to the traditional effective dose when every organ is irradiated uniformly and when the dependence of risk on organ volumes is ignored. Standard relative-risk tissue weighting factors can be used with the MED approach until more refined data become available. Results: The MED is intended for clinical and phantom dosimetry, and it provides an estimate of overall relative radiogenic stochastic risk for any given dose distribution. A result of the MED derivation is that the stochastic risk may increase with the volume of tissue (i.e., the number of cells) irradiated, a feature that can be activated when forthcoming radiobiological research warrants it. In this regard, the MED resembles neither the standard effective dose (E) nor the CT dose index (CTDI), but it is somewhat like the CT dose-length product (DLP). Conclusions: The MED is a novel, probabilistically and biologically based means of estimating stochastic-risk-weighted doses associated with medical imaging. Built in, ab initio, is the ability to link radiogenic risk to organ volume and other clinical factors. It is straightforward to implement when medical dose distributions are available, provided that one is content, for the time being, to accept the relative tissue weighting factors published by the International Commission of Radiological Protection (ICRP). It requires no new radiobiological data and avoids major problems encountered by the E, CTDI, and CT-E formalisms. It makes possible relative inter-patient dosimetry, and also realistic intercomparisons of stochastic risks from different protocols that yield images of comparable quality.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2012

Calculation of conversion factors for effective dose for various interventional radiology procedures

Purpose: To provide dose-area-product (DAP) to effective dose (E) conversion factors for complete interventional procedures, based on in-the-field clinical measurements of DAP values and using tabulated E/DAP conversion factors for single projections available from the literature. Methods: Nine types of interventional procedures were performed on 84 patients with two angiographic systems. Different calibration curves (with and without patient table attenuation) were calculated for each DAP meter. Clinical and dosimetric parameters were recorded in-the-field for each projection and for all patients, and a conversion factor linking DAP and effective doses was derived for each complete procedure making use of published, Monte Carlo calculated conversion factors for single static projections. Results: Fluoroscopy time and DAP values for the lowest-dose procedure (biliary drainage) were approximately 3-fold and 13-fold lower, respectively, than those for the highest-dose examination (transjugular intrahepatic portosystemic shunt, TIPS). Median E/DAP conversion factors from 0.12 (abdominal percutaneous transluminal angioplasty) to 0.25 (Nephrostomy) mSvGy{sup -1} cm{sup -2} were obtained and good correlations between E and DAP were found for all procedures, with R{sup 2} coefficients ranging from 0.80 (abdominal angiography) to 0.99 (biliary stent insertion, Nephrostomy and TIPS). The DAP values obtained in this study showed general consistency with the values provided in the literature and median E values ranged from 4.0 mSv (biliary drainage) to 49.6 mSv (TIPS). Conclusions: Values of E/DAP conversion factors were derived for each procedure from a comprehensive analysis of projection and dosimetric data: they could provide a good evaluation for the stochastic effects. These results can be obtained by means of a close cooperation between different interventional professionals involved in patient care and dose optimization.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2012

A characterization of robust radiation therapy treatment planning methods-from expected value to worst case optimization

Purpose: To characterize a class of optimization formulations used to handle systematic and random errors in radiation therapy, and to study the differences between the methods within this class. Methods: The class of robust methods that can be formulated as minimax stochastic programs is studied. This class generalizes many previously used methods, ranging between optimization of the expected and the worst case objective value. The robust methods are used to plan intensity-modulated proton therapy (IMPT) treatments for a case subject to systematic setup and range errors, random setup errors with and without uncertain probability distribution, and combinations thereof. As reference, plans resulting from a conventional method that uses a margin to account for errors are shown. Results: For all types of errors, target coverage robustness increased with the conservativeness of the method. For systematic errors, best case organ at risk (OAR) doses increased and worst case doses decreased with the conservativeness. Accounting for random errors of fixed probability distribution resulted in heterogeneous dose. The heterogeneities were reduced when uncertainty in the probability distribution was accounted for. Doing so, the OAR doses decreased with the conservativeness. All robust methods studied resulted in more robust target coverage and lower OAR doses than the conventional method. Conclusions: Accounting for uncertainties is essential to ensure plan quality in complex radiation therapy such as IMPT. The utilization of more information than conventional in the optimization can lead to robust target coverage and low OAR doses. Increased target coverage robustness can be achieved by more conservative methods.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Evaluation of an analytic linear Boltzmann transport equation solver for high-density inhomogeneities

Purpose: Acuros external beam (Acuros XB) is a novel dose calculation algorithm implemented through the ECLIPSE treatment planning system. The algorithm finds a deterministic solution to the linear Boltzmann transport equation, the same equation commonly solved stochastically by Monte Carlo methods. This work is an evaluation of Acuros XB, by comparison with Monte Carlo, for dose calculation applications involving high-density materials. Existing non-Monte Carlo clinical dose calculation algorithms, such as the analytic anisotropic algorithm (AAA), do not accurately model dose perturbations due to increased electron scatter within high-density volumes. Methods: Acuros XB, AAA, and EGSnrc based Monte Carlo are used to calculate dose distributions from 18 MV and 6 MV photon beams delivered to a cubic water phantom containing a rectangular high density (4.0-8.0 g/cm{sup 3}) volume at its center. The algorithms are also used to recalculate a clinical prostate treatment plan involving a unilateral hip prosthesis, originally evaluated using AAA. These results are compared graphically and numerically using gamma-index analysis. Radio-chromic film measurements are presented to augment Monte Carlo and Acuros XB dose perturbation data. Results: Using a 2% and 1 mm gamma-analysis, between 91.3% and 96.8% of Acuros XB dose voxels containing greater than 50% the normalized dose were in agreement with Monte Carlo data for virtual phantoms involving 18 MV and 6 MV photons, stainless steel and titanium alloy implants and for on-axis and oblique field delivery. A similar gamma-analysis of AAA against Monte Carlo data showed between 80.8% and 87.3% agreement. Comparing Acuros XB and AAA evaluations of a clinical prostate patient plan involving a unilateral hip prosthesis, Acuros XB showed good overall agreement with Monte Carlo while AAA underestimated dose on the upstream medial surface of the prosthesis due to electron scatter from the high-density material. Film measurements support the dose perturbations demonstrated by Monte Carlo and Acuros XB data. Conclusions: Acuros XB is shown to perform as well as Monte Carlo methods and better than existing clinical algorithms for dose calculations involving high-density volumes.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

A stochastic model of cell survival for high-Z nanoparticle radiotherapy

Purpose: The authors present a stochastic framework for the assessment of cell survival in gold nanoparticle radiotherapy. Methods: The authors derive the equations for the effective macroscopic dose enhancement for a population of cells with nonideal distribution of gold nanoparticles (GNP), allowing different number of GNP per cell and different distances with respect to the cellular target. They use the mixed Poisson distribution formalism to model the impact of the aforementioned physical factors on the effective dose enhancement. Results: The authors show relatively large differences in the estimation of cell survival arising from using approximated formulae. They predict degeneration of the cell killing capacity due to different number of GNP per cell and different distances with respect to the cellular target. Conclusions: The presented stochastic framework can be used in interpretation of experimental cell survival or tumor control probability studies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-161: Inverse Problems for Optical Parameters in Laser Induced Thermal Therapy

Purpose: Magnetic resonance-guided laser-induced thermal therapy (MRgLITT) is investigated as a neurosurgical intervention for oncological applications throughout the body by active post market studies. Real-time MR temperature imaging is used to monitor ablative thermal delivery in the clinic. Additionally, brain MRgLITT could improve through effective planning for laser fiber's placement. Mathematical bioheat models have been extensively investigated but require reliable patient specific physical parameter data, e.g. optical parameters. This abstract applies an inverse problem algorithm to characterize optical parameter data obtained from previous MRgLITT interventions. Methods: The implemented inverse problem has three primary components: a parameter-space search algorithm, a physics model, and training data. First, the parameter-space search algorithm uses a gradient-based quasi-Newton method to optimize the effective optical attenuation coefficient, μ-eff. A parameter reduction reduces the amount of optical parameter-space the algorithm must search. Second, the physics model is a simplified bioheat model for homogeneous tissue where closed-form Green's functions represent the exact solution. Third, the training data was temperature imaging data from 23 MRgLITT oncological brain ablations (980 nm wavelength) from seven different patients. Results: To three significant figures, the descriptive statistics for μ-eff were 1470 m{sup −1} mean, 1360 m{sup −1} median, 369 m{sup −1} standard deviation, 933 m{sup −1} minimum and 2260 m{sup −1} maximum. The standard deviation normalized by the mean was 25.0%. The inverse problem took <30 minutes to optimize all 23 datasets. Conclusion: As expected, the inferred average is biased by underlying physics model. However, the standard deviation normalized by the mean is smaller than literature values and indicates an increased precision in the characterization of the optical parameters needed to plan MRgLITT procedures. This investigation demonstrates the potential for the optimization and validation of more sophisticated bioheat models that incorporate the uncertainty of the data into the predictions, e.g. stochastic finite element methods.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-E-T-146: Effects of Uncertainties of Radiation Sensitivity of Biological Modelling for Treatment Planning

Purpose: The aim of this study was to evaluate the distribution of uncertainty of cell survival by radiation, and assesses the usefulness of stochastic biological model applying for gaussian distribution. Methods: For single cell experiments, exponentially growing cells were harvested from the standard cell culture dishes by trypsinization, and suspended in test tubes containing 1 ml of MEM(2x10{sup 6} cells/ml). The hypoxic cultures were treated with 95% N{sub 2}−5% CO{sub 2} gas for 30 minutes. In vitro radiosensitization was also measured in EMT6/KU single cells to add radiosensitizer under hypoxic conditions. X-ray irradiation was carried out by using an Xray unit (Hitachi X-ray unit, model MBR-1505R3) with 0.5 mm Al/1.0 mm Cu filter, 150 kV, 4 Gy/min). In vitro assay, cells on the dish were irradiated with 1 Gy to 24 Gy, respectively. After irradiation, colony formation assays were performed. Variations of biological parameters were investigated at standard cell culture(n=16), hypoxic cell culture(n=45) and hypoxic cell culture(n=21) with radiosensitizers, respectively. The data were obtained by separate schedule to take account for the variation of radiation sensitivity of cell cycle. Results: At standard cell culture, hypoxic cell culture and hypoxic cell culture with radiosensitizers, median and standard deviation of alpha/beta ratio were 37.1±73.4 Gy, 9.8±23.7 Gy, 20.7±21.9 Gy, respectively. Average and standard deviation of D{sub 50} were 2.5±2.5 Gy, 6.1±2.2 Gy, 3.6±1.3 Gy, respectively. Conclusion: In this study, we have challenged to apply these uncertainties of parameters for the biological model. The variation of alpha values, beta values, D{sub 50} as well as cell culture might have highly affected by probability of cell death. Further research is in progress for precise prediction of the cell death as well as tumor control probability for treatment planning.

60 APPLIED LIFE SCIENCES↗