Simulation of the Pulse Electrochemical Machining (PECM) process can provide information on system design and guidelines for practical use. The pulses that are applied to the PECM system have to be described on a time scale that can be orders of magnitude smaller than the physical time scales in the system. If the full detail of the applied pulses has to be taken into account, the time accurate calculation of the variable distribution evolutions in PECM can become a computationally very expensive procedure. In previous work of the authors, approximate techniques were introduced: the hybrid calculation and the Quasi Steady State Shortcut (QSSSC). In other previous work of the authors a model for PECM of steel in NaNO 3 was introduced. This model contains a changing polarization behaviour of the double layer as a function of the metal ion surface concentration, which brings a strong non-linearity in the system. In this paper a technique is introduced to integrate the non-linear model into the approximate methods. To achieve this, the strategy of the approximate methods is extended. For the QSSSC, the nonlinearity is handled using an extra convergence level. For the hybrid calculation, live averaging is used to take care of the non-linear effects. Performing this, the timesteps used during the high level calculations are no longer dictated by the pulse characteristics. Using this approach, computationally very cheap, yet satisfying results can be obtained. The technique is very general and very powerful and can be used in any multi-timescale system.
List of symbols aPolarization parameter 1 (S m -2 ); A Electrode surface (m 2 ); bPolarization parameter 2 (A m -2 ); c Concentration (mol m -3 );Thermal conductivity (W m -1 K -1 ) P dl Heat produced, in the double layer (W m -2 ) P bulk Heat produced in the bulk (W m -3 ) Pr t Turbulent Prandtl number (-) rGeneral location vector (m) ReReynoldsVelocity (m s -1 ) wWater depletion factor (-) x Distance (m) z Valence (-)Greek symbols a Duty cycle (-) g Overpotential (V) H Temperature (K) lDynamic viscosity (kg m -1 s -1 )