Adaptive Petal Reflector: In-Lab Software Configurable Optical Testing System Metrology and Modal Wavefront Reconstruction

This paper addresses two aspects of the metrology of spherical, petal polymer reflectors which are part of an effort by the European Space Agency (ESA) to develop actively controlled foldable reflectors, enabling larger apertures on CubeSats and small satellites. The first problem is that of measuring the surface figure error of the spherical reflector alone during the development phase, and to assess the quality before assembling the telescope (large stroke, low accuracy). The SCOTS (Software Configurable Optical Testing System) appears to provide a fast and satisfactory solution to this problem. The second problem is the wavefront error reconstruction when the petal reflector is mounted on the telescope, because parts of the petals are obscured by the secondary mirror, in such a way that the petals appear completely disconnected, making the gradient-based metrology impossible. Using the fact that the petals have common mechanical boundary conditions at the central support ring, the problem is solved by using a set of orthogonal modes satisfying the same boundary conditions. The vibration modes are used for this purpose; the modal amplitudes are reconstructed from slope data outside the obstruction, allowing for wavefront error reconstruction over the entire surface.


Introduction
This study is motivated by the ongoing, ESA-funded development of foldable, polymer, actively controlled petal reflectors for CubeSat applications, as seen in Figure 1 [1][2][3].As the reflector is foldable, larger apertures are enabled on small satellite platforms.Upon unfolding, active control is attained by a thin layer of piezoelectric PVDF-TrFE between aluminum electrodes, allowing for an accurate surface to be maintained.Two problems are addressed in the development and operation of the active reflector, namely:

•
The developer of the reflector must have access to a metrology system capable of a rapid evaluation of the manufacturing process, and to verify, without the availability of the whole telescope, that the surface figure aberration will remain within the envelope that can be corrected subsequently by the control system.• When the primary reflector is mounted on the telescope, part of the petals are in the shade of the secondary mirror and they appear as segments (i.e., independent of each other).However, unlike segmented mirrors which require edge sensors, all the petals have a common mechanical boundary condition at the inner ring of the reflector.This paper investigates the use of this property to implement a wavefront reconstruction based on the measurement of the x-y slopes at a set of grid points (Shack-Hartmann-type sensor).
Other concepts for deployable space telescopes for CubeSats include rigid segmented (hinged) reflectors [4] and photon sieve designs [5].In contrast, the folding of the reflector described herein is attained by its own flexibility.Further ESA R&D activities in optical design and control of space telescopes are reviewed in [6,7].
Figure 1.Illustration of the system design of the active petal reflector in folded and deployed states.The reflector (red) is coated with a piezoelectric material, allowing for active control in closed loop given a metrology system.Since the secondary will obstruct the primary mirror, the recovered sensor data will be segmented, and special wavefront reconstruction is needed to recover the error accurately.

In-Lab Metrology
For the metrology in a lab environment, a technique is needed which is able to provide a few-micron accuracy, sub-micron repeatable test, which at the same time has sufficient dynamic range to capture aberrations larger than 100 µm Peak-to-Valley (PV).While interferometric techniques would provide accuracy much better than required, the dynamic range would be challenging.Conversely, techniques such as photogrammetry would have sufficient range, but the accuracy would be challenging.The Software Configurable Optical Testing System (SCOTS) was developed specifically to overcome similar challenges in testing solar concentrators [8].The SCOTS test has since been progressively refined, and some implementations are successfully measuring sub-nm errors on X-ray mirrors [9].The test is based on deflectometry and requires three components: a reflecting surface under test, a screen that projects a series of fringe patterns with sinusoidal intensity, and a camera recording the distorted fringes as reflected by the mirror.A main difficulty in obtaining high-precision results is the placement of the camera, screen, and mirror components, and high accuracy requires measurements of the relative positions with, e.g., laser trackers [10].For the purposes of our work, such strict geometric control will not be necessary due to the relatively relaxed accuracy requirements; the purpose is to test that the manufacturing error is within the actuator stroke.
One challenge with the SCOTS test is that results are obtained easily, but corroborating them with reality is more difficult.One approach is to use a "closed-loop" system in which the resulting geometry from the measurement is introduced into a ray-tracing software, allowing the user to see if the experimentally recorded raw image data correspond to simulated image data.This approach will be used in this paper, and a similar approach is reported in [11].Unique to this research is the measurement and reconstruction of petal geometries from measured slope data, where specialized algorithms for phase unwrapping and shape integration are needed to handle discontinuities in the reflector.The methods developed in this article would also be extensible to other topologies and have been made available as open source [12].

Vibration Modal Wavefront Reconstruction
Once mounted on the telescope, the adaptive petal reflector is intended to be controlled by an array of piezoelectric actuators.This assumes the availability of a wavefront sensor.Several approaches exist to wavefront sensing and reconstruction, for instance zonal [13] or modal (Zernike, Legendre or Chebyshev [14,15]) reconstruction techniques for Shack-Hartmann sensors, neural network-based reconstruction for pyramid wavefront sensors [16] or even wavefront sensorless adaptive optics [17].In this work, we focus on gradient data obtained from, e.g., a Shack-Hartmann sensor, and the reconstruction of the wavefront from this data.
Due to the petal geometry, the central part of the reflector is obscured by the shade of the secondary mirror (Figure 1), in such a way that the petals appear disconnected (as segments), making the gradient-based metrology impossible.However, since the petals have common mechanical boundary conditions at the inner ring of the reflector, this feature can be exploited to implement a modal reconstruction method based on a complete set of orthogonal modes satisfying the same boundary conditions.The normal vibration modes [18,19] appear as a natural choice.Having reconstructed the modal amplitudes from the slope data in the accessible part of the reflector, an approximation of the wavefront aberration over the entire reflector is reconstructed from the modal expansion.It is worth noting that the only strict requirement on the vibration modes is that they satisfy the geometric boundary conditions at the inner ring of the reflector.Two sets of modes will be compared below: the modes of the spherical petal shell constituting the reflector and the modes of a uniform plate with the same profile.

SCOTS
In a SCOTS setup, successive images of sinusoidal fringes (phase screens) are displayed on a monitor.The images of the phase screens are then reflected by the mirror and captured by the camera.This allows the shape of the mirror to be reconstructed from the acquired camera images, as the slope at a given location on the mirror can be recovered from the distorted images.The points on the mirror in which the slope is reconstructed are referred to as the Mirror Pixels.
To understand the test, consider a Mirror Pixel with location coordinates (x m , y m , z m ) and a pinhole camera with location coordinates (x c , y c , z c ). Tracing (in reverse) the ray of the camera through its reflection with the mirror leads to a single screen point (x s , y s , z s ).The goal of SCOTS is the accurate identification of the two geometric vectors for each mirror pixel: the Mirror to Camera vector ( m2c) and the Mirror to Screen vector ( m2s).An illustration of the SCOTS test can be seen in Figure 2. Given these vectors, the surface vector normal to the mirror, N, is obtained as the sum of the normalized vectors: From the surface normal vector N, the slope of the mirror (w x , w y ) can then be calculated: The surface of the mirror can subsequently be reconstructed from the slope data.The challenge lies in the accurate identification of these geometric parameters for each mirror pixel.The position of the camera aperture can be measured directly.Furthermore, the position of any pixel on the screen can be determined by choosing an origin and measuring the screen pixel scale.Likewise, under the assumption that the camera creates a flat image, the location of Mirror Pixels may be found by selecting an origin and measuring its location along with the image pixel scale.
In the simplest possible implementation, a single Screen Pixel may be lit at a time.In this case, all nine geometric parameters of Figure 2 are known by simple measurement, and the slope can be calculated as shown.However, for a screen with millions of pixels, this is tedious.Alternatively, the entire screen can be used to project a vertical and horizontal sinusoidal fringe map, such that each point on the screen has two associated unique phase values.These unique phase values must then be recovered from the recorded image, allowing every Mirror Pixel to have an associated screen location.The process of recovering the phase happens in three steps: • Zero-phase recovery : A single pixel is lit up in the (0, 0)-phase location of the screen and its position in the recorded image is saved (Figure 2, red dot).• Determination of phase: Successive horizontal and vertical sinusoidal fringes with varying phase offsets are displayed.For each individual Mirror Pixel, a sine function is fit on the data from the images to determine the phase value of this pixel from −π to π (Figure 2, blue dot).

•
Phase unwrapping: as the previous step happened on a per-pixel basis, there is no global information about the phase.The global phase is recovered by a numerical phase unwrapping technique and offset according to the previously recovered zerophase location.
With the phase values known in both the recorded image and the projected image, all nine geometric parameters can be inferred, and the slope may be calculated.The process of finding the geometric parameters is formalized in Appendix A. Finally, the slopes may be integrated to find the shape.This step usually happens by the integration of the difference between the recorded and expected slope, rather than directly from the slope.

Raytracing
Before using the SCOTS test to measure a new topology, its important to verify the test.One approach to this is to simulate images using raytracing in a known virtual environment.To do this, a 3D model of the petal with a known aberration may be loaded into a raytracing software.The camera can be modeled as a pinhole, and the screen is an emissive source of varying intensity.In the raytracing software, the camera launches rays "in reverse" toward the 3D model; they are then reflected and intersect the emissive source (screen).The value of the emissive source is sampled to generate the image.The software used is Raysect [20] in Python, which can be found from the Matlab environment.An example aberration of the petal reflector is illustrated in Figure 3 (left) in which a 100 µm tip displacement is applied.The resulting ray-traced image can be seen in Figure 3 (right).The tip displacement is generated by applying a fixed displacement in an FEM software, ensuring the displacement is representative of expected deformations on the reflector.Notable parameters are as follows: radius of curvature-2.5 m; diameter-20 cm; camera location (x c , y c , z c )-(−50 mm, 0, 0); screen origin (x s , y s , z s )-(50 mm, 0, 0); mirror origin (x m , y m , z m )-(0, 0,−1400 mm).A set of images with horizontal and vertical fringes is generated, and the phase information is captured and reconstructed using a weighted phase unwrapping technique [21,22] where the petal geometry is used as a mask.The shape is then reconstructed from the slope information using the discontinuity-preserving Mumford Shah integration technique [23,24].The reconstructed aberration can be seen in Figure 4 (left), and the error between reconstruction and nominal aberration can be seen in Figure 4 (right).The main component of the error shown is from the assumption that the reflector is flat from the point of view of the camera.This assumption leads to a coma error in the case of a full reflector, but the petal integration leads to a more complex error.With this error known, it may be possible to subtract it in the slope calculation for a specific setup, but the accuracy shown is sufficient for our purposes: characterizing different manufacturing techniques and ensuring the resulting errors are within the control stroke of the adaptive reflector.With the ray-tracing setup, geometric tolerances can also be simulated for expected deviations.Ordinary machine shop measurement tools are used to place the screen and camera with expected tolerances of 2 mm and 2 degrees.Simulations are then made in which the camera is moved in the ray-tracing software, but the input to the reconstruction remains the same.The results are shown in Table 1.Similar deviations are expected from moving the screen or mirror.The errors reported here are higher than in the previous literature [8], with various possible explanations such as the geometric positions of the test components or the influence of the tested aberration.

Experimental Results
By using the ray-tracing software and comparing the expected deformation with the reconstructed deformation, we confirmed our software implementation of SCOTS works as expected even for the petal geometry.The next step would be to use the implemented software to measure samples.In this case, it is also helpful to have a verification of the hardware setup.As illustrated in Figure 5, this can also be performed using raytracing.A series of images may be captured experimentally, and passed to the SCOTS implementation to reconstruct a 3D model of the measured geometry.The 3D model may then be passed to the ray-tracer, simulating an image.The simulated image may then be overlaid with the real image.If the images match accurately, then no major errors are present in the reconstruction.Figure 6 shows two demonstrators manufactured by MateriaNova in the context of the previous project (MATS).The demonstrators consist of a curved polymer shell with spray-coated PVDF-TrFE between vacuum-deposited patterned electrodes, allowing for active control (see also Figure 1).The SCOTS implementation is tested on a similar sample (aluminum-coated polymer only) with 20 cm diameter and 2.5 m radius of curvature.The simulated image and real image are compared in Figure 7.This confirms a close match between the experimental setup and the virtual experiment.Thus, we expect that the tolerances reported in Table 1 are representative of the actual accuracy of the test.With this, we have laid the necessary ground work for measuring deployable petal reflectors in the laboratory.

Modal Wavefront Reconstruction
Consider the petal reflector of Figure 8 provided with a micro-lens sensor array (Shack-Hartmann) providing the average x-y slopes at the control points.Due to the obscuration of the secondary mirror, the petals appear disconnected, suggesting that the primary reflector should be considered as segmented.However, the petals have an interesting property that segmented mirrors do not have: they share the same mechanical boundary conditions at the inner support ring of the reflector.This property will be used to reconstruct an approximation of the aberration over the whole reflector surface from slope data.

Vibration Modes
Let us consider the vibration modes φ i (Figure 9) of the petal reflector.They constitute a complete set of orthogonal functions, solutions of the eigenvalue problem where K and M are, respectively, the stiffness and the mass matrices of a finite element (FE) model of the structure (the spherical shell petal reflector in this case) and ω i are the natural frequencies.The mode shape amplitudes φ i contain the numerical values of the degrees of freedom (d.o.f.) of the finite element model.Using a Guyan reduction [18,19], the degrees of freedom can be selected to contain only the z-components of the nodal displacements of the FE model.Equation ( 5) does not define the amplitude, but only the shape of the modes; it is customary to normalize them in such a way that φ T i Mφ i = 1 (unit modal mass).Figure 9 shows a few examples of the vibration modes of the petal reflector; they all satisfy the mechanical boundary conditions; their shapes have increasing complexity as the order increases.They satisfy the orthogonality condition φ T i Mφ j = 0 if (i = j).Using the matrix notation Φ = (φ 1 , φ 2 , . . ., φ n ), the orthogonality condition reads Any surface figure error w may be expanded in the mode shapes The modal components a are easily obtained from the condition of orthogonality: The accuracy of the modal expansion increases with the number of modes, as illustrated in Figure 10; the assumed surface figure error is a tip displacement of 100 µm in the corner of one petal; the reflector has a diameter of 30 cm and a radius of curvature of 0.85 m.Triangular Mindlin elements in FEM solver SAMCEF are used, with 3066 nodes in the Guyan reduction.The figure shows the relative RMS residual error as a function of the number of modes n included in the expansion.One sees that the residual error reduces rapidly as n increases.
Note that the set of modes used in the modal expansion, as seen in Equation ( 7), does not have to be the modes of the reflector itself, but they must have the same mechanical boundary conditions at the inner support ring.To illustrate this, Figure 10 also shows the relative RMS residual error when the mode shapes of a petal flat plate (with the same inner ring support and outer contour) are used instead of the spherical shell.

Jacobian
The Shack-Hartmann (S-H) sensor is one of the most popular wavefront sensors; it measures the average slopes (s x , s y ) of the wavefront at an array of discrete points corresponding to the various sub-apertures of the lenslet array.The total number of sensor outputs is twice the number of active lenslets, n s .Such a sensor alone is not appropriate for segmented mirrors, because it cannot detect relative rigid body displacements between segments along z.However, for the petal reflector with central obscuration considered here, the restriction is alleviated due to the common support condition at the inner ring.
Taking the partial derivative of Equation ( 7) with respect to x and y at the Shack-Hartman sensor points, where Φ x and Φ y are the matrices of the partial derivatives of the mode shapes at the regular S-H sensor points.These matrices are constructed by interpolation from the FE mode shapes.The matrix is the Jacobian of the mode shapes.If the sensor measurements are s = (s T x , s T y ) T , the modal amplitudes are solutions of the equation If the number of active lenslets is n s and the number of modes is n, the size of the Jacobian is (2n s × n).Assuming that 2n s > n, the solution is given by the Moore-Penrose pseudo-inverse a = J + s (12) can be readily obtained by singular value decomposition (SVD) [25].If J is illconditioned, the situation may be resolved by one of the following ways: (i) truncating the singular values lower to a threshold in the pseudo-inverse expansion; (ii) using the damped least squares (DLS) method [26].
To illustrate the reconstruction method, consider a sensor array of 30 × 30 lenslets; the method is used to reconstruct the tip displacement aberration considered earlier.The reconstructed surface based on slope measurements and the residual error are shown in Figure 11.The RMS residual error is compared with that based on the orthogonality condition in Table 2 for n s = 604 active lenslets and n = 50, 100 modes.

Central Obstruction
A central obscuration from the secondary mirror is now considered (Figure 8); the central obscuration is defined by the ratio D S /D P between the diameter of the secondary mirror D S and that of the primary mirror D P .The petals appear as disconnected when D S /D P > 20%; an obstruction of 30% is assumed, so that the petals appear as disconnected.The reduced Jacobian is easily constructed from the full one by considering only the sensor points outside the obstructed area where the measurement vector s is available.Figure 12 shows the reconstructed surface and the residual error outside the obstructed area.The RMS residual error is that shown in the last column of Table 2. Once the modal expansion of the aberration has been estimated, the z displacement can be reconstructed anywhere on the reflector.In particular, reconstructing the displacements on adjacent sides of the petals can provide virtual edge sensors that can be combined with the slope measurements in the shape control of the reflector [1][2][3].Further studies will be conducted to evaluate the benefit of the procedure for various shapes of the aberration.In particular, the reflector is sensitive to the clamping at the boundary condition, and imperfect clamping leading to sharp transitions may be difficult to recover by the proposed method.

Conclusions
Two aspects of the metrology of spherical petal reflectors have been addressed:

•
The first aspect is that of measuring the surface figure error of a spherical reflector alone, during the development and manufacturing phase.The requested accuracy is modest, because the reflector is intended to be actively controlled once in operation.
The SCOTS approach has been found to be fast and satisfactory.Experimental results have been presented and their consistency with a ray-tracing virtual experiment has been assessed.

•
The second part of this paper is concerned with the surface figure error reconstruction from slope measurements of a petal reflector when a central part of the mirror is obscured by the secondary mirror of the telescope, making the petals appear as completely disconnected (like segments).Using the fact that all petals have the same mechanical boundary conditions, the deformed shape is expanded in a set of orthogonal modes having the same boundary conditions (the vibration modes).The modal amplitudes are reconstructed from slope data (Shack-Hartmann) and an approximation of the surface figure error is obtained.

Figure 2 .
Figure 2. Illustration of the SCOTS test.The two vectors m2c and m2s together describe the surface slope.The blue dot illustrates how the recorded image relates a position of the Mirror Pixel to its corresponding screen location.The red dot shows the zero-phase location.Several images are needed to accurately relate the projected and recorded images in terms of phase.

Figure 3 .
Figure 3. (left) A 100−micron tip displacement aberration simulated with FEA.(right) The resulting ray-traced image from the simulation.Note that, on the right side, the fringes do not line up exactly due to the tip displacement.

Figure 4 .
Figure 4. (left) Reconstructed aberration from simulated images of 100−micron tip displacement.The geometry is slightly eroded to avoid edge problems.This also reduces the apparent PV slightly.(right) Error between the nominal aberration of Figure3(left) and the measured aberration of Figure4(left).The main discrepancy arises from the assumption of zero curvature on the primary mirror.

Figure 5 .
Figure 5. Illustration of the different operating modes of the software.(a) The raytracing can be used to test the software implementation by comparing to a known aberration.(b) the ray-tracing can be used to verify the hardware implementation by reconstructing a virtual image from the measured surface.

Figure 6 .
Figure 6.Demonstrators of the previous project [3] with individual electrode control.

Figure 7 .
Figure 7. Image overlay of the real image (green) and ray-traced simulation of the measured geometry (red).The two images closely match, showing that the geometry reported by the SCOTS implementation corresponds closely to the real world.

Figure 8 .
Figure 8. Petal reflector, Shack-Hartmann sensor array, and secondary mirror obscuration.The petals appear disconnected, like on a segmented mirror.

Figure 9 .
Figure 9. Example of vibration modes of the petal reflector, which will be used in the reconstruction of the aberration.

Figure 10 .
Figure 10.Evolution of the relative RMS residual error as a function of the number of modes in the expansion.Comparison between the mode shapes of the petal spherical shell and the corresponding petal plate.The surface figure error is a tip displacement of 100 µm in the corner of one petal.

Figure 11 .
Figure 11.Reconstructed surface (left) and residual error (right) of a 100 micron tip displacement.

Table 1 .
Error between the nominal aberration and the measured aberration including geometric tolerances.

Table 2 .
Relative RMS error after modal reconstruction based on orthogonality condition (8) and slopes measurements, for 0% and 30% central obstruction with 50 and 100 modes in the reconstruction.