Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

Joint DIII-D/EAST research on the development of a high poloidal beta scenario for the steady state missions of ITER and CFETR

Experimental and modeling investigations on the DIII-D and EAST tokamaks show the attractive transport and stability properties of fully noninductive, high poloidal-beta (β P ) plasmas, and their suitability for steady-state operating scenarios in ITER and CFETR. A key feature of the high-β P regime is the large-radius (r>0.6) internal transport barrier (ITB), often observed in all channels (ne, Te, Ti, rotation), and responsible for both excellent energy confinement quality and excellent stability properties. Experiments on DIII-D have shown that, with a largeradius ITB, very high β N and β P values (both≥4) can be reached by taking advantage of the stabilizing effect of a nearby conducting wall. Synergistically, higher plasma pressure provides turbulence suppression by Shafranov shift, leading to ITB sustainment independent of the plasma rotation. Experiments on EAST have been used to assess the long pulse potential of the high-β P regime. Using RF-only heating and current drive, EAST achieved minute-long fully noninductive steady state H-mode operation with strike points on an ITER-like tungsten divertor. Improved confinement (relative to standard H-mode) and steady state ITB features are observed with a monotonic q-profile with q min ~1.5. Separately, experiments have shown that increasing the density in plasmas driven by lower hybrid wave broadens the q-profile, a technique that could enable a large radius ITB. These experimental results have been used to validate MHD, current drive, and turbulent transport models, and to project the high-β P regime to a burning plasma. These projections suggest the Shafranov shift alone will not suffice to provide improved confinement (over standard H-mode) without rotation and rotation shear. However, increasing the negative magnetic shear (higher q on axis) provides a similar turbulence suppression mechanism to Shafranov shift, and can help devices such as ITER and CFETR achieve their steady-state fusion goals.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of window temperature variations on ITER toroidal interferometer and polarimeter (TIP) measurements

Abstract The toroidal interferometer and polarimeter (TIP) is one of the primary electron density diagnostics on ITER. To meet measurement requirements over several thousand seconds, environmental factors such as changes in the window temperature, air temperature, and humidity, which can cause uncompensated phase drifts, must be minimized. This paper reports measurements of the phase shift induced by ZnSe and BaF₂ windows as their temperature is varied utilizing the TIP prototype. To accomplish these measurements, test pieces of BaF₂ and ZnSe are placed in a small oven that is located in one leg of the 10.59/5.22 micron two-color interferometer TIP prototype. The oven temperature was varied from room temperature to ∼90 °C, which is slightly higher than the expected window temperatures on ITER. The vibration compensated phase shifts measured for these two materials are dϕᵥ꜀ /dT/L = 0.0121 and −0.275 deg. /°C/cm −1 for BaF₂ and ZnSe, respectively. From these measurements, it was concluded that temperature variations of the environmental/secondary confinement barrier window and primary vacuum window have to be suppressed within approximately 1°C and 2°C, respectively, in order to passively meet the electron density measurement requirements in ITER.

Akiyama, T (ORCID:0000000198469795)↗

Rotation and neoclassical ripple transport in ITER

Neoclassical transport in the presence of non-axisymmetric magnetic fields causes a toroidal torque known as neoclassical toroidal viscosity (NTV). The toroidal symmetry of ITER will be broken by the finite number of toroidal field coils and by test blanket modules (TBMs). The addition of ferritic inserts (FIs) will decrease the magnitude of the toroidal field ripple. 3D magnetic equilibria in the presence of toroidal field ripple and ferromagnetic structures are calculated for an ITER steady-state scenario using the Variational Moments Equilibrium Code (VMEC). Furthermore, neoclassical transport quantities in the presence of these error fields are calculated using the Stellarator Fokker-Planck Iterative Neoclassical Conservative Solver (SFINCS).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Advances in the steady-state hybrid regime in DIII-D – a fully non-inductive, ELM-suppressed scenario for ITER

Here, the hybrid regime with beta, collisionality, safety factor and plasma shape relevant to the ITER steady-state mission has been successfully integrated with ELM suppression by applying an odd parity n=3 resonant magnetic perturbation (RMP). Fully non-inductive hybrids in the DIII-D tokamak with high beta (< β> ≤ 2.8%) and high confinement (H 98y2 ≤ 1.4) in the ITER similar shape have achieved zero surface loop voltage for up to two current relaxation times using efficient central current drive from ECCD and NBCD. The n=3 RMP causes surprisingly little increase in thermal transport during ELM suppression. Poloidal magnetic flux pumping in hybrid plasmas maintains q above 1 without loss of current drive efficiency, except that experiments show that extremely peaked ECCD profiles can create sawteeth. During ECCD, Alfvén eigenmode (AE) activity is replaced by a more benign fishbone-like mode, reducing anomalous beam ion diffusion by a factor of 2. While the electron and ion thermal diffusivities substantially increase with higher ECCD power, the loss of confinement can be offset by the decreased fast ion transport resulting from AE suppression. Extrapolations from DIII-D along a dimensionless parameter scaling path as well as those using self-consistent theory-based modeling show that these ELM-suppressed, fully non-inductive hybrids can achieve the Q fus = 5 ITER steady-state mission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Use of reconstructed 3D equilibria to determine onset conditions of helical cores in tokamaks for extrapolation to ITER

The threshold for spontaneous symmetry breaking, which results in growth of m/n = 1/1 helical cores, in tokamaks is determined using VMEC. This model is based on a DIII-D hybrid discharge with helical core, and it predicts ITER (15MA scenario) to operate far in the helical core formation regime. Helical cores can only exist in tokamak discharges with monotonic but low, or reversed q-shear and qmin 1 in the core. The helical core is a saturated internal kink mode; its onset threshold is proportional to (dp/d)/B 2 t around q = 1. Below the threshold, applied 3-D fields can drive a helical core to finite size, as in DIIID. Above it, a random 3-D kick causes a bifurcation from axisymmetry and excites a spontaneous helical core. Helical cores occur frequently in C-Mod due to impurities; modeling shows the axisymmetric equilibrium becomes unstable at the end of the current ramp up due to the slow current penetration and a resultant reversed shear q-profile, which lowers the threshold. Lastly, a comparison of the helical core onset threshold for discharges from DIII-D, C-Mod and ITER confirms that while DIII-D is marginally stable, C-Mod and especially ITER are highly susceptible to helical core formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessment of the baseline scenario at q 95 ~ 3 for ITER

The International Tokamak Physics Activity topical group on integrated operational scenarios has compiled a database of stationary H-mode discharges at q 95 ~ 3 from AUG, C-Mod, DIII-D, JET and JT-60U, for both carbon wall and high-Z metal wall experiments with ~3300 entries. The analyses focus on discharges that are stationary for ≥5 thermal energy confinement times to evaluate the baseline scenario proposed for ITER at 15 MA for achieving its goals of Q = 10, fusion power of 500 MW at normalized pressure, β N = 1.8 and normalized confinement as predicted by the standard H-mode scaling, H 98y2 = 1. With the data restricted to stationary H-modes at q 95 ~ 3, the database shows significant variation of thermal energy confinement compared to the standard H-mode scaling (IPB98(y,2)) in dimensionless form. The data show similar scaling with normalized gyro-radius, but more favorable scaling towards lower collision frequency and more favorable scaling with plasma beta. Using all the engineering variables employed in IPB98(y,2), results in an overfit due to correlations among the data. Furthermore, there are significant residual trends in the confinement for plasma current, device size, loss power, and in particular for the plasma density. Significant differences between results obtained for devices with a carbon wall and high-Z metal wall are observed in the data, with data from carbon wall devices providing a larger operating space, encompassing ITER parameters or even exceeding them. H-modes in high-Z metal wall devices have, so-far, not accessed conditions at low collision frequencies, have lower normalized confinement (H 98y2 ~ 0.8–0.9) at low input power or beta, achieving H 98y2 ~ 1.0 only at input powers two times the L- to H-mode transition scaling predictions and at β N ~ 2.0. Hence, only the best H-modes with high-Z metal walls reach ITER baseline performance requirements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear burn control in ITER using adaptive allocation of actuators with uncertain dynamics

Abstract ITER will be the first tokamak to sustain a fusion-producing, or burning, plasma. If the plasma temperature were to inadvertently rise in this burning regime, the positive correlation between temperature and the fusion reaction rate would establish a destabilizing positive feedback loop. Careful regulation of the plasma’s temperature and density, or burn control, is required to prevent these potentially reactor-damaging thermal excursions, neutralize disturbances and improve performance. In this work, a Lyapunov-based burn controller is designed using a full zero-dimensional nonlinear model. An adaptive estimator manages destabilizing uncertainties in the plasma confinement properties and the particle recycling conditions (caused by plasma–wall interactions). The controller regulates the plasma density with requests for deuterium and tritium particle injections. In ITER-like plasmas, the fusion-born alpha particles will primarily heat the plasma electrons, resulting in different electron and ion temperatures in the core. By considering separate response models for the electron and ion energies, the proposed controller can independently regulate the electron and ion temperatures by requesting that different amounts of auxiliary power be delivered to the electrons and ions. These two commands for a specific control effort (electron and ion heating) are sent to an actuator allocation module that optimally maps them to the heating actuators available to ITER: an electron cyclotron heating system (20 MW), an ion cyclotron heating system (20 MW), and two neutral beam injectors (16.5 MW each). Two different actuator allocators are presented in this work. The first actuator allocator finds the optimal mapping by solving a convex quadratic program that includes actuator saturation and rate limits. It is nonadaptive and assumes that the mapping between the commanded control efforts and the allocated actuators (i.e. the effector model) contains no uncertainties. The second actuator allocation module has an adaptive estimator to handle uncertainties in the effector model. This uncertainty includes actuator efficiencies, the fractions of neutral beam heating that are deposited into the plasma electrons and ions, and the tritium concentration of the fueling pellets. Furthermore, the adaptive allocator considers actuator dynamics (actuation lag) that contain uncertainty. This adaptive allocation algorithm is more computationally efficient than the aforementioned nonadaptive allocator because it is computed using dynamic update laws so that finding the solution to a static optimization problem is not required at every time step. A simulation study assesses the performance of the proposed adaptive burn controller augmented with each of the actuator allocation modules.

Physics↗

Toward holistic understanding of the ITER-like resonant magnetic perturbation (RMP) ELM control on KSTAR

Here, KSTAR has clarified a set of unresolved 3D physics issues utilizing the ITER-like in-vessel, three-row, resonant magnetic perturbation (RMP) configurations. Since RMP-driven, edge-localized-modes (ELMs)-crash control elevates the divertor heat flux peak through its impact on edge plasma parameters and transport, a series of intentionally misaligned RMP configurations (IMCs) have been explored to investigate the relationship between RMP ELM control and divertor heat fluxes, while searching for an ideal IMC that could be favorable in both aspects. First of all, the contrasting influence of kink vs anti-kink phasing on the ELM-crash suppression has been articulated, demonstrating the synergistic benefit of 'kink' phasing on ELM-crash-suppression. On the other hand, the three-row IMC in the anti-kink phasing becomes more insensitive to the ELM-crashes at the sub-marginal level of RMP, consistent with theory. Meanwhile, the divertor 'wetted' area of ELM-crash-suppression gets narrower than that of ELM-crash-mitigation, suggesting that ELM-crash-mitigation remains advantageous over ELM-crash-suppression in terms of time-averaged divertor thermal loading. In comparison, based on a set of two-row IMCs, no evidence of divertor heat flux broadening was found during ELM-crash-suppression, supporting a hypothesis that the dispersal of the divertor heat flux in three-row IMCs cannot be driven by helically structured two-row RMPs alone. Among ITER-like three-rows, lower two-row RMPs have been found to be much more effective in suppressing the ELM-crashes than upper two-row RMPs. Although it is quite preliminary, the up/down asymmetric dependence of RMP coupling may be generically attributed to lower-single-null plasmas. Such a holistic understanding of RMP-driven, ELM-crash-control in KSTAR is expected not only to elucidate various subtle points in the vicinity of ELM-crash-suppression, but also to clarify the relevant divertor thermal loading issues for ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical error analysis of SOLPS-ITER simulations of EAST

Abstract Plasma edge simulations with codes like SOLPS-ITER are widely employed to interpret fusion experiments. However, numerical errors appearing in such simulations are rarely investigated, despite their potential large impact on simulation results. These errors consist of the statistical error and the bias, both resulting from the finite number of employed EIRENE Monte Carlo particles and incomplete convergence, and the discretization error due to the finite resolution of the computational grids. In this contribution, the resulting numerical errors on simulations of pure deuterium and neon seeded H-mode EAST discharges are examined. The statistical error can be kept small compared to other numerical error contributions by averaging the plasma profiles. This allows investigating the bias and discretization errors using Richardson extrapolation. It is shown that grid refinement and the number of employed Monte Carlo particles have the largest influence on the result, in agreement with similar studies of an ITER deuterium case. For the first time, numerical error bars on the entire simulated target profiles are determined showing that the largest numerical error is 17.9%, mainly due to the plasma grid discretization. On top, also numerical errors on simulated neutral pressures are investigated in detail, for which the statistical error is dominant. The analysis demonstrates which setup is needed to keep numerical errors limited: the SOLPS-ITER averaging procedure should be employed including enough EIRENE particles, and the involved grids should be sufficiently refined to reduce discretization errors.

Boeyaert, Dieter (ORCID:0000000309208660)↗

Overview of T and D–T results in JET with ITER-like wall

In 2021 JET exploited its unique capabilities to operate with T and D–T fuel with an ITER-like Be/W wall (JET-ILW). This second major JET D–T campaign (DTE2), after DTE1 in 1997, represented the culmination of a series of JET enhancements—new fusion diagnostics, new T injection capabilities, refurbishment of the T plant, increased auxiliary heating, in-vessel calibration of 14 MeV neutron yield monitors—as well as significant advances in plasma theory and modelling in the fusion community. DTE2 was complemented by a sequence of isotope physics campaigns encompassing operation in pure tritium at high T-NBI power. Carefully conducted for safe operation with tritium, the new T and D–T experiments used 1 kg of T (vs 100 g in DTE1), yielding the most fusion reactor relevant D–T plasmas to date and expanding our understanding of isotopes and D–T mixture physics. Furthermore, since the JET T and DTE2 campaigns occurred almost 25 years after the last major D–T tokamak experiment, it was also a strategic goal of the European fusion programme to refresh operational experience of a nuclear tokamak to prepare staff for ITER operation. The key physics results of the JET T and DTE2 experiments, carried out within the EUROfusion JET1 work package, are reported in this paper. Progress in the technological exploitation of JET D–T operations, development and validation of nuclear codes, neutronic tools and techniques for ITER operations carried out by EUROfusion (started within the Horizon 2020 Framework Programme and continuing under the Horizon Europe FP) are reported in (Litaudon et al Nucl. Fusion accepted), while JET experience on T and D–T operations is presented in (King et al Nucl. Fusion submitted).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Divertor-safe nonlinear burn control based on a SOLPS parameterized core-edge model for ITER

Abstract For ITER operations, the range of desirable burning-plasma regimes with high fusion power output will be restricted by various operational constraints. These constraints include the saturation of ITER’s various heating and fueling actuators such as the neutral beam injectors, the ion and electron cyclotron heating systems, the gas puffing system, and the deuterium–tritium pellet injectors. In addition to these actuator constraints, the H-mode power threshold, divertor detachment, and the heat load on the divertor targets may apply limitations to ITER’s operational space. In this work, Plasma Operation Contour (POPCON) plots that map the aforementioned constraints to the temperature-density space are used to investigate which constraints are most limiting towards accessing regimes with high fusion power output. The presented POPCON plots are based on a control-oriented core-edge model that couples the nonlinear density and energy response models for the core-plasma region with SOLPS4.3 parameterizations for conditions in the edge-plasma regions (scrape-off-layer and divertor). Using this control-oriented core-edge model, a nonlinear burn controller, which aims to regulate the plasma temperature and density in the core-plasma region, is constructed in this work. This controller is augmented with an online optimization scheme that governs the control references such that the plasma can be guided towards regimes with high fusion powers while protecting the divertor targets from dangerously high heat loads. A closed-loop simulation study illustrates the capability of this burn control scheme.

Physics↗

Progress in ITER ECE diagnostic design and integration

The ITER electron cyclotron emission (ECE) diagnostic system has primary roles in providing measurements of the core electron temperature profile and the electron temperature fluctuation associated with the neoclassical tearing modes. The ITER ECE system includes a radial and oblique line-of-sight. Four 43-meter long low-loss transmission lines (TLs) are designed to transmit millimeter wave power in the frequency range of 70–1000 GHz in both X- and O-mode polarization from the port plug to the ECE instrumentation room in the diagnostic building. The measurement instrumentation includes two Fourier transform spectrometer (FTS) systems and two radiometer systems. The Indian Domestic Agency (IN-DA) and United States Domestic Agency share the responsibility. The IN-DA scope excluding instrumentation and control has passed its preliminary design review and is progressing towards the final design review (FDR). In parallel, the diagnostic integration in different areas is ongoing. Furthermore, several captive components for the TLs have passed FDR and will be manufactured for installation in the tokamak building soon. A peer review meeting has been held on the prototype hot calibration source, and its integration and new thermal analysis in the diagnostic shield module are continuing. A prototype TL is being tested. A prototype polarizing Martin-Puplett type FTS, operating in the frequency range 70–1000 GHz, features an in-vacuo fast scanning mechanism and a cryo-cooled dual-channel THz detector system. Its performance has been assessed in detail against ITER requirements.

47 OTHER INSTRUMENTATION↗

tomoCAM : fast model-based iterative reconstruction via GPU acceleration and non-uniform fast Fourier transforms

X-ray-based computed tomography is a well established technique for determining the three-dimensional structure of an object from its two-dimensional projections. In the past few decades, there have been significant advancements in the brightness and detector technology of tomography instruments at synchrotron sources. These advancements have led to the emergence of new observations and discoveries, with improved capabilities such as faster frame rates, larger fields of view, higher resolution and higher dimensionality. These advancements have enabled the material science community to expand the scope of tomographic measurements towards increasingly in situ and in operando measurements. In these new experiments, samples can be rapidly evolving, have complex geometries and restrictions on the field of view, limiting the number of projections that can be collected. In such cases, standard filtered back-projection often results in poor quality reconstructions. Iterative reconstruction algorithms, such as model-based iterative reconstructions (MBIR), have demonstrated considerable success in producing high-quality reconstructions under such restrictions, but typically require high-performance computing resources with hundreds of compute nodes to solve the problem in a reasonable time. Here, tomoCAM , is introduced, a new GPU-accelerated implementation of model-based iterative reconstruction that leverages non-uniform fast Fourier transforms to efficiently compute Radon and back-projection operators and asynchronous memory transfers to maximize the throughput to the GPU memory. The resulting code is significantly faster than traditional MBIR codes and delivers the reconstructive improvement offered by MBIR with affordable computing time and resources. tomoCAM has a Python front-end, allowing access from Jupyter -based frameworks, providing straightforward integration into existing workflows at synchrotron facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optimal Iterative Method for Network Utility Maximization with Intertemporal Constraints

Network utility maximization (NUM) has been widely used to solve the optimal resource allocation problem in wireless network and power systems. The network planner aims to allocate generation and demand over producers and consumers such that social welfare is maximized while supply and demand are balanced, and network constraints such as branch capacity constraints are satisfied. However, when taking into account dynamics of agents, such as ramp rate constraints on power output for suppliers and impact of cumulative prior consumption on current demand for consumers, a single-period formulation is not sufficient. We propose an iterative method for the multi-period NUM with intertemporal constraints based on dual decomposition. It features a price iteration scheme that achieves optimal social welfare, and the algorithm is guaranteed to converge if proper concavity condition holds. We illustrate the iterative algorithm by power system examples. It is shown that similar to locational marginal price (LMP) in power systems, time-varying congestion prices arise when there exists binding branch constraints. Simulation results also demonstrate that there is increase in social welfare compared to a single-period formulation.

Ma, Ke↗

Assessment of the dose reduction potential of a model-based iterative reconstruction algorithm using a task-based performance metrology

Purpose: Different computed tomography (CT) reconstruction techniques offer different image quality attributes of resolution and noise, challenging the ability to compare their dose reduction potential against each other. The purpose of this study was to evaluate and compare the task-based imaging performance of CT systems to enable the assessment of the dose performance of a model-based iterative reconstruction (MBIR) to that of an adaptive statistical iterative reconstruction (ASIR) and a filtered back projection (FBP) technique. Methods: The ACR CT phantom (model 464) was imaged across a wide range of mA setting on a 64-slice CT scanner (GE Discovery CT750 HD, Waukesha, WI). Based on previous work, the resolution was evaluated in terms of a task-based modulation transfer function (MTF) using a circular-edge technique and images from the contrast inserts located in the ACR phantom. Noise performance was assessed in terms of the noise-power spectrum (NPS) measured from the uniform section of the phantom. The task-based MTF and NPS were combined with a task function to yield a task-based estimate of imaging performance, the detectability index (d′). The detectability index was computed as a function of dose for two imaging tasks corresponding to the detection of a relatively small and a relatively large feature (1.5 and 25 mm, respectively). The performance of MBIR in terms of the d′ was compared with that of ASIR and FBP to assess its dose reduction potential. Results: Results indicated that MBIR exhibits a variability spatial resolution with respect to object contrast and noise while significantly reducing image noise. The NPS measurements for MBIR indicated a noise texture with a low-pass quality compared to the typical midpass noise found in FBP-based CT images. At comparable dose, the d′ for MBIR was higher than those of FBP and ASIR by at least 61% and 19% for the small feature and the large feature tasks, respectively. Compared to FBP and ASIR, MBIR indicated a 46%–84% dose reduction potential, depending on task, without compromising the modeled detection performance. Conclusions: The presented methodology based on ACR phantom measurements extends current possibilities for the assessment of CT image quality under the complex resolution and noise characteristics exhibited with statistical and iterative reconstruction algorithms. The findings further suggest that MBIR can potentially make better use of the projections data to reduce CT dose by approximately a factor of 2. Alternatively, if the dose held unchanged, it can improve image quality by different levels for different tasks.

60 APPLIED LIFE SCIENCES↗

Soft syndrome iterative decoding of quantum LDPC codes and hardware architectures

In practical quantum error correction implementations, the measurement of syndrome information is an unreliable step—typically modeled as a binary measurement outcome flipped with some probability. However, the measured syndrome is in fact a discretized value of the continuous voltage or current values obtained in the physical implementation of the syndrome extraction. In this paper, we use this “soft” or analog information to benefit iterative decoders for decoding quantum low-density parity-check (QLDPC) codes. Syndrome-based iterative belief propagation decoders are modified to utilize the soft syndrome to correct both data and syndrome errors simultaneously. We demonstrate the advantages of the proposed scheme not only in terms of comparison of thresholds and logical error rates for quasi-cyclic lifted-product QLDPC code families but also with faster convergence of iterative decoders. Additionally, we derive hardware (FPGA) architectures of these soft syndrome decoders and obtain similar performance in terms of error correction to the ideal models even with reduced precision in the soft information. The total latency of the hardware architectures is about 600 ns (for the QLDPC codes considered) in a 20 nm CMOS process FPGA device, and the area overhead is almost constant—less than 50% compared to min-sum decoders with noisy syndromes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Ranger-based Iterative Random Forest

Iterative Random Forest (iRF) is an improvement upon the classic Random Forest, using weighted iterations to distill the forests. Ranger is a C++ implementation of random forest or recursive partitioning, particularly suited for high dimensional data. Ranger-Based iRF uses Ranger as a base Random Forest codebase upon which iterative and scalability aspects have been added.

Jacobson, DanielA↗

RISE: robust iterative surface extension for sub-nanometer X-ray mirror fabrication

Precision optics have been widely required in many advanced technological applications. X-ray mirrors, as an example, serve as the key optical components at synchrotron radiation and free electron laser facilities. They are rectangular silicon or glass substrates where a rectangular Clear Aperture (CA) needs to be polished to sub-nanometer Root Mean Squared (RMS) to keep the imaging capability of the incoming X-ray wavefront at the diffraction limit. The convolutional polishing model requires a CA to be extended with extra data, from which the dwell time is calculated via deconvolution. However, since deconvolution is very sensitive to boundary errors and noise, the existing surface extension methods can hardly fulfill the sub-nanometer requirement. On one hand, the figure errors in a CA were improperly modeled during the extension, leading to continuity issues along the boundary. On the other hand, uncorrectable high-frequency errors and noise were also extended. In this study, we propose a novel Robust Iterative Surface Extension (RISE) method that resolves these problems with a data fitting strategy. RISE models the figure errors in a CA with orthogonal polynomials and ensures that only correctable errors are fit and extended. Combined with boundary conditions, an iterative refinement of dwell time is then proposed to compensate the errors brought by the extension and deconvolution, which drastically reduces the estimated figure error residuals in a CA while the increase of total dwell time is negligible. To our best knowledge, RISE is the first data fitting-based surface extension method and is the first to optimize dwell time based on iterative extension. An experimental verification of RISE is given by fabricating two elliptic cylinders (10 mm × 80 mm CAs) starting from a sphere with a radius of curvature around 173 m using ion beam figuring. The figure errors in the two CAs greatly improved from 204.96 nm RMS and 190.28 nm RMS to 0.62 nm RMS and 0.71 nm RMS, respectively, which proves that RISE is an effective method for sub-nanometer level X-ray mirror fabrication.

36 MATERIALS SCIENCE↗