Unsteady Simulation

Local time scales can vary significantly over the computational domain making unsteady simulations expensive. However, local time stepping can be used to reduce drastically convergence for stationary flows. 

The list is merely suggestive. Complexity of reactive flows may greatly expand the list of issues on which further research is required. Another area which deserves mention here is modeling of inherently unsteady flows. Most flows in engineering equipment are unsteady (gas-liquid flow in a bubble column reactor, gas-solid flow in a riser reactor and so on). However, for most engineering purposes, all the details of these unsteady flows are not required to be known. Further work is necessary to evolve adequate representation of such flows within the CFD framework without resorting to full, unsteady simulations. This development is especially necessary to simulate inherently unsteady flows in large industrial reactors where full, unsteady simulations may require unaffordable resources (and therefore, may not be cost effective). Different reactor types and different classes of multiphase flows will have different research requirements based on current and future applications under consideration. 

Keywords hydrate, unsteady simulation, natural gas network, Tehran city... 

Spreadsheet Applications. The types of appHcations handled with spreadsheets are a microcosm of the types of problems and situations handled with fuU-blown appHcation programs that are mn on microcomputers, minis, and mainframes and include engineering computations, process simulation, equipment design and rating, process optimization, reactor kinetics—design, cost estimation, feedback control, data analysis, and unsteady-state simulation (eg, batch distillation optimization). 

Dawes, W., 1995, A Simulation of the Unsteady Interaction of a Centrifugal Impeller with its Vaned Diffuser Flows Analysis, ASME Journal of Turbomachinery, Vol. 117, pp. 213-222. 

As outlined earlier, in multizone models, perfect mixing is assumed in the individual zone. The spatial distribution of velocities, contaminant concentrations, and air temperatures in a zone can be determined only by using CFD. On the other hand, wind effects are easily accounted for in multizone models, and unsteady-state simulation is normally performed. On the combined use of the two methods, see Schaelin et al.--... 

The ODEs governing the unsteady CSTR are obtained by adding accumulation terms to Equations (4.1). The simulation holds the volume constant, and... 

The simulation example DRY is based directly on the above treatment, whereas ENZDYN models the case of unsteady-state diffusion, when combined with chemical reaction. Unsteady-state heat conduction can be treated in an exactly analogous manner, though for cases of complex geometry, with multiple heat sources and sinks, the reader is referred to specialist texts, such as Carslaw and Jaeger (1959). 

Equation based programs in which the entire process is described by a set of differential equations, and the equations solved simultaneously not stepwise, as in the sequential approach. Equation based programs can simulate the unsteady-state operation of processes and equipment. 

The principal advantage of equation based, dynamic, simulators is their ability to model the unsteady-state conditions that occur at start-up and during fault conditions. Dynamic simulators are being increasingly used for safety studies and in the design of control systems. 

Most plant simulations have been steady-state simulations. This is to be expected, since just as a baby must learn to crawl before he can walk, so the simpler steady-state problems must be solved before the unsteady-state ones can be tackled. However, unsteady- state plant simulations are being attempted, and undoubtedly sometime in the future this will be a common tool for plant designers. 

It is noted that the virtual body force Fp depends not only on the unsteady fluid velocity, but also on the velocity and location of the particle surface, which is also a function of time. There are several ways to specify this boundary force, such as the feedback forcing scheme (Goldstein et al., 1993) and direct forcing scheme (Fadlun et al., 2000). In 3-D simulation, the direct forcing scheme can give higher stability and efficiency of calculation. In this scheme, the discretized momentum equation for the computational volume on the boundary is given as... 

Fig. 3a indicates that the bubble-rise velocity measured based on the displacement of the top surface of the bubble ( C/bt) quickly increases and approaches the terminal bubble rise velocity in 0.02 s. The small fluctuation of Ubt is caused by numerical instability. The bubble-rise velocity measured based on the displacement of the bottom surface of the bubble (Ubb) fluctuates significantly with time initially and converges to Ubt after 0.25 s. The overshooting of Ubb can reach 45-50 cm/s in Fig. 3a. The fluctuation of Ubb reflects the unsteady oscillation of the bubble due to the wake flow and shedding at the base of the bubble. Although the relative deviation between the simulation results of the 40 X 40 x 80 mesh and 100 x 100 x 200 mesh is notable, the deviation is insignificant between the results of the 80 x 80 x 160 mesh and those of the 100 X 100 x 200 mesh. The agreement with experiments at all resolutions is generally reasonable, although the simulated terminal bubble rise velocities ( 20 cm/s) are slightly lower than the experimental results (21 25 cm/s). A lower bubble-rise velocity obtained from the simulation is expected due to the no-slip condition imposed at the gas-liquid interface, and the finite thickness for the gas-liquid interface employed in the computational scheme.

An inherent property of the LES approach is that the simulated flow field is no longer steady, but exhibits a transient character due to the presence and motion of large-scale eddies. The LES methodology has proven to be a powerful tool for studying and visualizing stirred tank flows (Eggels, 1996 Derksen et al. 1999 Bakker et al., 2000 Derksen, 2001 Bakker and Oshinowo, 2004), as it inherently takes the unsteady and periodic behavior of the flow (around impeller and baffles) into account. 

Another comparison is due to Van Wageningen et al. (2004) who performed a similar study (in terms of the numerical scheme used) on unsteady laminar flow in a Kenics static mixer. They found that the LB code was 500-600 times faster than FLUENT in terms of simulation time per grid node per time step and that FLUENT used about 5 times more memory than LB. 

In those early days, when computer power was limited, often use was made of a symmetry assumption each quarter of the vessel containing one of the four baffles at the vessel wall was supposed to behave identically hence, a steady flow in the RANS approach was simulated in just a quarter vessel. Such strong simplifications are no longer in use. Precessing vortices moving around the vessel centerline contribute to flow unsteadiness and, therefore, exclude models that just assume flow steadiness or allow for domain reductions through geometrical symmetries. The most correct response to this flow unsteadiness is the concept of LES. 

Large eddy simulations explicitly resolves the inherently unsteady character of the turbulent flow in a stirred tank into account, including the periodic phenomena associated with the motion of the impeller and their interaction with... 

The chemical species were treated as passive scalar tracers in the unsteady LBM equations. The reaction was simulated as being mass-transfer limited at low Re — 166, with diffusivities in the ratios DA DB Dc— 1 3 2. The concentration fields shown in Fig. 16 are different for each species due to the different diffusivities. The slow-diffusing species A is transported mainly by convection and regions of high or low concentration correspond to features of the flow field. A more uniform field is seen for the concentration of faster... 

The research on the flow regimes in packed tubes suggests that laminar flow CFD simulations should be reasonable for Re <100 approximately, and turbulent simulations for Re >600, also approximately. Just as RANS models provide steady solutions that are regarded as time averages of the real time-dependent turbulent flow, it may be suggested that CFD simulations in the unsteady laminar inertial range 100 time-averaged picture of the flow field. As with wall functions, comparisons with experimental data and an improved assessment of what information is really needed from the simulations will inform us as to how to proceed in these areas. 

Example 14-7 can also be solved using the E-Z Solve software (file exl4-7.msp). In this simulation, the problem is solved using design equation 2.3-3, which includes the transient (accumulation) term in a CSTR. Thus, it is possible to explore the effect of cAo on transient behavior, and on the ultimate steady-state solution. To examine the stability of each steady-state, solution of the differential equation may be attempted using each of the three steady-state conditions determined above. Normally, if the unsteady-state design equation is used, only stable steady-states can be identified, and unstable... 

In the MC simulation, these equations are treated numerically using fractional time stepping. For unsteady flow, the simulation time step At is determined from the FV code by... 

At the end of the chemical-reaction step, all particle properties (w n>, X(n), fl(n>) have been advanced in time to t + At. Particle-field estimates of desired outputs can now be constructed, and the MC simulation is ready to perform the next time step. For a constant-density flow, the particle-field estimates are not used in the FV code. Thus, for stationary flow, the particle properties can be advanced without returning to the FV code. For unsteady or variable-density flow, the FV code will be called first to advance the turbulence fields before calling the PDF code (see Fig. 7.3). 

Multiple-scale perturbation analysis and numerical simulation of the unsteady-state IEM model. Chemical Engineering Science 45, 2857-2876. 

---