**OPTIMIZATION FRAMEWORK AND METHODS FOR ADAPTIVE EV CHARGING**

JP2014514909 | A charging method and charging equipment for battery charge |

JPS57155696 | DATA PROCESSOR FOR AUTOMOBILE |

JP2009265888 | MOBILE INFORMATION ACQUISITION SYSTEM AND METHOD |

LEE, Zachary, J. (1124 Cordova Street, Apartment 2Pasadena, CA Pasadena, Pasadena, US)

LOW, Steven, H. (200 Mariner Views Lane, La Cañada, CA, 91011, US)

CHANG, Daniel, C. (471 Northlake Drive, Unit #18San Jose, CA, 95117, US)

JIN, Cheng (12474 Woodside Drive, Saratoga, CA, 95070, US)

LEE, Rand, B. (1680 Bedford Road, San Marino, CA, 91108, US)

LEE, George, S. (1200 E. California Boulevard, MC 305-16Pasadena, CA, 91125, US)

*;*

**G06Q50/30***;*

**G06Q20/14**

**G07F15/00**What is claimed is; 1. An electric vehicle charging network, comprising; one or more centralized computing systems; a communications network; a plurality of electric vehicle node controllers for charing a pluralit of electric vehicles tEVs), where each electric vehicle node controller in the plurality of node controller comprises; a network interface; a processor; a memory containing; an adaptive charing application; a plurality of electric vehicle node parameters describing charging parameters of an electric vehicle in the electric vehicle charging network; wherein the process is configured by the adaptive charing application to: send electric vehicle node parameters to the one or more centralized computing systems; and charge an electric vehicle using a charging rate received fro the one or more centralized computing systems; wherein the one or more centralized computing systems is configured to: receive the electric vehicle node parameters from the plurality of electric vehicle node controllers; calculate a plurality of charging rates for the plurality of electric vehicle node controllers using quadratic programming (QP), wherein the quadratic programming computes the plurality of charging rates based on the electric vehicle node parameters, a plurality of adaptive charging parameters and a quadratic cost funetiomand distributing the charging rates to the plurality of electric vehicle node controllers. 2. Hie electric vehicle charging network of claim 1 , wherein the charging rates meet energy demands of the plurality of electric vehicles by a plurality of departure times. 3. The electric vehicle charging network of claim 1 , wherein the charging rales for t e plurality of electric vehicle nodes a a time series of timing rates provided to each electric vehicle node controller. 4. The electric vehicle charging network of claim 1. wherein the electric vehicle node parameters include an electric vehicle departure time, a remaining energy demand, and a maximum charging rate. 5. The electric vehicle charging network of claim 3, wherein the time series of timing rates can be evaluated by the one or more centralized computing systems using the following expression: min c{r) (42a) subject to ^G,I 6, The electric vehicle charging network of claim L wherein the charging rates for a first set of electric vehicle node controllers are a maximum charging rate; and the charging rates for a second set of electric vehicle node controllers are a minimum charging rate. 7. The electric vehicle charging network of claim 1, wherein calculating die plurality of charging rates is a least laxity first process. 8. The electric vehicle charing network of claim I . wherein a first electric vehicle node controller is assigned a first group that is guaranteed that, for each EV in the first group, a requested energy will be delivered, and a second electric vehicle node controller is assigned to a second group that is guaranteed, for each EV in the second group, a minimum energy. 9 The electric vehicle charging network of claim 8, wherein charging rates for the first group and second group are determined sequentially by: solving for charging rates using QP for EVs in the first group only; computing left-over capacities for Group 2; and solving for charging rates using QP for EVs in the second group only using the left-over capacities. 10. The electric vehicle charging network of claim 1 , further comprising reducing rate fluctuations across a time period by using a penalty term within the quadratic cost function. 1 1. Tiie electric vehicle charging network of claim 1 , wherein the one or more centralized computing systems is configured to: set a charging rate to he greater than zero for an EV with a remaining energy demand. 12. The electric vehicle charging network of claim 1, wherein the one or more centralized computing systems is configured to: receive a request for admission of an electric vehicle from an electric vehicle node controller; determine whether to admit the electric vehicle based on existing electric vehicle node parameters from the plurality of electric vehicle node controllers and existing capacity constraints. 13. The electric vehicle charging network of claim 1. wherein the one or more centralized computing systems is configure to: prioritize the plurality of electric vehicle node controllers when there is insufficient capacity to meet energy demands of the plurality of electric vehicles. Ί 14. The electric vehicle charing network of claim 1. wherein the one or more centralized computing systems is configure to schedule charging by the plurality of electric vehicle node controllers based on solar generation. 15. The electric vehicle charing network of claim 14, wherein the one or more centralized computing systems is configure to select charging rates that minimize the distance between a forecasted solar generation and a total net load, 16. The electric vehicle charging network of claim 1, wherein the one or more centralized computing systems is configured to cap a total site load over a time period. 17. The electric vehicle charing network of claim 1. wherein the plurality of electric vehicle node controllers are connected in delta configurations providing an unbalanced t hree -phase i n i rastr uc to re . 18. The electric vehicle charging network of claim 17, further comprising providing phase constrains and line constraints on currents along legs of the delta configurations. 19. The electric vehicle charging network of claim I , wherein the one or more centralized computing .systems is configure to; determining that a minimum energy demand cannot he rnet for an EV and post- processing, using the QP, the plurality of charging rates. 20. The electric vehicle charing network of claim 1 , wherein the plurality of adaptive charging parameters are quadratic capacity constraints. /a. |

STATEMENT OF FEDERALLY SPONSORED RESEARCH

[0001] This invention was made with government support under Grant No, IIP 16021 19 an CCF 1637598 awarded by the National Science Foundation The government has certain rights in the invention.

CROSS-REFERENCE TO RELATED APPLICATIONS

[0002] The present invention claims priority to U.S. Provisional Patent Application Serial No. 62/593,755 entitled‘Optimization Framework and Methods for Adaptive EV Charging” ied December 1 , 2017 and to U.S. Provisional Patent Application Serial No. 62/678,877 entitled‘Optimization Framework and Methods for Adaptive EV Charging” filed May 31 , 201 . The disclosures of IPS. Provisional Patent Application Serial No. 62/593.755 and U.S. Provisional Patent Appli cation Serial No. 62/678,877 are herein incorporated by reference in their entirety.

FIELD OF THE I VENTION

[0003} The present invention generally relates to electric vehicles and more specifically relates to adaptive charging station optimization for electric vehicles.

BACKGROUND

10004} An incredible amount of infrastructure is relied upon to transport electricity from power stations, where the majority of electricity is currently generated, to where it is consumed by individuals. Power stations can generate electricity in a number of ways including using fossil fuels or using renewable energy sources such as solar, wind, and hydroelectric sources. Substations typically do not generate electricity, but can change the voltage level of the electricity as well as provide protection to other grid infrastructure during faults and outages. From here, the electricity travels over distribution lines to bring electricity to locations where ir is consumed such as homes, businesses, and schools. The term‘'smart grid ^{' }" describes a new approach to power distribution which leverages advanced technology to track and manage the distribution of electricity, A smart grid applies upgrades to existing power grid infrastructure including the addition of more renewable energy sources, advanced smart meters that digitally record power usage in real time, and bidirectional energy flow that enables the generation and s torage of e nergy i n dditional places along the electric grid.

[ 0051 Electric vehicles (EVs), which include plug-in hybrid electric vehicles (PHEVs), can use an electric motor for propulsion, EV adoption has been spurred by federal, state, and local gov ernment policies providing various incentives (e.g. rebates, fast lanes, parking, etc.). Continued EV adoption is likely to have a significant i mpact on the future s mart grid due to the additional stress load that EVs add to the grid (an EV Y power demand can be many times that of an average residential house).

SUMMARY OF THE INVENTION

[0006] Adaptive charging networks in accordance with embodiments of the invention enable the optimization of electric design of charging networks for electric vehicles. One embodiment includes an electric vehicle charging network, including: one or more centralized computing systems, a communications network, a plurality of electric vehicle node controllers for charing a plurality of electric vehicles ( EVs), where each electric vehicle node controller in the plurality of node controller comprises: a network interface; a processor; a memory containing; an adaptive charing application; a pluralit of electric vehicle node parameters describing chargin parameters of an electric vehicle in the electric vehicle charging networks herein the process is configured by the adaptive charing application to; send electric vehicle node parameters to the erne or more centralized computing systems; and charge an electric vehicle using a charging rate received from the one or more centralized computing systems; wherein the one or more centralized computing systems is configured to; receive She electric vehicle node parameters from the plurality of electric vehicle node controllers; calculate a plurality of charging rates for the plurality of electric vehicle node controllers using quadratic programming i QP . wherein the quadratic programming computes the plurality of charging rates based on the electric vehicle node parameters, a plurality of adaptive charging parameters and a quadratic cost funetiorgand di tributing the charging rates to the plurality of electric vehicle node controllers.

[0007] In a further embodiment, the charging rates meet energy demands of the plurality of electric vehicles by a plurality of departure times.

[0008] In another embodiment, the charging rates for the plurality of electric vehicle nodes are a time series of timing rates provided to each electric vehicle node controller.

[0009] In a still further embodiment, the electric vehicle node parameters include an electric vehicle departure time, a remaining energy demand, and a maximum charging rate.

[0010] in still another embodiment, the time series of timing rates can be evaluated by the one or more centralized computing systems using the following expression;

min, r)

subject t

all resources l wherec{r)isae<m‘exqmdratkx0$tf unction, Cj is an energy demand, is a minimum energy that will be delivered, T is an optimization horizon, r _{t }(! ) is the calculated charging rate, iy is a maximum rate, P‘ ) denotes the capacities of resources / at time /, and ,¼ are the coefficients that describe how EV fs are constrained by resources /, wherein for each QP instance, / can be t 1 , 2, . , 7. .

f 0 11] in a yet further embodiment, the charging rates for a first set of electric vehicle «ode controllers are a maximum charging rate; and the charging rates for a second set of electric vehicle node controllers are a minimum charging rate,

[0012] In yet another embodiment, calculating the plurality of charging rates is a least laxity first process.

[0013] In a further embodiment again, a first electric vehicle node controller is assigned a first group that is guaranteed that, for each BV in the first group, a requested energy will be delivered, and a econd electric vehicle node controller is assigned to a second group that is guaranteed for each EV in the second group, a minimum energy

[0014] in another embodiment again, charging rates for the first g roup and second group are determined sequentially by: solving for charging rates using QP for B V in the first group only;

computing left-over capacities for Group 2; and

solving for charging rates usin QP for EVs in the second group only using the left-over capacities.

[0015] la a further additional embodiment, the one or more centralized computing systems is configure to reduce rate 11 actuations ac ross a ti me pe riod by us ing a pe natty te r wi thin the quadratic cost function,

[0016] In another additional embodiment, the one or more centralized computing systems is configured c o: set a charging rate to be greater than zero for an E V w ith a remaining e nerg demand

[0017] In a still yet further embodiment, the one or more centralized computing systems is configured t o: receive a request for a mission f a n e lectric vehicle from a n electric vehicle node controller: determine whether to admit the electric vehicle based on existing electric vehicle node parameters from the plurality of electric vehicle node controllers and existing capacity constraints [0018] In still yet another embodiment, the one or more centralized computing systems is configure to: prioritize the plurality of electric vehicle node controllers when there is insufficient capacity to meet energy demands of the plurality of electric vehicles,

[0019] in a still further embodiment again, the one or more centralized computing systems is configure to schedule charging by the plurality of electric vehicle node controllers based on solar generation.

[0020] In still another embodiment again, the one or more centralized computing systems is configure to select charging rates that minimize the distance between a forecasted solar generation and a total net load. [0021] In still another embodiment again, the one or more centralizer! computing systems is configured to cap a total site load over a time period.

[0022] In still another embodiment yet again, the plurality of electric vehicle node controllers are connected in delta configurations providing an unbalanced three-phase infrastructure.

[1)023] in still another embodiment yet again still, the electric vehicle charging network further includes providing phase constrains and line constraints on currents along legs of the delta configurations,

[0024 j In yet still another embodiment again, the one or more centralized computing systems is configure to: determining that a minimum energy demand cannot be met for an EV and postprocessing, using the QP, the plurality of charging rates,

[0025] In another embodiment again, the plurality of adaptive charging parameters are quadratic capacity constraints.

BRIEF DESCRIPTION OF THE DRAWINGS

[0020] Fig. 1 is a diagram conceptually illustrating a power distribution network in accordance with an embodiment of the invention.

[0027] Fig. 2 is a diagram conceptually illustrating electric vehicle nodes connected to a communications network in accordance with an embodiment of the invention.

[0028] Fig. 3 is a block diagram of an electric vehicle node controller in accordance with an embodiment of the invention.

[0029] Fig. 4 illustrates a quadratic programming (QP) framework for computing charging rates for an electric vehicle (EV) in accordance with an embodiment of the invention.

[0030] Fig. 5 illustrates an BVSB configuration lor a single phase or Y-configumicm of a three- phase system in accordance with an embodiment of the invention.

[0031 ] Fig. 6 illustrates a three-phase design whereby a three-phase voltage source is connected to loads in parallel in accordance with an embodiment of the invention.

[0032] lug. 7a illustrates a three-phase transformer in accordance with an embodiment of the invention.

[0033] Fig, 7b illustrates an electric panel in accordance with an embodiment of the invention. (0034) Ftg. 8 illustrates a reference design transformer with voltages 480V - 120/2O8V and power rating of 75kVA in accordance with an embodiment of the invention,

(0035) Fig, 9a illustrates a 5-wire three-phase transformer in Wye configuration in accordance with an embodiment of the invention.

(0036) Fig. 9b illustrates a 5-wire transformer in Delta configuration with one of the phases center-tapped to provide three voltage levels in accordance with an embodiment of the invention.

10037] Fig 10a shows a Wye-configured voltage source connected to a set of loads in Delta configuration in accordance with an embodiment of the invention.

(0038) Fig. 10b illustrates an electrical panel arrangement to connect loads to voltage sources in accordance with an embodiment of the invention.

(0039) Fig, 1 1 illustrates computing inequalities for enforcing line limits using magnitudes and phases of loads in accord acne with an embodiment of the invention,

(0040) Fig. 12 illustrates parameters for groups in accordance with an embodiment of the invention.

(0041) Fig, 13 i llustrates an algorithm for checking the laxity of a new EV request in accordance with an embodiment of the invention

(0042) Fig 14 illustrates an algorithm for a QP with minimum rates for Group 2 in accordance with an embodiment of the invention

(0043) Fig 15 illustrates an algorithm for allocating minimum rates for the t-subprobiem in accordance with an embodiment of the invention.

[0044] Fig, 16 illustrates an algorithm for handling mfeasibi!ity in accordance with an embodi ment of the invention.

(0045] Fig. 17 illustrates an algorithm for QP with minimum rates for Group I in accordance with an embodiment of the invention.

10046] fig, 18 illustrates an algorithm for admission control for group 1 in accordance with an embodiment of the invention.

[00471 Fig. 19 illustrates an algorithm for maintaining feasibility by scaling up durations in accordance with an embodiment of the invention.

[0048] Pig. 20 illustrates an algorithm for maintaining feasibility through additional dwell times it) accordance with an embodiment of the invention.

[0049] Pig. 21 illustrates a structure of a basic online LP for EV charging in accordance with an embodimen of the invention.

[0050] Fig. 22 illustrates a structure of online BP for infrastructure protection and load shifting in accordance with an embodiment of the invention.

[0051] Fig. 23 illustrates a structure of quadratic program for load tracking it) accordance with an embodiment of the invention.

[0 52] Fig. 24 illustrates a structure of online LP with demand charge mitigation in accordance with an embodiment of the invention.

[0053] Fig, 25 illustrates real-time pricing with price cap in accordance with an embodiment of the invention,

[0054] Pig. 26 illustrates building blocks of optimization-based product feature in accordance with an embodiment of the invention.

[00551 Fig- 27 illustrates HVSBs connected in a delta configuration in accordance with an embodiment of the invention.

DETAILED DESCRIPTION

[0056] Turning now to the drawings, electric vehicle supply equipment and methods for supplying power to a set of electric vehicles using electric vehicle supply equipment in accordance with various embodiments of the invention are illustrated. An adaptive charging station (ACS) is a smart version of electrical vehicle supply equipment ( ^{" }EYSE) having dynamic adjustment features, BVSE generally can he any device which brings power to and/or fills an BVs Battery, and ar an intermediate Between an EV and a power source. F.VSF ear) utilize a variety of parameters including (but not limited to) voltages, amperages, current type, charging speeds, and/or plug types. Level1 charging (L.1 ) is generall the slowest form or charging and can connect an EV to a standard 1 10V or 120V outlet. Level 2 charging (L2) can provide additional voltage (generally up to 240V) and as such can provide a faster charge compared to L.l . Level 3 charging (Ϊ..3) generally uses up to 480V and can provide even faster chargi g than 1, 1 or 1,2. in several embodiments of the invention, the SAE ,1 1772 standard can be used to define AC charging I eveis. It should be readily apparent that other standardised systems for EVSBs can be utilized as appropriate including (but not limited to) CHAdeMG, SAE Combined Charging Solutions, and/or Tesla charging format and that EVSEs can include both alternating current (AC) and/or direct current (DC). Furthermore, the development of additional charging standards involving a variety of AC and/or DC charging profiles is contemplated,

{00571 ACS can be grouped together into an adaptive charging network (AC ). ACNs can be specifically designed for large scale d eployments such as ( but not l imited to) college campuses, corporate offices, airports, movie theaters, and/or malls.

{0058] In many embodiments, an ACM computes charging rates for EYs over a time horizon. In certain embodiments, an EV model is used to construct and solve a linear program (LP) over a rolling time window in several embodiments, a quadratic programming (QP) framework may be utilized to compute charging rates for EV , where the objective function is quadratic and the constraints are linear in certain embodiments, the system is an unbalanced 3 -phase system, where the capacity constraints are quadratic. Accordingly, certain embodiments provide a quadratically constrained quadratic program (QCQP). In many embodiments, the QCQP is a convex program that is polynomial-time solvable.

{00591 In man embodiments, the objective function may have several components such drat if is a weighted sum of the component function. Each component can be designed to achieve a certain purpose, such as charging as fast as possible, tracking a given signal (e.g., onsite solar generation, demand-response signal, among others), implementing different priorities, reducing temporal fluctuations, ensuring unique optimal solutions, among numerous o thers. Accordingly, different product features may be systematically implemented within a same QP framework.

[0060] In several embodiments, product features may lie implemented through the introduction of new constraints. Described in detail below are various features, including demand charge mitigation, demand response: load shifting/tracking, priority charging, joint EV/soiar/storage optimization. reducing temporal fluctuations, e «forcing m inimun r ates, a nd h andling infea i billy, among numerous other features as appropriate to the requirements of specific applications in accordance with many embodiments of the invention.

[0061) In particular, many embodiments provide a QP that provides load shifting/tracking. In particular, a demand response event may be to maintain a total ACN site load to below or at a certain load for a certain time period (e.g.. below 1 MW from 1 PM to 4PM today). Furthermore, the QP may have the total site load track a given profile for a particular time period,

[0062] Several embodiments of the QP may provide for priority charging among EVs by «sing an appropriate choice of parameter values. In particular, a higher priority EV can be assigned a larger weight, a larger minimum energy, and/or a larger maximum charging rate. In certain embodiments, a driver of an EV vehicle may pay a different price for prioritized charging.

[0063] In several embodiments, an EV may be grouped into different groups providing different priorities. The priorities may be implemented using different techniques, including admission control whereby an EV may be admitted or denied admission for charging based on the existing EVs within the ACN, and strict priority where certain groups are given strict priority over others, among various other techniques as discussed in detail below, In particular, different groups may be specified, including, for example, Group 1 where a defining feature may be that the EV is guaranteed its requeste energy whereas a Group 2 whereby it only guarantees a minimum energy (but not necessarily i s requested energy). Accordingly. Group 1 EVs may have strict priority over Group 2. in certain embodiments, admission control may be performed before a new BV is admitted to Group 1 in order to guarantee the requested energy of ail existing Group 1 EVs and the minimum energy of all existing Group 2 EVs.

[0064] in certain embodiments, an ACN site may provide for EV charging, solar generation. and a battery (onsite energy storage) and the QP may schedule EV charging and battery operation to track solar generation. In particular, the charging rates and battery operation may lie chosen in order to minimize the distance between the forecast solar generation and the total net load (e,g., EV -r background + battery draw). Certain embodiments may solve a convex optimization problem (QP), as described in detail below

[0065] In several embodiments, a QP framework may reduce rate fluctuations across a time period. In particular, the cost function may reduce the temporal fluctuation of the solution to a single QP instance, which is different from reducing fluctuations across d ifferent QP s olutions. Certain embodiments may utilize a penalty term to reduce temporal fluctuations, which can make the QP cost function strictly convex and hence provide a unique optimal solution.

(0066] In several embodiments, the QP framework may enforce a minimum rate. In particular, it may be undesirable to set a charging rate to zero before an EV has finished charging because this may cause the mechanical contact in the charger to open, and a charging profile with many zero and nonzero rates before an EV has finished charging can incur excessive wear and ear. As such, certain embodiments of the QP framew ork may set charging rates to a rate that is greater than zero as long as an EV has an energy demand,

[0067] In certain embodiments, the QP framework may handle mfeasibility in a situation where a minimum energy demand cannot be met for some EV, or when a capacity constraint may be violated at some resource. In certain situations, an EV rate may be determined based on laxity.

[0068] Charging stations and processes utilized to perform energy discovery protocol processes and determine EV charging rates for a set of vehicles in accordance with various embodiments of the invention are discussed further below,

ELEC TRIC: VEHICLE POWER DISTRIBUTION NETWORKS

[0069] A power distribution network in accordance with an embodiment of the invention is shown in FIG. I. Electricity is generated at power generator 102. Power transmission Lines 104 can transmit electricity between the power generator and power substation 106. Power substation 106 additionally can connect to one or more large storage batteries 108, which temporarily store eiectxi city, as well as power distribution lines 1 10. The power distribution lines 1 10 can transmit electricity from the power substation to an adaptive charging network (ACN) 112. ACNs 1 12 can include a battery 1 14, solar panels .1 16, and adaptive charging stations (ACSs) 120. Electric vehicles 1 18 can connect to the ACN {and therefore to the power distribution network) at the ACSs. Some charging networks can only charge a small number of electric vehicles, while others have the necessary infrastructure to charge dozens or even hundreds of electric vehicles simultaneously.

[0070} The power generator 102 can represent a power source including (but not limited to) those using fossil fuels, nuclear, solar, wind, or hydroelectric power, Substation 106 changes the voltage of the electricity for more efficient p ower d istribution. S olar panels 1 16 are distributed power generation sources, and can generate power to supply electric charging stations as well as generate additional power for the power grid.

[0071] Although many different systems are described above with reference to FIG. 1 , any of a variety of power distribution networks including adaptive charging station controllers may be utilized as appropriate to the requirements of specific applications in accordance w th embodiments of the invention. Adaptive charging stations connected to an adaptive charging network in accordance with various embodiments of the invention are discussed below.

ADAPTIVE CHARGING STATION CONTROLLERS

[0072] ACSs connected to an ACN in accordance with an embodiment of the invention are shown in FIG. 2. EV nodes 202 can connect to ACSs 204 using a wired and/or wireless charging connection. ACSs can connect to the ACN 208 using wired and/or wireless connections 210. ACSs can communicate with each other over this wired and/or wireless connection. In some embodiments, ACSs can communicate via data flowing irectly n t he p ower 1 ine. The adaptive charging network may also be connected to one or more centralized computers 212 to monitor calculations made by or to send instructions to multiple EV nodes. Additionally in many embod iments, a database management system 2 14 can be connected to the network to track EV node data which, for example, may be used to historically track power usage at various locations or at various times of day over time. In many embodiments, adaptive charging stations can use adap tive charging station controllers, where linked adaptive charging stations can distribute available charging capacity between upstream and downstream ACSs in the ACN. in several embodiments of the invention, ACSs can be connected in various topologies including (but not limited to) hub and spoke, bus, tree, daisy chaining, point-to-point, star, ring, mesh, and/or hybrid topologies.

[0073] Although many systems are described above with reference to FIG. 2, any of a variety of systems can be utilized to implement a network of adaptive charging stations connecte in an adaptive charging network in a variety of topologies as appropriate to the requirements of specific applications in accordance with various embodiments of the invention. Adaptive charging station controllers in accordance with several embodiment.» of the invention are discussed below.

[0074] An adaptive charging station con {roller (ACS controller) in accordance with an embodiment of the invention is shown in FIG. 3 In various embodiments, the ACS controller 300 can perform calculations to distribute charging capacity between linked ACSs upstream and downstream within an ACN. In many embodiments, an ACS controller can make requests for available capacity from upstream nodes and/or transmit available capacity to downstream nodes.

[0075] In the illustrated embodiment, the ACS controller includes at least one processor 302, an I/O interface 304, and memory 306 in many embodiments, the memory·' includes software including BV charging application 308 as well as EV parameters 3 10. adaptive charging parameters 312, and energy discovery protocol parameters 3 14. An ACS can calculate charging parame ters by using a combination of its own electric vehicle parameters, adaptive charging parameters, and/or energy discovery protocol parameters received through the I/O interface. Adaptive charging parameters can include specific c harging process parameters a nd/or o ptimization constraint parameters. Additionally, adaptive charging parameters can include parameters specific to adaptive charging tations and/or adaptive charging networks. Energy discovery protocol parameters can include (but are not limited to) parameters specific to available c apacity requested capacity from upstream ACSs, and/or transmitted capacity to downstream ACSs. In a number of embodi- ments, the ACS controller and/or the ACS includes a touch screen display that enables the operator of an BV to provide information concerning the B.V connected to at) AC'S and/or information concerning desired charging requirements te.g, information indicative of a power re uirement and an associated charging time such as (but. not limited to) departure time and/or desired additional miles to add to range of EV}. As is discussed further below, the ACS controller and/or ACS can also be connected to one or more sensors that can detect a vehicle occupying a specific parking space associated with the ACS when the vehicle is not drawing current from the ACS. In this way, the sensors enable the ACS controller to provide reliable information concerning the availability of the ACS to controllers within the ACN and/or operators of EVs more generally via web services

[00761 Although a number of different ACS controller implementations are described above with reference to FIG. 3, any of a variet of computing systems can be utilized to control an ACS within an ACN as appropriate to the requirements of specific applications in accordance with various embodiments of the invention. Various configurations of ACSs within an ACN in accordance with many embodiments of the invention are discussed below.

SYSTEM MODEL.

[00771 in many embodiments, processes are utilized that attempt to optimize charging of EVs within an ACNs in accordance with at least one optimization criterion, A system model for computing a charging rate for an EV using a linear program (LP) in accordance with an embodiment is described below.

Optimization horizon. Fix a time horizon T { 1 , 2. , 7 ^{' } [. This defines a rolling time window over which EV charging rates can be optimized repeatedly, as in model-predictive control. Specifically, at time t, ACN:

1 . assumes there w ill be no future EV arrivals (this is the online aspect);

2. computes the charging rate for ever}·’ EV that is active (e.g., still needs charging), by optimizing the charging rates of ail EVs over the time horizon t -f 1. , t + 2 + T;

3. for the period t + 1. charge all the active EVs at the calculated rates r(f -j~ 1 ) (n(l: ÷ 1 }. tor all EV i);

4. updates the .remaining energy demand of each BY, and possibly other state variables at the end of time t t 1 ;

5. updates new EV arrivals in time 1 - 1. if any; and repeats the procedure at time i r 1.

(0078] For example, if the charging rates are updated every minute and they are optimized over 10 hours, then r is in unit of minute and T— 600 minutes. Suppose the current time is i— 0 and there are 30 EVs i G 11.30]. Where the notation 6 j l _{< }«! and T\ denotes pectively. Then an ACN can compute the charging rates r{ / } :— For time period t ~~ L the ACN charges the 30 EVs at rates ri 1 ) : Then the ACN can update the remaining energy demand of these EVs and any new EV arrivals, and repeats the computation at time t ~ 1 to compute the rates

(0079] An electric vehicle (EV) model used to compute charging rates using a linear program in accordance with several embodiments of the invention is now described.

f¾080] EV mo el Many embodiments provide an EV model that is user! to construct and solve a LP at each time t to compute the charging rates r(r -f~ l ), r(/ 4- 2), + I ) over the rolling window j? 4 1 , T 4- 1), where rEvj is a vector of charging rates for time s. At the next time t 4 1 , a new I.P can be constructed for the computation of charging rates over the window (/ 4- 2. T 4- 2] In general, at each time .v, the ACN may compute the rates for the next T times in the window Is, T 4· $ - 1 1. lb differentiate each time t G is called the computation period.

[0081] Without loss of generality, the computation period .v— 0 may be focused on and the computation of rates may be considered over the time window \ K T\.

(0082] An ACN can describe an EV i using a vector (e;. ·. with the following interpre tation: • i.y > 0 is the remaining energy demand it is equal to the remaining energy demand at the beginning of the previous computation period minus the measured energy actually drawn by the EV in that period. 7 ^{' }] is the (expected) departure time ^{' }( ^{' }he initial departure time when EV i first arrives can be obtained from user input. For subsequent computation periods, di are the departure times in the previous periods minus 1.

[0083] Note, these parameters (e;. _{; }. /y i depend on the current computation period s ~~ 0 and will be updated at subsequent computation periods s ~ L .

[0084] Note that d; may be interpreted as the time when EV / is scheduled to departs and therefore the time available for charging is [Ld _{> } - I , not il . ,). This detail should not be forgotten because some of the design below may be modified if time available for charging is interpreted to be [1 , </,·], e.g , the cost function c;(r) in (23) should be modified to ci(r) - ( y -p i })/v(r) otherwise; see discussion below in section’’Charge as fast as possible” and Equation ( 18),

[0085] There may be two“corner eases” to handle. The first corner case is when an a new EV arrival has an expected departure time outside the window, _{<¾ }',· > T. One way is to assume T is large enough to satisfy any EV’s energy demand and set d; to be T. The second corner case is when an EV departs later than it originally specified (at computation periods after the departure time, d, will be negative). If an EV is not charged after its specified departure time, then d: · · 0 may be set after the specified departure time or the EV may be removed from the list of active EVs. This process has the nice feature that it discourages a driver from deliberately declaring a departure time that is earlier than the true departure time,

* ? (r _{f } ( i ) ? (1 , 71 ) is a vector of max charging rates. For the QP to be solved at computation period 0, this is a given vector it is determined based on a variety of factors such as EVSE limit, EV limit, or the ramp-down mechanism, e.g.. r; min ] EVSE limit, EV limit, ramp-down rate } By ^{« }setting r,(/ } -·· 0 for t > d he QP framework may not charge after EV fs specified departure time.

• r,· is a minimum charge rale. It is undesirable to set a charging rate to zero before the EV has finished charging because this will cause the mechanical contact in the charger to open. A charging profile n{t) with many zero and nonzero rates before the EV finishes charging can incur excessive wear and tear. Moreover, some EVs may not recover once its charging rate is set to zero, he., it will not resume charging until it is unplugged and plugged in again. Hence, many embodiments enforce minimum rates whenever possible.

• e, > 0 is a minimum energy that may be delivered, if <¾ ··· 0 then there is no minimum energy guarantee. This is for Group 2 EVs (as described in detail further below).

A NOTE ON R;(T)

in many embodiments, the AON uses an ac (alternating current) system, whereby the currents, voltages, and instantaneous power flows are sinusoidal functions of time { nstantaneous power has twice the frequency as currents or voltages). Accordingly, h(ή may be used to sometimes denote power, and .sometimes current. When h{() refers to current, it is the RMS value (which is equal to the magnitude of the complex current phasor). When r;(t) refers to power, it can either by the real power (in kWh) or the apparent power (in KVA). EV charging stations may draw power at unity power factor, reactive power is zero and apparent power equals the real power, and hence the distinction may not be important in these situations. In several embodiments where products provide reactive power compensation for grid services, then there may be a need to differentiate between real power and apparent power.

(00871 In many embodiments, most parts of the optimization processes may use current instead of power. If the voltages are fixed (both their magnitudes and phases), power is directly propor tional to current and either can be equivalently used. However, using currents may have several advantages including the following; • Line limits and breaker limits etc may all be specified in terms of current limits.

* Pilot signais/CT measurements can be in current.

• Whenever there may be a need to calculate power or energy, many embodiments may multiply available voltage with current (and time for energy), so if voltage changes, there may be no need to have to modify other parts of the process.

* Dealing with currents and voltages may stay in a linear system and can be simpler to handle.

Accordingly, processe in accordance with many embodiments of the invention can use n(t) to denote power or current as appropriate to the requirements of a given application.

QUADRATIC PROGRAM

[0088] Many embodiments use a quadratic program to compute the charging rates for a set of EVs over a time period. A QP framework for computing charging rates for an EV in accordance with several embodiments of the invention is illustrated in Fig. 4.

{<1089] Given (as illustrated in Pig 4):

• an optimization horizon (1 , 7 ^{' });

• a set of EVs specified the basic form of a quadratic program (QP) to compute the charging rates over [1 , 71 these BVs is: min c{r) ( l a) subject to

T ^{' } ,4 _{/i }qf/ ) < Q(/ ) _{? } t€ | L 7 h for all resources / ( i d) where c(r) is a convex quadratic cost function, / (/ ) denotes the capacities of resources / at time r, and , arc the coefficients that describe: bow EV i’s are constrained by resources / (as described in detail below). Note that for each QF instance, t can be not the real time but t 1 , 2 . .. .7\ . Although Fig 4 illustrate s a particular QF framework for computing charging rates for an EV, any of a variety of frameworks, including linear program with linear constraints, quadratic progrant with linear constraints, and quadratic program with quadratic constraints may be utilized as appropriate to the requirements of specific applications i n accordance with several embodiments of the invention. Capacity constraints for objective functions in accordance with numerous embodiments of invention is now described.

CAPACITY CONSTRAINTS

[0090J An objective function of a QF framework may compute charging rates for EVs based on the capacity constraints of a particular ACN. The capacities may vary based on the capacities of the cables connecting the chargers. Furthermore, the QF may use linear capacity constraints or quadratic capacity constraints when determining an optimal EV charging rate. Computation of capacity constraints may be based on the following.

[009 lj !t may be convenient to express the linear capacity constraints ( i d) in matrix form:

L p h < Pa ) where the matrix A has ,½ as its UJ) entry, r(t) (h(ή,ΐ ( Lnj) is the column vector of charging rates, and P(i) (/)(;) ,¥/) is the column vector of resource capacities

[0092] For example, if a set of 8 chargers are fed by a cable with capacity 80.4, the corresponding capacity constraint is

Q (Q

where B (?) ~ 80A for all t. This can be put in the form of ( 1 d) below, £^4 _{/ }pv(r) < P((t), by setting

1 if EV i shar es resource /

¾ - 0 otherwise

Sometimes, 5 _{/ } ma be used to denote the set of EVs that share resource l.

[0093 J For example (Single -phase or Y configuration), suppose a pane! is fed by a cable with capacit P _{\ } (constant for all time t), It serves a set H _{\ } of EVs and a subpanel with capacity /¾ that serves a set ,¾ of EVs, as shown in Fig, 5. Note, as illustrated in Fig. 5, each EVSE can be connected between a live wire and the neutral, e.g,, either in a single-phase configuration or F-configuration of a three-phase system. Suppose the capacity of the cable between the (sub)panel and an EVSE is not a bottleneck,

[0094] Then the capacity constraint on the subpane! is: and the capacity constraint on the main pane! is:

To write these two constraints in the form of ( I d) below, define if EV / Si

otherwise i†EV / £ S;

otherwise

Then the constraints (2) become

Let:

* > {() € ¾) denote the column vectors of charging rates for EVs in ,Vj and ¾ respectively; A ^{2 } (Ay, i€ .¾) denote the row vectors corresponding to EVs in S; and Sz respectively; then these two capacity constraints can be written in matrix form as:

[0095] Solving the above QP subject to the specified constraints can provide a time sequence of charging rates for each EV. However, additional constraints may be provided to further direct the manner in which available capacity is assigned to individual EVs by the ACN.

{Wm\ Although Fig. 5 illustrates various single-phase or Y configurations, any of a variety of configurations stray be specified as appropriate to the requirements of specific applications in accordance with various embodiments of the invention.Futthermore, solving the above QP subject to the specified constraints can provide a time sequence of charging rates for each EV. However, in many embodiments, additional constraints may be provided to further direct the manner in which available capacity is assigned to individual EVs by the ACN. The following section describes the manner in which the QP can be implemented in the context of a 3-phase system.

REFERENCE DES IGN 100971 In many embodiments, EVSEs are connected to three-phase circuits. Accordingly, many embodiments adapt the QP framework, as outlined above, within the context of a three-phase sys tem. Described now are QP frameworks within the context of three-phase systems in accordance with various embodiments of the invention. A Reference Design for a three-phase EVSE in accordance with many embodiments of the invention is illustrated in Fig 7a. An electric panel in accordance with several embodiments is illustrated in Fig. 7b. A circuit model in accordance with several embodiments is shown in Pig. 6 where the secondary’ side of a distribution transformer can be modeled as a three-phase voltage source on the left in Wye configuration that upplies a set of three-phase loads (charging stations) in parallel in Delta configuration. In p articular. Fig. 6 illustrates a three-phase design whereby a three-phase voltage source is connected to loads in parallel in accordance with various embodiments of the invention.

100981 Referring to the circuit model in Fig. 6, in many embodiments, the QP may include capacity limits on the line currents fro the secondary side of the distribution transformer (the three-phase voltage source) to the electric panel. Let I ^{mK } denote the line limit for each of the three lines, e.g„ require that the current magnitude be bounded by / ^{aK }:

The line limit / ^{!:<sa : } may depend on both the power rating of the transformer and the current rating of the wire connecting the .secondary side of the distribution transformer to the electric panel it is derived in detail below as (?):

100991 Certain embodiments may set bounds on the magnitudes of the charging currents k 1.2, . . .. that enforce the line limit (3). In certain embodiments, the circuit model illustrated in Fig. 6 may use a linear bound which is simple but conservative may take the following form (from (21 )):

[0100] The optimization variables rgf ) may denote the magnitudes of the charging currents 0101] A second embodiment of the circuit model illustrated in Fig. 6 may use a quadratic bound instead of linear bounds which can be more computationally expensive but may be tight can take the form (from ( 19)):

Although Fig 6 illustrates various three-phase designs where a three-phase voltage source is connected to loads i parallel any of a variety of designs may he specified as appropriate to the requirements of specific applications in accordance with various embodiments of the invention. Three-phase transformer designs in accordance with various embodiments are now described be low.

TRANSFORMERS

[0102] Three -phase transformers that can be utilized in the Reference Design illustrated in Fig 7a in accordance with various embodiments can be specified base upon limits on the line currents feeding into a three-phase transformer, in either Wye or Delta configuration. The three-phase transformer in the Reference De gn (as ill «strafed in for example. Fig. 7a) has the specification voltages : 480V - 120/208 V'

power rating : 75kVA

This means that the printary side may be three-phase arranged in Delta configuration with line- io-iine voltage magnitude of 480V, and the secondary side is three-phase arranged in Wye configuration with fine -to· neutral voltage magnitude of 120V and Hne-io-Hne voltage magnitude of

208V— y 3 · 120V. Fig. 8 illustrates a Reference Design transformer with voltages 480V ί 20/208 V and power rating of 75k VA in accordance with several embodiments of the invention. Sometimes this may be also more explicitly specified as 480L 208K / 120 with 208V being the iine-to-line voltage on the secondary side.

{0103] The power rating of a three-phase transformer typically refers to the three-phase power which is typically three times the single-phase power 5} _{h }ase (assuming balanced operation). The magnitude of the single-phase power for the Reference Design can therefore be

5 kVA

01 4] The power input on the primary side is typically equal to the power output on the sec ondary side under the assumption that the transformer loss is negligible. The power rating may be the maximum power the transformer can handle on the primary or secondary side. It implies a capacity limit on each phase wire on the primary as well as secondary side. In many embodiments, these limits are used as capacity constraints on the charging currents when determining charging rates for each EV, Various techniques for calculating these current Hunts in accordance with certain embodiments of the invention are described below. (01051 On the primary side, the configuration may be Delta nd hence the single-phase power

Stptew is given by (see Fig. 8}

_{* }^1 phase Kr// -7/ i Sii) where 1 ^{* },^ denotes the complex conjugate of the phase current _{/ },’. The phase current _{¾' } is related to the line current by (assuming balanced operation in positive sequence)

and hence

Substituting (5b) into (5a), results in

Hence the current limit on each line is

(0100] Since the power rating of a three-phase transformer may always specify the three-phase power ¾pjme instead of single-phase power 5} ! _{. }· «.- _{< } (6) >- often expressed as [0107] For the Reference Design, the current limit on each line is;

1 75 kVA

90 A

y34$0 V

[0108] In many embodiments, there is usually less concern with the current limit on the primary side because the goal of the system is to manage charging stations connected to the secondary side [0109] On the secondary side, the configuration may be Wye and hence the complex powers can be related as (assuming balanced operation)

[01101 since ~ ~ «, Tins and the power rating implies a limit on the line current on the seconds V side;

L 208.3 A

[01 11] This is the current limit due to transformer power rating.

[0112] The wire connecting the secondary side of the transformer and the electric panel may also have a current limit. The limit used in the capacity constraint of various optimization processes described: wire current limit} (7)

[0113] Although Fig, 7a illustrates various three-phase transformer designs, any of a variety of transformer designs may be utilized as appropriate to the requirements of specific applications in accord acne with various embodiments of the invention.

CURRENT LIMITS ON LOADS

[0114] Described in detail now are various embodiments of the QP framework that optimize a constraint for a current limit imposed on each charging station by distribution transformers. In particular, in various embodiments, the current limit in (7) on each line connecting a secondary side of a distribution transformer to an electric panel may impose a current limit on each HVSE. and the QP framework may optimize the constraints based on this current limit,

[01151 Figs. 9a and 9b show common distribution transformer configurations i n accordanee with several embodiments of the invention. In particular, Figs. 9a and 9b show the secondary sides of two common configurations of three-phase distribution transformers in accordance with various embodiments of the invention. Fig. 9a shows a 5·· wire three-phase transformer in Wye configura tion, as in the Reference Design, in accordance with various embodiments of the invention. Three phase wires (labeled U K C) and a neutral wire (labeled n) are shown. The fifth wire, not shown, is the earth ground wire, typically connected to neutral, A different voltage magnitude can be supplied to a load depending on how it is connected. The voltage magnitude between a phase wire and the neutral is 120V and that between a pair of phase wires is ! 20v 3 - 208V,

[0116] Fig. 9b shows a 5-wire transformer in Delta configuration with one of the phases center- tapped to provide three voltage levels in accordance with various embodiments of the invention. Four phase wires ( labeled aJh c, d) are shown but an earth ground wire is not shown. The voltage magnitude between wires ad or bd is 120V, whereas that between wire cd is 208V. The iine-to-ime voltage magnitude is 240V. Although Figs. 9a and 9b illustrate various distribution transformer configurations, any of a variet of transformer configurations may be utilized as appropriate to the requirements of speeifie applications in accordance with various embodiments of the invention.

[0117] Fig. 10a shows a Wye-configured voltage source c onnected to a set of loads in Delta configuration in accordance with various embodiments of the i nvention. The voltage source is the secondary side of a three-phase 120V /20SV ^{' } transformer shown in Fig 9a. The voltage magnitude across each load is the Hne-i.o-l.ine voltage 208V. This is one of the configurations described above in the context of the Reference Design. Although Fig. 10a illustrates various Wye-coniigured voltage sources, any o a variety of configured voltage sources may be utilized as appropriate to the requirements of specific appl ications in accordance with various embodiments of the invention, £0118] Fig 10b shows an electric panel arrangement to connect the loads to the voltage source in accordance with various embodiments of the Invention. The dot in the first row indicates that the wires numbere 1 and 2 a connected to phase «, the dot in the secon row indicates that the wires numbered 3 and 4 are connected to phase />, the dot in the third ro indicates that the wires numbered 5 and 6 tire connected to phase c, and so on. Therefore the load connected between wires 1 and 3 can be connected between phase a and phase b lines (see the corresponding labels on the loads in Fig. 10a. Similarly for the load connected between wires 2 and 4, and other loads connected between different phases. Although Fig. 10b illustrates various electrical pane! arrangements to connect loads to voltage sources, any of a variety of electrical panel arrangements may be utilized as appropriate to the requirements of specific applications in accordance with various embodiments of the invention.

|0U9] I many embodiments, there is an interest in the currents in supplied by the three-phase source to the loads. Suppose the wires connecting the three-phase source to the loads are rated at / ^{m }‘ ^{,x } . e.g„ it may set that the current magnitudes be bounded by / ^{ax }: / ^{ma }* p ii . h . i ' ^{( }8) where P ^{m } Is derived in (7).

10120J Suppose the loads are not impedance loads, but constant current loads that draw specified currents. Let the current drawn by the load in Fig. 10a between wires 3 and 3 be 4p. _{>; } . that between wires 9 and 1 1 be _{(G( } , that between wires 5 and 7 be 4, _{« }, . In general, let the load currents in the kth three-phase load be

4,,

Bounds on the load currents may now be derived ( _{* } A ~ 3 , . . . . K) that enforce the line limits (8). Recall that the magnitudes j/ _{/w } [„ p · a,b.c _{> } of the current phasors can be the root-mean-square (RMS) values of the corresponding sinusoidal functions. In practice, load currents are usually specified in terms of RMS values, and hence the bounds of interest depend on the magnitudes 1 \ of the current phasors, not their phases.

[0I21 | Applying KCL at nodes ( i .b i , Cj ) provides

where ( _{|ίίί }., ^ r _{¾ ., ίί }· _{ί,; : } } _{> } k ·--- 0, . . , . K ~ 1 , are the line currents from stage k to stage k -i- 1, in general, it has

. ···.·,. L/* 4 .4...; · k _{^ } O . K 1

Hence the total supply currents are given by when there are K three-phase constant current loads. Note that this expression does not require that the loads are balanced./ Note however that the line limit l ^{miK } in (?) is derived assuming balanced operation and is therefore only an approximate limit.) In particular, if a load (say) .^ is absent, then set ,.,¾ - 0 in (9)

[OI22j Let the total load current in each leg of the Delta configuration he denoted by

Then (9) can be written in terms of the total load currents as:

The line lirails (8) are ther fore

Enforcing line limits hence may need one to know not just the magnitudes j of the load currents, hut also their phases in order to compute their sums,

{§1231 Fig, 1 1 illustrates I _{((i <· ™ } } ^{at> }— l ^{a> }. Hence by the cosine rule \ - the angle between t ^{it,} } and l ^{m } .

As explained with reference to Fig 1 1, these inequalities are equivalent to;

[0124] If the angles p ~ aj>,e, are know between the total load currents (i ^{iih }J ^{h } J ^{i }) in each leg of the Delta configuration, then ( 12) are convex quadratic constraints on the magnitudes o i (r ^{,h }J ^{lk } _{) }.

[0125] ACNs can use the processes described above to determine the charging rates to provide to a number of EVs, when the ACN is supplied by three phase power. In the real world, the grid is less than ideal and the three-phase power supplied to the ACN may be unbalanced. Accordingly, many embodiments solve the QP for three-phase unbalanced supplies.

UNBALANCED THREE-PHASE INFRASTRUCTURE CONSTRAINTS

{01261 In many embodiments, EVSEs are usually connected on three-phase circuits. The EVSEs may be connected in a delta configuration as shown in Fie 27 in accordance with various etnbod- iments of the invention.

[0127] Because of differences in demand the loads in this delta configuration are often imbal anced which may require for careful consideration of the infrastructure constraints to ensure safe operation. For modeling simplicity, certain embodiments may assume that the line impedances leading to each station are negligible, allowing to lump all EVSEs between common phases into a single load represented by current phasors Certan embodiments may also assume that each EVSB is modeled as a controllable current source with unity power factor. With this model, certain embodiments may consider two types of constraints:

* Phase Constraints: A phase constraint pertains to the current along one leg of the delta, i.e. ah, he, or ca. For example constraints on the aggregate current of the A V and CC pods are phase constraints in an ACN. Certain embodiments may tieat the EVSEs within the same phase as in parallel and all have unity power factor, and can express the constraint as Lex, nit ^{) } < Rj vj

* Line Constraints: A line constraint pertains to the current along each line i.e, o. k. or e. Examples of line constraints in an ACN include limits on ή\, ij ,, l~ and /; where p€

10128] Explained in detail now- are how to derive the constraints on these current magnitudes p€. {aJ.i.c} _{> } using the circuit diagram in Fig. 27. Certain embodiments can calculate the line currents /:) from the phase currents; : /

4 - > ^{·· } , where each variable is a phasor. From this point on, certain embodiments will only consider one phase/) me in the interest of space, but all other constraints follow from this derivation. To find certain embodiments may defirse the set of al l EVSEs connected between lines a and b to likewise for be and ca. Certain embodiments can then define the magnitude of the aggregate phase current for each leg of the delta as

101 9] Certain embodiments can treat each EVSE as a constant current load with unity power factor, so the phase of each current matches the phase of the corresponding voltage. Certain embodiments may assume that they are able to measuie/caleuiate the phase angle of the voltage across each leg of the delta configuration. Denote the phase angle of each phase as ø¾·.·, and f _{a< }· re spectively. If measurements of voltage phase angles are not available, certain embodiments may assume that voltage angles are balanced, he,, each phasor is spaced 120 apart.

1(11301 In any case, certain embodiments may emphasize that in the phasor only the magnitude variable and the phase

Is known,

10131 ] From ( 1 3), the current constraint |/|| < R _{$M } becomes a constraints on Note that this constraint is a second-order cone fSOC) constraint in the magnitudes j/f? !- To see this, notice

{ ; ί4Ί sincjw, - j/f] Ί sin<j> ) ^{2 }

|0132]} In order to account for constraints on / , /.] and I , certain emhodlmetns may consider the effect of the delta- wye transformer ti , Using circuit analysis, certain embodiments can relate f to the aggregated EVSE currents: where n is the turns ratio of the transformer which in our system is 4. Hence the constraint on l~ can be expressed in terms of EVSE current magnitudes as: where /¾. _{« } is expressed as a current constraint on the primary side of / _{s } rather than reflecting it to the secondary side.

101331 Finally, certain embodiments can obtain /J and ! from / by adding the currents drawn from the DC fast charging and the auxiliary garage loads. Hence its constraints are:

Like (16a), the constraints ( 16b), ( 16c) and (16d) are SOC constraints. These constraints translate into constraints on the charging rates r.(/l through ( 14). [0134] In some applications these SOC constraints arc too computationally expensive, however. Simpler but more conservative constraints can be derived by observing

Hence the constraints { 16) can be relaxed to:

SPECIAL CASES

[01351 Described below are several special cases and the derivation of simple bounds on the magnitudes of the individual load currents that will enforce ( 12),

[01361 Assumption l : Current phasors I _{(Hj } have the same, and known, phase angle 8«/> for all

and constraints ( 12a) become where cosf _{rtj. 7i } is known. Similarly for constraints 12b) and (12c). These are quadratic constraints in the magnitudes of the individual load currents that will enforce ( 12), given the angles _{I ;!> }, , p ~ a, b c, between the load currents in different legs of the Delta configuration

£0137] Assumption 2; in addition to Assumption 1, the angles § _{>i ·~· } 120°, for p -.·. a,b,e. Then cos f , _{ί>( } -- 1 /2 and f 18) becomes

Similarly for constraints (12b) and ( 12c),

[0138! Assumption 3 (balanced case): ,4/7 load currents have the same magnitude and the phases of currents on different legs of the Delta differ by 120°, Tha is, assuming positive sequence, for ail k 1 _ ,A”, this provides where / is the common magnitude of the loa currents and

, t

Then the constraint ( 19) reduces to 3K ^{2 }l ^{2 } ( {R ^{, ά' }'U or a bound on the common magnitude / of individual load currents

[0139] Linear hounds fn many embodiments, an application may operate in unbalanced conditions, e.g., adaptive electric vehicle charging where the magnitudes !/ _{<·>< } of the load currents are to he determined and generally different. In these cases there can he two difficulties with the line limits ( 18) and ( 19). First, the angles (0 _{« }*,0 _{& }·,0 _{« },·) may not be known. Second, even when these angles are known, the constraints may be quadratic which can be computationally expensive to implement in real time in inexpensive devices. In this case, many embodiments can impose linear constraints which are simpler but more conservative.

[0140] Take phase a as an example. Since 1 ^ j ~ \p ^{m } - f ^{a> } \ < j/" ^{ft }j -t \i ^{w }\. a simple Hmit on the load currents that enforces \I _{U } | < P ^{K } can call for e.g.. the sum of the magnitudes of the total load currents in legs ah and ca should he less than the current rating / ^{m:i }V From ( 10), this provides ~ | £* <A I < å _{k } ! pU· Hence a simple linear bound ott the load current magnitudes is:

The constraints on phases h and c can he similar.

101411 For a balanced system, many embodiments can easily assess how conservative the bound (21 ) is compared with the exact limit (20) on the load currents. In the balanced case the bound (21) reduces to

Hence it is y 3/2 ~ 87% of that in (20), be., it is conservative by 1 tor a balanced system.

MULTIPLE EV GROUPS

[0142] In many embodiments, the ACN differentiates between different EV groups to provide different charging guarantees depending upon the group to which a particular EV belongs. In several embodiments, an EV can he moved fro orse group to another based on network state as well as its own condition. The interpretation of each group can be revised over time. In addition the process for determining the charging rate to provide to each EV in a group can be updated independently of other groups by maintaining interfaces and assumptions on each group.

GROUPS AND PROPERT IES

[0143] 1» several embodiments, the ACN divides EVs into at least three groups. The basic assumptions on each group are as follows,

1 , Group 1; EVs admitted to this group will be guaranteed (subject to certain assumptions) that the requested energy e,· will be delivered before the departure time fo, When an EV requests to receive service as a Group 1 EV, tests will be performed (detailed below) to determine if the requested energy e _{> } can be guaranteed under normal anticipated conditions, and an admission decision will be made and communicated to the driver. In many embodiments, the ACN may charge a driver of an EV a high price to secure guaranteed energy delivery.

2. Group 2: EVs in this group will be served with best effort, i .e.. there is no guarantee that requested energ ty will he delivered before time _{< }i,·. They will, however, lie guaranteed a minimum energy <¾·, which can be zero or positive, before time d For instance, an EV is guaranteed at least 15 miles of energy or 15% of requested a,·. The determination of ¾ is not specified as a property of Group 2, but can ire. decided and updated separately by another module. Note that initially, even e, cannot be fully guaranteed if > 0. In order to guarantee nonzero the process to determine «¾ shoul include a check on feasibility and can be reduced until feasibility is maintained. In many embodiments, the ACN can choose internally nonzero e, but they are not revealed to the drivers and therefore EVs i are likely to receive at least c _{; } though the drivers are not guaranteed this,

3. Group 3: This group contains EVs that cannot be served normally. Example: the driver has not provided input or has provided inconsistent input, or the driver input implies a laxit Oi < 0 tinfeasibie), or requested energy e _{< } has been delivered but the EV can still take more charge, etc. The simplest policy is not to charge any BV in Group 3. Since any EV that does not conform to Groups 1 and 2 for any reason may be sent to Group 3. it is likely that some of them may be handled differently from others. Accordingly, certain embodiments may subdivide Group 3 into subgroups. This may allow for updates to the processes and its implementation for handling one subgroup independently of other subgroups.

101441 While specific groups of HV's are described above, ACNs in accordance with many embodiments can utilise any of a variety of processes that treat different groups of BVs in different ways as appropriate to the requirements of a given application. Use of group priority in determining charging rate in accordance with a number of embodiments of the invention is discussed below.

GROUP PRIORITY

{01451 in many embodiments, the defining feature of Group I may be that EV i will be guaranteed its requested energy e _{s ( }up to a measurement error margin). This is in contrast with Group 2 that may guarantee only a minimum energ ¾· > 0, but not the requested energy a. In order to guarantee _{<? },·, a key design decision is:

Group 1 has strict priority over Group 2 and Group 3.

This means that all resources will be devoted to satisfy Group 1 BV's energy requests e; first. Group 2 EVs will be allocated left-over capacities. Specifically the charging rates of Group 1 and Group 2 EVs are determined sequentially:

1. Solve QP {with equality energy constraints) for Group 1 EVs only. Denote their rates by has (>' ^{* }(/),! - Group 1./ { 1 , 7 ^{' }j ) .

2. Compute left-over capacities for Group 2: where i (t) are the original resource capacities 3. Solve QP (with inequality energy constraints) for Group 2 EVs only using the left-over capacities

[0146) initially, certain embodiments may serve only Group 2. It is expected that most EVs may be in Group 2 even when Group 1 service is offered in the future.

10147) In certain embodiments, initially. Group 3 EVs may not be served. EVs in Group 3 may be re-assigned to Group 1 or Group 2 for service (e.g., an EV without driver input may be assigned default parameters (o,,ί/,,r; ~ 0) by ACN and sent to Group 2.), may be allocated charging separately from Groups 1 and 2 using left-over capacities, or may not be served, base on its subgroup.

j0148] Some of the key parameters of these groups are summarized in Fig. 12. As can readily be appreciated, any of a variety of parameters can be utilized in determining charging rates for different groups as appropriate to the requirements of a given application

DRIVER INPUT

|0149) Described now are processes for obtaining driver inputs in accordance with various embodiments of the invention. The processes may check if the inputs are valid, determine parameters for rate optimization, and assign an input request to an appropriate group,

PROCESSING EV INPUTS

[0150] When a new EV arrives, its driver can be expected to input via a user interface on an EVSB or mobile application information (including but not limited to) their energy demand <¾, departure time <4 and/or tin when group 1 service is offered) which group they desires !n many embodiments, the inputs can be automatically provided by the EV and/or a remote service that uti lizes data analysis and/or machine learning to estimate specific parameters such as (but not limite to) estimated departure time irrespective of the manner in which the information is received, some sanity checks may be performed, including:

1. The requested energy <?,· is too close to zero (e.g., less than the minimum rate ri) or negative, ask the driver to re-enter e _{t } or assign the EV to Group 3 for further handling.

2. The departure time </,· exceeds optimization horizon T, .reset d- ··-·.·. 7\ This means that an EV may not be charged for a duration longer than T even if the EV stays beyond T.

3. The laxity is negative, i.e., impossible to meet the requested energy even if the EV is charged at its peak rate f . throughout its duration i l - d; lj See Laxity and Group Assignment discussion below.

[01511 in addition to sanity cheeks, an ACN may also determine for each new EV request:

1. Its peak charging rate L, e.g., based on driver subscription information, EV make and model or ACN real-time testing on the E V.

2. Its minimum rate r,, e.g., based on EVSE· limit, or EV make and model

3. Its minimum minimum energy de aud _{<? },·, e.g., based on site policy. ACN policy, or user input. This parameter might lie changed in the middle of the charging session, e.g., if an EV is re - assigned to a different group.

[0152 j While specific examples of ways in which data received by the ACN can he validated to assist with the development of viable charging schedules, any of a variety of processes cat) be utilized to validate data received by the ACN as appropriate to the requirements of a given application in accordance with various embodiments of the invention.

LAXITY AND GROUP ASSIGNMENT

{(1153) For each new EV 0 request (<¾,h¾), the input module should first check its laxity and decide if it should be sent to Group 3, using a process such as (but not limited to) Algorithm 1 illustrated In Fig. 13. This group assignment may be performed on each new EV request. Other group re-assignments may also be needed, e.g., when the remaining energ demand c, has dropped to (approximately) zero as EV i is served (part of Group 3 design). In many systems, the ACN does not provide any guarantees and there is no Group 1 as described above. Accordingly, the following discussion starts by examining cost functions that can be utilized by AC N s in accordance with a number of embodiments of the i nvention to allocate charging capacity in the scenario corresponding to Group 2 EVs described above. A variety of cost functions that accommodate different objectives from priority to switching frequency can be utilized. The manner in which an ACN can provide specific guarantees with respect to charging of specific Group 1 EVs is then discussed.

GROUP 2

f 0154 j This section discusses specializations and modifications o f t his b asic form ( I ) to implement various features. In several embodiments, the first few features may be implemented by appropriate implementation of the cost function dr) utilized by the ACN to determine charging rates. The other features may be implemented mainly in the constraints. However, in certain embodiments, cost functions and constraints may not be chosen independently of each other, and it can therefore be important to understand their interactions.

{0155J The capacity constraints {Prft ) A i l , 7h,VY) may denote left-over capacities if Group I service is offered

CHARGE AS CLOSE TO £ _{/ } AS POSSIBLE

{0156 J In several embodiments, there may be two ways to charge as dose as possible to the specified energy demand ty, The first is to enforce it as a hard constraint:

101571 This approach may have the disadvantage that the QP can be infeasible, when EVs do not have sufficient laxity or the infrastructure does not have sufficient capacity, This can be used for Group 1 EVs, |0158| A second approach, as in ( 1), relaxes the equality constraint into an inequality constraint:

10159] This approach may be used for Group 2 EVs. With this inequality constraint, it may be important that the cost function c(r) is decreasing in

to drive £,L J / ^{* },·{>} to e,. for example, use for some constant > 0, subject to the inequality constraint ( lc).

CH ARGE AS FAST AS POSSIBLE

101601 t Jsing a cost function that is increasing in time / may encourage charging as fast as possible. An example is to modify ej (r) in (22) into:

Re-iterating that for each QP instance, t is not the real time but /— 1 , 2, . . . , T.

10161] This cost function may have the following properties.

1 . For t < e/;. {t— d )rii ) is decreasing in h{() and hence it attempts to charge as close as possible to e, with inequality constraint, as discussed hi detail above.

2. For i > </,·, ( t -~ d _{j })n(t) is increasing in />·(*) and hence it drives r _{; }(() towards zero at. times after the specified departure time /,·.

3. At time t ·· - d _{! } eWr) is independent of n(d·). At times of congestion, this will drive /y(c( ) to zero. Otherwise, its value may depend oft the QP solver and initial condition. For instance, if the process always initialize h(ΐ) · · 0 for all i.J before calling the QP solver, then it is likely that the solution will have r _{} }{d _{} } ) 0.

4. The cost function coi/) is increasing in t and hence it encourages charging as fast as possible.

10162) !n several embodiments, a cost coefficient should not be t - ai . This factor does not necessarily encourage charging as fast as possible since t - a,· > 0. It has the effect of giving priority to BVs that arrive earlier.

j0163 j An alternative to the cost function in (23) is

Since / T < 0 and is increasing .in /, use of this cost function by an ACN encourages charging as fast as possible and as close to «y as possible. Note the use of the cost coefficient t ^{. } d in

(23) instead of / ···· T may not have the effect of prioritizing BVs with earlier (or later) departure times implementation using the t T approach may seem simpler though the use of / · ύ· might offer certain advantages (e.g.. drives h(ί) ~ 0 for / > d·}. Unless complication arises, certain embodiments may use t -- d _{; , } As is discussed below, ACN in accordance with many embodiments of the invention ca further prioritize the charging of specific vehicles by explicitly adding priority constraints to the process for determining charging rates for individual BVs.

INCORPORATING PRIORITY

[0164) Many embodiments may add priority among EVs, e.g., to prioritize BVs that arrive earlier, or EVs that have lower laxity. Certain embodiments may assign different priorities based on a driver paying different amounts for charging a vehicle depending upon the charging rate /guarantees provided by the ACN. For example, a driver may pay extra in order to charge at a higher charging rate. Suppose each EV i has a parameter a,· arid, everything being equal, one wishes to charge EV i faster than EV j if 0 < m < a · . For example, a,· may be its arrival time, or a, may be its laxity

Suppose the departure time is greater than the current time, dt > t, then a < i . Moreover a-€· (0, 1 ] if and only if it is possible to deliver the requested energy t within the available time d, - i, assuming the infrastructure is not constrained. Otherwise, if «,· < 0 then it is impossible.

{0165] To prioritize EVs with smaller a„ the cost function in (23) may be modified to {provided

> 0):

f0166] In this way. the cost t can be weighted by a decreasing function of parameter m. As can readily be appreciated this weighting cast be achieved in different ways. The following example illustrates why this may prioritize smal l a, over large

EXAMPLE. OF UTILIZATION OF CHARGING PRIORIT Y

[0167] Consider charging two EVs I and 2 each with an energy demand of e* ~ to— 1 unit, over two time periods t ~ 1 .2, i,e„ </,· T 3 (not 2). Suppose m < <¾, i.e, want to prioritize EV I over EV 2. Suppose the capacity is I unit. The cost is: and the constraints are: ry ( 1 ) 4· r _{/ }(2) - I for i ~ 1 ,2

rj (?) f n(/ ) < 1 for t 1.2 Substituting h(2} I— r.( l ) into f>) to obtain:

Therefore to minimize cost 03 (r), we set. r, ( i ) ~ I and r ^{* }{ 1 )— 0 since 0 < a _{\ } < at _{* } ant! r (l ) 4· The optimal charging profile is r|( ! } ^{~ } 1 , r _{j } (2) ~ 0; and r _{2 { } i } ~ 0, r _{2 }(2) ~ .1 e.g., EV are charged as fast as possible and EV I enjoys a higher priority over EV 2.

fOlbSJ instead of l /o _{· }. many embodiments may also use other decreasing functions of _{< }¾ in tire cost function. For example, ACNs in accordance with many etnbodiments of the invention utilize the following cost function where the priority function /{«) is positive and .strictly decreasing in a. For example, to prioritize early arrival times a., we may use i— a _{t } instead of 1 /«*. If <¾ denotes laxity then a necessary condition for feasibility is m€ (0. 1.) and hence, instead of l /u,, one may use I The next example illustrates that any decreasing priority / prioritizes EV’s with smaller u _{t }.

ADDITIONAL. EXAMPLE OF PRIORITIZED CHARGING

f 0169] ACNs in accordance with a number of embodiments of the invention utilize the following moat generalized cost function in scenarios similar to those outlined with respect to the example above. where the priority function /in) > 0 is strictly decreasing in a. One can argue that an optimal solution may need to satisfy r _{\ } (t) 4 n _{\ }t ) ~ 1 for t ~ .1 ,2. Substitute mf?) ~ I ···· n (/) into coir), provides

The last two terms may be constant. If Oj < <¾ then /{<:?( ) - fiat) > 0 since f is decreasing, and hence the first term is increasing in /. This drives be as large as possible for small /, be., ~ 0. Hence the optimal charging profile remains

[0170] In various embodiments, to avoid numerical issues, in the cost function c y, the unit of ά > 0 should be chosen so that the cost coefficients (/ ···· «¾)/« _{/ } are not exceeding ;ly large or exceedingly small relative to each other, or within a reasonable range, e.g,, the range o f 10 ^{4 } and S O ^{4 } if possible. In this way, an AC can use a constant weight a > 0 to bring it )/ a within a reasonable range: mm While many of the· ACNs described above utilize cost functions that enable the utilization of charging rates that provide priority to particular EYs, ACNs in accordance with several embodiments of the invention utilize processes that accommodate a variety of other objectives.

REDUCING TEMPORAL FLUCTUATION

[0172] in many embodiments, it may be desirable to reduce rate fluctuations over a charging time period. This may be to reduce the wear and tear on a vehicle. Accordingly, it may be desirable to reduce rate fluctuations r;(/ } across / . In several embodiments, this can be achieved by penalizing the squared differences ffyff ) - r,(/ - 1 ) P of the rates by uti lizing a cost function ff ) that replaces the cost function in (23 ) with: or replaces ed r) in (24) with: where ?, 0} is the actual charging rate used in the last QP instance, and a > 0 and > 0 are weights.

[0173] In many embodiments, a QP at each time $ may solve for the entire charging profile ,,v 4- T. Only the first rates r*(s) of an optima! solution

U of the QP may be used for charging the EVs, An updated QP may then be solved at time s 4 1 and the cycle repeats it may be beneficial to that the first rates r*(s) do not fluctuate widely across different QP instances that are solved at different times. The cost function (24) attempts to reduce the temporal fluctuation of the solution to a single QP instance. Even though this is different from reducin fluctuations across different. QP solutions, it may achieve the desired effect

[0174] An alternative to the penalty term £,-.<(*)(/) ryff - .1 ) H is £,o ffvff }) ^{2 }. Both may reduce temporal fluctuation in r _{<, } The penalty term ff _{/ }ff ) ) ^{2 } has the advantage of making the QP cost function tq (r) strictly convex (and hence has a unique optimal solution) it may work well with the equality constraint on energy: but not with the inequality constraint:

Hence many embodiments .may use , (/yi/ } ···· r _{t }{t ···· 3 } r. Although the above describes a variety of mechanisms for modifying a cost function to reduce temporal fluctuations, any of a variety of mechanisms may be utilized as appropriate to the requirements of specific applications in accor dance with various embodiments of the invention. Certain other features that may be implemented using cost functions are described below.

ENFORCING MINIMUM CHARGING R TES

[0175] In general, it may be undesirable to set a charging rate to zero before an EV has finished charging because this may cause the mechanical contact in the charger to open. A charging profile r,(t) with many zero and nonzero rates before the EV finishes charging may incur excessive wear and tear. To avoid this, several embodiments may set charging rates h(? } > r, for some o greater than zero, as long as the energy demand ei > 0. For example, r, ···· 8A for constant charge (CC) chargers or tj 6A in J1772 standard. Note that, if the unit of charging rate r,{t) is in W, then 1 _{} } -- 6r· W at voltage level v V. To simplify exposition, certain embodiments often assume r;{t) is in A.

|0176] Since each QP may compute for every EV / an entire charging profile r* { r. (/ ) , / 1 , 71 } but only the first rate r ( .l ) may be executed, if may be conservative to impose the constraint dmt all future rates r,·!/ , and also increases the chance that the QP becomes infeasible. Several embodiments impose only /y( l ) > /y and pern-process the future rates, as described in Algorithm 2 below. 101771 A comparison of Algorithm 2 illustrated in Fig. 14 with the basic QP in ( 1 ) is described below.

1. The constraint ( lb) on the charging rates in the basic QP is separated into two constraints (14A-2.K14A-3) in Step 1 of Algorithm 2 illustrated in Fig. 14, the first on rates r( l) :~ (/?{ ! },/€ [ 1 , /?! } at time / 1 that enforces the minimum rate n and the second on future rates r(t }j > 1. As explained above, many embodiments do not impose the minimum rate on future rates.

2. Step 2 of Algorithm 2 post-processes the QP solution r* by choosing a rate r ^{* }* that is as close to r ^{* } as possible, subject to the constraints that ail rates r ^{* }* (/), t ~ 1 , 2, . . . , 7 ^{' }, are equal or higher than their minimum rates and that capacity limits are not exceeded. Notably, the energy requirements are relaxed, It means that the energy delivered by the charging profile r* ^{* } may not be as close to e< as possible or may even exceed ry, but one expects this effect to be negligible.

IMPLEMENTATION

{0178} Tiie problem ( 14B > is separable in t, and hence the subproblem for each / can lie solved separately in parallel:

Even if these (much smaller) /-subprobfems are solved sequentially, the computation seems to be faster than solving the overall problem ( 14B). Moreover, solving these r-subproblems separately helps handling infeasibility. Therefore, many embodiments solve them separately.

01791 Discrete rates. If the charging rate for EV i can only take values in a discrete set e.g.. CC chargers can charge at rates in A'. ( 64 , 12,4. 184 , 24,4, 30,4, 344 }. The final charging rates can be obtained from the solution of Algorithm 2 by rounding down to the nearest value in R,:

Simply rounding down may underutilize capacity if the gap between the discrete values is large compared with the capacity in that case, a more sophisticated algorithm can he used to map >· ^{* }*(/ } to a value in A',

HANDLING INFKASIBILITY

{0180J Certain embodiments may assume throughout. that 0 < r,· < A, ¾ < cy, and 0 < ¾ < e·, so that (I4A-2) (14A 4) can always he satisfied individually. Furthermore, assume f;(<A— 1 ) > <p, i.e., if each EV i is charged at its peak rate /y over its entire charging horizon, then its minimum energy < _{?, }· demand is always met.

{0181] Described now is how an ACN can handle the infeasibi!ity of QP ( 14A) in Algorithm 2 illustrated in Fig. 14, and then the infeasibility of the post-processing problem ( 14B _{) }. informally, the QP { 14 A } in Algorithm 2 can be infeasible when the minimum energy demand e,· in ( 14A-4) cannot be met for some EV i _{f } or when some capacity constraint in ( 14A-5) is violated at sonic resource /. Certain embodiments deal with this by setting in (!4A-2)and ( 14A-4) respectively. The constraints ( l 4A-2)-( l4A-5) can then be satisfied jointly (e.g,, setting >y(r)— 0 for all id). Note that this may not guarantee minimum energy delivered for this QP instance. In certain embodiments, this cart be refined by allowing some e, to remain nonzero, e.g., to guarantee all EVs 15 miles of charge or 15% of their requested energy. It also may not guarantee ATI ) > r _{( }- but this is handled in the post-processing step.

[0182] The post-processing problem can be solved by solving each /-subproblem separately. Recall that it is assumed .½ £ (0. 1 , ! } in general if ½ > 0, then a way a r-subprohiem can be infeasible is when even setting ryfr ) & for all / will violate (Ί4A-5) at some resource /, i.e,, where \S;\ is the number of BVs sharing resource /. In this case, some EV i may be assigned a rate lower than r, in order not to violate the capacity constraint. This means that some EVs may be charged at rate r, and others at zero. Certain embodiments may prioritize based on laxity by calling Algorithm 3, as illustrated in Pie. 15.

IΌ183] Even though Algorithm 3 illustrated in Fig, 15 is motivated by the case where Am— 0 or 1. it can ensure that the capacity constraints may be satished even when An can take value— 1. indeed, providing

It may be simple, but conservative since *Q can be chosen based on max _{/ } r, and independently of the group of EVs that share a constraint. It. also does not exploit the fact that t may exceed t; in which case — 0. In certain embodiments, if this scenario is only rarely encountered, a simple algorithm should be sufficient. Otherwise, the algorithm can be tightened to be «tore efficient.

101841 For some EVs, once their rates are set to zero, they may not recover (i.e., they ay not resume charging until they are unplugged and plugged in again). Since this may negatively impact user experience, many embodiments may avoid the situation where a resource constraint might be violated even if all EVs are charged at their minimum rates whena deployment is sized and if it participates in a demand response program

10185J In summary, certain embodiments may modify Algorithm 2 illustrated in Fig. 14 above by adding a mechanism to handle infeasibiiity and this is given as Algorithm 4 illustrated in Fig, 16.

PROVIDING GUARANTEES (SUPPORTING GROUP 1 EVS) 101861 In many embodiments. EVs may be groupe into different groups that provide different charging rates based on a variety of factors. For example a driver may pay to have a higher charging rale (in order to reduce ihe lime needed to charge their vehicle). Other factors may include reducin _{, } the charging rate for a vehicle that is expected to be parked for a longer duration of time and thus provide more charge to other vehicles that are expected to be parked for a shorter time period. In several embodiments, the groups may include a prioritized group (e.g. a Group 1 ), where the defining feature of the prioritized group is that EV i is guaranteed its requested energy e (up to a measurement error margin) under normal conditions. This is why Group 1 ma have strict priority over Group 2 and/or lower priority groups and why it is expected most EVs may be in Group 2 that guarantees only minimum energy > 0, but not requested energy <? _{/, }

[0187 j Algorithms that provide guarantees for certain EVs can be characterized in that they: i . Enforce equality energy constraint for each EV / in Group I : instead of the inequality constraints for Group 2,

2. In certain embodiments, in order to avoid infeasibility, perform admission control before a new EV / can be admitted to Group 1 . The purposes of the admission control may be to

(a) the requested energy <?/ of all exi ting Group I EVs j,

ib) the minimum energy £,· > 0 of all existing Group 2 EVs j.

Described below are various; charging algorithms that can be utilized by ACNs to determine charging rates for Group i EVs in accordance with certain embodiments, and then various processes that ACNs can utilize for admitting EVs into Group L

CHARGING ALGORITHM [0188] A charging algorithm for Group 1 (only ) may be QP followed by post- processing to enforce minimum rates ry. A charging algorithm for Group 1 in accordance with various embodiments of the invention is given in Algorithm 5 illustrated in Fig, 17. It can be the same as Algorithm 2 for Group 2 illustrated in fug. 14, except that the inequality energy constraints (I4A- 4} may replaced by equality energy constraints

[0189] An important assumption for the Group I algorithm is that, throughout the ACN execution, ail existing Group 1 and Group 2 EV requests can lie guaranteed, as long as resource capac ities i (/} (and other ACN conditions) remain unchanged, in that case, both the QP (17A) and the past-processing problem ( 17B) may always be feasible. Furthermore, conditions can change in unexpected ways, and therefore certain embodiments add mechanisms to Algorithm 5 illustrated in Fig. 17 to handle i «feasibility (c.f. Algorithm 4).

ADMISS ION CONT ROL

[0190] Certain embodiments may control admission of an BV to an ACN based on the current profile o f t he ACN, i ncluding i ts e xisting EYs charging demands, the ACNN c apacity, a mong a variety of other factors. Assuming the laxity of the new PV 0 request { ,<·/ ) has already been checked by Algorithm 1 illustrate in Fig. 13, before the admission to Group ! is determined. The admission control may be implemented by Algorithm 6 illustrated in Fig. I S and Fig. 19. Certain embodiments may handle differently new EV requests to Group 1 that have been rejected by the admission control. Certain embodiments may repeatedly ask the driver to re-enter their parameter (<¾,<:¾) until it is admitted to Group 1 , or a maximum number of tries are reached. In certain embodiments, the ACN may suggest a departure time of so that (eo.c/ _{j } can be admitted, or suggest an energy level so that uf, d(>) can be admitted. Note an alternative to admission control is not to have drivers specify both <?,- and </ , but only one of them and ACN may decide the other parameter so that the new request can always be admitted to Group 1. In certain embodiments, an ACN may assign the request (<¾,«:¼>) to Group 2 with appropriate parameter ¾ > 0, In several embodiments, the ACN may assign it to Group 3 and not serve it until subgroups are defined and implemented. As can readily he appreciated, any of a variety of admission control processes can be utilized by an ACN as appropriate to the requirements of a given application in accordance with various embodiments of the invention.

HANDLING !NFEASIBIUTY FOR QP IMPLEMENTATIONS

101915 in many embodiments, ACNs can utilize a variety of processes to determine charging rates for EVs when requested charging parameters provided with respect to the EVs would other wise yield an infeasible QP. In several of these embodiments the basic features of a QP implementation may include the following:

1. Each subscriber (EV driver) may input energy demands (e.g., must-have miles) and departure time, among others.

2. During registration, each subscriber may also set a willingness to pay for nice-to-have miles.

After the rnust-have miles have been delivered to the subscriber, energy may he delivered to the EV when the current energy price is lower or equal to the willingness to pay.

3. A price ($/kWh) may be computed dynamically in real time (e.g,, every minute). A subscriber may pay the energy delivered multiplied by the price, for both rnust-have miles and nice-to-have miles, regardless of her wilhngness-to-pay.

4. Charging rates may be computed using the QP with a cost function that is the weighted stmt of three components; where wy > 0 are nonnegative weights and

* The first component cq (r) is a measure of performance, e.g., to charge as fast as possi ble, to minimize expected energy payment, or to track forecast sealed solar generation * The second component is to reduce temporal fluctuation in charging rates.

* ^{' }The third component is to snake the cost function strictly convex so that the optimal charging rates {i) are unique. EVs are divided into two groups:

* Group 1 : EVs whose must-have miles have not been fully delivered,

* Group 2; EVs whose must-have miles have been fully delivered, and may receive niee- to-have miles when the dynamic prices do not exceed their willingness to pay.

Group I EVs have strict priority over Group 2 EVs At time i when QP is computed: fa) All Group I EVs will be included in QP computation at time ?, regardless of the current price pit) .

(b) A Group 2 EV will be included in QP computation if and only if its wiiiingness-to- pay for niee-to-have miles is higher or equal to the current price p{t) (and the QP is feasible; see below).

(c) All EVs that are scheduled for charging at time i wil l pay at price p{t), regardless of their willingness to pay for Group 2 EVs. The QP computation can he implemented as follows:

(a) Solve QP with:

* equality energy constraints for ail Group I EVs.

• inequality energy constraints for all Group 2 EVs that are included In scheduling.

(b) If the QP is feasible, then it will produce charging rates for the next control perio (e.g., next minute) for all Group 1 EVs and those Group 2 EVs that are included in scheduling for that control period. (c) Handling mfeasibility. If the QP is infeasible, then solve QP again with only Group

1 EVs an inequality energy constraints. Notify Group I EVs whose energy demands are forecast not to be met fully.

HANDLING INFEASIBILITY: DEADLINE EXTENSION

|0192] When QP is infeasible, certain embodiments may only schedule Group i EVs, Instead of changing all their energy constraints from equalities to inequalities many embodiments can extend their deadl ines (charging durations) to ma e QP feasible, and notify the drivers that the may have to stay longer to receive fully their must-have miles, or receive less energy at their original departure times.

{0193] Described now are two processes utilized in many ACNs that can be used to determine new deadlines for EV charging in the face of infeaslbiiity of originally requested deadlines. The first algorithm scales up all EV durations by a common factor it does not involve discrete variables but requires binary search over QF's. The second algorithm adds dwell times to the EVs in a way that minimizes weighted sum of the additional times, The weights allow the implementation of different priorities among these EVs. It involves binary variables (which may be computation- ally expensive) and thus certain embodiments may provide an algorithm that relaxes the integer constraints.

SCALE ALL DURATIONS BY A COMMON FACTOR

[9194] Recall that the departure time d _{t }- is incorporated into the maximum charging rate 7y{/), by setting r- i t ) ------ 0 for t > d . Then the constraints can ensure that the QP will not charge EVs i after their specified departure times. Here, / ··· 1 denotes the current time when the QP is called, and d- - I denotes the duration over which EV ϊ is available bar charging from time t ~ I ("current time) to time t— d,· - 1. |0195] When QP is infeasible, certain embodiments may extend the durations of all EYs by a common scaling factor y > 1 so that QP becomes feasible with equality energy constraints. The following algorithm attempts to determine a minimum such y.

[0196] Let d (h _{. }·L7) denote the vector of EV durations. Let yd denote the vector of durations scaled by y > 1 but tipper bounded by / . re..

(Consider the QP that include only Group 1 EYs with equality energy constraints, and suppose it. takes the form;

Here, .'¾("/</) denotes the feasible set for the charging rates r, as a function of the scaled durations d. Therefore the QP is a function of yd, and this dependence is emphasized by the notation QP( d). The algorithm will do a binary search (or other appropriate search) on y until QPfy ) becomes feasible.

[0197] Assumption: The QP when the durations of all TVs are T is feasible with equality e nergy eo ns trai ms .

{01981 The algorithm is described in Algorithm 6 illustrated in Figs 17 and 18 and described below.

ADDITION OF HETEROGENOUS DWELL TIMES TO COST FUNCTION

ff)!99] In certain embodiments, the (original) depasture times and the maximum charging rate conditions are combined into a single set of constraints (27). To specify the second algorithm to extend durations, (27) may be replaced by two separate set of constraints: where f- ) are the upper bounds on charging rates of EVs i at time t, e.g. due to EVSB or EV limits. Note that h(h are given constants. Suppose the QP takes the form:

Here r iR denote all QP constraints tha do not involve the durations , including the equality energy constraints. The deadline constraints are given by (29).

10200] To add to the durations d _{t }, we introduce binary’ variables xi(i), one for each time t - d· _{. } _ _ T, and each EV i. This adds a total of ); / ^{' } - -f 1) binary variables. If t > is an additional dwell time (i.e., before the new deadline), then Xj(t) ~~ I ; otherwise. ay(r) 0. Hence if Xjit )— 0 (i.e., t is after the new deadline), then x-{s) ~ 0 for ail subsequent times $ > /. The constraints (29) are replaced by ( note by our convention, .q( G)— () even though this is not enforced in ( 1 ). )

10201 } Condi if on (3 l b) says that if / > d- is an additional dwell time (.¾·,(/) ~ 1 ) then the rate /-,(/ ) is upper bounded by /)(/ }; otherwise (.v;(r) 0). nit) is set to zero. Condition (3.1c) says that if / > d, is not an additional dwell time (Le., ··· 0) then ah subsequent times will not be additional dwell times. This ensures that all rates are set to zero after the new deadline. Note that it is possible to choose y (Q _{· } · 1 and n(t) ~~ 0, i.e„ time t is before the new deadline even though QP chooses not to charge EV / at time t.

ALGORITHM

[0202] When the QP (30) with original deadlines d is infeasible, certain embodiments may solve the following problem (with only Group I EV and equality energy constraints):

Note that the optimization is over both the charging rates r and the binary variables

o(i) represents the additional dwell time added to EV rs duration. The optimization (32) therefore minimizes the weighted sum of additional dwell times subject to the same set of constraints as the original QP as well as constraints (31 ) discussed above. An optimal solution may provide both feasible charging rates r and the new durations ( deadlines) f ^{cW } for Group 1 EV :

[0203) A bigger weight w; for EV r means its deadline will be extended by a smaller amount. Hence the weights H) can be chosen to implement priorities, e.g., an EV that arrived earlier or has a greater laxity can be assigned a bigger ny and hence will incur a smaller deadline extension.

[0204] Since both the cost function and the constraints in (32) are linear, the optimization is a mixed integer linear program (Mil.PT Certain embodiments, can relax the binary constraints (3 Id) to jO. 1 {-interval constraints: instead of (32), certain embodiments solve the following linear program

Given an optimal solution (>, ) of (35), certain embodiments may need to discretize .vg/) to either

0 or 1 in order lo calculate the new durations <f ^{iew } using (33).

102051 In summary, the ACN ^{! } can solve the U1P .R (32 } if if is computationally manageable. Otherwise, the ACN can use Algorithm 7 illustrated in fag. 20 and described below.

[0200] Even though the charging rates produced by (35) in Step ! of Algorithm 7 are feasible for the new durations certain embodiments may re-compute QP( ^{v }) in Step 4 to obtain a potentially better set of charging rates.

ONLINE LINEAR PROGRAM

10207 j Described now is an LP that is to be constructed and solved at computation time 0 to compute charging rates it) over time window / _{< }£· [ 1 , 71.

1 , The decision variable is the vector of charging rates for ail active EVs / over the time window (1 , G|. This is an «T -dimensional vector.

2. The simplest linear cost function is

Here, (c, - (/)) is the unmet energy demand of EV r by the end of the time window j 1 , T\. Hence this cost function minimizes the total unmet demand of all active EVs.

Note that, in minimizing this cost, it is possible that the optimal charging rates r deliver more energy than the EVs have requested, i.e., (r optimality if this is undesirable, certain embodiments can include the constraints:

Y nin

^{' }This will ensure that no more than the requested energy is deliverer

Since ,·*? _{( }· is a constant, it can be dropped from the cost function, i.e . we can equivalently minimise the cost or, maximizing the total energy delivered to all active BVs. There may be several constraints on the charging rates r

* Rates do not exceed their ax rates:

By setting ?,·(/ } 0 for t > constraint (36.) also enforces it) 0 for t > i.e„ do not charge alter BV fs declared departure time

* Capacity constraints of the form ry(/) < P _{j }{t) for all resources

In summary, the simplest version of LP that an AC can solve at time 0 is:

1112 81 This is illustrated in Figure 21. which illustrates a structure of an online LP for EV charging.

1021)9] Certain modifications to two aspects from the processes described above have been made. First, weights uy are included in the objective (37a), This can be used to implement priority when it is not possible to satisfy every EV’s energy demand by their deadlines; online LP favors EVs with larger weights uy. When it is possible to satisfy every EV’s energy demand, the weights ny have no effect a all EVs i,

10210] Second, a lower bound tc on the energy delivered to EV ί is included in (37c). This is a service-level guarantee (minimum energy delivered). If e,— 0 for all EV i, then there is no guarantee and the LP (37) will also be feasible, e.g., 0 for all /, ? is a feasible solution.

Otherwise, It is possible that (37) is infeasible, but when online LP is feasible in ever computation period, then every EV i will get a minimum amount ¾ of energy by the end of time 7 ^{" }. In many embodiments, the st function utilized by an ACN can vary based upon a variety of factors including (but not limited to) time of day, cost of energy, actual and/or anticipated available capacity, number of EVs, and charging parameters requested with respect to specific EVs.

PRODUCT FEATURES AND LP EXTENSIONS

1021 I f Tire online LP framework of many embodiments is very versatile. Different product features can be implemented within this framework, e.g., by choosing different cost functions, different minimum energy guarantees ¾, different constraints /)(/), and/or different max charging rates Tjit ). Described belo are some of these features in accordance with various embodiments of the invention.

I NFRASTRUCTUR E PROTECTION

[0212] Consider the scenario:

• The infrastructure capacity of the site is ta constant) P.

• an ACN shares this capacity with some background (non-E-V) load.

• an ACN measures in real time the background load at the computation period r— 0; denote it by L().

[0213] Many embodiments ensure that the total background load plus the total EV load do not exceed P.

[0214] To construct the frannvework in the computation period 0, the ain task is to forecast the background load over the (future) time window ( I , T\ . Let Lit ) denote the forecast background load for In the simplest case, certain embodiments assume future background load remains the same as the current background load:

Lit ) : ····- /< . / T i Jj if historical data on background load is available, then set the forecast Lit) to be its historical value (e.g., the average value over the last week, perhaps depending on day of the week).

[021.5] Using the forecast Lit), certain embodiments can prevent overloading the infrastructure by adding to (or replacing some of) the constraints (37d) (by) the following constraint;

This ensures that the EV load plus the forecast background load do not exceed the infrastructure capacity. This is illustrated in Fig. 21 , which illustrates a structure of online LP for infrastructure protection and load shifting in accordance with various embodiments of the invention. Although Fig 21 illustrates a particular structure of an online LP for EV shifting, any of a variety of structures may be utilized as appropriate to the requirements of specific applications i accordance with various embodiments of the invention.

DEMAND RESPONSE: LOAD SHIFTING/TRACKING

10216} Consider a same scenario as described above where an AC shares the same capacity with some background load, and Uj ) denotes the forecast background load in the future based on historical data and/or real-time measurement of the current background load,

LOAD SHIFTING:

10217} Cap the total site load to D(l) over t€ (b d?]· For example, an DR event is to maintain the total site load below or at 1 MW from 1pm to 4pm today. This feature can be implemented in online LP, as illustrated in Fig. 22, by adding to (or replacing some of) the constraints (37d) ^{( }by) the following constraint:

This is illustrated in Fig. 22. which illustrates the structure of an online quadratic program for load tracking. Although Fig. 22 illustrates a particular structure of an online QP for load tracking any of a variety of structures may be utilized as appropriate to the requirements of specific applications in accordance with various embodiments of the invention

LOAD TRACKING;

10218} Have the total site load track a given profile />{?) over / e i i - d · This feature can be implemented by replacing the cost (37a) by the following cost (or adding to the cost):

[(12191 In this case, the constraint (38) should be included in (37d) to ensure the background load and the EV load does not exceed the infrastructure capacity. Note that in this ease, the objective is no longer linear and hence the problem is not an LP It is a simple (convex) quadratic progrant and can still be efficiently solved. This is illustrated in Fig. 23, which illustrates a structure of online quadratic program for load tracking in accordance with various embodiments of the invention, Although Fig. 23 illustrates a particular structure of an online quadratic program for load tracking, any of a variety of structures may be utilized as appropriate to the requirements of specific applications in accordance with various embodiments of the invention

PRIORITY CHARGING

[0220] Priority charging can be important, e.g., to coordinate between DCFC and L2 chargers, or if drivers pay different prices. Priority among the EVs can be implemented by the appropriate choice of parameter values in online LP (37); a higher-priority EV l can be assigned

1 , a larger weight ny in (37a), and/or

2 a larger min energy ¾ in (37d). and/or

3 a larger max charging rate ?, (/) in (37c),

[(12211 When the EVs have enough laxity and the infrastructure has enough capacity, these parameters have no effect on the final energy delivered, i,e., AON will meet the energy demand of every EV before their departure in that case. Priority makes a difference only when there is insufficient capacity or insufficient laxity.

[02221 Priority can be used to deal with a driver who does not provide input. In this case, the EV is assigned the lowest priority.

DEMAND CHARGE MITIGATION

10223] Certain embodiments include a demand charge mitigation feature assuming the peak demand used for demand charge determination is the max energy consumed in each time period t. More likely, the peak demand used for demand charge determination is the max energy consumed in (say) 15 mins which can span multiple time periods i. The following design can be easily extended to the more general case.

[0224] Let 0{ 0) be the peak demand (EV load + background load) observed from the beginning of the current Fulling cycle (e.g., from beginning of the mouth.) to the current time t ~ 0. To minimize demand charge, certain embodiments can add to (or replace some of) the constraints (37d) (by ) the following constraint: where L{t) are the forecast background loads over t e ( 1 , ^{' }F J , This guarantees that the total charging rates do not exceed the previous peak demand D{ 0} (nor the infrastructure capacity P). This is illustrated in Fag. 24, which illustrates a structure of an online LP with demand charge mitigation in accordance with various embodiments of the invention. Although Fig. 24 illustrates a particular structure of an online LP with demand charge mitigation, any of a variety of structures may be utilized as appropriate to the requirements of specific applications in accordance with various embodiments of the invention,

|0225] Just like the basic online LP (37), the LP with demand charge mitigation (39) will always be feasible if the minimum energy guarantees . 0 for all EV L Otherwise, the LP may be infeasible, in that case, there are two possible policies. In the first case, certain embodiments set all e, ~ 0 so that LP becomes always feasible, but may forgo minimum energy guarantees. In the second case certain embodiments remove the limit due to the current peak load D(r) and require only (note that this reverts to the formulation (38) for infrastructure protection and therefore may still be infeasible if g, > 0, but the infeasibility is not due to demand charge mitigation):

{0226] After an optimal charging rate vector is computed, the rates r ^{*: }( l ) : · · (r ^{* }( l ) J £. i Ϊ , /¾] ) is used to set the pilot signals for time / ·· 1. It is possible the actual total load exceed the current peak load Dip) . Then the current peak load Dip) is replaced by a measured new (higher) value D( I ) for the online LP problem at time t ~ 1.

REAL-TIME PRICE ADAPTATION

102271 in many embodiments, there can be multiple EV charging policies based on real-time prices

PRICE CAP

[0228] Each user specifies a cap PA ) on the real-time electricity prices with the interpretation that if the real-time price pit) at computation period t is less or equal to the cap />, { ) of EV i, then EV i will be included in the charging rate computation. Otherwise, its rate is set to zero and excluded in the compulation. This is illustrated in Fig. 25, which illustrates real-time pricing with price cap in accordance with various embodiments of the invention. Although Fig. 25 illustrates real-time pricing with price cap, any of a variety of mechanisms for price capping may be utilized as appropriate to the requirements of specific applications in accordance with various embodiments of the invention

PAYMENT CAP

[0229] Each user i specifies a budget cap so that the total charing payment is no more than b;. Let /)(/ ) be the published or forecast real-time prices over time window

[0230] Certain embodiments may minimize the expected electricity payment in each computation period, and update the remaining budget after charging in each period. Let b,(ϋ) be the remaining budget of EV / at computation period p>(0). Specifically, we solve the following online LP:

After an optimal charging rate vector r* : i/ } <4 is computed, the rates used to set the pilot signals for time i ·-- 1 . This incurs an electricity cost pi 1 )ή ^{( } 1 ). assuming the />( i ) = p{ 1 ) is the true real-time price for time i ~ I. The remaining budget is updated to: and the cycle repeats for computation period / 1.

TIERED PRICES/PAYMENTS

[0231] In certain embodiments, each user may specify two energy levels (ej ,ej) and three caps on prices { p- (t).pf {i) ,pj{t i ) of on payments (b , b?, b·' ) , with the following interpretation: if energy delivered < e , then charge EV / with a price (demand) cap of pj (t) f† ); if e- < energy delivered < f _{< } then charge EV i with a price (demand) cap of (r) ίb~); if energy delivered > e , then charge EV i with a price (demand) cap of pf it) (b;').

JOINT EV/SOL-AR/STORAGE OPTIMIZATION 10232] In many embodiments an AC includes EV charging solar generation, and at least one battery*. In certain embodiments, the OP .framework may schedule EV charging and battery (onsite energy storage) operation to track solar generation. This can be formulated as a convex program. Consider again the problem at time t 0.

Let be the forecast background load over time window t G ( 1 , 71, based on real-time measurement /,(0) and historical data.

* Sit } be the forecast solar generation over time window t€ (1 , 71, based on real-time measurement 5(0) and historical data.

* b(i) be the SoC of the battery at time t . Use a simple battery model (without loss) with b{t -r 1 } - / (/ ) · ! «(/} where u(t) is the amount of energy in time period t drawn by the battery if ait ) is positive. and the amount of energy released by the battery if u(t ) is negative. 0 3 J In certain embodiments, the QP framework may choose charging rates r and battery operation u(t } so as to minimize the distance between the forecast solar: generation S{i ! and the total net load (EV + background + battery draw). This can be accomplished by solving the following convex optimization problem (quadratic program):

Y rj(t) < Pt(j ) t |U 7 j* for all / (4 I d) [0234] Here the new constraints (41e) represent a simple (lossless) battery model that describes how SoC evolves according to the control w(/ j and that the SoC must He between empty and full (B being the battery capacity of the onsite energy storage facility).

[0235] Accordingly, many embodiments provide a unified algorithmic framework to guide the design of a dean, flexible a n e vo!vabie a rchitecture t o i mplement v arious optimization- based product features. Fig. 26 illustrates building blocks of optimization-based product features in accordance with various embodiments of the invention. In particular, these may include;

* User account and management, e.g., for priority or payment purposes;

* Background load /,(/) real-time measurement and forecast from historical data;

* Interface for DR commands !>(/ ) from site operator or utility;

* User priority determination;

* Current peak load measurement /. ^{) }(0};

* Interface for real-time prices />(/) and forecast from historical data;

* Solar generation S(l } real-time measurement and forecast from historical data;

[0236] Although Fig 26 illustrates various budding blocks for optimization-based product fea tures, any of a variety of building blocks may be utilized as appropriate to the requirements of specific applications in accordance with various embodiments of the invention

[0237] Although the present invention has been described in certain specific aspects, many additional modifications and variations would be apparent to those skilled in the an. It is therefore to be understood that the present invention may be practiced otherwise than specifically described as appropriate to the requirements of a given application. Thus, embodiments of the present invention should bo considered in all respects as illustrative and not restrictive.

**Previous Patent:**SYSTEMS AND METHODS OF PILOT ASSIST FOR SUBSEA VEHICLES

**Next Patent: FIBER-BASED COMPOSITE WITH FRACTURE-INDUCED MECHANO-ELECTRICAL SENSITIVITY**