Computed fluxes

Compute flux values and final zone temperatures ... 

Case (a) Compute Flux Density Using Exact Values of the WSGG Constants... 

The computed flux densities are nearly equal because there is asingle sink zone Aj. 

Fig. 32. Computed fluxes of hydrogen atoms in flame of Fig. 25. (a) Convective flux, Afiun (see eqn. (64)) (b) ordinary diffusions flux, j (c) thermal diffusional flux,yjf (d) overall flux, MGyf.

One of the most popular VOF methods is that due to Hirt and Nichols (1981). This method uses an approximate interface reconstruction that forces the interface to align with one of the co-ordinate axis, depending on the prevailing direction of the interface normal. A schematic diagram of reconstruction of a two-dimensional interface is shown in Fig. 7.9. To compute fluxes in a direction parallel to the reconstructed interface, upwind fluxes are used. Fluxes in a direction perpendicular to the reconstructed interface are estimated using a donor-acceptor method. In a donor-acceptor method, a computational cell is identified as a donor of some amount of fluid from one phase and another neighbor cell is identified as the acceptor of that donated amount of fluid. The amount of fluid from one phase that can be convected (donated) across a cell boundary is limited by the minimum of the filled volume in the donor cell or the free volume in the acceptor cell. This minimizes numerical diffusion at the interface. 

The computed flux densities are nearly equal because there is asingle sink zone Ai-(This example was developed as a MATHCAD 14 worksheet. Mathcad is a registered trademark of Parametric Technology Corporation.)... 

Subsequent field studies by Rutler et al. (11) showed extensive negative saturation anomalies for methylchloroform in surface waters of the Pacific Ocean. CFC-11, an unreactive gas that also dissolves in the upper layers of the sea, did not exhibit such anomalies. The methylchloroform data were used to compute fluxes of methylchloroform for various locations. As shown in Figure 4, the methylchloroform fluxes were variable but almost uniformly negative with a mean flux of -3.1 nmol/m2 per day, close to the computed value for hydrolysis (67). Figure 4, however, shows that the negative fluxes were much greater than the mean at a number of locations. 

Figure 4. Computed fluxes of CH3CCl3 in the Pacific Ocean as a function of latitude. Butler et al. reported that the most negative fluxes corresponded to upwelling regions. (Reproduced with permission from reference 11. Copyright 1991 American Geophysical Union.)...

The method of solution adopted in this case was to search for the value of the mass flux of methane that makes co= 0. Thus, instead of computing fluxes from the rate equations knowing the mass fractions in the bulk and at the interface, we solve the linear system... 

We calculated the flux vectors Jjfe/) = 0,(z.OVi(z,/), where v is the mean flow velocity at time 1 and depth z, and the index i denote macropore or micropore, respectively. Figures 4-6 and 4 7 show the computed flux fields for one of the numerical experiments. During rainlall, infiltration is dominated by the macropore. I ateral infiltration into the matrix occurs as water advances within the macropore (see the direction ol vectors in Fig. 4 6). As soon as the nearby walls are saturated. 

When the computed fluxes are inserted in equation (B.40), the roots (m ) and the discretization coefficients (Ajj) are computed by JCOBI and DFORP by considering the surface as the interpolation point where the concentrations at the surface are known and the centre is excluded. The roots and discretization coefficients (m, Aj ) for determination of are of course different from ( , Aj,) used for the determination of the concentration This is because the flux is specified at one boundary and the concentration at the other. The collocation points (Uj) for A are shifted one position from the collocation points Uj of the concentration profile, i.e. +1 = Uj. This is illustrated by a simple diagram in Figure B.l. 

Steady-state numerical simulations of fluid flow and cupric ion transport within an electrochemical fountain plating system are presented. Specifically, the diffusion-limit is determined directly from the computed flux of cupric ions to the wafer under the assumption of complete surface consumption. This maximum flux, in turn, determines the maximum ionic current that can be passed through the electrolyte to the wafer, which is called the limiting current. The goal of the present study is to predict variations in the limiting current density for different electrolyte volumetric flow rates and wafer (cathode) rotation rates. The efficacy of different computational models, including one-dimensional, two-dimensional axisymmetric, and three-dimensional approximations, are assessed via comparisons of numerical predictions with experimental data. 

Eddy correlation measurements require fast-response instrumentation to resolve the turbulent fluctuations that contribute primarily to the vertical flux. These requirements are particularly severe under stable conditions where response times on the order of 0.2 s or less may be required. In practice, it is often possible to use somewhat slower instruments and apply various corrections to the computed fluxes as compensation. The eddy correlation technique has been used in aircraft (Pearson and Steadman 1980 Lenschow et al. 1982) as well as with tower-mounted instruments. 

Figure 4 An approximately 6 year history of surface water He isotope ratio anomalies (A) and computed flux to the atmosphere near Bermuda (B).

The general solution of the differential equation (25) is a linear combination of the linearly independent solutions, where the constants of combination are determined by the initial conditions. In the special case considered below, from three to five terms of the asymptotic expansions in (26) and (27) are needed to compute fluxes to an accuracy of four decimal places. 

It follows from Table II that, for a fixed value of the parameter e, the quantity P (0) = P (x = 0 e, A(), which is the computed normalized flux for x = 0, converges to the quantity P(0 s) as N increases. However, the computed flux P (0) does not converge e-uniformly. In fact, for any given value of N one can find a value of the parameter e = e(N) such that the error Q(s, N) is no less than some positive constant, for example, e N) = const A ". ... 

Thus, in the case of the boundary value problem (1.16), (1.24), the use of the scheme (1.17), (1.18), (1.25) leads, for small values of the parameter s, to sharp underestimation of the computed normalized flux (and also of the solution gradient) on the boundary. Therefore, even qualitatively, the normalized flux cannot be approximated by the computed flux e-uniformly. 

It follows from the results given in Tables XXI-XXIII that the largest error in the computed flux arises from the approximation error for the singular part of the problem solution. We see that, as N and IVq increase, the computed diffusion flux P x, t )) approximates the real flux... 

P x, t M(3 7)( )) for a fixed value of the parameter. However, the computed flux does not converge c-uniformly. Nevertheless, the computed flux qualitatively well approximates the actual flux e-uniformly. 

A useful measure of the flux distortion is provided by the concept of the effective neutron temperature Tn- We define this temperature as that number which when used in (4.197) gives the best least-squares fit of a Maxwell-Boltzmann distribution to the computed flux (solid-line curve) in the range 0 < x < 35. The ratio of the effective neutron temperature to the moderator temperature Tn/Ts is indicated in the figure for the first two cases. The ratio has been omitted from the last case (A = 9, = 2 ) because the flux was so severely distorted from a... 

This concept of alternating reflective layers with nonconducting fibrous spacers has been known since 1951. This arrangement results in heat impedance excelling by far the performance of the best dewar vessels by reducing the apparent k factor to below 0.5 xw/cm- C. In spite of this spectacular reduction of the heat flux, the reduction in heat flow is considerably lower than what one would expect from the Stefan-Boltzmann equation for the applied number of shields and their emissivity. The ratio of observed heat flow to the theoretically computed flux can be referred to as the efficiency of the radiation shields (see Table 1). 

---