Numerical Issues

There are many numerical issues that we must discuss before proceeding to applications of FDM when solving polymer processing problems. The most important are, 

As we can imagine, most of these issues are directly related to the order of the approximation used in the finite difference representation. In fact, the truncation error (as shown in Chapter 7) is the difference between the PDE and the FD representation, which is represented by the terms collapsed in 0(Axn). For problems represented by PDEs with more than one independent variable the truncation error will be the sum of the truncation error for each FD representation. For example, for a transient one dimensional PDE, where we use a first order approximation for the time derivative and a second order for the spatial derivative, we will have a truncation error that is O(At) + 0 Ax2), which can also be written as 0(At, Ax2). 

The consistency of a finite difference approximation is the behavior of this representation when the mesh is refined. In a one dimensional case, for example, the mesh will indicate the value for Ax, which, as we discussed above, dictates the value of the truncation error. Thus, a finite difference representation of a PDE is said to be consistent if the truncation error goes to zero as the grid size (or Ax) goes to zero. 

The stability of a FD representation deals with the behavior of the truncation error as the calculation proceeds in time or marches in space, typically, transient problems, and problems with convection-convection derivatives. A stable FD scheme will not allow the errors to grow as the solution proceeds in time or space. The issue of stability for transient problems will be analyzed in depth later in this chapter. 

Unlike continuous distillation, batch distillation is inherently an unsteady state process. Dynamics in continuous distillation are usually in the form of relatively small upsets from steady state operation, whereas in batch distillation individual species can completely disappear from the column, first from the reboiler (in the case of CBD columns) and then from the entire column. Therefore the model describing a batch column is always dynamic in nature and results in a system of Ordinary Differential Equations (ODEs) or a coupled system of Differential and Algebraic Equations (DAEs) (model types III, IV and V). 

---

As alluded to previously, numerical issues actually create a more complex situation than that just been described. For starters, the density of states is almost never calculated directly, as it typically spans many orders of magnitudes. This, in turn, would quickly overwhelm standard double-precision calculations in personal computers. This is easily remedied by working instead with the dimensionless entropy / = In Q, which for the purposes of this chapter will inherit all of the same notation used for the density of states in Chap. 1 - subscripts tot, ex, etc. 

Despite the many advances, the users of CFD codes must keep in mind that the underlying transport equations are based on models, which may or may not be valid for a particular application. This fact often escapes the minds of newcomers to the field who are typically overwhelmed by the numerical issues associated with convergence, grid-independence, and post-processing. Even... 

This book presents the current state of the art in computational models for turbulent reacting flows, and analyzes carefully the strengths and weaknesses of the various techniques described. The focus is on formulation of practical models as opposed to numerical issues arising from their solution. 

The PDF codes presented in this chapter can be (and have been) extended to include additional random variables. The most obvious extensions are to include the turbulence frequency, the scalar dissipation rate, or velocity acceleration. However, transported PDF methods can also be applied to treat multi-phase flows such as gas-solid turbulent transport. Regardless of the flow under consideration, the numerical issues involved in the accurate treatment of particle convection and coupling with the FV code are essentially identical to those outlined in this chapter. For non-orthogonal grids, the accurate implementation of the particle-convection algorithm is even more critical in determining the success of the PDF simulation. 

The use of different electromigration-based separation techniques hyphenated with MS has become a standard technique in modem pharmaceutical analysis. With the possibility of several commercially available instraments, including various interfaces, the coupling between CE and MS can now be easily achieved. Therefore, numerous issues can be resolved according to the wide choice of operation techniques afforded in CE and the various ionization modes and/or analyzers. CE—MS has emerged as a good alternative for trace... 

The thickness of pharmaceutical tablet coatings was predicted using target factor analysis (TFA) applied to Raman spectra collected with a 532-mn laser, where the samples were photobleached in a controlled manner before spectra were acquired. The authors acknowledge numerous issues that limit the direct applicability of this approach to process control. These include potential damage or alteration of the samples from photobleaching, laser wavelength selection, and data acquisition time. However, most of the issues raised relate to the hardware selected for a particular implementation and do not diminish the demonstration [286]. 

There are numerous issues relating to science and society. In the Spotlight is a special feature of Conceptual Chemistry that highlights these issues and prods you to consider implications. These features can also sewe as a centerpiece foT potentially heated discussions with youT classmates and others. To inform youTself further on the issues, be sure to explore the related web sites listed on the last page of each preceding chapter. 

Complications linked to the farm controls are both in the frame and core system, the crucial factor for the development of organic agriculture. For example, in order to import tropical organic products into the European markets, smallholders have to consider numerous issues concerning the compliance with the EEC 2092/91 ... 

There are numerous issues that influence the current and future state of flame-retardant plastics. This is particularly evident in Europe, where regulatory and... 

When choosing a laboratory reactor for the measurement of reaction rate data, numerous issues must be resolved. The choice of the reactor is based on the characteristics of the reaction and for all practical matters by the availability of resources (i.e., reactors, analytical equipment, money, etc.). A good example of the issues involved in selecting a laboratory reactor and how they influence the ultimate choice is presented by Weekman [AIChE J., 20 (1974) 833]. Methods for obtaining reaction rate data from laboratory reactors that approximate the ideal reactors listed in Table 3.5.1 are now discussed. 

Imparato and Peliti " show how the work distribution for a given protocol can be obtained using a joint probability distribution for the work and the initial state, and apply it to system represented by a coordinate in a mean field. They show that it leads to the JE, and demonstrate some numerical issues when applied to large systems where it becomes difficult to observe fluctuations. They also provide an analytical treatment of a driven Brownian particle, obtaining the probability distribution for the work and FR for this system. ... 

With the field of hydrogen-bonded electron transfer in its infancy, new frontiers await exploration and numerous issues remain to be addressed. The D—[H]—A supramolecules of the preceding sections define two basic PCET reaction mechanisms ... 

Drug therapy involves considerations beyond the strictly scientific pharmacological aspects of medicines.These include numerous issues relating to prescribers themselves and to patients. 

Various issues in the development of a flow model and its numerical simulation have been already discussed in the previous section. It will be useful to make a few comments on the validation of the simulated results and their use in reactor engineering. More details are discussed in Part III and Part IV. Even before validation, it is necessary to carry out a systematic error analysis of the generated computer simulations. The influence of numerical issues on the predicted results and errors in integral balances must be checked to ensure that they are within the acceptable tolerances. The simulated results must be examined and analyzed using the available post-processing tools. The results must be checked to verify whether the model has captured the major qualitative features of the flow such as shear layers and trailing vortices. 

The recent progress in experimental techniques and applications of DNS and LES for turbulent multiphase flows may lead to new insights necessary to develop better computational models to simulate dispersed multiphase flows with wide particle size distribution in turbulent regimes. Until then, the simulations of such complex turbulent multiphase flow processes have to be accompanied by careful validation (to assess errors due to modeling) and error estimation (due to numerical issues) exercise. Applications of these models to simulate multiphase stirred reactors, bubble column reactors and fluidized bed reactors, are discussed in Part IV of this book. 

In this approach, the finite volume methods discussed in the previous chapter can be applied to simulate the continuous fluid (in a Eulerian framework). Various algorithms for treating pressure-velocity coupling, and the discussion on other numerical issues like discretization schemes are applicable. The usual interpolation practices (discussed in the previous chapter) can be used. When solving equations of motion for a continuous fluid in the presence of the dispersed phase, the major differences will be (1) consideration of phase volume fraction in calculation of convective and diffusive terms, and (2) calculation of additional source terms due to the presence of dispersed phase particles. For the calculation of phase volume fraction and additional source terms due to dispersed phase particles, it is necessary to calculate trajectories of the dispersed phase particles, in addition to solving the equations of motion of the continuous phase. 

By 1996, a number of European countries and organizations started to pay much more attention to this issue than before. The United States experience offered a number of lessons for the Europeans to consider as they discussed and considered numerous issues. 

---