Digital simulations boundary conditions

To model this method using digital simulation techniques, one need only change the electrode boundary conditions after some predetermined number of time iterations (representing tf) have taken place. The electrode boundary conditions become... 

One of the main uses of digital simulation - for some workers, the only application - is for linear sweep (LSV) or cyclic voltammetry (CV). This is more demanding than simulation of step methods, for which the simulation usually spans one observation time unit, whereas in LSV or CV, the characteristic time r used to normalise time with is the time taken to sweep through one dimensionless potential unit (see Sect. 2.4.3) and typically, a sweep traverses around 24 of these units and a cyclic voltammogram twice that many. Thus, the explicit method is not very suitable, requiring rather many steps per unit, but will serve as a simple introduction. Also, the groundwork for the handling of boundary conditions for multispecies simulations is laid here. 

Equations (6) to (10) and (3) form a system of second order partial differential equations with initial [Eq. (8)] and boundary conditions [Eqs. (3), (9) and (10)]. The presence of the time dependence in Eq. (3) prevents a closed form solution of this system. Several approaches have been used to calculate the concentrations c as a function of X and t, including Laplace transformation [7, 12] and digital simulation [25], all relying on numerical integration steps. Chapter 1.3 provides further information on digital simulation. 

The solution of this partial differential second-order equation depends on the initial and boundary conditions of the particular experiment, giving rise to a multitude of techniques. In Chapter 2.2, the digital simulation of voltammetry under stagnant and hydrodynamic conditions is described. By changing the electrode potential one can modify the boundary conditions and transient effects arise until a new steady state is reached. 

However, it is not easy to calculate the exact relation between the probe current and the probe-sample distance in the feedback mode due to the geometric complexity in the experimental system. For quantitative analyses, digital simulation has been used in the feedback mode SECM both for conductor and insulator substrates [ 18]. Digital simulation is a powerful technique for analysis of relatively complicated electrochemical systems such as SECM. In the simulation, the space between the probe and the substrate are divided into small volimie elements. Mass transfer based on diffusion for each volume element is calculated. If one chooses suitable initial and boundary conditions, the experimental situation including the mass transfer and various chemical reactions can be simulated. Erom the simulation, one can determine the heterogeneous rate... 

Simple analytical solution of Eqs. 1 and 2 is impossible even for the simplest initial and boundary conditions due to hyperbolic function of the enzymatic rate dependence on substrate concentration [6, 7]. Therefore, the analytical simulation of biosensor action was performed in the simplest cases for which analytical solutions still exist. This approach was used widely, especially at beginning of development of biosensors, to recognize principles of biosensor actimi. The approximal analytical solution gives information about critical cases. They are useful also to test correctness of digital calculations found at initial and boundary limiting conditions. 

This has also been solved for the dropping mercury electrode/ by Koryta (1953). If adsorption is not sufficiently strong to justify the assumption Cq = 0/ then Cq will/ at any instant/ be determined by the adsorption isotherm (Eq. 2.53). This boundary condition leads to mathematical problems the integral equation 2.52 then becomes a Volterra equation. This has been solved for only some very simple isotherms (for example the Henry isotherm/ Eq. 2.42/ Delahay and Trachtenberg 1957) and in most realistic cases, digital simulation is needed. The special problems then have to do with boundary values (F or Q, c ), and will be dealt with in Chapt. 4. There/ the opposite case of control by the rate of the adsorption step itself will also be discussed. 

Britz D, da Silva BM, Avaca LA, Gonzales ER (1990) The Saul yev method of digital simulation under derivative boundary conditions. Anal Chim Acta 239 87-93... 

In this expository article, the basic mathematical model of some simple electrochemical processes was discussed. The model is based on the concept of conservation of charge within the electrolyte. The boundary conditions, on the other hand, are problem-specific. The subject of electrode kinetics is central to the proper specification of the boundary conditions. In their most general form, the conditions are nonlinear, leading to a nonlinear boundary value problem. This is closely tied to the nonlinear polarization curves. The analytical solution of the mathematical model is formidable and for moderately simple two-dimensional regions is impossible to obtain. The only feasible approach is numerical simulation. The use of high-speed digital computers is an essential tool in solving such problems. 

Wang and Sun (2001) developed another numerical method to simulate textile processes and to determine the micro-geometry of textile fabrics. They called it a digital-element model. It models yams by pin-connected digital-rod-element chains. As the element length approaches zero, the chain becomes fully flexible, imitating the physical behavior of the yams. The interactions of adjacent yarns are modeled by contact elements. If the distance between two nodes on different yarns approaches the yam diameter, contact occurs between them. The yarn microstructure inside the fabric is determined by process mechanics, such as yarn tension and interyam friction and compression. The textile process is modeled as a nonlinear solid mechanics problem with boundary displacement (or motion) conditions. This numerical approach was identified as digital-element simulation rather than as finite element simulation because of a special yam discretization process. With the conventional finite element method, the element preserves... 

---