Models description

Dynamic models solve the differential material and energy balances for the reactor both before and during relief. A relief system size is assumed, and the model calculates the pressure versus time and temperature versus time histories. Examination of these can then determine whether the maximum accumulated pressure for the reactor would be exceeded with the assumed relief size. Multiple runs are required to find the optimum relief size which yields a maximum pressure which just equals the maximum accumulated pressure (see 5.2.1). 

Many companies have developed their own in-house codes. An increasing number of commercially available models are being developed. These typically cost several thousand to tens of thousands of pounds and licence arrangements vary. 

Some examples of commercially available dynamic models are  

SuperChems Expert 161 is a code developed by Arthur D Little Inc. for risk assessment consequence analysis, which also has a relief system sizing option. The code has a physical properties package that can handle highly non-ideal properties. It can also consider the effect of chemical reaction in the relief system piping. The code uses the DIERS drift flux methods for level swell and has the option of a rigorous two-phase slip model for the. relief system capacity. 

A version of the AspenPlus process flowsheeting package has been developed which allows dynamic simulation of a runaway reaction in order to size the relief system171. 

It will become evident in the mathematical derivation that the feed location in a total reflux column is immaterial. Hence, although the model is developed for 

The other assumption in the model relates to the vapor-liquid equilibrium coef-hcients, or / -values. The /f-values at a given pressure are assumed to be a function of temperature only, and not of composition. It is further assumed that the temperature dependence of the A -values for the different components is similar, that is, the ratio of the /f-values of any pair of components is independent of temperature. Thus, the relative volatilities, defined as the ratios of A -values of any two components, are assumed constant throughout the column. 

The prevailing heat and mass transfer processes in the periodically operated packed beds have been investigated with a pseudo-homogeneous, one-dimensional plug flow model with superimposed axial dispersion. 

The mass and energy conservation equations have been listed in Table 2.1. The constitutive equations for the transport parameters and the mass deposition rate have been summarized in Tables 2.2 and 2.3, respectively. The gas phase (mixture) properties have been computed according to Reid et al. [14], using the pure component data supplied by Daubert and Danner [15]. Uniform initial temperature profiles are taken without any mass deposited 

Standard Newton-Raphson technique for the linearly implicit treatment of the source terms. Moreover, time-step adaptation and local grid refinement procedures have been implemented, making effective use of the WENO smoothness indicators and interpolation polynomials [16], The steep temperature and mass deposition gradients in combination with the strongly non-linear subhmation kinetics require a very efficient and stable numerical implementation using higher order implicit schemes. 

From this scheme it can be seen that the copolymer composition is determined by the values of four monomer reactivity ratios. 

Fukuda et al.lj were the first to recognize that a further two radical reactivity ratios were required to completely define the polymerization kinetics. 

Tlic reactivity ratios taab, / bab. J bua and r aba are sometimes abbreviated to r=,, a ba, rBB and tab or to v., rh rR, rd respectively. The notation used (rAAR, rBAR, rBRA and rAR( ) is preferred since i allows discussion of situations involving more than two monomers. 

In traditional treatments of copolymerizaiion kinetics, the values of the ratios sA and % are implicitly set equal to unity (Section 7.3.1.2.2). Since they contain no terms from cross propagation, these parameters have no direct influence on either the overall copolymer composition or the monomer sequence distribution they only influence the rate of polymerization. 

The instantaneous copolymer composition is described by the following equation (eq. 12)  

The transport properties of the membrane were measured by using a membrane permeation unit as described in our previous papers.32 42 Two gas mixtures were used as the feed gas for the gas permeation tests one consisting of 20% C02, 40% H2, and 40% N2, and the other consisting of 1% CO, 17% C02, 45% H2, and 37% N2 (both on dry basis). The second composition was used to simulate the composition of the synthesis gas from autothermal reforming (ATR) of gasoline with air. Argon was used as the sweep gas for the ease of gas chromatography analysis. Unless otherwise stated, the feed pressure was maintained at about 2.0 atm, while the permeate pressure was set at approximately 1.0 atm. 

As one of the two common types of membrane modules, the hollow-fiber membrane module has shown excellent mass transfer performance due to its large surface area per unit volume (about 1000-3000ft2/ft3 for gas separation). In the modeling work, the WGS membrane reactor was configured to be a hollow-fiber membrane module with catalyst particles packed inside the fibers. 

The catalyst packed was assumed to be the commercial Cu/ZnO catalyst for lower-temperature WGS reaction. A number of studies on the reaction kinetics of the commercial WGS catalyst, Cu0/Zn0/Al203, have been published.43-48 Based on the experimental data of the commercial catalyst (ICI 52-1), Keiski et al.47 suggested two reaction rates for the low-temperature WGS reaction in the temperature range 160-250 °C. The first was dependent only on CO concentration and gave an activation energy of 46.2kJ/mol. The second reaction rate was dependent on CO and steam concentrations with a lower activation energy of 42.6kJ/mol. Because of the proximity of our operation conditions to theirs and the fact that steam is in excess in most of the membrane reactors, Keiski and coworkers first reaction rate expression was chosen for this work. The reaction rate is given in Equation 9.5, 

A uniform and constant Darcy velocity of 15 myear-1 and effective porosity of 0.3 were used along the entire cross-section. The tailings fluid has relatively high chloride concentrations (0.016 mol L 1) compared with the background water has only 0.0007 molL-1. It is widely believed that Cl- acts as a conservative solute in most aquifer systems, and thus its distribution can be used to retrieve dispersivity for the aquifer. By trial and error, a longitudinal dispersivity between 10 to 15 m appears to fit the concentration differences best in monitor wells sampled in September, 1994. It 

A flux (Cauchy or third type) boundary conditions are used for both ends of the onedimensional strip. To represent the reclamation conditions, the incoming fluid has the chemistry of tailings pore fluid for the first five years and of uncontaminated upgradient groundwater thereafter (Table 10.3). 

As discussed in previous sections, coupled reactive transport models generate numerous data and the predicted mass transport is complex. Here, we only give some basic information from this modeling exercise. Interested readers are referred to Zhu and Burden (2001) and Zhu et al. (2001a, 2002) for details. 

Results from this multi-component, multi-species transport model show that the transport of some components are in inter-locked step fashion. For example, the transport of sulfate is closely tied to calcite dissolution and precipitation. 

TITLE Transport Bear Creek example in Phreeqc 2.0 SOLUTION 0 incoming solution units ppm 

To investigate the process of nucleation in a ternary system, we have derived the dependence of Gibbs potential change AG at formation of a new phase nucleus 

To get the expression for AG from Equation 9.57, where the concentration composition of the parent phase and nucleus size enter exphcitly, some approximations such as those given below, have been made. 

A nucleus has a shape of spheroid (ellipsoid having two axes of equal length, making it a surface of revolution). 

The concentration profiles in the region of formation of the new phase nucleus are approximated by the linear functions 

The expression for the change of Ag (Gibbs ptotential per atom of the nucleus), after applying the parallel tangent rule, will be as follows  

The generation of the assay takes places in two stages. The first one is the generation of preliminary TBP and physical properties distribution through 

From the operational point of view, only a minimum set of crude oil bulk properties will be available to generate the updated assay. This is the less favorable case for the application of the model, and minimum required properties must be carefully selected. 

Although crude oil density, sulfur and viscosity are the minimum typical bulk properties required by the model, sometimes optional laboratory analysis may be required to improve the accuracy of the predictions. The methodology is independent of the available information for the crude oil, and partial crude oil assays can be used. The output information is always a complete updated assay. 

Property balances with bulk values are checked and consistency analysis of property values in the fractions is performed until the updated assay is validated. When the validation process is finished, the reliability of the updated assay, matching the quality of the received cargo, will be improved for LP optimization. 

Since the level can vary, a mass balance is relevant and can be written as  

Becanse of the last two assnmptions, the mass balance can be rewritten as  

Additional equations are required to complete the model description. The first one is the relationship between the outlet flow and the pressure, taking assumption four into accoimt  

Since there is only one component present in the tank (boiling pure hquid), there is also a relationship between the pressure in the tank and the temperature of the vapor (and liqnid, which is the same). This relationship can be well described by the law of Clansins-Clapeyron 

As can be seen, Eqn. (15.3) contains a double derivative, it is good practice to change it to a single derivative by combining Eqn. (15.3) with Eqn. (15.2), which results in  

---

The logic that leads us to this last result also limits the applicability of the ensuing derivation. Applying the fraction of total lattice sites vacant to the immediate vicinity of the first segment makes the model descriptive of a relatively concentrated solution. This is somewhat novel in itself, since theories of solutions more commonly assume dilute conditions. More to the point, the model is unrealistic for dilute solutions where the site occupancy within the domain of a dissolved polymer coil is greater than that for the solution as a whole. We shall return to a model more appropriate for dilute solutions below. For now we continue with the case of the more concentrated solution, realizing... 

Perrv, S. G., Bums, D. J., Adams, L. A., Paine, R. J., Dennis, M. G., Mills, M. T., Strimaitis, D. G., Yamartino, R. J., and Insley, E. M., "User s Guide to the Complex Terrain Dispersion Model plus Algorithms for Unstable Conditions (CTDMPLUS)," Vol. I "Model Description and User Instructions," EPA/600/8-89/041, U.S. Environmental Protection Agency, Research Triangle Park, NC, 1989. 

IDLH). For a combustible release, the code gives an estimate of the mass of vapor within the flammable limits. Model Description References Model description and references I... 

Siegel, A.I. et al., Maintenance Personnel Performance Simulation (MAPPS) Model Description of Model Content, Structure and Sensitivity Testing, 1984. 

Weiss, J. M., Morgan, P. H., Lutz, M. W., and Kenakin, T. P. (1996a). The cubic ternary complex receptor-occupancy model. I. Model description. J. Theroet. Biol. 178 151-167. 

Model Description of the Bulk Polymerization of MMA with Chain Transfer to Monomer... 

Heimann, M. and Keeling, C. D. (1989). A three-dimensional model of atmospheric CO2 transport based on observed winds 2. Model description and simulated tracer experiments. In "Aspects of Climate Variability in the Pacific and Western Americas," Geophys. Monogr. Ser. vol. 55 (D. H. Peterson, ed.), pp. 237-275, AGU, Washington, DC. 

Taguchi, S. (1996). A three-dimensional model of atmospheric CO2 transport based on analyzed winds Model description and simulation results for TRANSCOM, /. Geophys. Res. 101, 15099-15109. 

FIGURE 29.15 Flattened model description of the cross-section of a partially fiUed mixing chamber. 

In order to design a zeoHte membrane-based process a good model description of the multicomponent mass transport properties is required. Moreover, this will reduce the amount of practical work required in the development of zeolite membranes and MRs. Concerning intracrystaUine mass transport, a decent continuum approach is available within a Maxwell-Stefan framework for mass transport [98-100]. The well-defined geometry of zeoHtes, however, gives rise to microscopic effects, like specific adsorption sites and nonisotropic diffusion, which become manifested at the macroscale. It remains challenging to incorporate these microscopic effects into a generalized model and to obtain an accurate multicomponent prediction of a real membrane. 

The key attribute of flows in micro devices is their laminar character, which stands in contrast to the mostly turbulent flows in macroscopic process equipment. Owing to this feature, micro flows are a priori much more accessible to a model description than macro flows and can be described by first-principle approaches without any further assumptions. In contrast, for the simulation of turbulent flows usually a number of semi-heuristic models are applied, and in many situations it is not clear which description is most adequate for the problem under investigation. As a result, it stands to reason to assume that a rational design of micro reactors... 

This example is based on the model description of Sec. 3.3.4, and involves a multicomponent, semi-batch system, with both heating and boiling periods. The compositions and boiling point temperatures will change with time. The water phase will accumulate in the boiler. The system simulated is based on a mixture of n-octane and n-decane, which for simplicity will be assumed to be ideal but which has been simulated using detailed activity coefficient relations by Prenosil (1976). 

There are various parameters and assumptions defining radionuclide behavior that are frequently part of model descriptions that require constraints. While these must generally be determined for each particular site, laboratory experiments must also be conducted to further define the range of possibilities and the operation of particular mechanisms. These include the reversibility of adsorption, the relative rates of radionuclide leaching, the rates of irreversible incorporation of sorbed nuclides, and the rates of precipitation when concentrations are above Th or U mineral solubility limits. A key issue is whether the recoil rates of radionuclides can be clearly related to the release rates of Rn the models are most useful for providing precise values for parameters such as retardation factors, and many values rely on a reliable value for the recoil fluxes, and this is always obtained from Rn groundwater activities. These values are only as well constrained as this assumption, which therefore must be bolstered by clearer evidence. 

ICRP (1989) Americium Biokinetics Model Description of the model. 

VII. TRANSCELLULAR DIFFUSION AND METABOLISM A. Biophysical Model Description... 

---