Biofluid behaviour in 3D microchannel systems: Numerical analysis and design development of 3D microchannel biochip separators

This paper describes the design and development cycle of a 3D biochip separator and the modelling analysis of flow behaviour in the biochip microchannel features. The focus is on identifying the difference between 2D and 3D implementations as well as developing basic forms of 3D microfluidic separators. Five variants, based around the device are proposed and analysed. These include three variations of the branch channels (circular, rectangular, disc) and two variations of the main channel (solid and concentric). Ignoring the initial transient behaviour and assuming steady state flow has been established, the efficiencies of the flow between the main and side channels for the different designs are analysed and compared with regard to relevant biomicrofluidic laws or effects (bifurcation law, Fahraeus effect, cell-free phenomenon, bending channel effect and laminar flow behaviour). The modelling results identify flow features in microchannels, a constriction and bifurcations and show detailed differences in flow fields between the various designs. The manufacturing process using injection moulding for the initial base case design is also presented and discussed. The work reported here is supported as part of the UK funded 3D-MINTEGRATION project. © 2010 IEEE.

Numerical analysis of AC loss reduction in HTS superconducting coils using magnetic materials to divert flux

In this paper, the use of magnetic materials to divert flux in high-temperature superconductor superconducting coils and reduce transport ac loss is investigated. This particular technique is preferred over other techniques, such as striation, Roebel transposition, and twisted wires because it does not require modification to the conductor itself, which can be detrimental to the properties of the superconductor. The technique can also be implemented for existing coils. The analysis is carried out using a coil model based on the H formulation and implemented in comsol multiphysics. Both weakly and strongly magnetic materials are investigated, and it is shown that the use of such materials can divert flux and achieve a reduction in transport ac loss, which, in some cases, is quite significant. This analysis acts to provide a foundation for further optimization and experimental work in the future. © 2011 IEEE.

Experimental and numerical study of a YBCO pancake coil with a magnetic substrate

A finite element model for a YBCO pancake coil with a magnetic substrate is developed in this paper. An axial symmetrical H formulation and the E-J power law are used to construct the model, with the magnetic substrate considered by introducing an extra time-dependent term in the formula. A pancake coil is made and tested. The measurement of critical current and transport loss is compared to the model result, showing good consistency. The influence of magnetic substrate in the condition of AC and DC current is studied. The AC loss decreases without a magnetic substrate. It is observed that when the applied DC current approaches the critical current the coil turn loss profile changes completely in the presence of magnetic substrate due to the change of magnetic field distribution. © 2012 IOP Publishing Ltd.

Effects of turbulent reynolds number on the performance of algebraic flame surface density models for large eddy simulation in the thin reaction zones regime: A direct numerical simulation analysis

A direct numerical simulation (DNS) database of freely propagating statistically planar turbulent premixed flames with a range of different turbulent Reynolds numbers has been used to assess the performance of algebraic flame surface density (FSD) models based on a fractal representation of the flame wrinkling factor. The turbulent Reynolds number Ret has been varied by modifying the Karlovitz number Ka and the Damköhler number Da independently of each other in such a way that the flames remain within the thin reaction zones regime. It has been found that the turbulent Reynolds number and the Karlovitz number both have a significant influence on the fractal dimension, which is found to increase with increasing Ret and Ka before reaching an asymptotic value for large values of Ret and Ka. A parameterisation of the fractal dimension is presented in which the effects of the Reynolds and the Karlovitz numbers are explicitly taken into account. By contrast, the inner cut-off scale normalised by the Zel'dovich flame thickness ηi/δz does not exhibit any significant dependence on Ret for the cases considered here. The performance of several algebraic FSD models has been assessed based on various criteria. Most of the algebraic models show a deterioration in performance with increasing the LES filter width. © 2012 Mohit Katragadda et al.

Numerical analysis of thermally actuated magnets for magnetization of superconductors

Superconductors, such as YBCO bulks, have extremely high potential magnetic flux densities, comparing to rare earth magnets. Therefore, the magnetization of superconductors has attracted broad attention and contribution from both academic research and industry. In this paper, a novel technique is proposed to magnetize superconductors. Unusually, instead of using high magnetic fields and pulses, repeatedly magnetic waves with strength of as low as rare earth magnets are applied. These magnetic waves, generated by thermally controlling a Gadolinium (Gd) bulk with a rare earth magnet underneath, travel over the flat surface of a YBCO bulk and get trapped little by little. Thus, a very small magnetic field can be used to build up a very large magnetic field. In this paper, the modelling results of thermally actuated magnetic waves are presented showing how to transfer sequentially applied thermal pulses into magnetic waves. The experiment results of the magnetization of YBCO bulk are also presented to demonstrate how superconductors are progressively magnetized by small magnetic field © 2010 IOP Publishing Ltd.

Boundary integral technique for the numerical modelling of the air flow for the Aaberg exhaust system

A boundary integral technique has been developed for the numerical simulation of the air flow for the Aaberg exhaust system. For the steady, ideal, irrotational air flow induced by a jet, the air velocity is an analytical function. The solution of the problem is formulated in the form of a boundary integral equation by seeking the solution of a mixed boundary-value problem of an analytical function based on the Riemann-Hilbert technique. The boundary integral equation is numerically solved by converting it into a system of linear algebraic equations, which are solved by the process of the Gaussian elimination. The air velocity vector at any point in the solution domain is then computed from the air velocity on the boundary of the solution domains.

Experimental and numerical investigation of fracture mechanics of bitumen beams

This work presents an experimental and numerical investigation to characterise the fracture properties of pure bitumen (the binder in asphalt paving materials). The paper is divided into two parts. The first part describes an experimental study of fracture characterisation parameters of pure bitumen as determined by three-point bend tests. The second part deals with modelling of fracture and failure of bitumen by Finite Element analysis. Fracture mechanics parameters, stress intensity factor, KIC, fracture energy, GIC, and J-integral, JIC, are used for evaluation of bitumen's fracture properties. The material constitutive model developed by Ossa et al. [40,41] which was implemented into a FE code by Costanzi [18] is combined with cohesive zone models (CZM) to simulate the fracture behaviour of pure bitumen. Experimental and numerical results are presented in the form of failure mechanism maps where ductile, brittle and brittle-ductile transition regimes of fracture behaviour are classified. The FE predictions of fracture behaviour match well with experimental results. © 2012 Elsevier Ltd.

Numerical simulation of single-stream jets from a serrated nozzle

Hybrid large-eddy type simulations for cold jet flows from a serrated nozzle are performed at an acoustic Mach number Ma ac = 0.9 and Re = 1.03×10 6. Since the solver being used tends towards having dissipative qualities, the subgrid scale (SGS) model is omitted, giving a numerical type LES (NLES) or implicit LES (ILES) reminiscent procedure. To overcome near wall streak resolution problems a near wall RANS (Reynolds averaged Navier-Stokes) model is smoothly blended to the LES making a hybrid RANS-ILES. The geometric complexity of the serrated nozzle is fully considered without simplification or emulation. An improved but still modest hexahedral multi-block grid with circa 20 million grid points (with respect to 12.5 million in Xia et al.; Int J Heat Fluid Flow 30:1067-1079, 2009) is used. Despite the modest grid size, encouraging and improved results are obtained. Directly resolved mean and second-order fluctuating quantities along the jet centerline and in the jet shear layer compare favorably with measurements. The radiated far-field sound predicted using the Ffowcs Williams and Hawkings (FW-H) surface integral method shows good agreement with the measurements in directivity and sound spectra. © 2011 Springer Science+Business Media B.V.

Numerical computation of laminar flow in a heated rotating cavity with an axial throughflow of air

A heated rotating cavity with an axial throughflow of cooling air is used as a model for the flow in the cylindrical cavities between adjacent discs of a high-pressure gas-turbine compressor. In an engine the flow is expected to be turbulent, the limitations of this laminar study are fully realised but it is considered an essential step to understand the fundamental nature of the flow. The three-dimensional, time-dependent governing equations are solved using a code based on the finite volume technique and a multigrid algorithm. The computed flow structure shows that flow enters the cavity in one or more radial arms and then forms regions of cyclonic and anticyclonic circulation. This basic flow structure is consistent with existing experimental evidence obtained from flow visualization. The flow structure also undergoes cyclic changes with time. For example, a single radial arm, and pair of recirculation regions can commute to two radial arms and two pairs of recirculation regions and then revert back to one. The flow structure inside the cavity is found to be heavily influenced by the radial distribution of surface temperature imposed on the discs. As the radial location of the maximum disc temperature moves radially outward, this appears to increase the number of radial arms and pairs of recirculation regions (from one to three for the distributions considered here). If the peripheral shroud is also heated there appear to be many radial arms which exchange fluid with a strong cyclonic flow adjacent to the shroud. One surface temperature distribution is studied in detail and profiles of the relative tangential and radial velocities are presented. The disc heat transfer is also found to be influenced by the disc surface temperature distribution. It is also found that the computed Nusselt numbers are in reasonable accord over most of the disc surface with a correlation found from previous experimental measurements. © 1994, MCB UP Limited.

Numerical study of soil heterogeneity effects on contaminant transport in unsaturated soil (model development and validation)

The movement of chemicals through soil to groundwater is a major cause of degradation of water resources. In many cases, serious human and stock health implications are associated with this form of pollution. The study of the effects of different factors involved in transport phenomena can provide valuable information to find the best remediation approaches. Numerical models are increasingly being used for predicting or analyzing solute transport processes in soils and groundwater. This article presents the development of a stochastic finite element model for the simulation of contaminant transport through soils with the main focus being on the incorporation of the effects of soil heterogeneity in the model. The governing equations of contaminant transport are presented. The mathematical framework and the numerical implementation of the model are described. The comparison of the results obtained from the developed stochastic model with those obtained from a deterministic method and some experimental results shows that the stochastic model is capable of predicting the transport of solutes in unsaturated soil with higher accuracy than deterministic one. The importance of the consideration of the effects of soil heterogeneity on contaminant fate is highlighted through a sensitivity analysis regarding the variance of saturated hydraulic conductivity as an index of soil heterogeneity. © 2011 John Wiley & Sons, Ltd.