The Berlin Infrared Optical System satellite, which is scheduled for launch in 2016, will carry onboard a picosatellite and release it through a spring mechanism. After separation, it will perform proximity maneuvers in formation with the picosatellite solely based on optical navigation. Therefore, it is necessary to keep the distance of the two spacecraft within certain boundaries. This is especially challenging because the employed standard spring mechanism is designed to impart a separation velocity to the picosatellite. A maneuver strategy is developed in the framework of relative orbital elements. The goal is to prevent loss of formation while mitigating collision risk. The main design driver is the performance uncertainty of the release mechanism. The analyzed strategy consists of two maneuvers: the separation itself, and a drift-reduction maneuver of the Berlin Infrared Optical System satellite after 1.5 revolutions. The selected maneuver parameters are validated in a Monte Carlo simulation. It is demonstrated that both the risk of formation evaporation (separation of more than 50 km) as well as the eventuality of a residual drift toward the carrier are below 0.1%. In the latter case, formation safety is guaranteed by a passive safety achieved through a proper relative eccentricity/inclination vector separation. Notation a = semimajor axis of the carrier satellite, m e = eccentricity of the carrier satellite e i = delta-v error on the ith tangential component, m∕s f i = simulated performance factor of maneuver i i = inclination of the carrier satellite, deg k = odd natural number n = mean angular motion of the carrier satellite, rad∕s p = real number u = mean argument of latitude of the carrier satellite,°v = absolute value of velocity vector, m∕s ΔB = differential ballistic coefficient, m 2 ∕kg Δ• = finite variation of a quantity δα = nondimensional relative orbital elements set δe = nondimensional relative eccentricity vector δi = nondimensional relative inclination vector δλ = nondimensional relative longitude δv R , δv T , δv N = instantaneous velocity changes in local radial-tangential-normal frame, m∕s δ• = relative quantity ϵ xi , ϵ yi , ϵ zi = simulated errors of the carrier attitude control during maneuver i, rad θ = argument of latitude of the relative ascending node,°ξ = phase change of the relative eccentricity vector, rad ρ = atmospheric density, g∕km 3 φ = argument of latitude of the relative perigee,°χ = function of the difference between the arguments of latitude of the two maneuvers Ω = right ascension of ascending node of the carrier satellite, deg ω = argument of perigee of the carrier satellite, deg