Model semi-implicit method

EX 52 5.2 Solution of the Qreganator model by semi-implicit method M14,M15,M72... 

The pneumatic drying model was solved numerically for the drying processes of sand particles. The numerical procedure includes discretization of the calculation domain into torus-shaped final volumes, and solving the model equations by implementation of the semi-implicit method for pressure-linked equations (SIMPLE) algorithm [16]. The numerical procedure also implemented the Interphase Slip Algorithm (IPSA) of [17] in order to account the various coupling between the phases. The simulation stopped when the moisture content of a particle falls to a predefined value or when the flow reaches the exit of the pneumatic dryer. 

This indicates that after an initial overhead of 0.319 model runs to set up the algorithm, an additional 0.07 of a model-run was required for the computation of the sensitivity coefficients for each additional parameter. This is about 14 times less compared to the one additional model-run required by the standard implementation of the Gauss-Newton method. Obviously these numbers serve only as a guideline however, the computational savings realized through the efficient integration of the sensitivity ODEs are expected to be very significant whenever an implicit or semi-implicit reservoir simulator is involved. 

The same numerical methods as those used to solve the homogeneous reactor models (PFR, BR, and stirred tank reactor) as well as the heterogeneous catalytic packed bed reactor models are used for gas-Uquid reactor problems. For the solution of a countercurrent column reactor, an iterative procedure must be applied in case the initial value solvers are used (Adams-Moulton, BD, explicit, or semi-implicit Runge-Kutta). A better alternative is to solve the problem as a true boundary value problem and to take advantage of a suitable method such as orthogonal collocation. If it is impossible to obtain an analytical solution for the liquid film diffusion Equation 7.52, it can be solved numerically as a boundary value problem. This increases the numerical complexity considerably. For coupled reactions, it is known that no analytical solutions exist for Equation 7.52 and, therefore, the bulk-phase mass balances and Equation 7.52 must be solved numerically. 

A similar equation can be set up for the pressure drop. The combined model now contains K+1 coupled parabolic partial differential equations and one ordinary first order differential equation. They are solved by discretization in the radial direction by use of the orthogonal collocation method, and integration of the resulting set of coupled first order differential equations by use of a semi-implicit Runge-Kutta method. With this model and the used solution method, one can now concentrate on the effective transport properties given in PeH, Pcm and the wall heat transfer coefficient, with the latter being the most important parameter for design. 

Research into the modeling of polymer self-assembly resides in two areas. The first is the formation of accurate and stable numerical methods. Among many studies in the numerical method for the Cahn-Hilliard Equation [5, 6], an elegant approach was reported by David J. Eyre [7] and L. Q. Chen etc [11]. The application of a semi-implicit Fourier-spectral method was extensively studied in this issue. This scheme is shown to be unconditionally gradient stable and solvable for all time steps. The semi-implicit Fourier-spectral method has been employed to model the phase separation as it s more efficient thus allowing us to simulate large systems for a longer time [2, 3]. 

Conceptually, the self-consistent reaction field (SCRF) model is the simplest method for inclusion of environment implicitly in the semi-empirical Hamiltonian24, and has been the subject of several detailed reviews24,25,66. In SCRF calculations, the QM system of interest (solute) is placed into a cavity within a polarizable medium of dielectric constant e (Fig. 2.2). For ease of computation, the cavity is assumed to be spherical and have a radius ro, although expressions similar to those outlined below have been developed for ellipsoidal cavities67. Using ideas from classical electrostatics, we can show that the interaction potential can be expressed as a function of the charge and multipole moments of the solute. For ease... 

The most common approach to solvation studies using an implicit solvent is to add a self-consistent reaction field (SCRF) term to an ab initio (or semi-empirical) calculation. One of the problems with SCRF methods is the number of different possible approaches. Orozco and Luque28 and Colominas et al27 found that 6-31G ab initio calculations with the polarizable continuum model (PCM) method of Miertius, Scrocco, and Tomasi (referred to in these papers as the MST method)45 gave results in reasonable agreement with the MD-FEP results, but the AM1-AMSOL method differed by a number of kJ/mol, and sometimes gave qualitatively wrong results. 

The simple orbital basis expansion method which is used in the implementation of most models of molecular electronic structure consists of expanding each R as a linear combination of determinants of a set of (usually) atom-centred functions of one or two standard forms. In particular most qualitative and semi-quantitative theories restrict the terms in this expansion to consist of the (approximate) occupied atomic orbitals of the constituent atoms of the molecule. There are two types of symmetry constraint implicit in this technique. 

The most expensive part of a simulation of a system with explicit solvent is the computation of the long-range interactions because this scales as Consequently, a model that represents the solvent properties implicitly will considerably reduce the number of degrees of freedom of the system and thus also the computational cost. A variety of implicit water models has been developed for molecular simulations [56-60]. Explicit solvent can be replaced by a dipole-lattice model representation [60] or a continuum Poisson-Boltzmann approach [61], or less accurately, by a generalised Bom (GB) method [62] or semi-empirical model based on solvent accessible surface area [59]. Thermodynamic properties can often be well represented by such models, but dynamic properties suffer from the implicit representation. The molecular nature of the first hydration shell is important for some systems, and consequently, mixed models have been proposed, in which the solute is immersed in an explicit solvent sphere or shell surrounded by an implicit solvent continuum. A boundary potential is added that takes into account the influence of the van der Waals and the electrostatic interactions [63-67]. 

---