Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimization methods”

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

TH-CD-209-05: Impact of Spot Size and Spacing On the Quality of Robustly-Optimized Intensity-Modulated Proton Therapy Plans for Lung Cancer

Purpose: To investigate how spot size and spacing affect plan quality, especially, plan robustness and the impact of interplay effect, of robustly-optimized intensity-modulated proton therapy (IMPT) plans for lung cancer. Methods: Two robustly-optimized IMPT plans were created for 10 lung cancer patients: (1) one for a proton beam with in-air energy dependent large spot size at isocenter (σ: 5–15 mm) and spacing (1.53σ); (2) the other for a proton beam with small spot size (σ: 2–6 mm) and spacing (5 mm). Both plans were generated on the average CTs with internal-gross-tumor-volume density overridden to irradiate internal target volume (ITV). The root-mean-square-dose volume histograms (RVH) measured the sensitivity of the dose to uncertainties, and the areas under RVH curves were used to evaluate plan robustness. Dose evaluation software was developed to model time-dependent spot delivery to incorporate interplay effect with randomized starting phases of each field per fraction. Patient anatomy voxels were mapped from phase to phase via deformable image registration to score doses. Dose-volume-histogram indices including ITV coverage, homogeneity, and organs-at-risk (OAR) sparing were compared using Student-t test. Results: Compared to large spots, small spots resulted in significantly better OAR sparing with comparable ITV coverage and homogeneity in the nominal plan. Plan robustness was comparable for ITV and most OARs. With interplay effect considered, significantly better OAR sparing with comparable ITV coverage and homogeneity is observed using smaller spots. Conclusion: Robust optimization with smaller spots significantly improves OAR sparing with comparable plan robustness and similar impact of interplay effect compare to larger spots. Small spot size requires the use of larger number of spots, which gives optimizer more freedom to render a plan more robust. The ratio between spot size and spacing was found to be more relevant to determine plan robustness and the impact of interplay effect than spot size alone. This research was supported by the National Cancer Institute Career Developmental Award K25CA168984, by the Fraternal Order of Eagles Cancer Research Fund Career Development Award, by The Lawrence W. and Marilyn W. Matteson Fund for Cancer Research, by Mayo Arizona State University Seed Grant, and by The Kemper Marley Foundation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Farm Deployable Microbial Bioreactor for Fuel Ethanol Production

Research was conducted to develop a farm and field deployable microbial bioreactor for bioethanol production from biomass. Experiments were conducted to select the most efficient microorganisms for conversion of plant fiber to sugars for fermentation to ethanol. Mixtures of biomass and surface soil samples were collected from selected sites in Alabama black belt counties (Macon, Sumter, Choctaw, Dallas, Montgomery, Lowndes) and other areas within the state of Alabama. Experiments were conducted to determine the effects of culture parameters on key biomass saccharifying enzymes (cellulase, beta-glucosidase, xylanase and beta-xylosidase). A wide-scale sampling of locally-grown fruits in Central Alabama was embarked to isolate potential xylose fermenting microorganisms. Yeast isolates were evaluated for xylose fermentation. Selected microorganisms were characterized by DNA based methods. Factors affecting enzyme production and biomass saccharification were examined and optimized in the laboratory. Methods of biomass pretreatment were compared. Co-production of amylolytic enzymes with celluloytic-xylanolytic enzymes was evaluated; and co-saccharification of a combination of biomass, and starch-rich materials was examined. Simultaneous saccharification and fermentation with and without pre-saccharifcation was studied. Whole culture broth and filtered culture broth simultaneous saccahrifcation and fermentation were compared. A bioreactor system was designed and constructed to employ laboratory results for scale up of biomass saccharification.

09 BIOMASS FUELS↗
DOE OSTI2024

Evaluation of DED and LPBF Fe-based Alloys Process Application Envelopes based on Performance, Process Economics, Supply Chain Risks, and Reactor-specific Targeted Components

The U.S. Department of Energy (DOE), Office of Nuclear Energy (NE), Advanced Materials and Manufacturing Technologies (AMMT) program aims to develop extreme-environment materials solutions for use in the deployment of advanced nuclear reactors and the sustainment of the current fleet. To achieve this objective, a combination of experiment, a computational tool, and machine learning (ML) for the design of materials is adopted for the maturation of materials for nuclear technology. Through advanced manufacturing techniques such as laser powder bed fusion (LPBF) and laser powder direct energy deposition (LP-DED), components with complex geometries can be fabricated with reduced time and effort. Such advanced manufacturing methods can also provide the opportunity to improve materials performance through optimized microstructures and mechanical properties. However, existing engineering alloys are not always well suited for fabrication with additive manufacturing (AM), as their compositions have been tuned to optimize fabrication via conventional methods. Thus, similar alloys with modified compositions that are better suited for AM can be studied for improved performance. Over the past three years, the AMMT teams from Argonne National Laboratory (ANL) and Pacific Northwest National Laboratory (PNNL) studied various known Fe-based alloys by evaluating their initial printability using LPBF, and an AMMT-developed down-selection and decision matrix reduced the number of alloys to be studied from six to three in fiscal year (FY) 2024. Additionally, in FY 2024, for parallel evaluation, these three alloys were studied using LPDED. While LPBF is better for small- to medium-sized components with high detail and internal features, LP-DED combines a material feed system to place the powder onto the exact spot where the laser will melt the material. This AM method can be easily scaled to extremely large components and provides high build rate speeds compared to those of conventional LPBF systems. Additionally, DED is a better choice for complex geometries and compositional gradients.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Feedback-based quantum algorithm inspired by counterdiabatic driving

In recent quantum algorithmic developments, a feedback-based approach has shown promise for preparing quantum many-body system ground states and solving combinatorial optimization problems. This method utilizes quantum Lyapunov control to iteratively construct quantum circuits. Here, we propose a substantial enhancement by implementing a protocol that uses ideas from quantum Lyapunov control and the counterdiabatic driving protocol, a key concept from quantum adiabaticity. Our approach introduces an additional control field inspired by counterdiabatic driving. We apply our algorithm to prepare ground states in one-dimensional quantum Ising spin chains. Comprehensive simulations demonstrate a remarkable acceleration in population transfer to low-energy states within a significantly reduced time frame compared to conventional feedback-based quantum algorithms. This acceleration translates to a reduced quantum circuit depth, a critical metric for potential quantum computer implementation. We validate our algorithm on the IBM cloud computer, highlighting its efficacy in expediting quantum computations for many-body systems and combinatorial optimization problems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

WE-G-204-08: Optimized Digital Radiographic Technique for Lost Surgical Devices/Needle Identification

Purpose: Small surgical devices and needles are used in many surgical procedures. Conventionally, an x-ray film is taken to identify missing devices/needles if post procedure count is incorrect. There is no data to indicate smallest surgical devices/needles that can be identified with digital radiography (DR), and its optimized acquisition technique. Methods: In this study, the DR equipment used is a Canon RadPro mobile with CXDI-70c wireless DR plate, and the same DR plate on a fixed Siemens Multix unit. Small surgical devices and needles tested include Rubber Shod, Bulldog, Fogarty Hydrogrip, and needles with sizes 3-0 C-T1 through 8-0 BV175-6. They are imaged with PMMA block phantoms with thickness of 2–8 inch, and an abdomen phantom. Various DR techniques are used. Images are reviewed on the portable x-ray acquisition display, a clinical workstation, and a diagnostic workstation. Results: all small surgical devices and needles are visible in portable DR images with 2–8 inch of PMMA. However, when they are imaged with the abdomen phantom plus 2 inch of PMMA, needles smaller than 9.3 mm length can not be visualized at the optimized technique of 81 kV and 16 mAs. There is no significant difference in visualization with various techniques, or between mobile and fixed radiography unit. However, there is noticeable difference in visualizing the smallest needle on a diagnostic reading workstation compared to the acquisition display on a portable x-ray unit. Conclusion: DR images should be reviewed on a diagnostic reading workstation. Using optimized DR techniques, the smallest needle that can be identified on all phantom studies is 9.3 mm. Sample DR images of various small surgical devices/needles available on diagnostic workstation for comparison may improve their identification. Further in vivo study is needed to confirm the optimized digital radiography technique for identification of lost small surgical devices and needles.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Maximizing the Benefits of Volt/VAR Optimization in the Presence of Community Energy Storage

This paper reports on the need for coordinating smart grid volt/volt-ampere reactive (VAR) control (VVC) and community energy storage (CES) to maximize the benefits of both approaches. It addresses a key application of volt/VAR optimization (VVO): the conservation of voltage reduction (CVR). The lack of coordination among network controls, such as VVO, and asset controls at the community battery storage system can lead to detrimental effects on grid operation. A smart grid-enabled VVC uses the VVO to perform the CVR operation. VVO can perform numerous tasks for advanced distribution management systems with coordination of the system operator. To achieve the optimal VVC solution, we use a particle swarm optimization driven VVO method. The test system has been simulated for two scenarios-during light loading and high loading condition-with and without the CES system. In both scenarios, the operation of CVR and CES have been performed considering different case studies. The outcomes of the simulated results demonstrate that significant improvement in peak load-demand reduction is achieved by the proposed VVC and CES system with reduced system losses.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Performance Evaluation of an Offshore Wave Measurement Buoy in Monochromatic Waves

The accurate measurement of waves underpins marine energy resource characterization, device design, and project development. Datawell wave buoys are widely deployed and have long served as a trusted standard for wave measurements. We quantify the measurement performance, including wave elevation and energy flux estimation, of a Datawell DWR-MkIII buoy using prescribed monochromatic heave motions on a large-amplitude six-degree-of-freedom motion platform at the National Laboratory of the Rockies, assuming the buoy behaves as an ideal wave follower. Commanded motions were validated with an optical motion tracking system while buoy elevation and raw acceleration were recorded. Wave elevations were propagated to wave energy flux estimation using four methods, including one frequency-domain method and three time-domain methods. The Bayesian optimization was applied for design of experiments, and records from three test sites were also applied and evaluated in the present study. Results show two error regions within the nominal period range of 1.6 s to 30 s. For wave periods between 5 s and 25 s, the buoy provides accurate wave height measurements. For short periods less than 5 s, the 1.28 Hz sampling frequency induces sub-Nyquist artifacts that bias elevation and can drive maximum energy flux estimation errors above 100%. For long periods exceeding 25 s, the buoy reported elevation is underpredicted with error depending on period but relatively independent of wave height, with maximum wave height and wave energy flux errors reaching 64% and 87%, respectively. Furthermore, analysis of three field-derived cases shows that frequency-domain estimates at 1.28 Hz agree within 2% of the corresponding 100 Hz estimates, while larger method-dependent differences are observed for the Hilbert method.

16 TIDAL AND WAVE POWER↗
DOE OSTI2018

Matlab Code for Multiscale Topology Optimization

This Matlab script implements a simple topology optimization algorithm that minimizes the compliance of a structure within a 3D prismatic design domain. It uses a simple linear elastic finite element solver to solve the equilibrium equations, and the optimality criterion (00) method to update the design iteratively. This functionality is identical to, or the simplest possible extension from 2D to 3D of, a 99-line Matlab code released in 2001 in a journal article by Sigmund. What is new in this code is our implementation of simple surrogate models for several truss micro-architectures, which converts this code to a multiscale design tool, whereas the original code is single-scale. The development of these surrogate models and their implementation is described in a journal article we intend to submit. We plan to include the full source code in an Appendix in this paper. The code’s primary function is as a pedagogical tool rather than as a finished software product. Sigmund’s 2001 paper describes the implementation line-by-line. We are including this updated version of the code to allow the reader to understand exactly what changes must be made to a “standard” topology optimization code to implement the surrogate models and thus multiscale design optimization.

Watts, Seth↗
DOE OSTI2015

PROBING THE ISOTROPY OF COSMIC ACCELERATION TRACED BY TYPE Ia SUPERNOVAE

We present a method to test the isotropy of the magnitude–redshift relation of Type Ia Supernovae (SNe Ia) and single out the most discrepant direction (in terms of the signal-to-noise ratio) with respect to the all-sky data. Our technique accounts for possible directional variations of the corrections for SNe Ia and yields all-sky maps of the best-fit cosmological parameters with arbitrary angular resolution. To show its potential, we apply our method to the high-redshift SNe Ia from the recent Union2.1 compilation, building maps with three different angular resolutions. We use a Monte Carlo method to estimate the statistical significance with which we could reject the null hypothesis that the magnitude–redshift relation is isotropic based on the properties of the observed most discrepant directions. We find that, based on pure signal-to-noise arguments, the null hypothesis cannot be rejected at any meaningful confidence level. However, if we also consider that the strongest deviations in the Union2.1 sample closely align with the dipole temperature anisotropy of the cosmic microwave background, we find that the null hypothesis should be rejected at the 95%–99% confidence level, slightly depending on the angular resolution of the study. If this result is not due to a statistical fluke, it might either indicate that the SN data have not been cleaned from all possible systematics or even point toward new physics. We finally discuss future perspectives in the field for achieving larger and more uniform data sets that will vastly improve the quality of the results and optimally exploit our method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Development of an Uncertainty Quantification Predictive Chemical Reaction Model for Syngas Combustion

An automated data-centric infrastructure, Process Informatics Model (PrIMe), was applied to validation and optimization of a syngas combustion model. The Bound-to-Bound Data Collaboration (B2BDC) module of PrIMe was employed to discover the limits of parameter modifications based on uncertainty quantification (UQ) and consistency analysis of the model–data system and experimental data, including shock-tube ignition delay times and laminar flame speeds. Existing syngas reaction models are reviewed, and the selected kinetic data are described in detail. Empirical rules were developed and applied to evaluate the uncertainty bounds of the literature experimental data. Here, the initial H 2 /CO reaction model, assembled from 73 reactions and 17 species, was subjected to a B2BDC analysis. For this purpose, a dataset was constructed that included a total of 167 experimental targets and 55 active model parameters. Consistency analysis of the composed dataset revealed disagreement between models and data. Further analysis suggested that removing 45 experimental targets, 8 of which were self-inconsistent, would lead to a consistent dataset. This dataset was subjected to a correlation analysis, which highlights possible directions for parameter modification and model improvement. Additionally, several methods of parameter optimization were applied, some of them unique to the B2BDC framework. The optimized models demonstrated improved agreement with experiments compared to the initially assembled model, and their predictions for experiments not included in the initial dataset (i.e., a blind prediction) were investigated. The results demonstrate benefits of applying the B2BDC methodology for developing predictive kinetic models.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2016

TREAT Neutronics Analysis and Design Support, Part I: Multi-SERTTA

Experiment vehicle design is necessary in preparation for Transient Reactor Test (TREAT) facility restart and the resumption of transient testing to support Accident Tolerant Fuel (ATF) characterization and other future fuels testing requirements. Currently the most mature vehicle design is the Multi-SERTTA (Static Environments Rodlet Transient Test Apparatuses), which can accommodate up to four concurrent rodlet-sized specimens under separate environmental conditions. Robust test vehicle design requires neutronics analyses to support design development, optimization of the power coupling factor (PCF) to efficiently maximize energy generation in the test fuel rodlets, and experiment safety analyses. Calculations were performed to support analysis of a near-final design of the Multi-SERTTA vehicle, the design process for future TREAT test vehicles, and establish analytical practices for upcoming transient test experiments. Models of the Multi-SERTTA vehicle containing typical PWR-fuel rodlets were prepared and neutronics calculations were performed using MCNP6.1 with ENDF/B-VII.1 nuclear data libraries. Calculation of the PCF for reference conditions of a PWR fuel rodlet in clean water at operational temperature and pressure provided results between 1.10 and 1.74 W/g-MW depending on the location of the four Multi-SERTTA units with the stack. Basic changes to the Multi-SERTTA secondary vessel containment and support have minimal impact on PCF; using materials with less neutron absorption can improve expected PCF values, especially in the primary containment. An optimized balance is needed between structural integrity, experiment safety, and energy deposition in the experiment. Type of medium and environmental conditions within the primary vessel surrounding the fuel rodlet can also have a significant impact on resultant PCF values. The estimated reactivity insertion worth into the TREAT core is impacted more by the primary and secondary Multi-SERTTA vehicle structure with the experiment content and contained environment having a near negligible impact on overall system reactivity. Additional calculations were performed to evaluate the peak-to-average assembly powers throughout the TREAT core, as well as the nuclear heat generation for the various structural components of the Multi-SERTTA assembly. Future efforts include the evaluation of flux collars to shape the PCF for individual Multi-SERTTA units during an experiment such as to achieve uniformity in test unit environmental conditions impacted by the non-uniform axial flux/power profile of TREAT. Upon resumption of transient testing, experimental results from both the Multi-SERTTA and Multi-SERTTA-CAL will be compared against calculational results and methods for further optimization and design strategies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

TREAT neutronics analysis and design support, Part I: Multi-SERTTA

Experiment vehicle design is necessary in preparation for Transient Reactor Test (TREAT) facility restart and the resumption of transient testing to support Accident Tolerant Fuel (ATF) characterization and other future fuels testing requirements. Currently the most mature vehicle design is the Multi-SERTTA (Static Environments Rodlet Transient Test Apparatuses), which can accommodate up to four concurrent rodlet-sized specimens under separate environmental conditions. Robust test vehicle design requires neutronics analyses to support design development, optimization of the power coupling factor (PCF) to efficiently maximize energy generation in the test fuel rodlets, and experiment safety analyses. Calculations were performed to support analysis of a near-final design of the Multi-SERTTA vehicle, the design process for future TREAT test vehicles, and establish analytical practices for upcoming transient test experiments. Models of the Multi-SERTTA vehicle containing typical PWR-fuel rodlets were prepared and neutronics calculations were performed using MCNP6.1 with ENDF/B-VII.1 nuclear data libraries. Calculation of the PCF for reference conditions of a PWR fuel rodlet in clean water at operational temperature and pressure provided results between 1.10 and 1.73 W/g-MW depending on the location of the four Multi-SERTTA units with the stack. Basic changes to the Multi-SERTTA secondary vessel containment and support have minimal impact on PCF; using materials with less neutron absorption can improve expected PCF values, especially in the primary containment. An optimized balance is needed between structural integrity, experiment safety, and energy deposition in the experiment. Type of medium and environmental conditions within the primary vessel surrounding the fuel rodlet can also have a significant impact on resultant PCF values. The estimated reactivity insertion worth into the TREAT core is impacted more by the primary and secondary Multi-SERTTA vehicle structure with the experiment content and contained environment having a near negligible impact on overall system reactivity. Additional calculations were performed to evaluate the peak-to-average assembly powers throughout the TREAT core, as well as the nuclear heat generation for the various structural components of the Multi-SERTTA assembly. Future efforts include the evaluation of flux collars to shape the PCF for individual Multi-SERTTA units during an experiment such as to achieve uniformity in test unit environmental conditions impacted by the non-uniform axial flux/power profile of TREAT. Upon resumption of transient testing, experimental results from both the Multi-SERTTA and Multi-SERTTA-CAL will be compared against calculational results and methods for further optimization and design strategies. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

In Operando Monitoring of Dynamic Recovery in 3D-Printed Thermoset Nanocomposites by XPCS

Extrusion-based additive manufacturing methods, such as direct-write of carbon fiber-reinforced epoxy inks, have become an attractive route toward development of structural composites in recent years, because of emerging techniques such as big area additive manufacturing. The development of improved materials for these methods has been a major focus area; however, an understanding of the effects of the printing process on the structural and dynamic recovery in printed materials remains largely unexplored. Here, our goal is to capture multiscale and temporal morphology and dynamics within thermosetting composite inks to determine the parameters during the printing process that influence the recovery of the printed material. Herein, we use X-ray photon correlation spectroscopy in small-angle scattering geometry to reveal both morphology and recovery dynamics of a nanoparticle (layered-silicate Cloisite 30B) in a thermoset epoxy resin (EPON 826) during the printing process in real time. Our results show that the dynamics of the layered silicate particles during recovery are anisotropic and slow down to behavior which is characteristic of aging in colloidal clay suspensions around tage ≈ 12 s. The dynamics and alignment of the particles during recovery were tempo-spatially mapped, and the recovery post printing was shown to be strongly influenced by the deposition onto the build plate in addition to the extrusion through the print head. Finally, our in operando results provide insight into the parameters that must be considered when optimizing materials and methods for precisely tailored local properties during 3D printing.

36 MATERIALS SCIENCE↗
DOE OSTI2016

TH-CD-209-06: LET-Based Adjustment of IMPT Plans Using Prioritized Optimization

Purpose: In-vitro experiments suggest an increase in proton relative biological effectiveness (RBE) towards the end of range. However, proton treatment planning and dose reporting for clinical outcome assessment has been based on physical dose and constant RBE. Therefore, treatment planning for intensity-modulated proton therapy (IMPT) is unlikely to transition radically to pure RBE-based planning. We suggest a hybrid approach where treatment plans are initially created based on physical dose constraints and prescriptions, and are subsequently altered to avoid high linear energy transfer (LET) in critical structures while limiting the degradation of the physical dose distribution. Methods: To allow fast optimization based on dose and LET we extended a GPU-based Monte-Carlo code towards providing dose-averaged LET in addition to dose for all pencil beams. After optimizing an initial IMPT plan based on physical dose, a prioritized optimization scheme is used to modify the LET distribution while constraining the physical dose objectives to values close to the initial plan. The LET optimization step is performed based on objective functions evaluated for the product of physical dose and LET (LETxD). To first approximation, LETxD represents a measure of the additional biological dose that is caused by high LET. Regarding optimization techniques, LETxD has the advantage of being a linear function of the pencil beam intensities. Results: The method is applicable to treatments where serial critical structures with maximum dose constraint are located in or near the target. We studied intra-cranial tumors (high-grade meningiomas, base-of-skull chordomas) where the target (CTV) overlaps with the brainstem and optic structures. Often, high LETxD in critical structures can be avoided while minimally compromising physical dose planning objectives. Conclusion: LET-based re-optimization of IMPT plans represents a pragmatic approach to bridge the gap between purely physical dose-based and RBE-based planning. It can be realized with fast GPU-based Monte-Carlo dose calculation. The project was in part supported by NIH grants U19 CA 021239-36 and C06 CA 059267.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

In Situ Observation of Directed Nanoparticle Aggregation During the Synthesis of Ordered Nanoporous Metal in Soft Templates

The prevalent approach to developing new nanomaterials is a trial and error process of iteratively altering synthesis procedures and then characterizing the resulting nanostructures. This is fundamentally limited in that the growth processes that occur during synthesis can only be inferred from the final synthetic structure. Directly observing real-time nanomaterial growth provides unprecedented insight into the relationship between synthesis conditions and product evolution, and facilitates a mechanistic approach to nanomaterial development. Here we use in situ liquid stage scanning transmission electron microscopy to observe the growth of mesoporous palladium in a solvated block copolymer (BCP) template under various synthesis conditions, and ultimately determine a refined synthesis procedure that yields ordered pores. We find that at low organic solvent (tetrahydrofuran, THF) content, the BCP assembles into a rigid, cylindrical micelle array with a high degree of short-range order, but poor long-range order. Upon slowing the THF evaporation rate using a solvent-vapor anneal step, the long-range order is greatly improved. The electron beam induces nucleation of small particles in the aqueous phase around the micelles. The small particles then flocculate and grow into denser structures that surround the micelles, forming an ordered mesoporous structure. The microscope observations revealed that template disorder can be addressed prior to reaction, and is not invariably induced by the growth process itself, allowing us to more quickly optimize the synthetic method. This work was conducted in the William R. Wiley Environmental Molecular Sciences Laboratory (EMSL), a national scientific user facility sponsored by DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory. The Pacific Northwest National Laboratory is operated by Battelle for the U.S. Department of Energy under contract DE-AC05-76RL01830. This research was funded in part by: the Presidential Early Career Award for Scientist and Engineers for I.A., the University of California Academic Senate and the University of California Laboratory fee research grant, the Laboratory-Directed Research and Development program at Sandia National Laboratories, and the Chemical Imaging Initiative at Pacific Northwest National Laboratory. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

cylindrical micelle↗
DOE OSTI2014

Solving security constrained optimal power flow problems by a structure exploiting interior point method

In this paper we present a new approach to solve the DC (n - 1) security constrained optimal power flow (SCOPF) problem by a structure exploiting interior point solver. Our approach is based on a reformulation of the linearised SCOPF model, in which most matrices that need to be factorized are constant. Hence, most factorizations and a large number of back-solve operations only need to be performed once. However, assembling the Schur complement matrix remains expensive in this scheme. To reduce the effort, we suggest using a preconditioned iterative method to solve the corresponding linear system. We suggest two main schemes to pick a good and robust preconditioner based on combining different “active” contingency scenarios of the SCOPF model. These new schemes are implemented within the object-oriented parallel solver (OOPS), a structure-exploiting primal-dual interior-point implementation. We give results on several SCOPF test problems. The largest example contains 500 buses. We compare the results from the original interior point method (IPM) implementation in OOPS and our new reformulation.

30 DIRECT ENERGY CONVERSION↗