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,333 records · Page 74

Strengthening the SDP Relaxation of AC Power Flows with Convex Envelopes, Bound Tightening, and Valid Inequalities

Here this work revisits the Semidefine Programming (SDP) relaxation of the AC power flow equations in light of recent results illustrating the benefits of bounds propagation, valid inequalities, and the Convex Quadratic (QC) relaxation. By integrating all of these results into the SDP model a new hybrid relaxation is proposed, which combines the benefits from all of these recent works. This strengthened SDP formulation is evaluated on 71 AC Optimal Power Flow test cases from the NESTA archive and is shown to have an optimality gap of less than 1% on 63 cases. This new hybrid relaxation closes 50% of the open cases considered, leaving only 8 for future investigation.

42 ENGINEERING↗

Optimizing DER Participation in Inertial and Primary-Frequency Response

This paper develops an approach to enable the optimal participation of distributed energy resources (DERs) in inertial and primary-frequency response alongside conventional synchronous generators. Leveraging a reduced-order model description of frequency dynamics, DERs' synthetic inertias and droop coefficients are designed to meet time-domain performance objectives of frequency overshoot and steady-state regulation. Furthermore, an optimization-based method centered around classical economic dispatch is developed to ensure that DERs share the power injections for inertial- and primary-frequency response in proportion to their power ratings. Simulations for a modified New England test-case system composed of ten synchronous generators and six instances of the IEEE 37-node test feeder with frequency-responsive DERs validate the design strategy.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Global Solution Strategies for the Network-Constrained Unit Commitment Problem with AC Transmission Constraints

Here, we propose a novel global solution algorithm for the network-constrained unit commitment problem that incorporates a nonlinear alternating current (AC) model of the transmission network, which is a nonconvex mixed-integer nonlinear programming (MINLP) problem. Our algorithm is based on the multi-tree global optimization methodology, which iterates between a mixed-integer lower-bounding problem and a nonlinear upper-bounding problem. We exploit the mathematical structure of the unit commitment problem with AC power flow constraints (UC-AC) and leverage second-order cone relaxations, piecewise outer approximations, and optimization-based bounds tightening to guarantee a globally optimal solution at convergence. Numerical results on four benchmark problems illustrate the effectiveness of our algorithm, both in terms of convergence rate and solution quality.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Multi-operator Production Cost Modeling

As contiguous power systems are often operated by several regional system operators, optimization models have traditionally been designed to generate optimal system-wide scheduling decisions. This report introduces an approach to reflect the existence of multiple system operators making scheduling decisions within an interconnected power system. The approach is explicitly designed to reflect each operator's limited knowledge about other regions when making scheduling decisions. Detailed analysis using a modified IEEE Reliability Test System reveals that the approach mimics real system operations by providing imperfect information to each operator about external operations and external constraints. The approach scales to a 78,000-bus model of the Eastern Interconnection, where a regression analysis provides quantitative insight into which constraints cause inefficiencies when improperly managed in a multi-operator simulation. Because the approach generates distinct regional scheduling problems, unit commitment and economic dispatch simulations can be applied to individual regions and solved in parallel, enabling production cost simulations at a previously intractable scale.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Failure Probability Constrained AC Optimal Power Flow

Despite cascading failures being the central cause of blackouts in power transmission systems, existing operational and planning decisions are made largely by ignoring their underlying cascade potential. This paper posits a reliability-aware AC Optimal Power Flow formulation that seeks to design a dispatch point which has a low operator-specified likelihood of triggering a cascade starting from any single component outage. By exploiting a recently developed analytical model of the probability of component failure, our Failure Probability-constrained ACOPF (FP-ACOPF) utilizes the system's expected first failure time as a smoothly tunable and interpretable signature of cascade risk. Here, we use techniques from bilevel optimization and numerical linear algebra to efficiently formulate and solve the FP-ACOPF using off-the-shelf solvers. Extensive simulations on the IEEE 118-bus case show that, when compared to the unconstrained and N-1 security-constrained ACOPF, our probability-constrained dispatch points can significantly lower the probabilities of long severe cascades and of large demand losses, while incurring only minor increases in total generation costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Robust Recovery of Temporally Smooth Signals From Under-Determined Multiple Measurements

In this study, we consider the problem of recovering jointly sparse vectors from underdetermined measurements that are corrupted by both additive noise and outliers. This can be viewed as the robust extension of the Multiple Measurement Vector (MMV) problem. To solve this problem, we propose two general approaches. As a benchmark, the first approach preprocesses the input for outlier removal and then employs state-of-the-art technologies for signal recovery. The second approach, as the main contribution of this paper, is based on formulation of an innovative regularized fitting problem. By solving the regularized fitting problem, we jointly remove outliers and recover the sparse vectors. Furthermore, by exploiting temporal smoothness among the sparse vectors, we improve noise robustness of the proposed approach and avoid the problem of over-fitting. Extensive numerical results are provided to illustrate the excellent performance of the proposed approach.

47 OTHER INSTRUMENTATION↗

Route-Based Online Energy Management of a PHEV and Sensitivity to Trip Prediction

In this paper, we present a method of optimizing the energy management of a plug-in hybrid electric vehicle (PHEV) using GIS-assisted stochastic trip prediction. A process was developed to synthesize speed profiles through a combination of Markov chains and information from a geographical information system (GIS) about the future route. In a potential real-world scenario, the future trip (speed, grade, stops, etc.) can be estimated, but not deterministically known. The stochastic trip prediction process models such uncertainty. The route-based energy management presented in this paper uses the Pontryagin Minimum Principle (PMP). A PMP strategy was implemented in a Simulink controller for a model of Prius-like PHEV and compared to a baseline strategy using Autonomie, an automotive modeling environment. An itinerary was defined, and several speed profiles were synthesized. It was then possible to evaluate the sensitivity of PMP tuning to the speed profile, providing insights about the applicability of PMP control in real-world situations.

energy mangement↗

Experiment on the treatment of waste extraction solvent from the molybdenum-99 process

In the Mo-99 (Molybdenum-99) isotope extraction test process for radiopharmaceutical applications, organic solvent is used to extract Mo-99 from an irradiated UO{sub 2} dissolution. The extraction solvent was stored when the test work was stopped. A total of about 120 liters of waste solvent was stored at INER (Institute of Nuclear Energy Research, Taiwan). The extraction solvent consisted of 5% di-(2-ethylhexyl)-phosphoric acid (D2EHPA) and kerosene. The radionuclides found in the waste solvent include Cs-137, Am-241, Tc-99, and Sr-90, which give off gross alpha and beta radioactivity of 1898 and 471 Bq/ml, respectively. This study aims to remove radionuclides from the waste solvent using sodium carbonate and sodium hydroxide solutions in different concentrations. After mixing the waste solvent with the alkaline solution followed by settling, a third phase other than organic and aqueous phase appeared which is expected due to the saponification reaction. The experimental results showed that increasing the number of washing and the alkaline solution concentration could enhance the radionuclides removal rate. An optimal removal method was proposed using 2 M Na{sub 2}CO{sub 3} solution twice followed by 1 M NaOH solution one time for the third phase generated early in the mixing stages. The remaining gross alpha and beta radioactivity of the treated organic solvent was 2 and 3 Bq/ml, respectively. The treated solvent could be stabilized by ashing at 500 deg. C and then immobilized. The alkaline solution would be neutralized by hydrochloric or nitric acid and then treated using a variety of adsorbents or bone char via adsorption to remove nuclides to meet the wastewater discharge limitation. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Training data selection for event classification in a highly variable environment

A problem of interest for nuclear nonproliferation is monitoring activities at nuclear facilities, where proliferation events may only take place a few times and often under variable conditions. Machine learning has revolutionized data analytics by enabling the use of measurable signatures to generate predictive models of facility operations. However, traditional methods for training these models require large, reliable data sets with labeled observations, a challenge for nonproliferation. Highly variable conditions further complicate this as events from training data may have occurred in conditions quite different from the event of interest. Our hypothesis is that when events occur in a highly variable environment, careful training data selection for each test event could outperform the standard approach of using all available training data. We developed a method to optimize training data selection for the given test event and applied it to predicting the power level of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory. In this study, the reactor startup exhibits variability between occurrences due to natural variability in environmental conditions and operational procedures. Using a combination of analysis techniques, a similitude assessment was performed on data collected from HFIR to isolate clusters that were optimal for training a predictive model. Concepts such as dynamic time warping and Jaccard similarity were used in conjunction with clustering analysis. In order to validate this approach, the model was trained on every combination of unique training events and the predictive performance was compared to the performance using a subset of the training data selected by isolated clusters found through the similitude assessment.

Iyer, A↗

Towards local progression estimation of pulmonary emphysema using CT

Purpose: Whole lung densitometry on chest CT images is an accepted method for measuring tissue destruction in patients with pulmonary emphysema in clinical trials. Progression measurement is required for evaluation of change in health condition and the effect of drug treatment. Information about the location of emphysema progression within the lung may be important for the correct interpretation of drug efficacy, or for determining a treatment plan. The purpose of this study is therefore to develop and validate methods that enable the local measurement of lung density changes, which requires proper modeling of the effect of respiration on density. Methods: Four methods, all based on registration of baseline and follow-up chest CT scans, are compared. The first naïve method subtracts registered images. The second employs the so-called dry sponge model, where volume correction is performed using the determinant of the Jacobian of the transformation. The third and the fourth introduce a novel adaptation of the dry sponge model that circumvents its constant-mass assumption, which is shown to be invalid. The latter two methods require a third CT scan at a different inspiration level to estimate the patient-specific density-volume slope, where one method employs a global and the other a local slope. The methods were validated on CT scans of a phantom mimicking the lung, where mass and volume could be controlled. In addition, validation was performed on data of 21 patients with pulmonary emphysema. Results: The image registration method was optimized leaving a registration error below half the slice increment (median 1.0 mm). The phantom study showed that the locally adapted slope model most accurately measured local progression. The systematic error in estimating progression, as measured on the phantom data, was below 2 gr/l for a 70 ml (6%) volume difference, and 5 gr/l for a 210 ml (19%) difference, if volume correction was applied. On the patient data an underlying linearity assumption relating lung volume change with density change was shown to hold (fitR{sup 2} = 0.94), and globalized versions of the local models are consistent with global results (R{sup 2} of 0.865 and 0.882 for the two adapted slope models, respectively). Conclusions: In conclusion, image matching and subsequent analysis of differences according to the proposed lung models (i) has good local registration accuracy on patient data, (ii) effectively eliminates a dependency on inspiration level at acquisition time, (iii) accurately predicts progression in phantom data, and (iv) is reasonably consistent with global results in patient data. It is therefore a potential future tool for assessing local emphysema progression in drug evaluation trials and in clinical practice.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-E-T-542: Comparison of Stereotactic Radiosurgery (SRS) of Brain Lesions Using Gamma Knife, VMAT, IMRT, and Conformal Arcs

Purpose: To compare dose distributions in stereotactic radiation surgery of brain lesions using gamma Knife, VMAT, conformal arcs, and IMRT in order to provide an optimal treatment. Methods: Dose distributions from single shot of 4C model of Gamma Knife at the helmet collimation sizes of 4, 8, 14, and 18 mm in diameter were compared with full arcs with the square shapes of 4×4 (or 5×5), 8×8 (or 10×10), and spherical shapes of 16 or 20 mm in diameter using EDR3 films in the same gamma knife QA phantom. Plans for ten SRS cases with single and multiple lesions were created in gamma knife plans and Pinnacle plans. The external beam plans had enlarged field size by 2-mm and used single conformal full circle arc for solitary lesion and none coplanar arcs/beams for multiple lesions. Coverage, conformity index, dose to critical organs, and integral dose to the brain and nearby critical structures were compared on all plans. Structures and dose matrices were registered in a Velocity deformable image registration system. Results: Single full circle arc from Elekta beam-modulate MLC (4-mm leaf thickness) and agility MLC (5-mm leaf thickness) have larger penumbra and less flatness than that of Gamma Knife single shot. None-coplanar arcs or beams were required to achieve similar dose distribution. In general, Gamma Knife plans provided significant less integral dose than that of linac-based plans. Benefits of IMRT and VMAT versus gamma Knife and conformal arcs were not significant. Conclusion: Our dose measurement and treatment planning evaluation clearly demonstrated dose distribution differences amount current popular SRS modalities for small solitary and multiple brain lesions. The trend of using MLC shape beams or arcs to replace conventional cones should be revisited in order to keep lower integral dose if the late correlates with some radiation-induced side effects. Pilot grant from Elekta LLC.

60 APPLIED LIFE SCIENCES↗

SU-E-T-192: Commissioning of a Commercial 3D Dose Calculation Program

Purpose: To commission a commercial software package (CSP) that is used as secondary dose calculation check. The CSP uses an independent golden data beam model. However, some parameters can be modified to generate a customer specific model. Plan comparisons and point dose measurements were performed to test if and to what extent the beam model needed adjustment to optimize results. Methods: Beam parameter configurations were compared between the CSP and both TPS. Twelve phantom test plans ranging from simple to complex were generated in two treatment planning systems (TPS). Tests included small field, off axis, EDW, IMRT and VMAT plans. For each plan a point dose was measured to establish ground truth. Lastly, patient plans were compared for both TPS systems and the CSP. Results: Beam parameters agreed within 2%. The output factors for small fields were changed for the 15 MV beam by 2 and 1.5 % for the 1 cm and 2 cm field sizes, respectively. For the 6 MV beam output factors were adjusted by 3−0.8% for field sizes ranging from 1 to 5 cm. The MLC dynamic leaf gap was adjusted by 1.5 mm for 18 MV beam. Differences between the CSP and the TPS were noted in the built-up region. These differences affected the gamma pass rate in the surface region, however this effect is reduced with increasing number of beam angles and does not affect point dose calculations at depth. All IMRT and VMAT plans agreed with the CSP using a gamma pass rate of 95% (3%, 3mm). Conclusion: The CSP is used to verify point doses for all 3D plans generated in our clinic for the last 6 months. No point dose mismatches were encountered since the CSP was implemented. Next, the CSP will be adapted for secondary checks of all IMRT plans. KL had a beta tester agreement with Mobius Medical for an in-kind equipment and software loan.

60 APPLIED LIFE SCIENCES↗

SU-D-207A-04: Use of Gradient Echo Plural Contrast Imaging (GEPCI) in MR-Guided Radiation Therapy: A Feasibility Study Targeting Brain Treatment

Purpose: To implement the Gradient Echo Plural Contrast Imaging(GEPCI) technique in MRI-simulation for radiation therapy and assess the feasibility of using GEPCI images with advanced inhomogeneity correction in MRI-guided radiotherapy for brain treatment. Methods: An optimized multigradient-echo GRE sequence (TR=50ms;TE1=4ms;delta-TE=4ms;flip angle=300,11 Echoes) was developed to generate both structural (T1w and T2*w) and functional MRIs (field and susceptibility maps) from a single acquisition. One healthy subject (Subject1) and one post-surgical brain cancer patient (Subject2) were scanned on a Philips Ingenia 1.5T MRI used for radiation therapy simulation. Another healthy subject (Subject3) was scanned on a 0.35T MRI-guided radiotherapy (MR-IGRT) system (ViewRay). A voxel spread function (VSF) was used to correct the B0 inhomogeneities caused by surgical cavities and edema for Subject2. GEPCI images and standard radiotherapy planning MRIs for this patient were compared focusing the delineation of radiotherapy target region. Results: GEPCI brain images were successfully derived from all three subjects with scan times of <7 minutes. The images derived for Subjects1&2 demonstrated that GEPCI can be applied and combined into radiotherapy MRI simulation. Despite low field, T1-weighted and R2* images were successfully reconstructed for Subject3 and were satisfactory for contour and target delineation. The R2* distribution of grey matter (center=12,FWHM=4.5) and white matter (center=14.6, FWHM=2) demonstrated the feasibility for tissue segmentation and quantification. The voxel spread function(VSF) corrected surgical site related inhomogeneities for Subject2. R2* and quantitative susceptibility map(QSM) images for Subject2 can be used to quantitatively assess the brain structure response to radiation over the treatment course. Conclusion: We implemented the GEPCI technique in MRI-simulation and in MR-IGRT system for radiation therapy. The images demonstrated that it is feasible to adopt this technique in radiotherapy for structural delineation. The preliminary data also enable the opportunity for quantitative assessment of radiation response of the target region and normal tissue.

60 APPLIED LIFE SCIENCES↗

SU-F-J-08: Quantitative SPECT Imaging of Ra-223 in a Phantom

Purpose: Ra-223 therapy of prostate cancer bone metastases is being used to treat patients routinely. However, the absorbed dose distribution at the macroscopic and microscopic scales remains elusive, due to the inability to image the small activities injected. Accurate activity quantification through imaging is essential to calculate the absorbed dose in organs and sub-units in radiopharmaceutical therapy, enabling personalized absorbed dose-based treatment planning methodologies and more effective and optimal treatments. Methods: A 22 cm diameter by 20 cm long cylindrical phantom, containing a 3.52 cm diameter sphere, was used. A total of 2.01 MBq of Ra-223 was placed in the phantom with 177.6 kBq in the sphere. Images were acquired on a dual-head Siemens Symbia T16 gamma camera using three 20% full-width energy windows and centered at 84, 154, and 269 keV (120 projections, 360° rotation, 45 s per view). We have implemented reconstruction of Ra-223 SPECT projections using OS-EM (up to 20 iterations of 10 subsets) with compensation for attenuation using CT-based attenuation maps, collimator-detector response (CDR) (including septal penetration, scatter and Pb x-ray modeling), and scatter in the patient using the effective source scatter estimation (ESSE) method. The CDR functions and scatter kernels required for ESSE were computed using the SIMIND MC simulation code. All Ra-223 photon emissions as well as gamma rays from the daughters Rn-219 and Bi-211 were modeled. Results: The sensitivity of the camera in the three combined windows was 107.3 cps/MBq. The visual quality of the SPECT images was reasonably good and the activity in the sphere was 27% smaller than the true activity. This underestimation is likely due to partial volume effect. Conclusion: Absolute quantitative Ra-223 SPECT imaging is achievable with careful attention to compensate for image degrading factors and system calibration.

60 APPLIED LIFE SCIENCES↗

SU-G-IeP3-07: High-Resolution, High-Sensitivity Imaging and Quantification of Intratumoral Distributions of Gold Nanoparticles Using a Benchtop L-Shell XRF Imaging System

Purpose: To demonstrate the ability to perform high-resolution imaging and quantification of sparse distributions of gold nanoparticles (GNPs) within ex vivo tumor samples using a highly-sensitive benchtop L-shell x-ray fluorescence (XRF) imaging system. Methods: An optimized L-shell XRF imaging system was assembled using a tungsten-target x-ray source (operated at 62 kVp and 45 mA). The x-rays were filtered (copper: 0.08 mm & aluminum: 0.04 mm) and collimated (lead: 5 cm thickness, 3 cm aperture diameter) into a cone-beam in order to irradiate small samples or objects. A collimated (stainless steel: 4 cm thickness, 2 mm aperture diameter) silicon drift detector, capable of 2D translation, was placed at 90° with respect to the beam to acquire XRF/scatter spectra from regions of interest. Spectral processing involved extracting XRF signal from background, followed by attenuation correction using a Compton scatter-based normalization algorithm. Calibration phantoms with water/GNPs (0 and 0.00001–10 mg/cm{sup 3}) were used to determine the detection limit of the system at a 10-second acquisition time. The system was then used to map the distribution of GNPs within a 12×11×2 mm{sup 3} slice excised from the center of a GNP-loaded ex vivo murine tumor sample; a total of 110 voxels (2.65×10{sup −3} cm{sup 3}) were imaged with 1.3-mm spatial resolution. Results: The detection limit of the current cone-beam benchtop L-shell XRF system was 0.003 mg/cm{sup 3} (3 ppm). Intratumoral GNP concentrations ranging from 0.003 mg/cm{sup 3} (3 ppm) to a maximum of 0.055 mg/cm{sup 3} (55 ppm) and average of 0.0093 mg/cm{sup 3} (9.3 ppm) were imaged successfully within the ex vivo tumor slice. Conclusion: The developed cone-beam benchtop L-shell XRF imaging system can immediately be used for imaging of ex vivo tumor samples containing low concentrations of GNPs. With minor finetuning/optimization, the system can be directly adapted for performing routine preclinical in vivo imaging tasks. Supported by NIH/NCI grant R01CA155446 This investigation was supported by NIH/NCI grant R01CA155446.

60 APPLIED LIFE SCIENCES↗

Formation of correlated states and tunneling for a low energy and controlled pulsed action on particles

We consider a method for optimizing the tunnel effect for low-energy particles by using coherent correlated states formed under controllable pulsed action on these particles. Typical examples of such actions are the effect of a pulsed magnetic field on charged particles in a gas or plasma. Coherent correlated states are characterized most comprehensively by the correlation coefficient r(t); an increase of this factor elevates the probability of particle tunneling through a high potential barrier by several orders of magnitude without an appreciable increase in their energy. It is shown for the first time that the formation of coherent correlated states, as well as maximal |r(t)|{sub max} and time-averaged 〈|r(t)|〉 amplitudes of the correlation coefficient and the corresponding tunneling probability are characterized by a nonmonotonic (oscillating) dependence on the forming pulse duration and amplitude. This result makes it possible to optimize experiments on the realization of low-energy nuclear fusion and demonstrates the incorrectness of the intuitive idea that the tunneling probability always increases with the amplitude of an external action on a particle. Our conclusions can be used, in particular, for explaining random (unpredictable and low-repeatability) experimental results on optimization of energy release from nuclear reactions occurring under a pulsed action with fluctuations of the amplitude and duration. We also consider physical premises for the observed dependences and obtain optimal relations between the aforementioned parameters, which ensure the formation of an optimal coherent correlated state and optimal low-energy tunneling in various physical systems with allowance for the dephasing action of a random force. The results of theoretical analysis are compared with the data of successful experiments on the generation of neutrons and alpha particles in an electric discharge in air and gaseous deuterium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Properties and parameters of the electron beam injected into the mirror magnetic trap of a plasma accelerator

The parameters of the injector of an axial plasma beam injected into a plasma accelerator operating on the basis of gyroresonance acceleration of electrons in the reverse magnetic field are determined. The trapping of the beam electrons into the regime of gyroresonance acceleration is numerically simulated by the particle- in-cell method. The optimal time of axial injection of the beam into a magnetic mirror trap is determined. The beam parameters satisfying the condition of efficient particle trapping into the gyromagnetic autoresonance regime are found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗