Numerical modelling computer simulation

The time dependent solvation funetion S(t) is a directly observed quantity as well as a convenient tool for numerical simulation studies. The corresponding linear response approximation C(t) is also easily eomputed from numerical simulations, and can also be studied using suitable theoretical models. Computer simulations are very valuable both in exploring the validity of such theoretical calculations, as well as the validity of linear response theory itself (by comparing S(t) to C(t)). Furthermore they can be used for direct visualization of the solute and solvent motions that dominate the solvation process. Many such simulations were published in the past decade, using different models for solvents such as water, alcohols and acetonitrile. Two remarkable outcomes of these studies are first, the close qualitative similarity between the time evolution of solvation in different simple solvents, and second, the marked deviation from the simple exponential relaxation predicted by the Debye relaxation model (cf Eq. [4.3.18]). At least two distinct relaxation modes are... 

In this chapter we consider the derivation of unsteady-state models of chemical processes from physical and chemical principles. Unsteady-state models are also referred to as dynamic models. We first consider the rationale for dynamic models and then present a general strategy for deriving them from first principles such as conservation laws. Then dynamic models are developed for several representative processes. Finally, we describe how dynamic models that consist of sets of ordinary differential equations and algebraic relations can be solved numerically using computer simulation. 

A numerical model for simulation of the global heat transfer and the melt flow in the Czochralski growth of large silicon crystals is presented. The key model features are an extended 3D domain for the 2D/3D computations and a hybrid LES/RANS approach to turbulence modeling. It is shown that use of parallel computations on affordable multiprocessor systems assembled from the COTS hardware could reduce the turn-around time of simulation by an order of magnitude. The model validation using the experimental data on the growth of 100-mm and 300-mm silicon crystals in the industrial pullers Ekz-1300 and Ekz-2405 has proved its predictive power. 

The complexity of polymeric systems make tire development of an analytical model to predict tlieir stmctural and dynamical properties difficult. Therefore, numerical computer simulations of polymers are widely used to bridge tire gap between tire tlieoretical concepts and the experimental results. Computer simulations can also help tire prediction of material properties and provide detailed insights into tire behaviour of polymer systems. A simulation is based on two elements a more or less detailed model of tire polymer and a related force field which allows tire calculation of tire energy and tire motion of tire system using molecular mechanisms, molecular dynamics, or Monte Carlo teclmiques 1631. 

Computational fluid dynamics (CFD) is the numerical analysis of systems involving transport processes and solution by computer simulation. An early application of CFD (FLUENT) to predict flow within cooling crystallizers was made by Brown and Boysan (1987). Elementary equations that describe the conservation of mass, momentum and energy for fluid flow or heat transfer are solved for a number of sub regions of the flow field (Versteeg and Malalase-kera, 1995). Various commercial concerns provide ready-to-use CFD codes to perform this task and usually offer a choice of solution methods, model equations (for example turbulence models of turbulent flow) and visualization tools, as reviewed by Zauner (1999) below. 

To illustrate the relationship between the microscopic structure and experimentally accessible information, we compute pseudo-experimental solvation-force curves F h)/R [see Eq. (22)] as they would be determined in SEA experiments from computer-simulation data for T z [see Eqs. (93), (94), (97)]. Numerical values indicated by an asterisk are given in the customary dimensionless (i.e., reduced) units (see [33,75,78] for definitions in various model systems). Results are correlated with the microscopic structure of a thin film confined between plane parallel substrates separated by a distance = h. Here the focus is specifically on a simple fluid in which the interaction between a pair of film molecules is governed by the Lennard-Jones (12,6) potential [33,58,59,77,79-84]. A confined simple fluid serves as a suitable model for approximately spherical OMCTS molecules confined... 

In the light of the above questions, it is tempting to refer to the results emerging from numerous theoretical and computer simulation studies [40,41,85-88,129-131] of the random field Ising model, and we shall do so, but only after completing the present discussion. 

But a computer simulation is more than a few clever data structures. We need algorithms to manipulate our system. In some way, we have to invent ways to let the big computer in our hands do things with the model that is useful for our needs. There are a number of ways for such a time evolution of the system the most prominent is the Monte Carlo procedure that follows an appropriate random path through configuration space in order to investigate equilibrium properties. Then there is molecular dynamics, which follows classical mechanical trajectories. There is a variety of dissipative dynamical methods, such as Brownian dynamics. All these techniques operate on the fundamental degrees of freedom of what we define to be our model. This is the common feature of computer simulations as opposed to other numerical approaches. 

A complete set of intermolecular potential functions has been developed for use in computer simulations of proteins in their native environment. Parameters have been reported for 25 peptide residues as well as the common neutral and charged terminal groups. The potential functions have the simple Coulomb plus Lennard-Jones form and are compatible with the widely used models for water, TIP4P, TIP3P and SPC. The parameters were obtained and tested primarily in conjunction with Monte Carlo statistical mechanics simulations of 36 pure organic liquids and numerous aqueous solutions of organic ions representative of subunits in the side chains and backbones of proteins... 

The numerical accuracy of simulations performed using this model is affected by several factors. These include a) the degree of triangulation, b) the number of marching steps taken along the flow direction and c) the order of the polynomial basis function. Numerical accuracy improves as a, b and c increase, however the computational time can become excessive. Therefore, it was necessary to quantitatively determine the effects of these variables on numerical accuracy. 

Beyond the clusters, to microscopically model a reaction in solution, we need to include a very big number of solvent molecules in the system to represent the bulk. The problem stems from the fact that it is computationally impossible, with our current capabilities, to locate the transition state structure of the reaction on the complete quantum mechanical potential energy hypersurface, if all the degrees of freedom are explicitly included. Moreover, the effect of thermal statistical averaging should be incorporated. Then, classical mechanical computer simulation techniques (Monte Carlo or Molecular Dynamics) appear to be the most suitable procedures to attack the above problems. In short, and applied to the computer simulation of chemical reactions in solution, the Monte Carlo [18-21] technique is a numerical method in the frame of the classical Statistical Mechanics, which allows to generate a set of system configurations... 

There have been books on droplet-related processes. However, the present book is probably the first one that encompasses the fundamental phenomena, principles and processes of discrete droplets of both normal liquids and melts. The author has attempted to correlate many diverse mechanisms and effects in a single and common framework in an effort to provide the reader with a new perspective of the identical basic physics and the inherent relationship between normal liquid and melt droplet processes. Another distinct and unique feature of this book is the comprehensive review of the empirical correlations, analytical and numerical models and computer simulations of droplet processes. These not only provide practical and handy approaches for engineering calculations, analyses and designs, but also form a useful basis for future in-depth research. Therefore, the present book covers the fundamental aspects of engineering applications and scientific research in the area. 

---