The effect of thermal radiation on the natural convection boundary layer flow over a wavy horizontal surface

In this article, natural convection boundary layer flow is investigated over a semi-infinite horizontal wavy surface. Such an irregular (wavy) surface is used to exchange heat with an external radiating fluid which obeys Rosseland diffusion approximation. The boundary layer equations are cast into dimensionless form by introducing appropriate scaling. Primitive variable formulations (PVF) and stream function formulations (SFF) are independently used to transform the boundary layer equations into convenient form. The equations obtained from the former formulations are integrated numerically via implicit finite difference iterative scheme whereas equations obtained from lateral formulations are simulated through Keller-box scheme. To validate the results, solutions produced by above two methods are compared graphically. The main parameters: thermal radiation parameter and amplitude of the wavy surface are discussed categorically in terms of shear stress and rate of heat transfer. It is found that wavy surface increases heat transfer rate compared to the smooth wall. Thus optimum heat transfer is accomplished when irregular surface is considered. It is also established that high amplitude of the wavy surface in the boundary layer leads to separation of fluid from the plate.

Moffatt vortices induced by the motion of a contact line

A recent hydrodynamic theory of liquid slippage on a solid substrate (Kirkinis & Davis, Phys. Rev. Lett., vol. 110, 2013, 234503) gives rise to a sequence of eddies (Moffatt vortices) that co-move with a moving contact line (CL) in a liquid wedge. The presence of these vortices is established through secular equations that depend on the dynamic contact angle α and capillary number Ca. The limiting case α→O is associated with the appearance of such vortices in a channel. The vortices are generated by the relative motion of the interfaces, which in turn is due to the motion of the CL. This effect has yet to be observed in experiment.

Effect of water management practice on pesticide behavior in paddy water

The fate and transport of three herbicides commonly used in rice production in Japan were compared using two water management practices. The herbicides were simetryn, thiobencarb and mefenacet. The first management practice was an intermittent irrigation scheme using an automatic irrigation system (AI) with a high drainage gate and the second one was a continuous irrigation and overflow drainage scheme (CI) in experimental paddy fields. Dissipation of the herbicides appeared to follow first order kinetics with the half-lives (DT₅₀) of 1.6-3.4 days and the DT₉₀ (90% dissipation) of 7.4-9.8 days. The AI scheme had little drainage even during large rainfall events thus resulting in losses of less than 4% of each applied herbicide through runoff. Meanwhile the CI scheme resulted in losses of about 37%, 12% and 35% of the applied masses of simetryn, thiobencarb and mefenacet, respectively. The intermittent irrigation scheme using an automatic irrigation system with a high drainage gate saved irrigation water and prevented herbicide runoff whereas the continuous irrigation and overflow scheme resulted in significant losses of water as well as the herbicides. Maintaining the excess water storage is important for preventing paddy water runoff during significant rainfall events. The organic carbon partition coefficient K_oc seems to be a strong indicator of the aquatic fate of the herbicide as compared to the water solubility (S_W). However, further investigations are required to understand the relation between K_oc and the agricultural practices upon the pesticide fate and transport. An extension of the water holding period up to 10 days after herbicide application based on the DT₉₀ from the currently specified period of 3-4 days in Japan is recommended to be a good agricultural practice for controlling the herbicide runoff from paddy fields. Also, the best water management practice, which can be recommended for use during the water holding period, is the intermittent irrigation scheme using an automatic irrigation system with a high drainage gate. © 2006 Elsevier B.V. All rights reserved.

Rattling modes and the intrinsic vibrational spectrum of beetle-type scanning tunneling microscopes

It is known that the vibrational spectra of beetle-type scanning tunneling microscopes with a total mass of ≈3–4 g contain extrinsic ‘rattling’ modes in the frequency range extending from 500 to 1700 Hz that interfere with image acquisition. These modes lie below the lowest calculated eigenfrequency of the beetle and it has been suggested that they arise from the inertial sliding of the beetle between surface asperities on the raceway. In this paper we describe some cross-coupling measurements that were performed on three home-built beetle-type STMs of two different designs. We provide evidence that suggests that for beetles with total masses of 12–15 g all the modes in the rattling range are intrinsic. This provides additional support for the notion that the vibrational properties of beetle-type scanning tunneling microscopes can be improved by increasing the contact pressure between the feet of the beetle and the raceway.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter = 9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 10(8). Under these conditions the convection is turbulent, and the time-averaged velocity at any point is `zero'. The Reynolds number based on the Taylor microscale, Re-lambda, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6-7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as (RaSc1/2)-Sc-1/2, and the Reynolds number would scale as (RaSc-1/2)-Sc-1/2. The velocity and the flux measurements appear to be consistent with the Ra-1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are Much higher compared to what would be obtained in Rayleigh-Benard (R-B) convection for similar density differences.

Effects of the filling process on the evolution of the mushy zone and macrosegregation in alloy casting

A numerical model of the entire casting process starting from the mould filling stage to complete solidification is presented. The model takes into consideration any phase change taking place during the filling process. A volume of fluid method is used for tracking the metal–air interface during filling and an enthalpy based macro-scale solidification model is used for the phase change process. The model is demonstrated for the case of filling and solidification of Pb–15 wt%Sn alloy in a side-cooled two-dimensional rectangular cavity, and the resulting evolution of a mushy region and macrosegregation are studied. The effects of process parameters related to filling, namely degree of melt superheat and filling velocity on macrosegregation in the cavity, are also investigated. Results show significant differences in the progress of the mushy zone and macrosegregation pattern between this analysis and conventional analysis without the filling effect.

Purification and functional characterization of type II DNA topoisomerase from rat testis and comparison with topoisomerase II from liver

A number of studies in yeast have shown that DNA topoisomerase TI is essential for chromosome condensation and disjunction during mitosis at the metaphase/anaphase transition and meiosis I. Accordingly, kinetic and mechanistic studies have implied a role for topoisomerase rr in chromosome disjunction. As a step toward understanding the nature and role of topoisomerase II in a mammalian germline in vivo, we have purified topoisomerase II from rat testis to homogeneity and ascertained several of its catalytic activities in conjunction with that of the purified enzyme from liver. The purified enzymes appeared to be monomers under denaturing conditions; however, they differed in their relative molecular mass. Topoisomerase II from testis and liver have apparent molecular masses of 150 +/- 10 kDa and 160 +/- 10 kDa, respectively. The native molecular mass of testis topoisomerase II as assayed by immunoblot analysis of cell-foe extracts, prepared in the presence of SDS and a number of protease inhibitors, corroborated with the size of the purified enzyme. Both enzymes are able to promote decatenation and relax supercoiled DNA substrates in an ATP and Mg2+-dependent manner. However, quantitative comparison of catalytic properties of topoisomerase II from testis with that of the enzyme from liver displayed significant differences in their efficiencies. Optimal pH values for testis enzyme are 6.5 to 8.5 while they are 6 to 7.5 for the liver enzyme. Intriguingly, the relaxation activity of liver topoisomerase II was inhibited by potassium glutamate at 1 M, whereas testis enzyme required about half its concentration. These findings argue that topoisomerase II from rat testis is structurally distinct from that of its somatic form and the functional differences between the two enzymes parallels with the physiological environment that is unique to these two tissues.

Process Analytical Technology Approach on Fluid Bed Granulation and Drying : Identifying Critical Relationships and Constructing the Design Space

Fluid bed granulation is a key pharmaceutical process which improves many of the powder properties for tablet compression. Dry mixing, wetting and drying phases are included in the fluid bed granulation process. Granules of high quality can be obtained by understanding and controlling the critical process parameters by timely measurements. Physical process measurements and particle size data of a fluid bed granulator that are analysed in an integrated manner are included in process analytical technologies (PAT). Recent regulatory guidelines strongly encourage the pharmaceutical industry to apply scientific and risk management approaches to the development of a product and its manufacturing process. The aim of this study was to utilise PAT tools to increase the process understanding of fluid bed granulation and drying. Inlet air humidity levels and granulation liquid feed affect powder moisture during fluid bed granulation. Moisture influences on many process, granule and tablet qualities. The approach in this thesis was to identify sources of variation that are mainly related to moisture. The aim was to determine correlations and relationships, and utilise the PAT and design space concepts for the fluid bed granulation and drying. Monitoring the material behaviour in a fluidised bed has traditionally relied on the observational ability and experience of an operator. There has been a lack of good criteria for characterising material behaviour during spraying and drying phases, even though the entire performance of a process and end product quality are dependent on it. The granules were produced in an instrumented bench-scale Glatt WSG5 fluid bed granulator. The effect of inlet air humidity and granulation liquid feed on the temperature measurements at different locations of a fluid bed granulator system were determined. This revealed dynamic changes in the measurements and enabled finding the most optimal sites for process control. The moisture originating from the granulation liquid and inlet air affected the temperature of the mass and pressure difference over granules. Moreover, the effects of inlet air humidity and granulation liquid feed rate on granule size were evaluated and compensatory techniques used to optimize particle size. Various end-point indication techniques of drying were compared. The ∆T method, which is based on thermodynamic principles, eliminated the effects of humidity variations and resulted in the most precise estimation of the drying end-point. The influence of fluidisation behaviour on drying end-point detection was determined. The feasibility of the ∆T method and thus the similarities of end-point moisture contents were found to be dependent on the variation in fluidisation between manufacturing batches. A novel parameter that describes behaviour of material in a fluid bed was developed. Flow rate of the process air and turbine fan speed were used to calculate this parameter and it was compared to the fluidisation behaviour and the particle size results. The design space process trajectories for smooth fluidisation based on the fluidisation parameters were determined. With this design space it is possible to avoid excessive fluidisation and improper fluidisation and bed collapse. Furthermore, various process phenomena and failure modes were observed with the in-line particle size analyser. Both rapid increase and a decrease in granule size could be monitored in a timely manner. The fluidisation parameter and the pressure difference over filters were also discovered to express particle size when the granules had been formed. The various physical parameters evaluated in this thesis give valuable information of fluid bed process performance and increase the process understanding.

Numerical studies on columnar-to-equiaxed transition in directional solidification of binary alloys

A numerical study on columnar-to-equiaxed transition (CET) during directional solidification of binary alloys is presented using a macroscopic solidification model. The position of CET is predicted numerically using a critical cooling rate criterion reported in literature. The macroscopic solidification model takes into account movement of solid phase due to buoyancy, and drag effect on the moving solid phase because of fluid motion. The model is applied to simulate the solidification process for binary alloys (Sn-Pb) and to estimate solidification parameters such as position of the liquidus, velocity of the liquidus isotherm, temperature gradient ahead of the liquidus, and cooling rate at the liquidus. Solidification phenomena under two cooling configurations are studied: one without melt convection and the other involvin thermosolutal convection. The numerically predicted positions of CET compare well with those of experiments reported in literature. Melt convection results in higher cooling rate, higher liquidus isotherm velocities, and stimulation of occurrence of CET in comparison to the nonconvecting case. The movement of solid phase aids further the process of CET. With a fixed solid phase, the occurrence of CET based on the same critical cooling rate is delayed and it occurs at a greater distance from the chill.

Mathematical modelling of the impaction and spreading of spray droplets on leaves

This thesis concerns the development of mathematical models to describe the interactions that occur between spray droplets and leaves. Models are presented that not only provide a contribution to mathematical knowledge in the field of fluid dynamics, but are also of utility within the agrichemical industry. The thesis is presented in two parts. First, thin film models are implemented with efficient numerical schemes in order to simulate droplets on virtual leaf surfaces. Then the interception event is considered, whereby energy balance techniques are employed to instantaneously predict whether an impacting droplet will bounce, splash, or adhere to a leaf.

On the origin of the inflectional instability of a laminar separation bubble

This is an experimental and theoretical Study of a laminar separation bubble and the associated linear stability mechanisms. Experiments were performed over a flat plate kept in a wind tunnel, with an imposed pressure gradient typical of an aerofoil that would involve a laminar separation bubble. The separation bubble was characterized by measurement of surface-pressure distribution and streamwise velocity using hot-wire anemometry. Single component hot-wire anemometry was also used for a detailed study of the transition dynamics. It was foundthat the so-called dead-air region in the front portion of the bubble corresponded to a region of small disturbance amplitudes, with the amplitude reaching a maximum value close to the reattachment point. An exponential growth rate of the disturbance was seen in the region upstream of the mean maximum height of the bubble, and this was indicative of a linear instability mechanism at work. An infinitesimal disturbance was impulsively introduced into the boundary layer upstream of separation location, and the wave packet was tracked (in an ensemble-averaged sense) while it was getting advected downstream. The disturbance was found to be convective in nature. Linear stability analyses (both the Orr-Sommerfeld and Rayleigh calculations) were performed for mean velocity profiles, starting from an attached adverse-pressure-gradient boundary layer all the way up to the front portion of the separation-bubble region (i.e. up to the end of the dead-air region in which linear evolution of the disturbance could be expected). The conclusion from the present work is that the primary instability mechanism in a separation bubble is inflectional in nature, and its origin can be traced back to upstream of the separation location. In other words, the inviscid inflectional instability of the separated shear layer should be logically seen as an extension of the instability of the upstream attached adverse-pressure-gradient boundary layer. This modifies the traditional view that pegs the origin of the instability in a separation bubble to the detached shear layer Outside the bubble, with its associated Kelvin-Helmholtz mechanism. We contendthat only when the separated shear layer has moved considerably away from the wall (and this happens near the maximum-height location of the mean bubble), a description by the Kelvin-Helmholtz instability paradigm, with its associated scaling principles, Could become relevant. We also propose a new scaling for the most amplified frequency for a wall-bounded shear layer in terms of the inflection-point height and the vorticity thickness and show it to be universal.

Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh

A numerical scheme is presented for accurate simulation of fluid flow using the lattice Boltzmann equation (LBE) on unstructured mesh. A finite volume approach is adopted to discretize the LBE on a cell-centered, arbitrary shaped, triangular tessellation. The formulation includes a formal, second order discretization using a Total Variation Diminishing (TVD) scheme for the terms representing advection of the distribution function in physical space, due to microscopic particle motion. The advantage of the LBE approach is exploited by implementing the scheme in a new computer code to run on a parallel computing system. Performance of the new formulation is systematically investigated by simulating four benchmark flows of increasing complexity, namely (1) flow in a plane channel, (2) unsteady Couette flow, (3) flow caused by a moving lid over a 2D square cavity and (4) flow over a circular cylinder. For each of these flows, the present scheme is validated with the results from Navier-Stokes computations as well as lattice Boltzmann simulations on regular mesh. It is shown that the scheme is robust and accurate for the different test problems studied.

Structural analyses of (1->3),(1->4)-β-D-glucan of oats and barley

The structures of (1→3),(1→4)-β-D-glucans of oat bran, whole-grain oats and barley and processed foods were analysed. Various methods of hydrolysis of β-glucan, the content of insoluble fibre of whole grains of oats and barley and the solution behaviour of oat and barley β-glucans were studied. The isolated soluble β-glucans of oat bran and whole-grain oats and barley were hydrolysed with lichenase, an enzyme specific for (1→3),(1→4)-β-D-β-glucans. The amounts of oligosaccharides produced from bran were analysed with capillary electrophoresis and those from whole-grains with high-performance anion-exchange chromatography with pulse-amperometric detection. The main products were 3-O-β-cellobiosyl-D-glucose and 3-O-β-cellotriosyl-D-glucose, the oligosaccharides which have a degree of polymerisation denoted by DP3 and DP4. Small differences were detected between soluble and insoluble β-glucans and also between β-glucans of oats and barley. These differences can only be seen in the DP3:DP4 ratio which was higher for barley than for oat and also higher for insoluble than for soluble β-glucan. A greater proportion of barley β-glucan remained insoluble than of oat β-glucan. The molar masses of soluble β-glucans of oats and barley were the same as were those of insoluble β-glucans of oats and barley. To analyse the effects of cooking, baking, fermentation and drying, β-glucan was isolated from porridge, bread and fermentate and also from their starting materials. More β-glucan was released after cooking and less after baking. Drying decreased the extractability for bread and fermentate but increased it for porridge. Different hydrolysis methods of β-glucan were compared. Acid hydrolysis and the modified AOAC method gave similar results. The results of hydrolysis with lichenase gave higher recoveries than the other two. The combination of lichenase hydrolysis and high-performance anion-exchange chromatography with pulse-amperometric detection was found best for the analysis of β-glucan content. The content of insoluble fibre was higher for barley than for oats and the amount of β-glucan in the insoluble fibre fraction was higher for oats than for barley. The flow properties of both water and aqueous cuoxam solutions of oat and barley β-glucans were studied. Shear thinning was stronger for the water solutions of oat β-glucan than for barley β-glucan. In aqueous cuoxam shear thinning was not observed at the same concentration as in water but only with high concentration solutions. Then the viscosity of barley β-glucan was slightly higher than that of oat β-glucan. The oscillatory measurements showed that the crossover point of the G´ and G´´ curves was much lower for barley β-glucan than for oat β-glucan indicating a higher tendency towards solid-like behaviour for barley β-glucan than for oat β-glucan.

Dynamics of dense sheared granular flows. Part 1. Structure and diffusion

Shear flows of inelastic spheres in three dimensions in the Volume fraction range 0.4-0.64 are analysed using event-driven simulations.Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to the line Joining the centres). Here, we have considered both e(t) = +1 and e(t) = e(n) (rough particles) and e(t) =-1 (smooth particles), and the normal coefficient of restitution e(n) was varied in the range 0.6-0.98. Care was taken to avoid inelastic collapse and ensure there are no particle overlaps during the simulation. First, we studied the ordering in the system by examining the icosahedral order parameter Q(6) in three dimensions and the planar order parameter q(6) in the plane perpendicular to the gradient direction. It was found that for shear flows of sufficiently large size, the system Continues to be in the random state, with Q(6) and q(6) close to 0, even for volume fractions between phi = 0.5 and phi = 0.6; in contrast, for a system of elastic particles in the absence of shear, the system orders (crystallizes) at phi = 0.49. This indicates that the shear flow prevents ordering in a system of sufficiently large size. In a shear flow of inelastic particles, the strain rate and the temperature are related through the energy balance equation, and all time scales can be non-dimensionalized by the inverse of the strain rate. Therefore, the dynamics of the system are determined only by the volume fraction and the coefficients of restitution. The variation of the collision frequency with volume fraction and coefficient of estitution was examined. It was found, by plotting the inverse of the collision frequency as a function of volume fraction, that the collision frequency at constant strain rate diverges at a volume fraction phi(ad) (volume fraction for arrested dynamics) which is lower than the random close-packing Volume fraction 0.64 in the absence of shear. The volume fraction phi(ad) decreases as the coefficient of restitution is decreased from e(n) = 1; phi(ad) has a minimum of about 0.585 for coefficient of restitution e(n) in the range 0.6-0.8 for rough particles and is slightly larger for smooth particles. It is found that the dissipation rate and all components of the stress diverge proportional to the collision frequency in the close-packing limit. The qualitative behaviour of the increase in the stress and dissipation rate are well Captured by results derived from kinetic theory, but the quantitative agreement is lacking even if the collision frequency obtained from simulations is used to calculate the pair correlation function used In the theory.

Dynamics of dense sheared granular flows. Part II. The relative velocity distributions

The distribution of relative velocities between colliding particles in shear flows of inelastic spheres is analysed in the Volume fraction range 0.4-0.64. Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to line joining the centres). The distribution or pre-collisional normal relative velocities (along the line Joining the centres of the particles) is Found to be an exponential distribution for particles with low normal coefficient of restitution in the range 0.6-0.7. This is in contrast to the Gaussian distribution for the normal relative velocity in all elastic fluid in the absence of shear. A composite distribution function, which consists of an exponential and a Gaussian component, is proposed to span the range of inelasticities considered here. In the case of roughd particles, the relative velocity tangential to the surfaces at contact is also evaluated, and it is found to be close to a Gaussian distribution even for highly inelastic particles.Empirical relations are formulated for the relative velocity distribution. These are used to calculate the collisional contributions to the pressure, shear stress and the energy dissipation rate in a shear flow. The results of the calculation were round to be in quantitative agreement with simulation results, even for low coefficients of restitution for which the predictions obtained using the Enskog approximation are in error by an order of magnitude. The results are also applied to the flow down an inclined plane, to predict the angle of repose and the variation of the volume fraction with angle of inclination. These results are also found to be in quantitative agreement with previous simulations.