Mathematical modelling and numerical simulation of phosphine flow during grain fumigation in leaky cylindrical silos

The phosphine distribution in a cylindrical silo containing grain is predicted. A three-dimensional mathematical model, which accounts for multicomponent gas phase transport and the sorption of phosphine into the grain kernel is developed. In addition, a simple model is presented to describe the death of insects within the grain as a function of their exposure to phosphine gas. The proposed model is solved using the commercially available computational fluid dynamics (CFD) software, FLUENT, together with our own C code to customize the solver in order to incorporate the models for sorption and insect extinction. Two types of fumigation delivery are studied, namely, fan- forced from the base of the silo and tablet from the top of the silo. An analysis of the predicted phosphine distribution shows that during fan forced fumigation, the position of the leaky area is very important to the development of the gas flow field and the phosphine distribution in the silo. If the leak is in the lower section of the silo, insects that exist near the top of the silo may not be eradicated. However, the position of a leak does not affect phosphine distribution during tablet fumigation. For such fumigation in a typical silo configuration, phosphine concentrations remain low near the base of the silo. Furthermore, we find that half-life pressure test readings are not an indicator of phosphine distribution during tablet fumigation.

Evaluation of k-epsilon and RNG turbulence models for confined swirling and non-swirling flow

In a very recent study [1] the Renormalisation Group (RNG) turbulence model was used to obtain flow predictions in a strongly swirling quarl burner, and was found to perform well in predicting certain features that are not well captured using less sophisticated models of turbulence. The implication is that the RNG approach should provide an economical and reliable tool for the prediction of swirling flows in combustor and furnace geometries commonly encountered in technological applications. To test this hypothesis the present work considers flow in a model furnace for which experimental data is available [2]. The essential features of the flow which differentiate it from the previous study [1] are that the annular air jet entry is relatively narrow and the base wall of the cylindrical furnace is at 90 degrees to the inlet pipe. For swirl numbers of order 1 the resulting flow is highly complex with significant inner and outer recirculation regions. The RNG and standard k-epsilon models are used to model the flow for both swirling and non-swirling entry jets and the results compared with experimental data [2]. Near wall viscous effects are accounted for in both models via the standard wall function formulation [3]. For the RNG model, additional computations with grid placement extending well inside the near wall viscous-affected sublayer are performed in order to assess the low Reynolds number capabilities of the model.

Changes in testicular function following specific deprivation of LH in the adult male rabbit

Sexually mature male rabbits actively immunized against highly purified ovine LH (oLH) were used as a model system to study the effects of endogenous LH deprivation (and therefore testosterone) on spermatogenesis as well as pituitary FSH secretion. Immunization against oLH generated antibody titres capable of cross-reacting and neutralizing rabbit LH and this resulted in a significant reduction (P<0.01) in serum testosterone levels by 2-4 weeks of immunization. A significant increase in circulating FSH concentration (from a basal level of similar to 1 ng to 60-100 ng/ml; P<0.01) was observed within 4-6 weeks of immunization, perhaps a consequence of the negative feedback effect of the lack of testosterone. The effect of LH deprivation on spermatogenesis assessed by DNA flow cytometry and histological analyses of testicular biopsy tissue revealed that lack of testosterone primarily results in a rapid reduction and complete absence of round (1C) and elongated (HC) spermatids. The immediate effect of LH/testosterone deprivation thus appears to be at the step of meiotic transformation of primary spermatocytes (4C) to 1C. A significant reduction (>80%; P<0.01) in the 4C population and a relative accumulation (>90%; P<0.01) in spermatogonia (2C) was also observed, suggesting a need for testosterone during the transformation of 2C to 1C. In all but one of the rabbits, both qualitative and quantitative recovery in spermatogenesis occurred during the recovery phase, even at a time when only a marginal increase in serum testosterone (compared with the preimmunization) levels was observed as a result of a rapid decline in the cross-reactive antibody titres. These results clearly show that LH/testosterone deprivation in addition to primarily affecting the meiotic step also regulates the conversion of 2C to 4C during spermatogenesis.

Unsteady Natural Convection Flow over a Heated Cylinder Buried in a Fluid Saturated Porous Medium

Transient natural convection flow on a heated cylinder buried in a semi-infinite liquid-saturated porous medium has been studied. The unsteadiness in the problem arises due to the cylinder which is heated (cooled) suddenly and then maintained at that temperature. The coupled partial differential equations governing the flow and heat transfer are cast into stream function-temperature formulation, and the solutions are obtained from the initial time to the time when steady state is reached. The heat transfer is found to change significantly with increasing time in a small time interval immediately after the start of the impulsive change, and steady state is reached after some time. The average Nusselt number is found to increase with Rayleigh number When the surface of the cylinder is suddenly cooled, there is a change in the direction of the heat transfer in a small time interval immediately after the start of the impulsive change in the surface temperature;however when the surface temperature is suddenly increased, no such phenomenon is observed.

Carotid Stenosis, Endarterectomy, and the Brain : Brain microcirculation, diffusion, and cognitive function before and after carotid endarterectomy in patients with a high-grade carotid stenosis

Carotid atherosclerotic disease is a major cause of stroke, but it may remain clinically asymptomatic. The factors that turn the asymptomatic plaque into a symptomatic one are not fully understood, neither are the subtle effects that a high-grade carotid stenosis may have on the brain. The purpose of this study was to evaluate brain microcirculation, diffusion, and cognitive performance in patients with a high-grade stenosis in carotid artery, clinically either symptomatic or asymptomatic, undergoing carotid endarterectomy (CEA). We wanted to find out whether the stenoses are associated with diffusion or perfusion abnormalities of the brain or variation in the cognitive functioning of the patients, and to what extent the potential findings are affected by CEA, and compare the clinically symptomatic and asymptomatic subjects as well as strictly healthy controls. Coagulation and fibrinolytic parameters were compared with the rate microembolic signals (MES) in transcranial Doppler (TCD) and the macroscopic appearance of stenosing plaques in surgery. Patients (n=92) underwent CEA within the study. Blood samples pertaining to coagulation and fibrinolysis were collected before CEA, and the subjects underwent repeated TCD monitoring for MES. A subpopulation (n= 46) underwent MR imaging and repeated neuropsychological examination (preoperative, as well 4 and 100 days after CEA). In MRI, the average apparent diffusion coefficients were higher in the ipsilateral white matter (WM), and altough the interhemispheric difference was abolished by CEA, the levels remained higher than in controls. Symptomatic stenoses were associated with more sluggish perfusion especially in WM, and lower pulsatility of flow in TCD. All patients had poorer cognitive performance than healthy controls. Cognitive functions improved as expected by learning effect despite transient postoperative worsening in a few subjects. Improvement was greater in patients with deepest hypoperfusion, primarily in executive functions. Symptomatic stenoses were associated with higher hematocrit and tissue plasminogen activator antigen levels, as well as higher rate of MES and ulcerated plaques, and better postoperative improvement of vasoreactivity and pulsatility. In light of the findings, carotid stenosis is associated with differences in brain diffusion, perfusion, and cognition. The effect on diffusion in the ipsilateral WM, partially reversible by CEA, may be associated with WM degeneration. Asymptomatic and symptomatic subpopulations differ from each other in terms of hemodynamic adaptation and in their vascular physiological response to removal of stenosis. Although CEA may be associated with a transient cognitive decline, a true improvement of cognitive performance by CEA is possible in patients with the most pronounced perfusion deficits. Mediators of fibrinolysis and unfavourable hemorheology may contribute to the development of a symptomatic disease in patients with a high-grade stenosis.

Ischemia-reperfusion injury in human liver transplantation : Mechanisms and effects on graft function

Liver transplantation is an established therapy for both acute and chronic liver failure. Despite excellent long-term outcome, graft dysfunction remains a problem affecting up to 15-30% of the recipients. The etiology of dysfunction is multifactorial, with ischemia-reperfusion injury regarded as one of the most important contributors. This thesis focuses on the inflammatory response during graft procurement and reperfusion in liver transplantation in adults. Activation of protein C was examined as a potential endogenous anti-inflammatory mechanism. The effects of inflammatory responses on graft function and outcome were investigated. Seventy adult patients undergoing liver transplantation in Helsinki University Central Hospital, and 50 multiorgan donors, were studied. Blood samples from the portal and the hepatic veins were drawn before graft procurement and at several time points during graft reperfusion to assess changes within the liver. Liver biopsies were taken before graft preservation and after reperfusion. Neutrophil and monocyte CD11b and L-selectin expression were analysed by flow cytometry. Plasma TNF-α, IL-6, IL-8, sICAM-1, and HMGB1 were determined by ELISA and Western-blotting. HMGB1 immunohistochemistry was performed on liver tissue specimens. Plasma protein C and activated protein C were determined by an enzyme-capture assay. Hepatic IL-8 release during graft procurement was associated with subsequent graft dysfunction, biliary in particular, in the recipient. Biliary marker levels increased only 5 7 days after transplantation. Thus, donor inflammatory response appears to influence recipient liver function with relatively long-lasting effects. Hepatic phagocyte activation and sequestration, with concomitant HMGB1 release, occurred during reperfusion. Neither phagocyte activation nor plasma cytokines correlated with postoperative graft function. Thus, activation of the inflammatory responses within the liver during reperfusion may be of minor clinical significance. However, HMGB1 was released from hepatocytes and were also correlated with postoperative transaminase levels. Accordingly, HMGB1 appears to be a marker of hepatocellular injury.

Optimal multicommodity flow through the complete graph with random edge capacities

We consider a multicommodity flow problem on a complete graph whose edges have random, independent, and identically distributed capacities. We show that, as the number of nodes tends to infinity, the maximumutility, given by the average of a concave function of each commodity How, has an almost-sure limit. Furthermore, the asymptotically optimal flow uses only direct and two-hop paths, and can be obtained in a distributed manner.

Heat transfer in rarefied couette flow on the basis of discrete ordinate method

The method of discrete ordinates, in conjunction with the modified "half-range" quadrature, is applied to the study of heat transfer in rarefied gas flows. Analytic expressions for the reduced distribution function, the macroscopic temperature profile and the heat flux are obtained in the general n-th approximation. The results for temperature profile and heat flux are in sufficiently good accord both with the results of the previous investigators and with the experimental data.

Potential functions for hydrogen bond interactions - III. Empirical Potential Function for the Peptide N-H...O=C Hydrogen Bond

A careful comparison of the distribution in the (R, θ)-plane of all NH ... O hydrogen bonds with that for bonds between neutral NH and neutral C=O groups indicated that the latter has a larger mean R and a wider range of θ and that the distribution was also broader than for the average case. Therefore, the potential function developed earlier for an average NH ... O hydrogen bond was modified to suit the peptide case. A three-parameter expression of the form {Mathematical expression}, with △ = R - Rmin, was found to be satisfactory. By comparing the theoretically expected distribution in R and θ with observed data (although limited), the best values were found to be p1 = 25, p3 = - 2 and q1 = 1 × 10-3, with Rmin = 2·95 Å and Vmin = - 4·5 kcal/mole. The procedure for obtaining a smooth transition from Vhb to the non-bonded potential Vnb for large R and θ is described, along with a flow chart useful for programming the formulae. Calculated values of ΔH, the enthalpy of formation of the hydrogen bond, using this function are in reasonable agreement with observation. When the atoms involved in the hydrogen bond occur in a five-membered ring as in the sequence[Figure not available: see fulltext.] a different formula for the potential function is needed, which is of the form Vhb = Vmin +p1△2 +q1x2 where x = θ - 50° for θ ≥ 50°, with p1 = 15, q1 = 0·002, Rmin = 2· Å and Vmin = - 2·5 kcal/mole. © 1971 Indian Academy of Sciences.

Non-linear couette flow with heat transfer using the B-G-K model

In recent years a large number of investigators have devoted their efforts to the study of flow and heat transfer in rarefied gases, using the BGK [1] model or the Boltzmann kinetic equation. The velocity moment method which is based on an expansion of the distribution function as a series of orthogonal polynomials in velocity space, has been applied to the linearized problem of shear flow and heat transfer by Mott-Smith [2] and Wang Chang and Uhlenbeck [3]. Gross, Jackson and Ziering [4] have improved greatly upon this technique by expressing the distribution function in terms of half-range functions and it is this feature which leads to the rapid convergence of the method. The full-range moments method [4] has been modified by Bhatnagar [5] and then applied to plane Couette flow using the B-G-K model. Bhatnagar and Srivastava [6] have also studied the heat transfer in plane Couette flow using the linearized B-G-K equation. On the other hand, the half-range moments method has been applied by Gross and Ziering [7] to heat transfer between parallel plates using Boltzmann equation for hard sphere molecules and by Ziering [83 to shear and heat flow using Maxwell molecular model. Along different lines, a moment method has been applied by Lees and Liu [9] to heat transfer in Couette flow using Maxwell's transfer equation rather than the Boltzmann equation for distribution function. An iteration method has been developed by Willis [10] to apply it to non-linear heat transfer problems using the B-G-K model, with the zeroth iteration being taken as the solution of the collisionless kinetic equation. Krook [11] has also used the moment method to formulate the equivalent continuum equations and has pointed out that if the effects of molecular collisions are described by the B-G-K model, exact numerical solutions of many rarefied gas-dynamic problems can be obtained. Recently, these numerical solutions have been obtained by Anderson [12] for the non-linear heat transfer in Couette flow,

Anomalous behavior of hydrodynamic modes in the two dimensional shear flow of a granular material

The growth rates of the hydrodynamic modes in the homogeneous sheared state of a granular material are determined by solving the Boltzmann equation. The steady velocity distribution is considered to be the product of the Maxwell Boltzmann distribution and a Hermite polynomial expansion in the velocity components; this form is inserted into them Boltzmann equation and solved to obtain the coeificients of the terms in the expansion. The solution is obtained using an expansion in the parameter epsilon =(1 - e)(1/2), and terms correct to epsilon(4) are retained to obtain an approximate solution; the error due to the neglect of higher terms is estimated at about 5% for e = 0.7. A small perturbation is placed on the distribution function in the form of a Hermite polynomial expansion for the velocity variations and a Fourier expansion in the spatial coordinates: this is inserted into the Boltzmann equation and the growth rate of the Fourier modes is determined. It is found that in the hydrodynamic limit, the growth rates of the hydrodynamic modes in the flow direction have unusual characteristics. The growth rate of the momentum diffusion mode is positive, indicating that density variations are unstable in the limit k--> 0, and the growth rate increases proportional to kslash} k kslash}(2/3) in the limit k --> 0 (in contrast to the k(2) increase in elastic systems), where k is the wave vector in the flow direction. The real and imaginary parts of the growth rate corresponding to the propagating also increase proportional to kslash k kslash(2/3) (in contrast to the k(2) and k increase in elastic systems). The energy mode is damped due to inelastic collisions between particles. The scaling of the growth rates of the hydrodynamic modes with the wave vector I in the gradient direction is similar to that in elastic systems. (C) 2000 Elsevier Science B.V. All rights reserved.

The hard-particle model for a dense granular flow down an inclined plane

Perfectly hard particles are those which experience an infinite repulsive force when they overlap, and no force when they do not overlap. In the hard-particle model, the only static state is the isostatic state where the forces between particles are statically determinate. In the flowing state, the interactions between particles are instantaneous because the time of contact approaches zero in the limit of infinite particle stiffness. Here, we discuss the development of a hard particle model for a realistic granular flow down an inclined plane, and examine its utility for predicting the salient features both qualitatively and quantitatively. We first discuss Discrete Element simulations, that even very dense flows of sand or glass beads with volume fraction between 0.5 and 0.58 are in the rapid flow regime, due to the very high particle stiffness. An important length scale in the shear flow of inelastic particles is the `conduction length' delta = (d/(1 - e(2))(1/2)), where d is the particle diameter and e is the coefficient of restitution. When the macroscopic scale h (height of the flowing layer) is larger than the conduction length, the rates of shear production and inelastic dissipation are nearly equal in the bulk of the flow, while the rate of conduction of energy is O((delta/h)(2)) smaller than the rate of dissipation of energy. Energy conduction is important in boundary layers of thickness delta at the top and bottom. The flow in the boundary layer at the top and bottom is examined using asymptotic analysis. We derive an exact relationship showing that the a boundary layer solution exists only if the volume fraction in the bulk decreases as the angle of inclination is increased. In the opposite case, where the volume fraction increases as the angle of inclination is increased, there is no boundary layer solution. The boundary layer theory also provides us with a way of understanding the cessation of flow when at a given angle of inclination when the height of the layer is decreased below a value h(stop), which is a function of the angle of inclination. There is dissipation of energy due to particle collisions in the flow as well as due to particle collisions with the base, and the fraction of energy dissipation in the base increases as the thickness decreases. When the shear production in the flow cannot compensate for the additional energy drawn out of the flow due to the wall collisions, the temperature decreases to zero and the flow stops. Scaling relations can be derived for h(stop) as a function of angle of inclination.

Aspiration-based utility functions in a planning model for timber flow management.

The study presents a theory of utility models based on aspiration levels, as well as the application of this theory to the planning of timber flow economics. The first part of the study comprises a derivation of the utility-theoretic basis for the application of aspiration levels. Two basic models are dealt with: the additive and the multiplicative. Applied here solely for partial utility functions, aspiration and reservation levels are interpreted as defining piecewisely linear functions. The standpoint of the choices of the decision-maker is emphasized by the use of indifference curves. The second part of the study introduces a model for the management of timber flows. The model is based on the assumption that the decision-maker is willing to specify a shape of income flow which is different from that of the capital-theoretic optimum. The utility model comprises four aspiration-based compound utility functions. The theory and the flow model are tested numerically by computations covering three forest holdings. The results show that the additive model is sensitive even to slight changes in relative importances and aspiration levels. This applies particularly to nearly linear production possibility boundaries of monetary variables. The multiplicative model, on the other hand, is stable because it generates strictly convex indifference curves. Due to a higher marginal rate of substitution, the multiplicative model implies a stronger dependence on forest management than the additive function. For income trajectory optimization, a method utilizing an income trajectory index is more efficient than one based on the use of aspiration levels per management period. Smooth trajectories can be attained by squaring the deviations of the feasible trajectories from the desired one.

Nonsimilar boundary layers for non-Darcy mixed convection flow about a horizontal surface in a saturated porous medium

The nonsimilar non-Darcy mixed convection flow about a heated horizontal surface in a saturated porous medium has been studied when the surface temperature is a power function of distance (Tw = T∞ ± Axλ). The analysis is performed for the cases of parallel and stagnation flows with favourable induced pressure gradient. The partial differential equations governing the flow have been solved numerically using the Keller box method. The heat transfer is enhanced due to the buoyancy parameter and wall temperature, but the non-Darcy parameter reduces it. For non-Darcy flow, the similarity solution exists only for the case of parallel flow.

Criteria for prediction of flow instabilities and microstructural manifestations during warm working of AISI 304L stainless steel

The deformation characteristics of 304L stainless steel in compression in the temperature range 20–700°C and strain rate range 0·001–100 s−1 have been studied with the aim of characterising the .flow instabilities occurring in the microstructure. At higher temperatures and strain rates the stainless steel exhibits flow localisation, whereas at temperatures below 500°C and strain rates lower than 0·1 s−1 the flow instabilities are due to dynamic strain aging. Strain induced martensite formation is responsible for the flow instabilities at room temperature and low strain rates (0·01 s−1). In view of the occurrence of these instabilities, cold working is preferable to warm working to achieve dimensional tolerance and reproducible properties in the product. Among the different criteria tested to explain the occurrence of instabilities, the continuum criterion, developed on the basis of the principles of maximum rate of entropy production and separability of the dissipation function, predicts accurately all the above instability features.