Numerical investigation of the active role of the actin cytoskeleton in the compression resistance of cells.

Numerous in-vitro studies have established that cells react to their physical environment and to applied mechanical loading. However, the mechanisms underlying such phenomena are poorly understood. Previous modelling of cell compression considered the cell as a passive homogenous material, requiring an artificial increase in the stiffness of spread cells to replicate experimentally measured forces. In this study, we implement a fully 3D active constitutive formulation that predicts the distribution, remodelling, and contractile behaviour of the cytoskeleton. Simulations reveal that polarised and axisymmetric spread cells contain stress fibres which form dominant bundles that are stretched during compression. These dominant fibres exert tension; causing an increase in computed compression forces compared to round cells. In contrast, fewer stress fibres are computed for round cells and a lower resistance to compression is predicted. The effect of different levels of cellular contractility associated with different cell phenotypes is also investigated. Highly contractile cells form more dominant circumferential stress fibres and hence provide greater resistance to compression. Computed predictions correlate strongly with published experimentally observed trends of compression resistance as a function of cellular contractility and offer an insight into the link between cell geometry, stress fibre distribution and contractility, and cell deformability. Importantly, it is possible to capture the behaviour of both round and spread cells using a given, unchanged set of material parameters for each cell type. Finally, it is demonstrated that stress distributions in the cell cytoplasm and nucleus computed using the active formulation differ significantly from those computed using passive material models.

Experimental and numerical study of grout injections in silty soils

Compensation grouting is increasingly employed as a mitigation technique of settlements induced by tunnelling and its effectiveness both in clayey and sandy soils is reported in a wide number of case histories. However, the results are highly dependent on grout properties, injection characteristics and soil properties. An experimental study was conducted to investigate the parameters that control grout injections in silty soils. The results from one injection test in a large sample of silty soil show that the compensation efficiency, defined as the ratio of the volume of heave obtained at ground surface and the injected grout volume, is much lower than one and tends to decrease with time, while the initial volume of grout lost due to pressure filtration is small. Finally, results from finite elements back analyses of the laboratory test show that a good agreement with the experimental data can be obtained if the development of large strains is taken into account. © 2012 Taylor & Francis Group.

Numerical simulation of long-term twin-tunnel behaviour at St James's park

The twin-tunnel construction of the Jubilee Line Extension tunnels beneath St James's Park was simulated using coupled-consolidation finite-element analyses. The effect of defining different permeabilities for the final consolidation stage was investigated, and the performance of a fissure softening model was also evaluated. The analyses suggested an unexpectedly high permeability anisotropy for soil around the tunnel crown, possibly due to stress-induced permeability changes, or low-permeability laminations. Also, the permeability profile and lining conductivity were found to differ between the tunnels. Inclusion of the fissure model gave a narrower settlement trough, more alike that in the field, by preferentially softening simple shear behaviour. Long-term settlements at the site continue to increase at an unexpectedly high rate, suggesting the possibility of creep or unexpected soil softening during excavation. © 2012 Taylor & Francis Group.

Progress towards the design and numerical analysis of a 3D microchannel biochip separator

This paper reports the design and numerical analysis of a three-dimensional biochip plasma blood separator using computational fluid dynamics techniques. Based on the initial configuration of a two-dimensional (2D) separator, five three-dimensional (3D) microchannel biochip designs are categorically developed through axial and plenary symmetrical expansions. These include the geometric variations of three types of the branch side channels (circular, rectangular, disc) and two types of the main channel (solid and concentric). Ignoring the initial transient behaviour and assuming that steady-state flow has been established, the behaviour of the blood fluid in the devices is algebraically analysed and numerically modelled. The roles of the relevant microchannel mechanisms, i.e. bifurcation, constriction and bending channel, on promoting the separation process are analysed based on modelling results. The differences among the different 3D implementations are compared and discussed. The advantages of 3D over 2D separator in increasing separation volume and effectively depleting cell-free layer fluid from the whole cross section circumference are addressed and illustrated. © 2011 John Wiley & Sons, Ltd.

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.