Orbital flight of CubeSats at altitudes between 150 and 250 km has the potential to enable a new class of scientific, commercial, and defense-related missions. A study is presented to demonstrate the feasibility of extending the orbital lifetime of a CubeSat in a 210 km orbit. Propulsion consists of an electrospray thruster operating at a 2 W, 0.175 mN thrust, and an specific impulse (Isp) of 500 s. The mission consists of two phases. In phase 1, the CubeSat is deployed from a 414 km orbit and uses the thruster to deorbit to the target altitude of 210 km. In phase 2, the propulsion system is used to extend the mission lifetime until propellant is fully expended. A control algorithm based on maintaining a target orbital energy is presented that uses an extended Kalman filter to generate estimates of the orbital dynamic state, which are periodically updated by Global Positioning System measurements. For phase 1, the spacecraft requires 25.21 days to descend from 414 to 210 km, corresponding to a ΔV 96.25 m∕s and a propellant consumption of 77.8 g. Phase 2 lasts 57.83 days, corresponding to a ΔV 119.15 m∕s, during which the remaining 94.2 g of propellant are consumed. Nomenclature a = semimajor axis, km a d = acceleration due to atmospheric drag, km∕s 2 a ns = acceleration due to aspherical geopotential, km∕s 2 a T = control acceleration vector, km∕s 2 C = controller gain C d = drag coefficient E e = estimated energy per unit mass, km 2 ∕s 2 E SA = energy generated by solar arrays, J E t = target energy per unit mass, km 2 ∕s 2 F = Jacobian of nonlinear model G = state matrix of process noise H k = Jacobian of measurement model i = orbital inclination, deg K k = Kalman gain mt = spacecraft mass, kg m p;0 = initial propellant mass, kg _ mt = mass flow rate, mg∕s P EP = power available for electric propulsion system, W P k = covariance matrix of residual noise P S∕C = power available for spacecraft bus and payload, W Q = covariance of process noise R E = equatorial radius of Earth, km r a = actual position vector, km r e = estimated position vector, km r t = target position vector, km S = reference area, 0.01 m 2 T d = time spent in sunlight during one orbital period, s T e = time spent in eclipse during one orbital period, s T T = commanded thrust in the alongtrack direction, mN v a = actual velocity of the spacecraft, km∕s v e = estimated velocity of the spacecraft, km∕s v rel = velocity of the spacecraft relative to the atmosphere, km∕s v t = target velocity of the spacecraft, km∕s W = work done by atmospheric drag, J X d = power transfer efficiency: solar array to loads X e = power transfer efficiency: solar array to batteries x a = actual state, 6 × 1 x e = estimated state, 6 × 1 y k = Global Positioning System state measurement, 6 × 1 ΔE = energy per unit mass error, km 2 ∕s 2 Δt = integration time step, 0.1 s ΔV = change in velocity, m∕s δ = declination (latitude), rad δr = error between estimated and actual positions, m δ _ r = error between estimated and actual velocities, cm∕s ε = eccentricity ε GPS;k = Glob...