Discretisation boundary conditions

Computational fluid dynamics (CFD) is essentially a computer-based numerical analysis approach for fluid flow, heat transfer and related phenomena. CFD techniques typically consist of the following five subprocesses geometrical modelling, geometry discretisation, boundary condition definition, CFD-based problem solving, and post-processing for solution visualisation. 

In terms of computational implementation, the container that stores the concentration grid may still be rectangular one simply sets the initial concentration of all species at all points in this exclusion zone to zero. In addition, any coefficients for the Thomas algorithm (ogj, fij) that refer to spatial points inside this zone are also set to zero, and the discretised boundary conditions derived from (10.42) are applied in the appropriate places. The current is calculated in exactly the same manner as for a microdisc electrode. 

Irrespective of the type of extended system we are interested in we impose periodic boundary conditions in position space - "the large period" BK. Such conditions imply a discretisation of momentum and reciprocal space 27] which means that integrations are replaced by summations ... 

The coefficients are used to compute the derivatives, but can also be useful in the discretisation of derivative boundary conditions or in the setting up of discretisation matrices in some problems. 

In some cases, for example electrochemical pdes with derivative boundary conditions, the discretisation process for both the pde and the boundary conditions leads to a mix of a differential equation system and one or more plain algebraic equations. They might be, for example, equations of the form... 

The two methods are BDF and extrapolation. Both methods are used for the numerical solution of odes and are described in Chap. 4. The extension to the solution of pdes is most easily understood if the pde is semidiscretised that is, if we only discretise the right-hand side of the diffusion equation, thus producing a set of odes. This is the Method of Lines or MOL. Once we have such a set, as seen in (8.9), the methods for systems of odes can be applied, after adding boundary conditions. 

This subroutine is a general routine for computing the first or second derivative on a number n of points, referred to the ith one among that number, on an arbitrarily spaced grid of points. The derivatives are computed as a linear sum of terms, and the coefficients in that sum are also passed back, for use, for example, in the discretisation of boundary conditions or the spatial second derivative. The number of points is in principle unrestricted, but the routine will fail for values n > 12, where the accuracy abruptly drops. A value, in any case, exceeding about 8, is perhaps impractical. This routine can be used instead of the algebraic expressions shown in Chap. 7, or if n values greater than 4 are required. 

In the finite difference method, the derivatives dfl/dt, dfl/dx and d2d/dx2, which appear in the heat conduction equation and the boundary conditions, are replaced by difference quotients. This discretisation transforms the differential equation into a finite difference equation whose solution approximates the solution of the differential equation at discrete points which form a grid in space and time. A reduction in the mesh size increases the number of grid points and therefore the accuracy of the approximation, although this does of course increase the computation demands. Applying a finite difference method one has therefore to make a compromise between accuracy and computation time. 

The derivative is formed in the outward normal direction, and a and can be dependent on time. For the discretisation of (2.253) it is most convenient if the boundary coincides with a grid line, Fig. 2.45, as the boundary temperature which appears in (2.253) can immediately be used in the difference formula. The replacement of the derivative d d/dn by the central difference quotient requires grid points outside the body, namely the temperatures i9k or ()k,, which, in conjunction with the boundary condition, can be eliminated from the difference equations. 

The integration of the system 12.1 needs appropriate Initial and boundary conditions (see also Chapter 4). For simulation purposes, the spatial coordinates in distributed parameter models must be discretised, for example, by finite differences, finite elements, or other techniques. 

The Galerkin weighted residual method is employed to formulate the finite element discretisation. An implicit mid-interval backward difference algorithm is implemented to achieve temporal discretisation. With appropriate initial and boundary conditions the set of non-linear coupled governing differential equations can be solved. 

Discretisation and boundary conditions were taken from the simulations of Experiment 1 with the flow rates and influx concentrations being adapted to the new situation. The input concentrations for toluene and o-xylene for the 109 days of column operation are displayed in Figme 16.6. The only observations in Experiment 2 besides the input values were the o-xylene concentrations at the column outlet. 

For the one-dimensional system under consideration, the state of the system (which is represented in a discretised form by the set of concentrations Co,..., C i) at any moment in time depends on two things the spatial boundary conditions (the electrode surface and bulk solution boundaries), and the state of the system at the previous moment in time. We may discretise time in the same manner as space, with a constant interval between adjacent timesteps, AT. When modelling, we are interested in the state of the system from the starting time of the experiment, To, to its finishing time, Tmax, so the total number of timesteps, m, is given by... 

In this chapter, the case of first-order mechanisms will be considered, analysing the changes in the corresponding differential equations, boundary conditions and problem-solving methodology. Strategies for the resolution of the linear equation systems resulting from discretisation will be detailed for some representative mechanisms. The case of multiple-electron transfers will also be considered. This will enable us to simulate other common situations where the electroactive molecule transfers more than one electron. 

A more general and rigorous approach is the determination of the concentrations by employing a multidimensional root-finding algorithm. Indeed, the unknown concentrations can be viewed as the roots of a set of non-linear simultaneous equations to be zeroed, which are given by the discretised differential equations and boundary conditions ... 

Equilibrium equation (2) is used as reaction balance for H2S. After discretisation of the set of partial differential equations (6)-(10) and linearization by the Newton-Raphson technique the system is solved by an iterative numerical procedure. Boundary conditions and more detailed information are given elsewhere [4]. 

The first term on the right hand side represents heat transfer due to conduction and second term represents the heat released due to heterogeneous reactions within the electrodes which vanishes in the case of cathode. Two source terms, the radiative heat source term Qr, and the convective heat source term enters Eq. 4.48 as boundary condition at the interface between electrode and the flow channel, and the electrochemical heat source term enters as boundary condition at the interface between the anode and the electrolyte. The radiative heat transfer between the interconnect and the outer most discretised cell in the porous electrode is given by... 

Coelho et al. (1997) avoided an exphcit treatment of isolated aggregates by the variation of box size and subsequent data regression because this case requires a discretisation of large fluid domains or artificial boundary conditions when the... 

The equations for the catalyst bed combined with the boundary conditions can be solved using modem methods such as those available in CFD. A fast method uses discretisation in the radial direction by the orthogonal collocation method by Villadsen and Michelsen, where the dependent variable is only considered in the so-called collocation points [512]. This method reduces each partial differential equation to a corresponding number of coupled ordinary first-order differential equations with the axial distance as independent variable. The selected number of collocation points depends on the stiffiiess of the actual problem. Only with very steep gradients are more than five collocation points necessary to predict the extrapolated properties close to the reactor wall with sufficient accuracy to predict the risk of carbon formation. 

Secondly, the general balances in both phases and at the interface leads to a system of differential and algebraic equations with boundary conditions at each end. So a discretisation method is used and the resulting algebraic system is solved by a traditional Newton s method. This general balances use values resulting from the integration of the equations in the diffusional layer. 

---