Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PROPOSALS”

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 847 records · Page 47

Transactive Control Design for Commercial Buildings to Provide Demand Response

Transactive coordination and control is a new type of distributed strategy that is designed to engage various self-interested distributed energy resources such as controllable loads, distributed generators, and energy storage for providing demand response. In this paper, a novel transactive control design is proposed for commercial buildings to coordinate the electricity demand of the package rooftop air conditioning units (RTUs) for demand response. Such a coordination is realized by introducing a local relationship between the temperature setpoint and the market clearing price. With the proposed market-based control, the demand of individual RTUs is adjusted by changing their temperature setpoints as a response to the market clearing prices. Two dierent scenarios of demand response are considered herein. One is to limit the peak demand of commercial buildings, in which a new market clearing strategy based on proportionalintegral- derivative control with respect to the demand limit is proposed. The other one is to perform price responsive control, in which the total demand of commercial buildings will respond to the real-time pricing. The proposed transactive strategy is also tested and demonstrated on a real occupied building located at Pacic Northwest National Laboratory. The eld testing results conrm the eectiveness of the proposed transactive coordination and control for commercial building with RTUs to provide demand response.

Huang, Sen↗

Consensus Weighting of a Multi-Agent System for Load Shedding

An agent-based scheme is proposed for distributed underfrequency load shedding (UFLS). The key concept is a consensus weighting protocol (CWP) for agents to reach an agreement. Unlike the well adopted average-consensus protocol (ACP), the proposed CWP enables each agent to converge to its weighted portion of the sum of all initial values. Thus, the proposed CWP is more suitable for the UFLS application, as it allows load buses with a higher level of loading to shed more, rather than reducing by the same average value over all load buses as required by the ACP. Proof of convergence for the proposed CWP is derived and presented in this paper, together with the constraints. For implementation of the multi-agent system (MAS), the monitoring, estimation, and distribution steps are developed. Two study cases and the simulation results are provided to validate the performance of the proposed agent-based UFLS scheme.

Xie, Jing↗

DC Link Voltage Balancing of the Active Front-End for the Extreme Fast Charging Stations

One of the most important challenges of the extreme fast charging (XFC) stations is to balance the DC link voltages of the active front-end (AFE). The multilevel cascaded H-bridge (CHB) converter is a great topology for the AFE applications due to the high power density, low total harmonic distortion (THD), and a low number of active and passive components. In this paper, a technique is proposed for balancing DC-link capacitor voltages of a three-phase CHB converter, when a phase shiftPWM modulation technique is used on all cells of the converter. To reach this goal, the redundant states of the CHB converters are used for generating different AC voltage levels by sorting the voltages of the converters. The proposed voltage balancing technique uses the principles that are employed for the DC link voltage balancing of the low-frequency modulation techniques to balance the DC link voltages of the CHB converter. As shown in this paper, the DC link voltages of the CHB converter can accurately be balanced very fast without using PI controllers, which are used in the conventional DC link voltage balancing techniques. A controller also is proposed which will be used to control the dq currents of the CHB converter. To verify the advantages of the proposed technique, the simulation results of the proposed technique will be carried out in a three-phase 7-cell CHB AFE.

Grid-tied converters, cascaded H-bridge, DC link v↗

Hybrid Cyber-attack Detection in Photovoltaic Farms

Here, to address the cyber-physical security in PV farms, a hybrid cyber-attack detection is proposed in this manuscript. To secure PV farms, the proposed method integrates model-based and data-driven methods by fusing the detection score at the device and system levels. First, a model-based cyber-attack detection method is developed for each PV inverter. A residual between the estimation of the Kalman filter and measurement is calculated. By leveraging the calculated residual from all inverters, a squared Mahalanobis distance is developed for device detection score generation. At the system level, a convolutional neural network (CNN) is proposed to detect cyber-attack using the waveform data at the point of common coupling (PCC) in PV farms. To improve the CNN detection accuracy, a set of well-designed features are extracted from the raw waveform data. Finally, a weighted detection score fusion method is proposed to combine device and system detection scores by using their complementary strength. The feasibility and robustness of the proposed method are validated by testing cases and a comparative experiment.

14 SOLAR ENERGY↗

Distributed Energy Resources Coordination over Time-varying Directed Communication Networks

In this paper, we consider an optimal coordination problem for distributed energy resources (DERs) including distributed generators and energy storage devices. We first propose an algorithm based on the push-sum and gradient method to optimally coordinate storage devices and distributed generators in a distributed manner. In the proposed algorithm, each DER only maintains a set of variables and updates them through information exchange with a few neighboring DERs over a time-varying directed communication network. We show that the proposed distributed algorithm with appropriately chosen step-sizes solve the optimal DER coordination problem if the time-varying directed communication network is uniformly jointly strongly connected, which is a mild condition on the connectivity of communication topologies. The proposed distributed algorithm is illustrated and validated by numerical simulations.

Yang, Tao↗

Distance-preserving rigidity penalty on deformable image registration of multiple skeletal components in the neck

Purpose: This study aims at developing and testing a novel rigidity penalty suitable for the deformable registration of tightly located skeletal components in the head and neck from planning computed tomography (CT) and daily cone-beam CT (CBCT) scans of patients undergoing radiotherapy. Methods: The proposed rigidity penalty is designed to preserve intervoxel distances within each bony structure. This penalty was tested in the intensity-based B-spline deformable registration of five cervical vertebral bodies (C1–C5). The displacement vector fields (DVFs) from the registrations were compared to the DVFs generated by using rigid body motions of the cervical vertebrae, measured by the surface registration of vertebrae delineated on CT and CBCT images. Twenty five pairs of planning CT (reference) and treatment CBCTs (target) from five patients were aligned without and with the penalty. An existing penalty based on the orthonormality of the deformation gradient tensor was also tested and the effects of the penalties compared. Results: The mean magnitude of the maximum registration error with the proposed distance-preserving penalty was (0.86, 1.12, 1.33) mm compared to (2.11, 2.49, 2.46) without penalty and (1.53, 1.64, 1.64) with the existing orthonormality-based penalty. The improvement in the accuracy of the deformable image registration was also verified by comparing the Procrustes distance between the DVFs. With the proposed penalty, the average distance was 0.11 (σ 0.03 mm) which is smaller than 0.53 (0.1 mm) without penalty and 0.28 (0.04 mm) with the orthonormality-based penalty. Conclusions: The accuracy of aligning multiple bony elements was improved by using the proposed distance-preserving rigidity penalty. The voxel-based statistical analysis of the registration error shows that the proposed penalty improved the integrity of the DVFs within the vertebral bodies.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

The potential of positron emission tomography for intratreatment dynamic lung tumor tracking: A phantom study

Purpose: This study aims to evaluate the potential and feasibility of positron emission tomography for dynamic lung tumor tracking during radiation treatment. The authors propose a center of mass (CoM) tumor tracking algorithm using gated-PET images combined with a respiratory monitor and investigate the geometric accuracy of the proposed algorithm. Methods: The proposed PET dynamic lung tumor tracking algorithm estimated the target position information through the CoM of the segmented target volume on gated PET images reconstructed from accumulated coincidence events. The information was continuously updated throughout a scan based on the assumption that real-time processing was supported (actual processing time at each frame ≈10 s). External respiratory motion and list-mode PET data were acquired from a phantom programmed to move with measured respiratory traces (external respiratory motion and internal target motion) from human subjects, for which the ground truth target position was known as a function of time. The phantom was cylindrical with six hollow sphere targets (10, 13, 17, 22, 28, and 37 mm in diameter). The measured respiratory traces consisted of two sets: (1) 1D-measured motion from ten healthy volunteers and (2) 3D-measured motion from four lung cancer patients. The authors evaluated the geometric accuracy of the proposed algorithm by quantifying estimation errors (Euclidean distance) between the actual motion of targets (1D-motion and 3D-motion traces) and CoM trajectories estimated by the proposed algorithm as a function of time. Results: The time-averaged error of 1D-motion traces over all trajectories of all targets was 1.6 mm. The error trajectories decreased with time as coincidence events were accumulated. The overall error trajectory of 1D-motion traces converged to within 2 mm in approximately 90 s. As expected, more accurate results were obtained for larger targets. For example, for the 37 mm target, the average error over all 1D-motion traces was 1.1 mm; and for the 10 mm target, the average error over all 1D-motion traces was 2.8 mm. The overall time-averaged error of 3D-motion traces was 1.6 mm, which was comparable to that of the 1D-motion traces. There were small variations in the errors between the 3D-motion traces, although the motion trajectories were very different. The accuracy of the estimates was consistent for all targets except for the smallest. Conclusions: The authors developed an algorithm for dynamic lung tumor tracking using list-mode PET data and a respiratory motion signal, and demonstrated proof-of-principle for PET-guided lung tumor tracking. The overall tracking error in phantom studies is less than 2 mm.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Scatter and crosstalk corrections for {sup 99m}Tc/{sup 123}I dual-radionuclide imaging using a CZT SPECT system with pinhole collimators

Purpose: The energy spectrum for a cadmium zinc telluride (CZT) detector has a low energy tail due to incomplete charge collection and intercrystal scattering. Due to these solid-state detector effects, scatter would be overestimated if the conventional triple-energy window (TEW) method is used for scatter and crosstalk corrections in CZT-based imaging systems. The objective of this work is to develop a scatter and crosstalk correction method for {sup 99m}Tc/{sup 123}I dual-radionuclide imaging for a CZT-based dedicated cardiac SPECT system with pinhole collimators (GE Discovery NM 530c/570c). Methods: A tailing model was developed to account for the low energy tail effects of the CZT detector. The parameters of the model were obtained using {sup 99m}Tc and {sup 123}I point source measurements. A scatter model was defined to characterize the relationship between down-scatter and self-scatter projections. The parameters for this model were obtained from Monte Carlo simulation using SIMIND. The tailing and scatter models were further incorporated into a projection count model, and the primary and self-scatter projections of each radionuclide were determined with a maximum likelihood expectation maximization (MLEM) iterative estimation approach. The extracted scatter and crosstalk projections were then incorporated into MLEM image reconstruction as an additive term in forward projection to obtain scatter- and crosstalk-corrected images. The proposed method was validated using Monte Carlo simulation, line source experiment, anthropomorphic torso phantom studies, and patient studies. The performance of the proposed method was also compared to that obtained with the conventional TEW method. Results: Monte Carlo simulations and line source experiment demonstrated that the TEW method overestimated scatter while their proposed method provided more accurate scatter estimation by considering the low energy tail effect. In the phantom study, improved defect contrasts were observed with both correction methods compared to no correction, especially for the images of {sup 99m}Tc in dual-radionuclide imaging where there is heavy contamination from {sup 123}I. In this case, the nontransmural defect contrast was improved from 0.39 to 0.47 with the TEW method and to 0.51 with their proposed method and the transmural defect contrast was improved from 0.62 to 0.74 with the TEW method and to 0.73 with their proposed method. In the patient study, the proposed method provided higher myocardium-to-blood pool contrast than that of the TEW method. Similar to the phantom experiment, the improvement was the most substantial for the images of {sup 99m}Tc in dual-radionuclide imaging. In this case, the myocardium-to-blood pool ratio was improved from 7.0 to 38.3 with the TEW method and to 63.6 with their proposed method. Compared to the TEW method, the proposed method also provided higher count levels in the reconstructed images in both phantom and patient studies, indicating reduced overestimation of scatter. Using the proposed method, consistent reconstruction results were obtained for both single-radionuclide data with scatter correction and dual-radionuclide data with scatter and crosstalk corrections, in both phantom and human studies. Conclusions: The authors demonstrate that the TEW method leads to overestimation in scatter and crosstalk for the CZT-based imaging system while the proposed scatter and crosstalk correction method can provide more accurate self-scatter and down-scatter estimations for quantitative single-radionuclide and dual-radionuclide imaging.

60 APPLIED LIFE SCIENCES↗

Trust-region based return mapping algorithm for implicit integration of elastic-plastic constitutive models

A new method for the solution of the non-linear equations forming the core of constitutive model integration is proposed. Specifically, the trust-region method that has been developed in the numerical optimization community is successfully modified for use in implicit integration of elastic-plastic models. Although attention here is restricted to these rate-independent formulations, the proposed approach holds substantial promise for adoption with models incorporating complex physics, multiple inelastic mechanisms, and/or multiphysics. As a first step, the non-quadratic Hosford yield surface is used as a representative case to investigate computationally challenging constitutive models. The theory and implementation are presented, discussed, and compared to other common integration schemes. Multiple boundary value problems are studied and used to verify the proposed algorithm and demonstrate the capabilities of this approach over more common methodologies. Robustness and speed are then investigated and compared to existing algorithms. As a result through these efforts, it is shown that the utilization of a trust-region approach leads to superior performance versus a traditional closest-point projection Newton-Raphson method and comparable speed and robustness to a line search augmented scheme.

97 MATHEMATICS AND COMPUTING↗

150-kW Three-Port Custom-Core Transformer Design Methodology

The design of a 150-kW high-frequency three-port transformer (TPT) for a triple active bridge (TAB) converter application is the focus of this paper. A non-tradeoff inclusive design methodology for custom-core transformers is proposed. Classic transformer design methodologies based on complex trade-offs mechanisms require significant effort to adjust the design for commercially available pre-sized cores. The proposed design is based on the development of a set of design equations that allows to select the core dimensions and the turns numbers for the windings that best accomplish the transformer specifications. Experimental evidence verifying the proposed theoretical design methodology is presented in this work.

Christian, Shamar↗

Sizing Energy Storage System for Energy Arbitrage in Extreme Fast Charging Station

This paper proposes a non-linear programming (NLP) model to optimally size the energy storage system (ESS) and obtain an optimal energy management for energy arbitrage of an extreme fast charging station (XFCS) for electric vehicles (EVs), with minimized total cost of XFCS operation and ESS investment. Different from most reported work on sizing the ESS for EV charging stations, this paper proposes a pragmatic approach to model the ESS life degradation and accurately count the ESS cycles. Moreover, this work incorporates the peak demand charges in the operational cost of the charging station which are often overlooked in the literature. The proposed model is formulated and solved using AIMMS. Finally, a thorough sensitivity analysis is performed to offer insights into how different input parameters impact the ESS sizing and savings from the energy arbitrage perspective.

25 ENERGY STORAGE↗

Deadband Voltage Control and Power Buffering for Extreme Fast Charging Station

Voltage fluctuation is one of the most common challenges that electric vehicle charging station (EVCS) may introduce to the power grid. Local reactive power compensation (Q-compensation) capability of bi-directional electric vehicle (EV) chargers can mitigate the steady-state voltage violations caused by the EV charging itself or changes in the neighboring loads. Power buffering, using energy storage system (ESS), can be utilized to address the voltage transients (sags and swells) as a result of EV charging at the EVCS. To address PI controller’s ‘hunting’ issue, this paper proposes a Q-sign triggered deadband voltage control (V-control) method at the point of common coupling (PCC). In addition, to ensure the ramp rate of EV charging is within the allowable limits set forth by the grid code, a ramp rate control is proposed that uses the ESS as a ‘power buffer’. Lastly, different from most reported work in the literature where no explicit limit of the power electronic converters (PECs) is considered, this work considers a reasonable apparent power capacity limit of the PECs when achieving the V-control. This limit also affects the amount of active power that can be obtained from the grid, and subsequently may require ESS to function as ‘load sharing’ device to provide supplemental active power to satisfy EV load. A case study simulated in MATLAB (interfaced with PLECS) is presented to demonstrate the effectiveness of the proposed approaches for EVCS operation.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Combined iterative reconstruction and image-domain decomposition for dual energy CT using total-variation regularization

Purpose: Dual-energy CT (DECT) is being increasingly used for its capability of material decomposition and energy-selective imaging. A generic problem of DECT, however, is that the decomposition process is unstable in the sense that the relative magnitude of decomposed signals is reduced due to signal cancellation while the image noise is accumulating from the two CT images of independent scans. Direct image decomposition, therefore, leads to severe degradation of signal-to-noise ratio on the resultant images. Existing noise suppression techniques are typically implemented in DECT with the procedures of reconstruction and decomposition performed independently, which do not explore the statistical properties of decomposed images during the reconstruction for noise reduction. In this work, the authors propose an iterative approach that combines the reconstruction and the signal decomposition procedures to minimize the DECT image noise without noticeable loss of resolution. Methods: The proposed algorithm is formulated as an optimization problem, which balances the data fidelity and total variation of decomposed images in one framework, and the decomposition step is carried out iteratively together with reconstruction. The noise in the CT images from the proposed algorithm becomes well correlated even though the noise of the raw projections is independent on the two CT scans. Due to this feature, the proposed algorithm avoids noise accumulation during the decomposition process. The authors evaluate the method performance on noise suppression and spatial resolution using phantom studies and compare the algorithm with conventional denoising approaches as well as combined iterative reconstruction methods with different forms of regularization. Results: On the Catphan©600 phantom, the proposed method outperforms the existing denoising methods on preserving spatial resolution at the same level of noise suppression, i.e., a reduction of noise standard deviation by one order of magnitude. This improvement is mainly attributed to the high noise correlation in the CT images reconstructed by the proposed algorithm. Iterative reconstruction using different regularization, including quadratic orq-generalized Gaussian Markov random field regularization, achieves similar noise suppression from high noise correlation. However, the proposed TV regularization obtains a better edge preserving performance. Studies of electron density measurement also show that our method reduces the average estimation error from 9.5% to 7.1%. On the anthropomorphic head phantom, the proposed method suppresses the noise standard deviation of the decomposed images by a factor of ∼14 without blurring the fine structures in the sinus area. Conclusions: The authors propose a practical method for DECT imaging reconstruction, which combines the image reconstruction and material decomposition into one optimization framework. Compared to the existing approaches, our method achieves a superior performance on DECT imaging with respect to decomposition accuracy, noise reduction, and spatial resolution.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

An angle-dependent estimation of CT x-ray spectrum from rotational transmission measurements

Purpose: Computed tomography (CT) performance as well as dose and image quality is directly affected by the x-ray spectrum. However, the current assessment approaches of the CT x-ray spectrum require costly measurement equipment and complicated operational procedures, and are often limited to the spectrum corresponding to the center of rotation. In order to address these limitations, the authors propose an angle-dependent estimation technique, where the incident spectra across a wide range of angular trajectories can be estimated accurately with only a single phantom and a single axial scan in the absence of the knowledge of the bowtie filter. Methods: The proposed technique uses a uniform cylindrical phantom, made of ultra-high-molecular-weight polyethylene and positioned in an off-centered geometry. The projection data acquired with an axial scan have a twofold purpose. First, they serve as a reflection of the transmission measurements across different angular trajectories. Second, they are used to reconstruct the cross sectional image of the phantom, which is then utilized to compute the intersection length of each transmission measurement. With each CT detector element recording a range of transmission measurements for a single angular trajectory, the spectrum is estimated for that trajectory. A data conditioning procedure is used to combine information from hundreds of collected transmission measurements to accelerate the estimation speed, to reduce noise, and to improve estimation stability. The proposed spectral estimation technique was validated experimentally using a clinical scanner (Somatom Definition Flash, Siemens Healthcare, Germany) with spectra provided by the manufacturer serving as the comparison standard. Results obtained with the proposed technique were compared against those obtained from a second conventional transmission measurement technique with two materials (i.e., Cu and Al). After validation, the proposed technique was applied to measure spectra from the clinical system across a range of angular trajectories [−15°, 15°] and spectrum settings (80, 100, 120, 140 kVp). Results: At 140 kVp, the proposed technique was comparable to the conventional technique in terms of the mean energy difference (MED, −0.29 keV) and the normalized root mean square difference (NRMSD, 0.84%) from the comparison standard compared to 0.64 keV and 1.56%, respectively, with the conventional technique. The average absolute MEDs and NRMSDs across kVp settings and angular trajectories were less than 0.61 keV and 3.41%, respectively, which indicates a high level of estimation accuracy and stability. Conclusions: An angle-dependent estimation technique of CT x-ray spectra from rotational transmission measurements was proposed. Compared with the conventional technique, the proposed method simplifies the measurement procedures and enables incident spectral estimation for a wide range of angular trajectories. The proposed technique is suitable for rigorous research objectives as well as routine clinical quality control procedures.

61 RADIATION PROTECTION AND DOSIMETRY↗

Energy Efficient Clothes Dryer with IR Heating and Electrostatic Precipitator

The project goal was to develop a revolutionary energy saving technology for residential clothes drying. The team developed an IR (infrared) heating system and NESP (Nebulizer and Electro-Static Precipitator) for integration into a ventless clothes dryer. The proposed technology addresses two of the major inefficiencies in current electric vented dryers by providing effective energy transfer for the removal of the water and recapture of the vapor latent heat. The IR heaters operating in the mid wave (2.5-10um) are very efficient as they target the 3-micron peak absorption of the water molecule. This allows direct energy absorption, unlike conventional element heaters where heat is transferred by convection. The low power NESP removes water vapor from the exhausted stream and recaptures the latent heat in the ESP (Electro-Static Precipitator) exchanger section. This allows the warm dry air to be recirculated back into the drum for additional efficiency savings. The remaining majority of the dryer hardware stays the same. Summing the efficiency gain from the two subcomponents we anticipated the EF (Efficiency Factor) to exceed the goal of 4.04. EF is obtained by dividing the weight (lbs) of water removed by the energy (kWhr) used, where the test load size is 8.45 lbs of bone dry clothing wetted to 57.5% or 4.8lbs of water, and dried to a remaining moisture content of 2.5-5%. Additional benefits include not having to recondition (heat or cool) the large amounts of make-up air to replace the air exhausted by a vented dryer. It was anticipated that the NESP/heat exchanger would be the most challenging and highest risk element in the program. Therefore, the team focused their efforts during Phase 1 of the program on the design, construction, testing, and optimization of the NESP/heat exchanger. At the end Phase 1, the team compared the performance of the NESP/heat exchanger with the system level requirements and made a Go/No-Go decision on proceeding with the second portion of the program. Phase 2 of the program was structured to develop the IR heating system and then integrate it and the NESP/heat exchanger into a residential clothes dryer prototype for final testing. The proposed technology utilizes heat recovery which is known to have the biggest impact on dryer efficiency. The two current mainstream recovery approaches are air to air exchangers and heat pump condenser systems. Air to air exchanges can be very efficient but require large surface areas which are prone to fouling from uncaptured lint. Dryers based on heat pump condenser recovery systems have shown efficiency improvements of 20–60% and are commercially available. The issue with a heat pump condenser approach is the added cost, as typical prices are twice that of standard vented dryers and they are only available in small to medium capacities. The energy factor (EF) for these systems is 5.50 to 6.88 pounds/kWh compared to conventional dryers at 2.75 to 3.67 pounds/kWh. The efficiency improvements for the proposed technology come from the use of IR heating and the NESP. As the concept is in its infancy, and these improvements were difficult to predict without experimental data, assumptions were made based on available literature. IR radiant drying times, when compared to convection, are typically 30% less. This is a result of the fact that radiant energy heats directly and is absorbed at and below the surface, unlike convection heating, that must conduct the heat through the boundary film of air at the clothes surface and rely on wicking of the moisture to the surface. The second area of improvement comes from the NESP. The NESP operation is as follows: 1. Highly charged, micron sized, droplets of water are injected into the dryer exhaust by the Nebulizer. 2. These charged droplets attract water molecules and continue to grow in size, until losing their charge. During this process, latent heat is rejected back into the air stream. 3. The large droplets enter the ESP, where they are recharged and drawn to the ESP wall, to be extracted at the bottom of the ESP. The warm dry air is then recirculated back into the dryer. The proposed technology, at the time of the proposal submission was estimated to have an EF of 4.79. At program completion we have designed, built, tested, integrated and optimized the proposed technologies into a prototype “high efficiency” residential clothes dryer. Testing of the integrated prototype dryer provided insight into the technical effectiveness and economic feasibility of the proposed technologies. The program has the potential of greatly impacting energy savings. The predicted EF of the proposed technology is 4.79, yielding an energy savings of 42% when compared to conventional electric vented dryers. Given there are approximately 84 million dryers in the US consuming ~64 billion kWh per year, a net savings of 27 billion kWh per year or 0.092 quads could be realized. Since the dryer is not vented, make up air from the room is not needed, adding an additional savings of ~1kWh per load, or 23.8 billion kWh per year. In addition, it is envisioned that the proposed ESP, when successful, would potentially find applications in industrial and residential dehumidification.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

An Integrated Supercritical System for Efficient Produced Water Treatment and Power Generation

The main goal of this project was to evaluate the feasibility of an integrated supercritical (SC) water desalination and power generation system for treatment of produced waters from CO 2 sequestration, oil production, and coal-bed methane recovery. This system may offer a transformative approach to generating power from coal or natural gas (NG) and purifying saline water or produced water with high total dissolved solids (TDS) in a zero-liquid discharge plant. The closed-loop steam cycle of conventional power plants was replaced with an open-loop steam cycle that heated the pretreated produced water above the SC point of water (374 °C and 221 bar) to precipitate and separate the salts dissolved in the water. After separation of the precipitated salt, the SC steam required further treatment to reduce the salt content from ~100 mg/L (ppm) to <1 ppm to make it suitable for use in turbines for power generation. A SC membrane distillation unit was used to polish the SC steam, which could be sent through a series of turbines for power generation, and the steam exiting the low-pressure turbine could be condensed by cooling water to produce pure water. Project tasks included a process simulation and a techno-economic evaluation of the integrated system; design and assembly of SC salt precipitation and membrane distillation systems; development and characterization of advanced carbon membranes for SC membrane distillation; and pretreatment, desalination, and purification of different produced water samples with salt concentrations of 30,000–200,000 ppm to a high-purity (<1 ppm) level. Process simulations were performed for 550 MW SC coal-fired or NG baseline cases and various scenarios for the proposed integrated SC system using produced water as the water source and coal or NG as the fuel source. Simulation results indicated that the net plant efficiency loss of the proposed SC cogeneration process would be in the range of <1% to ~4%, in which NG and lower TDS produced water cases would have higher efficiencies. The proposed process for coal-fired cases could provide all water needs for the power generation process and generate additional pure water for beneficial use. A techno-economic analysis was performed to estimate the cost of produced water treatment (dollars /kgal of treated water) by conventional treatment processes compared with the proposed cogeneration process. The cost estimation for the proposed cogeneration system included the cost of pretreatment, the cost of power loss to integrate the water desalination process into the power generation system, and the cost of plant modification or retrofitting. The cost of produced water treatment was calculated by assuming three different scenarios in which converting a baseline power plant to the proposed cogeneration plant would increase the cost of the boiler and feed water system by 10%, 50%, and 100%. Results indicated that the cost of produced water treatment, even with a 100% cost increase, was in the range of $8.6–11.6 dollars per kgal, which is about half of the treatment cost by conventional methods. The cost of produced water treatment by the proposed cogeneration system could be further reduced or a considerable income could be generated if the values of saltwater and freshwater products are considered. Overall, the techno-economic analysis results indicated that the proposed integrated SC system could potentially treat produced water at a cost significantly lower than the cost of produced water treatment by conventional technologies. To evaluate the feasibility of the proposed treatment system at a laboratory scale, a SC salt precipitation and membrane distillation system was designed and fabricated to treat a simulated brine or pretreated produced water by precipitating salts and further polishing the SC steam containing ~100 ppm salts to a high-purity steam. Through baseline tests with NaCl solutions, with concentrations ranging from 30,000 to 200,000 ppm, it was demonstrated that concentrated feed solutions could be treated by precipitating salt at SC conditions to a TDS level of 130 ppm. Produced water samples were collected from a potential CO 2 sequestration site, a coal-bed methane site, and an oilfield in Illinois. Representative samples were characterized for various water quality parameters, including pH; turbidity; conductivity; and concentrations of selected anions, cations, alkalinity, dissolved organic carbon, total suspended solids, TDS, and total petroleum hydrocarbons. The produced water samples were first treated by conventional pretreatment methods to remove residual oil, dissolved organic matter, and scale-forming species. Pretreatment methods included coagulation, sedimentation, filtration, activated carbon adsorption, and ion-exchange adsorption. Pretreated produced waters were treated at SC conditions in which salts were precipitated. Analysis of treated samples revealed that the concentration of Na (the main cation) in produced waters was reduced dramatically from 14,000–47,400 ppm in the feed to 21–80 ppm in treated waters. Concentrations of other main cations (Ca, Mg, Sr, K) also decreased significantly from initial values of ~1,000–22,000 ppm to ~1 ppm or less. However, a significant increase in the concentrations of Fe, Ni, Cu, and Co was observed, mainly caused by leaching of these elements from the system (vessels, filters, tubing), which is manufactured from alloys containing Fe, Ni, Cu, and Co. To polish the SC steam containing a salt content of ~100 ppm to a high-purity level of <1 ppm, a SC membrane distillation treatment using carbon membranes was employed. Different types of self-standing or supported carbon membranes were developed from graphite, graphene, carbon nanotubes, carbon nanofibers, and pyrolytic carbon. Screening tests revealed that only a limited number of developed carbon membranes were suitable for the purification of steam at SC conditions. Membranes that were either fragile and lacked mechanical strength or had large pores and showed limited salt rejection (<20%) were found to be unsuitable for the purification of SC steam. Promising carbon membranes that were successfully tested for purifying SC steam included two graphite membranes (i.e., AF-5, a porous graphite material, and FGS, a flexible graphite sheet) and a carbon nanotube membrane grown on a Hastelloy substrate and then further modified using a proprietary method (i.e., membrane C). The AF-5 membrane exhibited a stable water flux of ~50 kg/m 2 ∙h and ~40% salt rejection. The FGS membrane had a high salt rejection of 95%–98% but a lower water flux of ~15 kg/m 2 ∙h. Membrane C showed an extremely high flux of 200-800 kg/m 2 ∙h, which is one order of magnitude higher than the water flux of conventional membranes. Membrane C also exhibited 87%–93% salt rejection. Overall, a two-stage or multistage polishing treatment by these membranes could lower the salt content of the SC steam from initial values of ~100–200 ppm to the steam quality of <1 ppm or to parts per billion levels required for the power cycle.

20 FOSSIL-FUELED POWER PLANTS↗

Disjunctive linear separation conditions and mixed-integer formulations for aircraft conflict resolution

In this paper, we address the aircraft conflict resolution problem in air traffic control. We introduce new mixed-integer programming formulations for aircraft conflict resolution with speed, heading and altitude control which are based on disjunctive linear separation conditions. We first examine the two-dimensional aircraft conflict resolution problem with speed and heading control represented as continuous decision variables. We show that the proposed disjunctive linear separation conditions are equivalent to the classical nonlinear conditions for aircraft separation. Further, we characterise conflict-free trajectories based on aircraft velocity bounds and propose a simple pre-processing algorithm to identify aircraft pairs which are either always conflict-free, or which cannot be separated using speed and heading control only. We then incorporate altitude control and propose a lexicographic optimisation formulation that aims to minimise the number of flight level changes before resolving outstanding conflicts via two-dimensional velocity control. The proposed mixed-integer programming formulations are nonconvex, and we propose convex relaxations, decomposition methods and constraint generation algorithms to solve the two-dimensional and lexicographic optimisation formulations to guaranteed optimality. Numerical experiments on four types of conflict resolution benchmarking instances are conducted to test the performance of the proposed mixed-integer formulations. Further, the proposed method is compared against two benchmarks based on state-of-the-art approaches for the aircraft conflict resolution problem. Our numerical results show that the proposed method largely outperforms both benchmarks in terms of runtime and is able to solve significantly more instances to global optimality.

97 MATHEMATICS AND COMPUTING↗

Precise ROCOF estimation algorithm for low inertia power grids

The precise estimation of Rate Of Change Of Frequencies (ROCOFs) in a generation trip event can be helpful on power system inertia estimation, fast system response, and accurate event capturing. However, the ROCOF estimations from the existing Synchronized Measurement Devices (SMDs) are usually simply calculated by the finite difference between two adjacent frequency measurement points. The ineluctable noises, disturbances, and spikes from real-world frequency measurements can bring large dynamics to the ROCOF estimation and thus can result in an inaccurate estimation of the initial ROCOF. This issue becomes more serious when the target power grid has low inertia, and a large amount of distributed energy sources are deployed. To address this issue, a precise ROCOF estimation algorithm is designed based on least square method with flexible window size. In addition, a median filter is also designed and applied on the frequency measurements before using the proposed algorithm. The window size and thresholds in the proposed algorithm are determined with historical event data analysis. The proposed algorithm is deployed in a low cost, flexible, and distribution level universal grid analyzer (UGA) platform. In conclusion, multiple experiments are conducted in both a laboratory and the Hawaiian Islands to verify the effectiveness of the proposed algorithm.

24 POWER TRANSMISSION AND DISTRIBUTION↗