Prediction of slurry transport in SAG mills using SPH fluid flow in a dynamic DEM based porous media

DEM modelling of the motion of coarse fractions of the charge inside SAG mills has now been well established for more than a decade. In these models the effect of slurry has broadly been ignored due to its complexity. Smoothed particle hydrodynamics (SPH) provides a particle based method for modelling complex free surface fluid flows and is well suited to modelling fluid flow in mills. Previous modelling has demonstrated the powerful ability of SPH to capture dynamic fluid flow effects such as lifters crashing into slurry pools, fluid draining from lifters, flow through grates and pulp lifter discharge. However, all these examples were limited by the ability to model only the slurry in the mill without the charge. In this paper, we represent the charge as a dynamic porous media through which the SPH fluid is then able to flow. The porous media properties (specifically the spatial distribution of porosity and velocity) are predicted by time averaging the mill charge predicted using a large scale DEM model. This allows prediction of transient and steady state slurry distributions in the mill and allows its variation with operating parameters, slurry viscosity and slurry volume, to be explored. (C) 2006 Published by Elsevier Ltd.

Transient temperature rise in a mouse due to low-frequency regional hyperthermia

A refined nonlinear heat transfer model of a mouse has been developed to simulate the transient temperature rise in a neoplastic tumour and neighbouring tissue during regional hyperthermia using a 150 kHz inductive coil. In this study, we incorporate various bio-energetic enhancements to the heat transfer equation and numerical validations based on experimental findings for the mouse, in terms of nonlinear metabolic heat production, homeothermy, blood perfusion parameters, thermoregulation, psychological and physiological effects. The discretized bio-heat transfer equation has been validated with the commercial software FEMLAB on a canonical multi-sphere object before applying the scheme to the inhomogeneous mouse voxel phantom. The time-dependent numerical results of regional hyperthermia of mouse thigh have been compared with the available experimental temperature results with only a few small disparities. During the first 20 min of local unfocused heating, the temperature in the tumour and the surrounding tissue increased by around 7.5 degrees C. The objective of this preliminary study was to develop a validated electrothermal numerical scheme for inductive hyperthermia of a small mammal with the intention of expanding the model into a complete numerical solution involving ferromagnetic nanoparticles for targeted heating of tumours at low frequencies. In addition, the numerical scheme herein could assist in optimizing and tailoring of focused electromagnetic fields for hyperthermia.

Modelling of transient heat transfer in annealing furnaces

This is a study of heat transfer in a lift-off furnace which is employed in the batch annealing of a stack of coils of steel strip. The objective of the project is to investigate the various factors which govern the furnace design and the heat transfer resistances, so as to reduce the time of the annealing cycle, and hence minimize the operating costs. The work involved mathematical modelling of patterns of gas flow and modes of heat transfer. These models are: Heat conduction and its conjectures in the steel coils;Convective heat transfer in the plates separating the coils in the stack and in other parts of the furnace; and Radiative and convective heat transfer in the furnace by using the long furnace model. An important part of the project is the development of numerical methods and computations to solve the transient models. A limited number of temperature measurements was available from experiments on a test coil in an industrial furnace. The mathematical model agreed well with these data. The model has been used to show the following characteristics of annealing furnaces, and to suggest further developments which would lead to significant savings: - The location of the limiting temperature in a coil is nearer to the hollow core than to the outer periphery. - Thermal expansion of the steel tends to open the coils, reduces their thermal conductivity in the radial direction, and hence prolongs the annealing cycle. Increasing the tension in the coils and/or heating from the core would overcome this heat transfer resistance. - The shape and dimensions of the convective channels in the plates have significant effect on heat convection in the stack. An optimal design of a channel is shown to be of a width-to-height ratio equal to 9. - Increasing the cooling rate, by using a fluidized bed instead of the normal shell and tube exchanger, would shorten the cooling time by about 15%, but increase the temperature differential in the stack. - For a specific charge weight, a stack of different-sized coils will have a shorter annealing cycle than one of equally-sized coils, provided that production constraints allow the stacking order to be optimal. - Recycle of hot flue gases to the firing zone of the furnace would produce a. decrease in the thermal efficiency up to 30% but decreases the heating time by about 26%.

Investigating viscous fluid flow in an internal mixer using computational fluid dynamics

This thesis presents an effective methodology for the generation of a simulation which can be used to increase the understanding of viscous fluid processing equipment and aid in their development, design and optimisation. The Hampden RAPRA Torque Rheometer internal batch twin rotor mixer has been simulated with a view to establishing model accuracies, limitations, practicalities and uses. As this research progressed, via the analyses several 'snap-shot' analysis of several rotor configurations using the commercial code Polyflow, it was evident that the model was of some worth and its predictions are in good agreement with the validation experiments, however, several major restrictions were identified. These included poor element form, high man-hour requirements for the construction of each geometry and the absence of the transient term in these models. All, or at least some, of these limitations apply to the numerous attempts to model internal mixes by other researchers and it was clear that there was no generally accepted methodology to provide a practical three-dimensional model which has been adequately validated. This research, unlike others, presents a full complex three-dimensional, transient, non-isothermal, generalised non-Newtonian simulation with wall slip which overcomes these limitations using unmatched ridding and sliding mesh technology adapted from CFX codes. This method yields good element form and, since only one geometry has to be constructed to represent the entire rotor cycle, is extremely beneficial for detailed flow field analysis when used in conjunction with user defined programmes and automatic geometry parameterisation (AGP), and improves accuracy for investigating equipment design and operation conditions. Model validation has been identified as an area which has been neglected by other researchers in this field, especially for time dependent geometries, and has been rigorously pursued in terms of qualitative and quantitative velocity vector analysis of the isothermal, full fill mixing of generalised non-Newtonian fluids, as well as torque comparison, with a relatively high degree of success. This indicates that CFD models of this type can be accurate and perhaps have not been validated to this extent previously because of the inherent difficulties arising from most real processes.

Rapid aquaporin translocation regulates cellular water flow:the mechanism of hypotonicity-induced sub-cellular localization of the aquaporin 1 water channel

The control of cellular water flow is mediated by the aquaporin (AQP) family of membrane proteins. The family's structural features and the mechanism of selective water passage through the AQP pore are established, but there remains a gap in our knowledge of how water transport is regulated. Two broad possibilities exist. One is controlling the passage of water through the AQP pore, but this has only been observed as a phenomenon in some plant and microbial AQPs. An alternative is controlling the number of AQPs in the cell membrane. Here we describe a novel pathway in mammalian cells whereby a hypotonic stimulus directly induces intracellular calcium elevations, through transient receptor potential channels, that trigger AQP1 translocation. This translocation, which has a direct role in cell volume regulation, occurs within 30s and is dependent on calmodulin activation and phosphorylation of AQP1 at two threonine residues by protein kinase C. This direct mechanism provides a rationale for the changes in water transport that are required in response to constantly-changing local cellular water availability. Moreover, since calcium is a pluripotent and ubiquitous second messenger in biological systems, the discovery of its role in the regulation of AQP translocation has ramifications for diverse physiological and pathophysiological processes, as well as providing an explanation for the rapid regulation of water flow that is necessary for cell homeostasis.

A method of fundamental solutions for transient heat conduction in layered materials

In this paper we investigate an application of the method of fundamental solutions (MFS) to transient heat conduction in layered materials, where the thermal diffusivity is piecewise constant. Recently, in Johansson and Lesnic [A method of fundamental solutions for transient heat conduction. Eng Anal Boundary Elem 2008;32:697–703], a MFS was proposed with the sources placed outside the space domain of interest, and we extend that technique to numerically approximate the heat flow in layered materials. Theoretical properties of the method, as well as numerical investigations are included.

Quantifying the effects of temperature and concentration on variable-density flow in numerical modeling of groundwater systems: Implications for predictive uncertainty and data collection

Groundwater systems of different densities are often mathematically modeled to understand and predict environmental behavior such as seawater intrusion or submarine groundwater discharge. Additional data collection may be justified if it will cost-effectively aid in reducing the uncertainty of a model's prediction. The collection of salinity, as well as, temperature data could aid in reducing predictive uncertainty in a variable-density model. However, before numerical models can be created, rigorous testing of the modeling code needs to be completed. This research documents the benchmark testing of a new modeling code, SEAWAT Version 4. The benchmark problems include various combinations of density-dependent flow resulting from variations in concentration and temperature. The verified code, SEAWAT, was then applied to two different hydrological analyses to explore the capacity of a variable-density model to guide data collection. ^ The first analysis tested a linear method to guide data collection by quantifying the contribution of different data types and locations toward reducing predictive uncertainty in a nonlinear variable-density flow and transport model. The relative contributions of temperature and concentration measurements, at different locations within a simulated carbonate platform, for predicting movement of the saltwater interface were assessed. Results from the method showed that concentration data had greater worth than temperature data in reducing predictive uncertainty in this case. Results also indicated that a linear method could be used to quantify data worth in a nonlinear model. ^ The second hydrological analysis utilized a model to identify the transient response of the salinity, temperature, age, and amount of submarine groundwater discharge to changes in tidal ocean stage, seasonal temperature variations, and different types of geology. The model was compared to multiple kinds of data to (1) calibrate and verify the model, and (2) explore the potential for the model to be used to guide the collection of data using techniques such as electromagnetic resistivity, thermal imagery, and seepage meters. Results indicated that the model can be used to give insight to submarine groundwater discharge and be used to guide data collection. ^

Lattice Boltzmann modeling of fluid flow and solute transport in karst aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Lattice Boltzmann Modeling of Fluid Flow and Solute Transport in Karst Aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Thirty minutes of handgrip exercise enhances brachial artery flow-mediated dilation in response to two different shear stress profiles.

The walls of blood vessels are lined with a single-cell layer of endothelial cells. As blood flows through the arteries, a frictional force known as shear stress is sensed by mechanosensitive structures on the endothelium. Short and long term changes in shear stress can have a significant influence on the regulation of endothelial function. Acutely, shear stress triggers a pathway that culminates in the release of vasodilatory molecules from the endothelium and subsequent vasodilation of the artery. This endothelial response is known as flow mediated dilation (FMD). FMD is used as an index of endothelial function and is commonly assessed using reactive hyperemia (RH)-FMD, a method which elicits a large, short lived increase in shear stress following the release of a brief (5 min) forearm occlusion. A recent study found that a short term exposure (30 min) to a sustained elevation in shear stress potentiates subsequent RH-FMD. FMD can also result from a more prolonged, sustained increase in shear stress elicited by handgrip exercise (HGEX-FMD). There is evidence to suggest that interventions and conditions impact FMD resulting from sustained and transient shear stress stimuli differently, indicating that HGEX-FMD and RH-FMD provide different information about endothelial function. It is unknown whether HGEX-FMD is improved by short term exposure to shear stress. Understanding how exercise induced FMD is regulated is important because it contributes to blood flow responses during exercise. The study purpose was therefore to assess the impact of a handgrip exercise (intervention) induced sustained elevation in shear stress on subsequent brachial artery (BA) HGEX-FMD. Twenty healthy male participants (22±3yrs) preformed a 30-minute HGEX intervention on two experimental days. BA-FMD was assessed using either an RH or HGEX shear stress stimulus at 3 time points: pre-intervention, 10 min post and 60 min post. FMD and shear stress magnitude were determined via ultrasound. Both HGEX and RH-FMD increased significantly from pre-intervention to 10 min-post (p<0.01). These findings indicate that FMD stimulated by exercise induced increases in shear stress is potentiated by short term shear stress exposure. These findings advance our understanding regarding the regulation of endothelial function by shear stress.

Does estrogen fluctuation throughout the menstrual cycle impact flow-mediated dilation during exercise in healthy premenopausal women?

The endothelium is the inner most layer of cells that lines all arteries. A primary function of endothelial cells is to regulate responses to increased blood flow and the resulting frictional forces or shear stress by producing factors such as nitric oxide that mediate arterial dilation (flow mediated dilation (FMD)). Menstrual cycle variations in estrogen (E2) have been shown to influence brachial artery (BA) FMD in response to transient increases in shear stress brought about by the release of a brief forearm occlusion (reactive hyperemia (RH)). FMD can also be assessed in response to a sustained shear stress stimulus such as that created with handgrip exercise (HGEX), and studies have shown that RH- and HGEX stimulated FMD provide unique information regarding endothelial function. However, the impact of menstrual phase on HGEX-FMD is unknown. Therefore, the purpose of this study was to determine the impact of cyclical changes in E2 levels on HGEX-FMD over two discrete phases of the menstrual cycle. FMD was assessed via ultrasound. 12 subjects (21 ± 2yrs) completed two experimental visits: (1) low estrogen phase (early follicular) and (2) High estrogen phase (late follicular). In each visit both RH- and HGEX-FMD (6 min handgrip exercise) were assessed. Results are mean ± SD. E2 increased from the low to the high estrogen phase of the menstrual cycle (low: 34 ± 8, high: 161 ± 113pg/mL, p = 0.004). There was no change in mean FMD between phases (RH-FMD: 7.7 ± 4.3% vs. 6.4 ± 3.1%, p = 0.139; HGEX-FMD: 4.8 ± 2.8% vs. 4.8 ± 2.3%, p = 0.979). The observation that both RH- and HGEX-FMD did not differ between phases indicates that menstrual cycle fluctuations in estrogen may not universally impact endothelial function in young, healthy premenopausal women. Further research is needed to improve our understanding of the mechanisms that underlie variability in the impact of menstrual phase on both transient and sustained FMD responses.

Parallelisation of ESAUNA within EUROPORT-1: A structured/unstructured aeronautical CFD flow code

General-purpose parallel processing for solving day-to-day industrial problems has been slow to develop, partly because of the lack of suitable hardware from well-established, mainstream computer manufacturers and suitably parallelized application software. The parallelization of a CFD-(computational fluid dynamics) flow solution code is known as ESAUNA. This code is part of SAUNA, a large CFD suite aimed at computing the flow around very complex aircraft configurations including complete aircraft. A novel feature of the SAUNA suite is that it is designed to use either block-structured hexahedral grids, unstructured tetrahedral grids, or a hybrid combination of both grid types. ESAUNA is designed to solve the Euler equations or the Navier-Stokes equations, the latter in conjunction with various turbulence models. Two fundamental parallelization concepts are used—namely, grid partitioning and encapsulation of communications. Grid partitioning is applied to both block-structured grid modules and unstructured grid modules. ESAUNA can also be coupled with other simulation codes for multidisciplinary computations such as flow simulations around an aircraft coupled with flutter prediction for transient flight simulations.

Characterization of Chimney Flue Gas Flows : Flow Rate Measurements with Averaging Pitot Probes

Performance testing methods of boilers in transient operating conditions (start, stop and combustion power modulation sequences) need the combustion rate quantified to allow for the emissions to be quantified. One way of quantifying the combustion rate of a boiler during transient operating conditions is by measuring the flue gas flow rate. The flow conditions in chimneys of single family house boilers pose a challenge however, mainly because of the low flow velocity. The main objectives of the work were to characterize the flow conditions in residential chimneys, to evaluate the use of the Pitot-static method and the averaging Pitot method, and to develop and test a calibration method for averaging Pitot probes for low

Beam-energy Dependence Of The Directed Flow Of Protons, Antiprotons, And Pions In Au+au Collisions.

Rapidity-odd directed flow (v1) measurements for charged pions, protons, and antiprotons near midrapidity (y=0) are reported in sNN=7.7, 11.5, 19.6, 27, 39, 62.4, and 200 GeV Au+Au collisions as recorded by the STAR detector at the Relativistic Heavy Ion Collider. At intermediate impact parameters, the proton and net-proton slope parameter dv1/dy|y=0 shows a minimum between 11.5 and 19.6 GeV. In addition, the net-proton dv1/dy|y=0 changes sign twice between 7.7 and 39 GeV. The proton and net-proton results qualitatively resemble predictions of a hydrodynamic model with a first-order phase transition from hadronic matter to deconfined matter, and differ from hadronic transport calculations.

Communication: Transient Anion States Of Phenol…(h₂o)n (n = 1, 2) Complexes: Search For Microsolvation Signatures.

We report on the shape resonance spectra of phenol-water clusters, as obtained from elastic electron scattering calculations. Our results, along with virtual orbital analysis, indicate that the well-known indirect mechanism for hydrogen elimination in the gas phase is significantly impacted on by microsolvation, due to the competition between vibronic couplings on the solute and solvent molecules. This fact suggests how relevant the solvation effects could be for the electron-driven damage of biomolecules and the biomass delignification [E. M. de Oliveira et al., Phys. Rev. A 86, 020701(R) (2012)]. We also discuss microsolvation signatures in the differential cross sections that could help to identify the solvated complexes and access the composition of gaseous admixtures of these species.