Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “loading pattern optimization”

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.

156 records · Page 9

Analyzing the Effects of Atmospheric Turbulent Fluctuations on the Wake Structure of Wind Turbines and Their Blade Vibrational Dynamics

In recent trends, a rising demand for renewable energy has driven wind turbines to larger proportions, where lighter blade designs are often adopted to reduce the costs associated with logistics and production. This causes modern utility-scale wind turbine blades to be inherently more flexible, and their amplified aeroelastic sensitivity results in complex multi-physics reactions to variant atmospheric conditions, including dynamic patterns of aerodynamic loading at the rotor and vortex structure evolutions within the wake. In this paper, we analyze the influence of inflow variance for wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex wake interactions for a range of flexible NRT blade variations, operating in differing conditions of variant atmospheric flow. CODEF solutions of turbine operation in Steady-In-The-Average (SITA) wind conditions are compared to SITA wind conditions featuring a controlled gust-like pulse overimposed, to isolate the effects of typical wind fluctuations. Finally, simulations of realistic time-varying wind conditions from SWiFT meteorological tower measurements are compared to the solutions of SITA wind conditions. These increasingly complex atmospheric inflow variations are tested to show the differing effects evoked by various patterns of spatiotemporal atmospheric flow fluctuations. An analysis is presented for solutions of wind turbine aeroelastic response and vortex wake evolution, to elucidate the consequences of variant inflow, which pertain to wind turbine dynamics at an individual and farm-collective scale. The comparisons of simulated farm flow for SITA and measured fluctuating wind conditions show that certain regions of the wake contain up to a 12% difference in normalized axial velocity, due to the introduction of wind fluctuations. The findings of this study prove valuable for practical applications in wind farm control and optimization strategies, with particular significance for modern utility-scale wind power plants operating in variant atmospheric conditions.

Farrell, Alayna (ORCID:000000023555720X)↗

Next Generation Hydrogen Storage Vessels Enabled by Carbon Fiber Infusion with a Low Viscosity, High Toughness System

Onboard hydrogen storage for fuel cell vehicles holds considerable commercial promise which could be accelerated by lower cost storage vessels. In past years, the use of lower cost reinforcing fiber had been studied as an intuitive avenue to reduce vessel costs. Although the matrix resin does not play a large role in bearing stress loads in many areas of a filament wound pressure vessel, the matrix is relied upon for stress transfer between fibers and appears to play an important role in certain regions of the vessel construction such as the “tangent region” where the cylinder section transitions to the dome. Also, the matrix resin can play a role in the damage tolerance of vessels after impact events. Furthermore, the effect of dry spots or voids in the laminate can affect stress transfer between fibers in a manner that is challenging to quantify. Therefore, excess fiber must be used to account for this variability in performance. Unfortunately, in the case hydrogen pressure vessels the cost of the vessel is strongly dependent on the amount of carbon-fiber used in the design and manufacturing. An extensive investigation on the use of non-traditional fabrication processes enabled by an ultra-low viscosity thermoset resin (<10 cP) with high fracture toughness was conducted to determine the opportunity to reduce carbon-fiber usage and the associated high cost. The hypothesis was based on enabling a lower variation in vessel performance from a reduction in void defects (often found with traditional wet-winding) and a more defect-tolerant vessel resulting from use of a matrix resin with high fracture toughness. These two features were anticipated to justify a more aggressive design approach to the vessel winding pattern. From the work in this project, the potential use of the vacuum infusion approach to produce high quality vessels was validated. However, the team recognized the need to re-optimize the design, the dry-process, and the vacuum infusion process at each scale-up stage. During the planning phase, the critical need to re-optimize at this frequency was under-estimated. Also, the team identified the lack of industry testing standards on subscale parts for providing predictive value of vessel performance. Therefore, significant effort is needed to produce statistically significant numbers of vessels. To produce such numbers of vessels the project needed in retrospect to partner closely with a series vessel manufactures to develop a proof of concept production line to make the vessels. However, without the appropriate data, a manufacturing partner was challenging to fully engage during this timeline. New technologies relying on changing raw materials and manufacturing processes carry inherent risks. This project did not reach its ultimate objectives of demonstrated a vessel with lower carbon-fiber. But some companies expressed interest in learning from the partial successes of this project for future challenges in pressure vessel applications.

08 HYDROGEN↗

Improving I/O Performance for Exascale Applications through Online Data Layout Reorganization

The applications being developed within the U.S. Exascale Computing Project (ECP) to run on imminent Exascale computers will generate scientific results with unprecedented fidelity and record turn-around time. Many of these codes are based on particle-mesh methods and use advanced algorithms, especially dynamic load-balancing and mesh-refinement, to achieve high performance on Exascale machines. Yet, as such algorithms improve parallel application efficiency, they raise new challenges for I/O logic due to their irregular and dynamic data distributions. Thus, while the enormous data rates of Exascale simulations already challenge existing file system write strategies, the need for efficient read and processing of generated data introduces additional constraints on the data layout strategies that can be used when writing data to secondary storage. We review these I/O challenges and introduce two online data layout reorganization approaches for achieving good tradeoffs between read and write performance. We demonstrate the benefits of using these two approaches for the ECP particle-in-cell simulation WarpX, which serves as a motif for a large class of important Exascale applications. Here, we show that by understanding application I/O patterns and carefully designing data layouts we can increase read performance by more than 80 percent.

97 MATHEMATICS AND COMPUTING↗

Coiled Tube Gas Heaters For Nuclear Gas-Brayton Power Conversion

This project developed an alternative design for heat exchangers for application to heating supercritical carbon dioxide (S-CO 2 ) or air for power conversion. We have identified an annular coiled tube bundle configuration–where hot sodium enters tubes from multiple vertical inlet manifold pipes, flows in a spiral pattern radially inward and downward, and then exits into an equal number of vertical outlet manifold pipes–as a potentially attractive option. The S-CO 2 gas or air flows radially outward through the tube bundle. Coiled tube gas heaters (CTGHs) are expected to have excellent thermal shock, long-term thermal creep, in-service inspection, and reparability characteristics, compared to alternative options. CTGHs have significant commonality with modern nuclear steam generators. Extensive experience exists with the design, manufacture, operation, in-service inspection and maintenance of nuclear steam generators. The U.S. Nuclear Regulatory Commission also has extensive experience with regulatory guidance documented in NUREG 0800. CTGHs leverage this experience and manufacturing capability. The most important difference between steam generators and gas-Brayton cycles such as the S-CO 2 cycle is that the heat exchangers must operate with counter flow with high effectiveness to minimize the pinch-point temperature difference between the hot liquid coolant and the heated gas. S-CO 2 -cycle gas heaters also operate at sufficiently elevated temperatures that time dependent creep is important and allowable stresses are relatively low. Designing heat exchangers to operate in this regime requires configurations that minimize stresses and stress concentrations. The cylindrical tubes and cylindrical manifold pipes used in CTGHs are particularly effective geometries. The first major goal of this research project was to develop and experimentally validate a detailed, 3-D multi-phase (gas-solid-liquid) heat transport model for CTGHs, using methods similar to earlier UCB multi-scale models for PCHEs, which will enable optimization of CTGH designs with respect to the number of manifold pipes, S-CO 2 and sodium circulating power, and other design parameters. CTGHs are relatively novel because the tubes loaded in compression as well as the manifold pipes which form the tube sheets. The second major goal of the experimental and modeling effort was to optimize the assembly and heat treatment of the tube-to-tubesheet joints, using a novel tapered joint configuration, and to develop diffusion-bonded joints with high shear strength and creep resistance, as has been done in other applications for tapered plugs. The project would also conduct creep testing of these tubesheet joints in pressurized furnaces under prototypical conditions of stress and temperature, however this was not feasible due to numerous hinderances. Detailed microstructural investigations was conducted on the diffusion bond interface, while the bond strength was evaluated by tensile testing. The mechanical testing quantified the overall strength of the bond and measured the pull out force as well as long term behavior, which complemented the microstructural analysis and gave a complete understanding of the diffusion bond.

30 DIRECT ENERGY CONVERSION↗

Final Report DE-FE0031785 Mohsen Ahmadian, Ph.D. Demonstration of Proof of Concept of a Multiphysics Approach for Real-Time Remote Monitoring of Dynamic Changes in Pressure and Salinity in Hydraulically Fractured Networks

Hydraulic fracturing has evolved into a multistep process with varying flow rates, carrier fluids (e.g., gel or slickwater), proppant loadings, and proppant grain sizes. As a result, primary recovery from a hydraulically fractured tight-oil reservoir is often a tiny fraction of the original oil in place, ranging between 5 and 10%. As stated in the FOA1990, “part of this problem is due to the inability of current well completion processes to effectively stimulate the entire reservoir volume in contact with the wellbore. Innovative technologies are needed that can help improve the effectiveness of reservoir completion methods, maximize stimulated reservoir volumes, and optimize recovery over the entire producing life span of a well”. We first need to enhance the current fracture diagnostic techniques to improve a well-completion design. However, detecting and delineating a subsurface hydraulic fracture is extremely difficult because the induced fracture network is only fractionally propped, and these propped fractures are generally very thin. Microseismic and tiltmeter monitoring techniques can provide information on the fracture extent but provide little or no information on the movement and final distribution of proppant or production fluids. On the other hand, electromagnetic (EM) imaging has shown the capability to monitor proppant distribution throughout the fracture area, especially in the presence of Electrically Active Proppants (EAPs). A previous EM survey of hydraulic fracturing at the Devine Fracture Pilot Site (DFPS) and subsequent EM code developments demonstrated this survey as a robust technique to remotely interrogate the extent of the EAP-filled hydraulic fracture during its propagation. The objectives of the project were threefold: (1) to capitalize on the material properties of an EAP to demonstrate remote monitoring of relative changes in pressure, pressure, and flow that are commonly encountered during production from a hydraulically fractured reservoir; (2) to evaluate EM imaging tools, to achieve Objective 1 in near real-time; and (3) to develop a multi-physics joint inversion approach to precisely predict flow patterns and physiochemical changes within an EAP-filled fracture network. This research project was built upon our previous work at the Devine Test Site managed by the Bureau of Economic Geology (BEG) at The University of Texas at Austin (UT-Austin). It also leveraged a significant investment from the Advanced Energy Consortium (AEC) to address the DOE's interest in subsurface flow, containment, and characterization by multiple signals. This three-year and three-month project succeeded in demonstrating the feasibility of a real-time dynamic fluid flow mapping technique at Technology Readiness Level 5 (TRL5) by utilizing a commercially available surface-based Controlled-Source Electromagnetic (CSEM) method (Objectives 1, 2). We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, the observed electric field in our study is influenced by bottomhole pressure, flow rate, and salinity, which is demonstrated by comparing these parameters with the electrical field potentials. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce the whole measured electric field magnitudes. Preliminary estimates showed that including Streaming Potential (SP) in our geophysical model is likely needed to reduce the simulation misfit.

02 PETROLEUM↗

Sparse chronology strategy for integrating seasonal energy storage in capacity expansion models

Here, this study develops the sparse chronology method to enhance the representative period framework in capacity expansion models, enabling the effective integration of long-duration energy storage modeling. Traditional representative period methods cannot capture the state of charge of seasonal energy storage systems because they do not establish effective inter-day linkages to connect the state of charge between periods. The sparse chronology approach addresses this limitation by establishing inter-day linkages that allow state of charge to shift inter-seasonally. At the same time, it groups identical representative days into partitions, applying constraints sparsely and implicitly to reduce computational load further. Validation results demonstrate that this method successfully simulates long-duration energy storage patterns, achieving close alignment with a continuous yearly benchmark model, with seasonal trends and state of charge cycles clearly represented. The computational load analysis reveals that the sparse chronology method efficiently applies constraints on maximum and minimum state of charge limits within the representative day framework, eliminating the need for detailed constraints on each individual day. By partitioning representative days and constraining only the start and end of each partition, the method significantly decreases computational requirements. Simulation results show that sparse chronology closely approximates the continuous yearly method's accuracy, even with as few as 20 representative days, achieving correlation values with the benchmark of nearly 0.9 in state of charge plots. Furthermore, it maintains computational efficiency, requiring only 4 % of the solver time compared to the continuous yearly method with 20 representative days. This approach allows capacity expansion models to incorporate long-duration energy storage with high temporal, spatial, and technological resolution, enabling more detailed modeling for large-scale power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Machine Learning Approach to Improve Air Traffic Management Initiatives

Collaborating closely with commercial air carriers and related organizations, the Federal Aviation Administration(FAA) regulates air traffic and ensures the safety and efficiency of air operations. Air traffic controllers make strategic decisions, such as delaying, rerouting, or canceling flights, partly based on guidance provided by the FAA’s Air TrafficControl System Command Center (ATCSCC). The guidance includes, among other things, control measures known asTraffic Management Initiatives (TMIs) designed to enhance safety and improve operational efficiency. TMIs play a crucial role in managing the demand and capacity within the U.S. National Airspace System (NAS). Two major TMIs that are routinely used (primarily to mitigate the adverse effects of bad weather) are Ground Delay Programs (GDPs) andGround Stops (GSs). In a GDP, flights destined for airports facing thunderstorm activity experience delays at their origin airports. This proactive approach minimizes the risk of routing aircraft through hazardous weather conditions and also replaces (fuel burning) airborne delays with ground delays. In a GS, a temporary restriction is imposed on the departure or arrival of aircraft at a specific airport or within a designated airspace. Although other TMIs (e.g., miles-in-trail) are also implemented as part of (air) traffic flow management in the NAS, the focus of this work is on GDPs and GSs. Since TMIs, by design, lead to flight delays or cancellations, it is crucial to put in place the right set of parameters(e.g., scope and duration of the GDP). For example, when the end time of a GDP extends beyond what is necessary, it imposes unnecessary delays on departing flights. This situation could occur as a result of inaccurate prediction of the(required) duration of the GDP based on the weather forecast. On the other hand, if a GDP ends prematurely before the underlying capacity constraints are resolved at the destination airport, it may result in airborne holding. The delicate balance lies in matching the termination of the GDP precisely with the resolution of capacity constraints, avoiding both the imposition of unnecessary ground delays and the need for airborne holding due to premature program termination.Failing to specify the right parameters for TMIs also leads to flight delays, creating a significant obstacle in managing the increasing traffic volumes causing increased work load for the controllers. To address this issue, we propose the integration of Machine Learning (ML) models in the traffic flow management(TFM) pipeline. In current operations, decisions are made by human experts based on extensive training, historical patterns, available traffic and weather data. Since we have an abundance of data from past events that tell us the likely impact of various TMIs, by ingesting historical data, properly trained ML models can offer valuable insights and aid human decision-making. With the FAA increasingly exploring advanced analytics, ML emerges as a focal point for enhancing TFM within the National Airspace System (NAS). As a first step, this study aims to provide traffic controllers with decision-making support for the issuance and adjustment of TMIs. Data analytics and machine learning have been previously employed to address some of the challenges associated with TMIs. Numerous studies have concentrated on various facets of TMI issuance, exploring factors influencing TMI parameters, including arrival rate, airport capacity, and delay prediction. For example, using weather forecasts, several statistical methods were used to produce probabilistic capacity profiles which in conjunction with deterministic models provided insights into the GDP planning process [1–4]. The downside of using deterministic models is that they rely on fixed inputs and predetermined rules, which lack the ability to account for the inherent uncertainty and variability present in real-world scenarios. In a separate series of studies, researchers aimed to predict the occurrences of GDPs and GSs. The majority of these studies utilized various supervised learning methods, including Decision Trees, Naive Bayes, Support VectorMachines, and Random Forests to analyze the influence of weather conditions and arrival demand on TMI incidents[5–8]. However, these studies primarily focused on predicting the incidence of TMIs without explicitly addressing the scope of TMIs, including their duration and their geographical coverage. Furthermore, the emphasis of these studies was largely on GDPs, given their higher frequency and longer duration when compared to GSs. A limited number of studies focused on predicting the parameters of TMIs, specifically addressing their duration and extent. In one such study focusing on optimizing the TMI parameters at San Francisco International Airport (SFO),the authors utilized a probabilistic forecast of fog [9]. They simulated various capacity scenarios based on the (fog)burn-off forecasts, selecting GDP parameters that minimized airborne and overall ground delays. However, this approach exclusively emphasizes stratus (fog) burn-off as the primary determinant of GDP and GS, neglecting other influential factors like severe weather events, runway closures, lower capacity than traffic demand, and other important variables. Given the complexity of predicting the TMI and determining its scope, we seek a more holistic approach. We aim to consider all significant factors that could impact TMIs and their parameters. What sets this research apart is the fusion of all data sources relevant to the issuance and adjustment of TMIs and it represents the first comprehensive attempt to optimize TMIs in this manner. Since this comprehensive solution involves various aspects, we break down the problem into smaller components and input all parameters into a unified model called the “TMI Adjuster”. Figure 1 shows the overall framework and the list of datasets used in each model. The objective of the TMI Adjuster module is to deliver reliable, consistent and expedited recommendations for the progression, adjustment, and termination of TMIs. The ML solution entails developing a pipeline capable of predicting the necessity of a TMI (e.g., GS or GDP) along with its various parameters. For example, in the case of a GS, this includes the scope of the GS either in terms of distance from the destination airport or based on pre-defined airspace sectors. Here, scope refers to those regions and departing airports that are subject to the GS. In this paper, we concentrate on the issuance of GSs in the three major airports in the New York area — LaGuardia(LGA), John F. Kennedy International (JFK), and Newark Liberty International (EWR). We fuse traffic, weather and other relevant aviation data from years 2017 to 2019 to train and validate the ML models. In particular, we use the following datasets: •Terminal Aerodrome Forecast (TAF): meteorological forecasts specific to each airport, issued four times a day, covering predefined time periods. •TMI data: includes all GSs and GDPs along with their respective parameters. •Aviation System Performance Metrics (ASPM): includes traffic related data such as aircraft delays, arrival, and departure rates. •Notices to Airmen (NOTAMs): utilized to extract runway closure data and manage interdependencies between terminals in close proximity. •Flight cancellation data •Airspace Flow Programs (AFP): includes information on flight airborne holdings caused by TMIs. The data preprocessing entails transforming ASPM, TMI, AFP, NOTAMs, and weather data into an hourly format and consolidating all datasets by merging them based on date and time as the primary key. The TMI Adjuster framework comprises two parallel models: one dedicated to GS and a second model focused on GDP. As previously mentioned, our specific focus is on the GS model as a multi-classification problem. In this framework, each data point of the GS model input summarizes ten hours of data. Specifically, the data loader for the GS model generates the input and output of the model as follows: at a given time step, the input includes the actual traffic, weather, and TMI data from the two-hour window before the time step, alongside the weather forecast and scheduled traffic for the next 8 hours starting from the time step. Based on this information, the output of the GS model for each time interval consists of three dimensions. The first dimension represents a binary decision on whether there should be a GS in place for the next hour or not. The second dimension is related to the scope of the GS in the United States, and the third dimension is related to the scope of the GS in Canada (i.e., to determine if the GS impacts airports in Canada).One of the challenges with TMI modeling is the sparsity of TMI events, particularly regarding its scope. To address this challenge in the scope of the GS model output, we implement grouping. The GS scope for the US region is defined based on a list of centers that should be included when the GS is in place. With 20 centers in the US, we utilized historical data to group them into 4 categories. In particular, we summarized our historical data in a graph format where nodes represent centers, and link weights are defined based on the co-occurrence of centers in the scope parameter ofTMIs. By identified strongly connected components in this graph, we were able to partition the centers into four groups. We consider two model structures for the GS Model. Firstly, a hierarchical classification model [10], where the human decision-making for a GS is of hierarchical nature. The decision-maker first decides whether there is a need fora GS, and if the answer is yes, determines the scope. A hierarchical classification model organizes the problem into a class hierarchy, typically a tree or a Directed Acyclic Graph (DAG) structure, and considers the dependency of the decision in the previous step to the next component [10]. Here, we employ the local classifier per level approach, which involves training one multi-class classifier for each level of the class hierarchy. The second structure is the independent structure. In this setting, as the name suggests, we do not consider the dependency of the decisions in the different dimensions of the output of the model. Instead, for each dimension, we train a multi-class classifier independently. Table 1 summarizes GS model statistics for training, validation and testing. The table documents the effect of limiting data to the time steps when there was actually a TMI in place or when a TMI had just terminated. This resulted in a more balanced distribution of the GS class(GS positive class)versus “No GS”(GS negative class), which might help the training process. While JFK and LGA follow very similar distributions, with 40% and 42% GS positive class respectively, EWR has proportionally fewer GS incidents at 28%. Our subsequent phase involves evaluating the performance of both hierarchical structure and independent structure using different state-of-the-art multi-class classifier models such as Random Forest, Decision Trees, K-nearest Neighbors, and Logistic Regression and forecast the duration and scope of the GSs.

Farzan Masrour Shalmani↗

Prefabricated High-Strength Rebar Systems with High-Performance Concrete for Accelerated Construction of Nuclear Concrete Structures

This project aimed to reduce the field construction times and fabrication costs of reinforced concrete (RC) nuclear structures through the use of: 1) high-strength steel deformed reinforcing bar (rebar); 2) prefabricated high-strength rebar (HSR) assemblies with headed anchorages; and 3) high-strength concrete (HSC). The focus of the project was not to develop new innovations on materials, but rather to fill a major knowledge gap on the effectiveness, code conformity, and viability (i.e., practicality, commercial availability) of existing high-strength materials in nuclear structures, especially in shear walls, their connections/joints, and around penetrations. These advances promote the Nuclear Energy Enabling Technologies Advanced Methods for Manufacturing (NEET-1) goals of accelerating deployment schedules by at least 6 months and reducing component fabrication costs by 20% or more. Toward this vision, the following objectives were identified and completed: Develop limit-benefit and cost-benefit analysis frameworks: The structural and economic limits and benefits of HSR and HSC in nuclear walls were evaluated. Importantly, it was found that modest increases in concrete strength up to 10.0 ksi and rebar yield strength up to 100 ksi result in the most economic designs and improved structural performance (Barbachyn et al., 2017a and Barbachyn et al., 2017b). Also, an industry survey was distributed to determine the cost benefits associated with prefabrication. Importantly, on-site worker-hours are significantly reduced, allowing for significantly accelerated construction schedules (Devine et al. 2018b). Develop an optimization methodology for structural design: An optimization study was deemed not necessary since the trivial solution with maximized use of high-strength materials and prefabrication always resulted in the largest cost and time savings. Instead, the benefits (and limits) of high-strength materials and prefabrication on construction costs and on-site construction time were investigated through a parametric study. Cost-reductions through the use of HSR and HSC were found to be greater than 40%. Also, on-site construction time of nuclear shear walls can be reduced by up to 83% through the use of prefabricated rebar assemblies. Conduct experimental evaluations of nuclear shear walls and wall-foundation joints: A large scale experimental program was completed on shear-critical RC cantilever deep beams. From these experiments, it was demonstrated that increasing rebar yield strength had a greater effect on increasing lateral strength than increasing the concrete compressive strength. Further, using HSR together with HSC resulted in the highest lateral strength and displacement ductility (Devine et al. 2018a). Additionally, large-scale reversed-cyclic lateral load tests were conducted on RC shear walls. It was demonstrated that through the utilization of HSR and HSC, a wall with 55% reduced vertical and horizontal rebar area was able to achieve 91% of the peak lateral strength of a conventional wall using state-of-practice materials with the same geometry. Further, trim reinforcement with headed anchorages around the penetrations proved to be effective. Develop validated numerical simulation models: Numerical modeling was completed in two phases during this project using 3 different finite element software packages. The first phase used existing measured data from previous experimental studies on stocky RC walls to validate simplified predictive models. Models using two different programs were validated based on a database of 38 stocky walls (Barbachyn et al., 2017a). The second phase occurred after experimental testing using a third and more detailed finite element package. Prediction improvements for the initial stiffness, lateral deformations, damage propagation, and failure mechanisms were targeted. The resulting detailed models more accurately predicted the peak lateral strength, initial and secant stiffnesses, lateral deformations, and cracking patterns of the wall test specimens subjected to reversed-cyclic lateral loading as compared to the simplified pretest models. Develop validated design procedures, tools, and criteria: This project evaluated the accuracy of existing closed-form design equations from U.S. building design standards and previous research (Barbachyn et al. 2017a and Devine et al. 2018a). It was shown that the current predictive equations for shear strength capacity need to be revisited as they provide significant scatter and gross overestimates, indicating the potential for unconservative designs. Further, for both the deep beam and shear wall specimens, the numerical models and flexural capacity predictions provided the most accurate strength predictions. Develop cost-effective field construction procedures: To promote the use of prefabricated rebar assemblies, the following construction conditions where prefabrication would be used instead of in-place rebar tying were identified from an industry survey: 1) to save on construction schedule; 2) to improve safety and/or quality control; 3) for areas with heavy rebar congestion; and 4) for structures with significant repetition in rebar layout. All of these conditions directly apply and benefit nuclear construction where construction costs and schedules are extensive, safety and quality are of utmost importance, the reinforcement is extremely dense, and the reinforcement layout is repetitive in many regions of the shear walls. Additionally, an experimental study was conducted to investigate the changes in bar spacing (with respect to code-required rebar spacing and placement tolerances) when a horizontally assembled prefabricated two-dimensional (2D) rebar mat or three-dimensional (3D) rebar cage is lifted (or tripped) to a vertical position. It was found that the bars involved in both the lateral translation of the rebar assemblies (i.e., movement of the assembly in its horizontal position) and the tripping process were most susceptible to large spacing changes, and thus, their positions within the assembly would need to be checked after the placement of the assembly in place within the formwork. Importantly, the spacings of all of the other bars in the assembly were not affected significantly by the lateral translation and tripping process, thereby eliminating the need to recheck their positions within the assembly (Devine et al. 2018b).

36 MATERIALS SCIENCE↗

Storage Enabled Flexibility of Conventional Generation Assets (StorFlex)

The power systems have faced progressively more demanding operational requirements over the last two decades. Several factors contribute to these challenging operating conditions, including load growth, aging infrastructure, increasing penetrations of distributed energy resources (DERs), electrification of the economy, and policy initiatives such as decarbonization. The power system and its components must provide high operational flexibility to mitigate these challenges. For example, the proliferation of intermittent DERs such as wind and solar has increased the need for conventional generation assets like hydropower plants to respond to sudden load-generation imbalances. The higher flexibility requirements for hydropower plants cause more wear and tear, potentially shortening the useful lifespan of hydropower turbines. To reduce the need for hydropower plants to follow sudden changes in the dispatch signal, we investigate their combined operation with the energy storage systems (ESSs; “ESS-based hybridization”). Our analyses focuses on improving the lifespan of hydropower plants through ESS-based hybridization. Wear and tear on hydropower turbines (particularly Francis turbines) is modeled using a loss-of-life concept that is based on damage experienced by the turbine due to various cycles of operation. Then, we show that using ESSs to offset some of the high variation increases the remaining life of the hydropower plants. To demonstrate this, a few modeling tools were developed for this work: (1) a dynamic model for various components of the turbine and its governor; (2) a control strategy that assigns a slow-varying dispatch signal to a hydropower unit versus a fastmoving signal to ESS, such that the overall power request remains the same; and (3) models for the financial analysis to quantify the economic merits of such a framework. We used the models we developed to analyze the dispatch pattern of an actual hydropower plant with a power output of 50 MW and a head height of 152 m. This work showed that ESS-based hybridization could extend the life of the hydropower plant by 5% on average. This extension in life was then used to estimate the economic benefit in terms of cost deferrals associated with hydropower plant maintenance and replacement: on average, $3.6 million. Sensitivity analysis with respect to the size of ESS and cost of turbines was performed to show the variation in benefits over the range of turbine costs and ESS sizes. Crucially, stacking damage reduction and lifetime extension with other ESS value streams such as providing ancillary services could substantially increase the financial benefits of ESS-based hybridization. The higher costs associated with ESS of appropriate size would make more financial sense when multiple value streams are stacked and co-optimized to extract the maximum benefit. This dimension will be explored in future work.

13 HYDRO ENERGY↗

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗

Data Locality Enhancement of Dynamic Simulations for Exascale Computing (Final Report)

The development of modern processors exhibits two trends that complicate the optimizations of modern software. The first is the increasing sensitivity of processors' throughput to irregularities in computation. With more processors produced through a massive integration of simple cores, future systems will increasingly favor regular data-level parallel computations, but deviate from the needs of applications with complex patterns. Some evidences are already shown on Graphic Processing Units (GPU): Irregular data accesses (e.g., indirect references A[D[i]]) and conditional branches are limiting many GPU applications' performance at a level an order of magnitude lower than the peak of GPU. The second hardware trend is the growing gap between memory bandwidth and the aggregate speed—that is, the sum of all cores' computing power—of a Chip Multiprocessor (CMP). Despite the capped growth of the peak CPU speed, the aggregate speed of a CMP keeps increasing as more cores get into a single chip. It is expected that by 2018, node concurrency in an exascale system will increase by hundreds of times, whereas, memory bandwidth will expand by only 10 to 20 times. Consequently, data movement and storage is expected to consume more than 70% of the total system power. Bridging this gap is difficult; the complexities of modern CMP memory hierarchy make it even harder: Data cache becomes shared among computing units, and the sharing is often non-uniform—whether two computing units share a cache depends on their proximity and the level of the cache. On the recent IBM Power7 architecture, for instance, four hardware contexts (or SMT threads) in a core share the entire memory hierarchy, all cores in one chip share an on-chip L3 cache, and cores across chips share L3 and main memory through off-chip connections. These two trends complicate the translation of computing power into performance, especially for a program with either intensive data accesses or complex patterns in data accesses or control flow paths. Unfortunately, both attributes present and will persist in a class of important applications. For instance, many scientific simulations deal with a large volume of data. And meanwhile, as most real-world processes are non-uniform and evolving (e.g., the evolution of a galaxy or the process of a drug injection), both the computations and data accesses of these programs tend to be irregular and dynamically changing. Currently, the lack of support to these applications on modern CMP severely limits their performance. On GPU, as our recent study shows and other studies echo, performance enhancement of a factor of integers is possible when memory accesses or control flows are streamlined for a set of GPU applications. On multicore CPU, our studies show that traditional locality enhancement, for being oblivious to the new features of multicore memory hierarchy, may even cause large slowdown to data-intensive dynamic applications. The severity of the issues is expected to worsen as the two hardware trends continue. Some recent studies try to match software with the trends, but in a limited scope or manner. For irregularities on GPU, most studies focus on irregularities analyzable through static analysis (e.g., data accesses in regular loops). Dynamic irregularities are harder to address because the needed analysis and transformations typically have to happen at run time. Some other research resorts to hardware extensions, an actual adoption of which is unclear for the entailed space cost and complexity. For data locality, recent years have seen some exploitations of the new memory hierarchy on multicore for performance, but most of them are on process or thread scheduling, rather than program transformations. Our recent study reveals that program-level transformations may magnify the scheduling benefits by a factor of seven, concluding that program-level transformation should play a central role for data locality enhancement on modern CMP. But research in this direction has been sparse, and most have focused on data layout or cache performance modeling, rather than program transformations to match with the new memory hierarchy features. Overall, it is still an open question how to bridge the gap between dynamic computations and the two prominent properties of modern processors. The goal of this project is to develop a set of techniques and software tools to enhance the matching between memory accesses in dynamic simulations and the prominent features of modern and future CMP, alleviating the memory performance issues for petascale and exascale computing. This report summarizes the discoveries and products produced throughout this project. It includes free launch, a new software approach to overcoming the shortcomings of both methods; coherence-free multiview, an approach that allows multiple views of a single data object to co-exist on GPU memory during a GPU kernel execution; algorithmic optimizations to data analytics problems, especially those that involve lots of distance calculations. Supporting dynamic parallelism is important for GPU to benefit a broad range of applications. There are currently two fundamental ways for programs to exploit dynamic parallelism on GPU: a software-based approach with software-managed worklists, and a hardware-based approach through dynamic subkernel launches. Neither is satisfactory. The former is complicated to program and is often subject to some load imbalance; the latter suffers large runtime overhead. In this work, we propose free launch, a new software approach to overcoming the shortcomings of both methods. It allows programmers to use subkernel launches to express dynamic parallelism. It employs a novel compiler-based code transformation named subkernel launch removal to replace the subkernel launches with the reuse of parent threads. Coupled with an adaptive task assignment mechanism, the transformation reassigns the tasks in the subkernels to the parent threads with a good load balance. The technique requires no hardware extensions, immediately deployable on existing GPUs. It keeps the programming convenience of the subkernel launch-based approach while avoiding its large runtime overhead. Meanwhile, its superior load balancing makes it outperform manual worklist-based techniques by 3X on average. The work was published at the 48th Annual IEEE/ACM International Symposium on Microarchitecture (Micro'2015). A Graphic Processing Unit (GPU) system is typically equipped with many types of memory (e.g., global, constant, texture, shared, cache). Data placement determines what data are placed on which type of memory, essential for GPU memory performance. Prior optimizations of data placement always require a single view of a data object on memory, which limits the optimization effectiveness. In this work, we propose coherence-free multiview, an approach that allows multiple views of a single data object to co-exist on GPU memory during a GPU kernel execution. We demonstrate that under certain conditions, the multiple views can remain incoherent while facilitating enhanced data placement. We present a theorem and some compiler support to ensure the soundness of the usage of coherence-free multiview. We further develop reference-discerning data placement, a new way to enhance data placements on GPU. It enables more flexible data placements by using coherence-free multiview to leverage the slack in coherence requirement of some GPU programs. Experiments on three types of GPU systems show that, with less than 200KB space cost, the new data placement technique can provide a 1.6X average (up to 4.27X) speedup. In addition, we have examined algorithmic optimizations to data analytics problems, especially those that involve lots of distance calculations. Computing distances among data points is an essential part of many important algorithms in data analytics, graph analysis, and other domains. In each of these domains, developers have spent significant manual effort optimizing algorithms, often through novel applications of the triangle equality, in order to minimize the number of distance computations in the algorithms. In this work, we observe that many algorithms across these domains can be generalized as an instance of a generic distance-related abstraction. Based on this abstraction, we derive seven principles for correctly applying the triangular inequality to optimize distance-related algorithms. Guided by the findings, we develop Triangular OPtimizer (TOP), the first software framework that is able to automatically produce optimized algorithms that either matches or outperforms manually designed algorithms for solving distance-related problems. TOP achieves up to 237x speedups and 2.5X on average. The work has been published at the 32nd International Conference on Machine Learning (ICML'15) and the 41st International Conference on Very Large Data Bases (VLDB'15).

42 ENGINEERING↗

Modular Controllable Transformers (MCT)

Large Power Transformers (LPTs) are at the heart of the electricity infrastructure and are critical pieces in today’s grid. Aging assets, long turn-around times, transportation delays and dependence on foreign manufacturing for these components have created conditions for points of failure across the system. Moreover, there is no available dynamic control on these passive sections of the grid. For all the automation embedded on the low voltage side of the low voltage side of the system there is absolutely no flexibility on the high voltage side. Changes in power flow patterns are brought about by altering dispatch patterns for generators which requires a coordinated effort and is often time consuming. Traditional FACTS based approaches such as UPFC or HVDC light involve massive modifications to the existing infrastructure. Recent instances of physical attacks on LPTs and the slow recovery process associated with fixing or replacing these damaged units has created an increased awareness for this issue. LPT designs are highly customized owing to specific field requirements. Long lead times coupled with transportation delays make procuring and installing LPTs a lengthy process. In an attempt to minimize the recovery time following an LPT failure, transformer manufacturers have developed mobile transformer units. These units typically consist of either a single-phase or three-phase transformer mounted on a truck and characterized by compact form factors. By reducing the assembly time and transport time these units offer a fast and temporary solution to LPT failures. Several efforts led by DOE, Edison Electric Company (EEI) and North American Electric Reliability Corporation (NERC) have led to the RecX program. RecX transformers consist of single phase transformers transported on specialized assemblies to enable fast restoration and minimal installation time. NERC and EEI have instituted transformer sharing programs like Spare Equipment Database (SED) as well as Spare Transformer Equipment Database (STEP) which encourage utilities to maintain and share a fleet of LPTs. Mobile transformers are typically limited to less than 100 MVA to enable transportation on a trailer. Also, the efficiency of mobile transformer units and RecX transformers is lower since they are temporary solutions. The ideal replacement solution would be one that could utilize a set of standardized transformers while incorporating flexibility to influence metrics like power flow, apparent impedance and voltage profiles. In view of these issues the team at Georgia Tech proposed a concept called Modular Controllable Transformers (MCTs). This project was commissioned as a part of DOE’s Transformer Resiliency and Advance Component (TRAC) program. The approach consists of splitting one large unit rated upwards of 100 MVA with smaller modular, standardized units. Each of these smaller units are augmented with a converter to add power flow, voltage and impedance control. Splitting the power level into multiple units has two important implications. By making smaller units, the transportation times can be reduced drastically. Moreover, even though these units have very specific impedance and other parameters, the augmented converters can make it emulate the unit being replaced. Thus, faster replacement times and added controls are achieved making this an extremely resilient approach. In the event of one modular units failure, the other two controllable units can operate in the most resilient fashion keeping the net system overloads to a minimum. Moreover, for the same probability of failure only a fraction of the system capacity is lost. The converter is augmented on to the tertiary winding of the transformer making it retrofittable in the field without massive modifications to the system. The material proposed in this report highlights design considerations to minimize recovery time as well as results displaying overall improvement in resiliency. The cost analysis of this approach shows significant savings compared to UPFC and HVDC lite based approaches. Extensive analysis was conducted on numerous aspects of this project over the course of this project. Simulation studies were conducted to prove the 5MVA converter design. Challenging BIL constraints and transient issues were addressed over the course of this project. Simulations on EMTP platform presented in consequent sections were used to finely tune the control for the same. A final design for a back-to-back 5MVA neutral-point-clamped (NPC) converter was completed in this project. Furthermore, switching components like high power IGBTs from IXYS Westcode Press-Pack series and ABB Stackpack series were identified for this application. These devices exhibit a blocking voltage of 6.5 kV and a nominal collector current of 900 A. Potential housing options as well as existing commercial option were identified for this converter design. A 56 MVA modular transformer unit was designed with Delta Star Inc. with consultation and inputs from all project members. The transformer cooling systems were chosen to be Oil Natural Air Forced (ONAF) to optimize ease of transport. A full Product Requirement Document (PRD) for the converter and the transformer was generated. The team also investigated key aspects of resiliency and how to quantify the same. A metric was developed to analyze resiliency and numerous simulation cases were presented. Increased (N-1) resiliency and outage resistance was showcased through simulation on the Texas Large scale power system. Dynamic control capabilities were also showcased on numerous system level simulations. An IEEE 30 bus system was used to show the ability of the MCT approach to reduce line congestion. Power flow control; a part of the MCTs abilities, showcased a 29% reduction in transformer loading on the same system. Full control over tie line flows was demonstrated on a modified IEEE 13 bus system. Dynamic LTC capabilities were shown through simulation and were compared to traditional LTCs. Furthermore, the ability of the MCTs to alter apparent transformer impedance was simulated to promote standardization of modular transformer units while retaining the ability to emulate a varied range of impedance characteristics. Thus granular control over power flow, voltages, impedances and was presented. Regulatory analysis conducted by Oak Ridge National Laboratory (ORNL) further bolstered the efficacy of the MCT approach and verified the added resiliency achieved due to this approach. As a conclusion, it was proved that the resilience improvement obtained from using the MCT approach is justified by just one failure in a 30-year period. This means that the benefits obtained from using the MCT approach far outweigh the cost of setting this approach up. As a part of the future efforts, the team is proposing building a (12.47-24) kV/400 kVA scaled down unit to control 5MVA of power flows. Emphasis would be laid on using commercial components to speed up the build of these units. Further the actual transport of the modular transformer unit would be a part of this operation and would provide insight into transportation issues. The approach shows much promise and adds key dynamic control on passive sub transmission and transmission sectors.

24 POWER TRANSMISSION AND DISTRIBUTION↗