A model for heat transfer in a honey bee swarm

A swarm is a temporary structure formed when several thousand honey bees leave their hive and settle on some object such as the branch of a tree. They remain in this position until a suitable site for a new home is located by the scout bees. A continuum model based on heat conduction and heat generation is used to predict temperature profiles in swarms. Since internal convection is neglected, the model is applicable only at low values of the ambient temperature T-a. Guided by the experimental observations of Heinrich (1981a-c, J. Exp. Biol. 91, 25-55; Science 212, 565-566; Sci. Am. 244, 147-160), the analysis is carried out mainly for non-spherical swarms. The effective thermal conductivity is estimated using the data of Heinrich (1981a, J. Exp. Biol. 91, 25-55) for dead bees. For T-a = 5 and 9 degrees C, results based on a modified version of the heat generation function due to Southwick (1991, The Behaviour and Physiology of Bees, PP 28-47. C.A.B. International, London) are in reasonable agreement with measurements. Results obtained with the heat generation function of Myerscough (1993, J. Theor. Biol. 162, 381-393) are qualitatively similar to those obtained with Southwick's function, but the error is more in the former case. The results suggest that the bees near the periphery generate more heat than those near the core, in accord with the conjecture of Heinrich (1981c, Sci. Am. 244, 147-160). On the other hand, for T-a = 5 degrees C, the heat generation function of Omholt and Lonvik (1986, J. Theor. Biol. 120, 447-456) leads to a trivial steady state where the entire swarm is at the ambient temperature. Therefore an acceptable heat generation function must result in a steady state which is both non-trivial and stable with respect to small perturbations. Omholt and Lonvik's function satisfies the first requirement, but not the second. For T-a = 15 degrees C, there is a considerable difference between predicted and measured values, probably due to the neglect of internal convection in the model.

Modeling Of Hydrolytic Polymerization In A Semibatch Nylon 6 Reactor

Hydrolytic polymerization of caprolactam to Nylon 6 in a semibatch reactor is carried out by heating a mixture of water and caprolactam. Evaporation of volatiles caused by heating results in a pressure build-up. After the pressure reaches a predetermined value, vapors are vented to keep the pressure constant for some time, and thereafter, to lower the pressure to a value slightly above atmospheric in a preprogrammed manner. The characteristics of the polymer are determined by the chemical reactions and the vaporization of water and caprolactam. The semibatch operation has been simulated and the predictions have been compared with industria data. The observed temperature and pressure histories were predicted with a fair degree of accuracy. It was found that the predictions of the degree of polymerization however are sensitive to the vapor-liquid equilibrium relations. A comparison with an earlier model, which neglected mass transfer resistance, indicates that simulation using the VLE data of Giori and Hayes and accounting for mass transfer resistance is more reliable.

Critical earthquake input power spectral density function models for engineering structures

The problem of determining optimal power spectral density models for earthquake excitation which satisfy constraints on total average power, zero crossing rate and which produce the highest response variance in a given linear system is considered. The solution to this problem is obtained using linear programming methods. The resulting solutions are shown to display a highly deterministic structure and, therefore, fail to capture the stochastic nature of the input. A modification to the definition of critical excitation is proposed which takes into account the entropy rate as a measure of uncertainty in the earthquake loads. The resulting problem is solved using calculus of variations and also within linear programming framework. Illustrative examples on specifying seismic inputs for a nuclear power plant and a tall earth dam are considered and the resulting solutions are shown to be realistic.

On the dual head-discharge characteristics of a modified chimney weir

This paper is concerned with the dual head-discharge characteristics of a modified chimney weir. It is shown by an optimization procedure that the modified chimney weir having an inward trapezoidal weir over an inverted V-notch can produce discharges proportional to both the logarithm of the head as well as the linear power of the head reckoned independently over two different reference planes, within a prefixed maximum allowable percentage error from the theoretical discharge. A new technique is adopted to optimize the weir parameters, in order to obtain the maximum ranges of measurement under logarithmic, linear as well as combined characteristics. In the case of linear weir it is shown that it is possible to enhance the linearity range of the chimney weir by more than 540% and for a weir with constant indication accuracy by more than 350%. In addition, about 86% of the overall depth of the designed linear weir and over 90% of the logarithmic weir is converted as the corresponding measurable ranges. Experiments with four typical weirs give consistent constant average coefficient of discharge for each weir confirming the theory. The practical application of the weir in minor irrigation, hydraulic and other engineering fields is highlighted.

Engineered dimer interface mutants of triosephosphate isomerase: the role of inter-subunit interactions in enzyme function and stability

The role of inter-subunit interactions in maintaining optimal catalytic activity in triosephosphate isomerase (TIM) has been probed, using the Plasmodium falciparum enzyme as a model. Examination of subunit interface contacts in the crystal structures suggests that residue 75 (Thr, conserved) and residue 13 (Cys, variable) make the largest number of inter-subunit contacts. The mutants Cys13Asp (C13D) and Cys13Glu (C13E) have been constructed and display significant reduction in catalytic activity when compared with wild-type (WT) enzyme (similar to 7.4-fold decrease in k(cat) for the C13D and similar to 3.3-fold for the C13E mutants). Analytical gel filtration demonstrates that the C13D mutant dissociates at concentrations < 1.25 mu M, whereas the WT and the C13E enzymes retain the dimeric structure. The order of stability of the mutants in the presence of chemical denaturants, like urea and guanidium chloride, is WT > Cys13Glu > Cys13Asp. Irreversible thermal precipitation temperatures follow the same order as well. Modeling studies establish that the Cys13Asp mutation is likely to cause a significantly greater structural perturbation than Cys13Glu. Analysis of sequence and structural data for TIMs from diverse sources suggests that residues 13 and 82 form a pair of proximal sites, in which a limited number of residue pairs may be accommodated.

Oxidation of Sulfur Dioxide in Aqueous Suspensions of Activated Carbon

Sulfur dioxide in aqueous solutions at low pH levels exists both in molecular SO2(aq) and in hydrolyzed ionic form HSO3-. Experiments indicate that only HSO3- is the reacting species in the oxidation catalyzed by activated carbon, while SO2(aq) deactivates by competing with HSO3 for the active sites of the catalyst particles. A mechanism is proposed and a rate model is developed that also accounts for the effect of sulfuric acid (product of the oxidation) on the solubility of sulfur dioxide. It predicts first order in HSO3-, half order in dissolved oxygen, and a linear deactivation effect of SO2(aq), which are confirmed by experimental data. The deactivation reaches a constant level corresponding to saturation of the active sites by SO2(aq). Activation energy for the oxidation is 93.55 kJ mol(-1) and for deactivation is 21.4 kJ mol(-1).

Some Basic Aspects of Reaction Engineering of Precipitation Processes

Analysis of precipitation reactions is extremely important in the technology of production of fine particles from the liquid phase. The control of composition and particle size in precipitation processes requires careful analysis of the several reactions that comprise the precipitation system. Since precipitation systems involve several, rapid ionic dissociation reactions among other slower ones, the faster reactions may be assumed to be nearly at equilibrium. However, the elimination of species, and the consequent reduction of the system of equations, is an aspect of analysis fraught with the possibility of subtle errors related to the violation of conservation principles. This paper shows how such errors may be avoided systematically by relying on the methods of linear algebra. Applications are demonstrated by analyzing the reactions leading to the precipitation of calcium carbonate in a stirred tank reactor as well as in a single emulsion drop. Sample calculations show that supersaturation dynamics can assume forms that can lead to subsequent dissolution of particles that have once been precipitated.

Gas-phase controlled mass transfer in a foam-bed reactor

Gas-phase controlled absorption of ammonia in foams made of solutions of sulphuric acid has been studied experimentally. Effects of gas-phase concentration of ammonia and type of surfactant on the performance of the foam-bed reactor are investigated. Gas-phase controlled absorption from a spherical bubble is anaylzed using the asymptotic value of Sherwood number (Sh = 6.58), for both negligible as well as significant changes in the volume of the bubble. The experimental data are shown to be in good agreement with the single-stage model of the foam-bed reactor using these asymptotic sub-models, as well as the diffusion-in-sphere analysis available in literature. Influence of effective diffusivity on the time dependence of fractional gas absorption has been found to be unimportant for foam columns with large times of contact. The asymptotic sub-models have been compared and use of the rigid-sphere asymptotic sub-model is recommended for foam columns of practical relevence.

Prediction of separation factor in foam separation of proteins

A phenomenological model has been developed for predicting separation factors obtained in concentrating protein solutions using batch-foam columns. The model considers the adsorption of surface active proteins onto the air-water interface of bubbles, and drainage of liquid from the foam, which are the two predominant processes responsible for separation in foam columns. The model has been verified with data collected on casein and bovine serum albumin (BSA) solutions, for which adsorption isotherms are available in the literature. It has been found that an increase in liquid pool height above the gas distributor and the time allowed for drainage result in a better separation. Further, taller foam columns yield poorer separation at constant time of drainage. The model successfully predicts the observed results. (C) 1997 Elsevier Science Ltd.

Transient analysis of mammalian cell growth in hollow fibre bioreactor

A mathematical model describing the dynamics of mammalian cell growth in hollow fibre bioreactor operated in closed shell mode is developed. Mammalian cells are assumed to grow as an expanding biofilm in the extra-capillary space surrounding the fibre. Diffusion is assumed to be the dominant process in the radial direction while axial convection dominates in the lumen of the bioreactor. The transient simulation results show that steep gradients in the cell number are possible under the condition of substrate limitation. The precise conditions which result in nonuniform growth of cells along the length of the bioreactor are delineated. The effect of various operating conditions, such as substrate feed rate, length of the bioreactor and diffusivity of substrate in different regions of the bioreactor, on the bioreactor performance are evaluated in terms of time required to attain the steady-state. The rime of growth is introduced as a measure of effectiveness factor for the bioreactor and is found to be dependent on two parameters, a modified Peclet number and a Thiele modulus. Diffusion, reaction and/or convection control regimes are identified based on these two parameters. The model is further extended to include dual substrate growth limitations, and the relative growth limiting characteristics of two substrates are evaluated. (C) 1997 Elsevier Science Ltd.

Enzymatic synthesis of flavors in supercritical carbon dioxide

Commercially important flavor esters of isoamyl alcohol, catalyzed by crude hog pancreas lipase (HPL), were synthesized under solvent-free conditions and in supercritical carbon dioxide. The esters synthesized were isoamyl acetate, isoamyl propionate, isoamyl butyrate, and isoamyl octanoate. Very low yields (3-4%) of isoamyl acetate were obtained, but high yields for the other three esters were obtained under both supercritical and solvent-free conditions. The yields of esters of the even-carbon acids, isoamyl acetate, butyrate, and octanoate, increased with increasing chain length, whereas the yield of isoamyl propionate was higher than that of isoamyl butyrate. The optimum temperature of the reaction was higher under supercritical conditions (45 degreesC) than under solvent-free conditions (35-40 degreesC). The effects of other parameters such as alcohol concentration, water concentration, and enzyme loading were investigated. An increase in the water concentration decreased the conversion significantly in supercritical carbon dioxide but not under solvent-free conditions. The optimum ratio of alcohol to acid was dependent on the extent of inhibition by the acid. Although providing a higher apparent yield by being run in a highly concentrated medium, the overall conversion under solvent-free conditions was lower than that under supercritical conditions for similar enzyme concentrations, indicating that the synthesis of esters in supercritical carbon dioxide might be a viable option.

Amino acid interaction preferences in helical membrane proteins

Membrane proteins are involved in a number of important biological functions. Yet, they are poorly understood from the structure and folding point of view. The external environment being drastically different from that of globular proteins, the intra-protein interactions in membrane proteins are also expected to be different. Hence, statistical potentials representing the features of inter-residue interactions based exclusively on the structures of membrane proteins are much needed. Currently, a reasonable number of structures are available, making it possible to undertake such an analysis on membrane proteins. In this study we have examined the inter-residue interaction propensities of amino acids in the membrane spanning regions of the alpha-helical membrane (HM) proteins. Recently we have shown that valuable information can be obtained on globular proteins by the evaluation of the pair-wise interactions of amino acids by classifying them into different structural environments, based on factors such as the secondary structure or the number of contacts that a residue can make. Here we have explored the possible ways of classifying the intra-protein environment of HM proteins and have developed scoring functions based on different classification schemes. On evaluation of different schemes, we find that the scheme which classifies amino acids to different intra-contact environment is the most promising one. Based on this classification scheme, we also redefine the hydrophobicity scale of amino acids in HM proteins.

Effect of Reagent Addition Rate and Temperature on Synthesis of Gold Nanoparticles in Microemulsion Route

Nanoparticle synthesis in a microemulsion route is typically controlled by changing the water to surfactant ratio, concentration of precursors, and/or concentration of micelles. The experiments carried out in this work with chloroauric acid and hydrazine hydrate as precursors in water/AOT-Brij30/isooctane microemulsions show that the reagent addition rate can also be used to tune the size of stable spherical gold nanoparticles to some extent. The particle size goes through a minimum with variation in feed addition rate. The increase in particle size with an increase in reaction temperature is in agreement with an earlier report. A population balance model is used to interpret the experimental findings. The reduced extent of nucleation at low feed addition rates and suppression of nucleation due to the finite rate of mixing at higher addition rates produce a minimum in particle size. The increase in particle size at higher reaction temperatures is explained through an increase in fusion efficiency of micelles which dissipates supersaturation; increase in solubility is shown to play an insignificant role. The moderate polydispersity of the synthesized particles is due to the continued nucleation and growth of particles. The polydispersity of micelle sizes by itself plays a minor role.

Kinetics of catalytic degradation of polycarbonate in benzene

The degradation kinetics of polycarbonate [poly(bisphenol A carbonate)] in benzene catalyzed by commercial (rutile) TiO2 (BET surface area = 11 m(2)/g), anatase TiO2 (156 m(2)/g), and 1 atom % Pt/TiO2 (111 m2/g), prepared by the solution combustion technique, was investigated at various temperatures (230-280 degreesC) and 50 atm. The time evolution of the molecular weight distribution (MWD) was determined by gel permeation chromatography (GPC) and modeled with continuous distribution kinetics to obtain the degradation rate coefficients. The rate coefficients for the catalytic degradation of polycarbonate increased by factors of 20, 3.5, and 1.3 compared to the rate coefficients for thermal degradation when catalyzed by nanosized TiO2 anatase, Pt/TiO2 anatase, and commercial TiO2, respectively, at 280 degreesC. The increased catalytic activity of combustion for synthesized TiO2 and 1% Pt/TiO2 might be due to the increased acidity and BET surface area. The activation energies, determined from the temperature dependencies of the rate coefficients, were 16.3, 21.5, and 39.1 kcal/mol for commercial TiO2, combustion-synthesized Pt/TiO2, and anatase TiO2, respectively.