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.

Porohyperelastic analysis to explore mechanical properties of chondrocytes using numerical modeling and experiments : a finite element study

There are many continuum mechanical models have been developed such as liquid drop models, solid models, and so on for single living cell biomechanics studies. However, these models do not give a fully approach to exhibit a clear understanding of the behaviour of single living cells such as swelling behaviour, drag effect, etc. Hence, the porohyperelastic (PHE) model which can capture those aspects would be a good candidature to study cells behaviour (e.g. chondrocytes in this study). In this research, an FEM model of single chondrocyte cell will be developed by using this PHE model to simulate Atomic Force Microscopy (AFM) experimental results with the variation of strain rate. This material model will be compared with viscoelastic model to demonstrate the advantages of PHE model. The results have shown that the maximum value of force applied of PHE model is lower at lower strain rates. This is because the mobile fluid does not have enough time to exude in case of very high strain rate and also due to the lower permeability of the membrane than that of the protoplasm of chondrocyte. This behavior is barely observed in viscoelastic model. Thus, PHE model is the better model for cell biomechanics studies.

Finite element analysis of the in-plane shear behaviour of masonry panels confined with reinforced grouted cores

A combined experimental and numerical program was conducted to study the in-plane shear behaviour of hollow concrete masonry panels containing reinforced grout cores. This paper is focused on the numerical program. A two dimensional macromodelling strategy was used to simulate the behaviour of the confined masonry (CM) shear panels. Both the unreinforced masonry and the confining element were modelled using macromasonry properties and the steel reinforcement was modelled as an embedded truss element located within the grout using perfectly bonded constraint. The FE model reproduced key behaviours observed in the experiments, including the shear strength, the deformation and the crack patterns of the unconfined and confined masonry panels. The predictions of the validated model were used to evaluate the existing in-plane shear expressions available in the national masonry standards and research publications.

Numerical modelling and analysis of the blast performance of laminated glass panels and the influence of material parameters

This paper presents a comprehensive numerical procedure to treat the blast response of laminated glass (LG) panels and studies the influence of important material parameters. Post-crack behaviour of the LG panel and the contribution of the interlayer towards blast resistance are treated. Modelling techniques are validated by comparing with existing experimental results. Findings indicate that the tensile strength of glass considerably influences the blast response of LG panels while the interlayer material properties have a major impact on the response under higher blast loads. Initially, glass panes absorb most of the blast energy, but after the glass breaks, interlayer deforms further and absorbs most of the blast energy. LG panels should be designed to fail by tearing of the interlayer rather than failure at the supports to achieve a desired level of protection. From this aspect, material properties of glass, interlayer and sealant joints play important roles, but unfortunately they are not accounted for in the current design standards. The new information generated in this paper will enhance the capabilities of engineers to better design LG panels under blast loads and use better materials to improve the blast response of LG panels.

Numerical analysis of geocell reinforced ballast overlying soft clay subgrade

Geotextiles and geogrids have been in use for several decades in variety of geo-structure applications including foundation of embankments, retaining walls, pavements. Geocells is one such variant in geosynthetic reinforcement of recent years, which provides a three dimensional confinement to the infill material. Although extensive research has been carried on geocell reinforced sand, clay and layered soil subgrades, limited research has been reported on the aggregates/ballast reinforced with geocells. This paper presents the behavior of a railway sleeper subjected to monotonic loading on geocell reinforced aggregates, of size ranging from 20 to 75 mm, overlying soft clay subgrades. Series of tests were conducted in a steel test tank of dimensions 700 mm x 300 mm x 700 mm. In addition to the laboratory model tests, numerical simulations were performed using a finite difference code to predict the behavior of geocell reinforced ballast. The results from numerical simulations were compared with the experimental data. The numerical and experimental results manifested the importance that the geocell reinforcement has a significant effect on the ballast behaviour. The results depicted that the stiffness of underlying soft clay subgrade has a significant influence on the behavior of the geocell-aggregate composite material in redistributing the loading system.

Three-dimensional numerical analysis of geocell-reinforced soft clay beds by considering the actual geometry of geocell pockets

Due to its complex honeycomb structure, the numerical modeling of the geocell has always been a big challenge. Generally, the equivalent composite approach is used to model the geocells. In the equivalent composite approach, the geocellsoil composite is treated as the soil layer with improved strength and stiffness values. Though this approach is very simple, it is unrealistic to model the geocells as the soil layer. This paper presents a more realistic approach of modeling the geocells in three-dimensional (3D) framework by considering the actual curvature of the geocell pocket. A square footing resting on geocell reinforced soft clay bed was modeled using the ``fast Lagrangian analysis of continua in 3D'' (FLAC(3D)) finite difference package. Three different material models, namely modified Cam-clay, Mohr-Coulomb, and linear elastic were used to simulate the behaviour of foundation soil, infill soil and the geocell, respectively. It was found that the geocells distribute the load laterally to the wider area below the footing as compared to the unreinforced case. More than 50% reduction in the stress was observed in the clay bed in the presence of geocells. In addition to geocells, two other cases, namely, only geogrid and geocell with additional basal geogrid cases were also simulated. The numerical model was systematically validated with the results of the physical model tests. Using the validated numerical model, parametric studies were conducted to evaluate the influence of various geocell properties on the performance of reinforced clay beds.

A theoretical analysis on stochastic behaviour of short cracks and fatigue damage of metals

In this paper, the closed form of solution to the stochastic differential equation for a fatigue crack evolution system is derived. and the relationship between metal fatigue damage and crack stochastic behaviour is investigated. It is found that the damage extent of metals is independent of crack stochastic behaviour ii the stochastic deviation of the crack growth rate is directly proportional to its mean value. The evolution of stochastic deviation of metal fatigue damage in the stage close to the transition point between short and long crack regimes is also discussed.

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.

An investigation into numerical analysis alternatives for predicting reingestion in turbine disc rim cavities

Reliable means of predicting ingestion in cavities adjacent to the main gas path are increasingly being sought by engineers involved in the design of gas turbines. In this paper, analysis is to be presented that results from an extended research programme, MAGPI, sponsored by the EU and several leading gas turbine manufactures and universities. Extensive use is made of CFD modelling techniques to understand the aerodynamic behaviour of a turbine stator well cavity, focusing on the interaction of cooling air supply with the main annulus gas. The objective of the study has been to benchmark a number of CFD codes and numerical techniques covering RANS and URANS calculations with different turbulence models in order to assess the suitability of the standard settings used in the industry for calculating the mechanics of the flow travelling between cavities in a turbine through the main gas path. The modelling methods employed have been compared making use of experimental data gathered from a dedicated two-stage turbine rig, running at engine representative conditions. Extensive measurements are available for a range of flow conditions and alternative cooling arrangements. The limitations of the numerical methods in calculating the interaction of the cooling flow egress and the main stream gas, and subsequent ingestion into downstream cavities in the engine (i.e. re-ingestion), have been exposed. This has been done without losing sight of the validation of the CFD for its use for predicting heat transfer, which was the main objective of the partners of the MAGPI Work- Package 1 consortium. Copyright © 2012 by ASME.

Modelling and analysis of modal behaviour in SOI slot waveguides

Mode behaviour for SOI slot waveguides is modelled and analysed using a numerical full vectorial method based on the film mode matching method (MMM). Only the quasi-TE mode is investigated. Waveguide heights and slot widths, as well as silicon widths are properly chosen with respect to the single mode behaviour in the slot region. Comparison between the effective index method and our side loss method shows that our single mode condition is creditable. The optical power confinement in slot region for the quasi-TE mode is also studied and presented. We demonstrate that the maximum achievable optical power confinement P-slot and the maximum normalized average optical intensity I-slot are 42% and 26 mu m(-2), respectively.

Numerical analysis of intralaminar failure mechanisms in composite structures, Part I: FE implementation

This paper presents a three-dimensional continuum damage mechanics-based material model which was implemented in an implicit finite element code to simulate the progressive intralaminar degradation of fibre reinforced laminates. The damage model is based on ply failure mechanisms and uses seven damage variables assigned to tensile, compressive and shear damage at a ply level. Non-linear behaviour and irreversibility were taken into account and modelled. Some issues on the numerical implementation of the damage model are discussed and solutions proposed. Applications of the methodology are presented in Part II

Numerical analysis of stiffener runout sections

The recent trend of incorporating more composite material in primary aircraft structures has highlighted the vulnerability of stiffened aerostructures to through-thickness stresses, which may lead to delamination and debonding at the skin-stiffener interface, leading to collapse. Stiffener runout regions are particularly susceptible to this problem and cannot be avoided due to the necessity to terminate stiffeners at rib intersections or at cutouts, interrupting the stiffener load path. In this paper, experimental tests relating to two different stiffener runout specimens are presented and the failure modes of both specimens are discussed in detail. A thinner-skinned specimen showed sudden and unstable crack propagation, while a thicker-skinned specimen showed initially unstable but subsequent stable crack growth. Detailed finite element models of the two specimens are developed, and it is shown how such models can explain and predict the behaviour and failure mode of stiffener runouts. The models contain continuum shell elements to model the skin and stiffener, while cohesive elements using a traction-separation law are placed at the skin-stiffener interface to effectively model the debonding which promotes structural failure.

Shear behaviour of reinforced cconcrete slab under concentrated load: an investigation through non-linear and sequentially linear analysis

English: The assessment of safety in existing bridges and viaducts led the Ministry of Public Works of the Netherlands to finance a specific campaing aimed at the study of the response of the elements of these infrastructures. Therefore, this activity is focused on the investigation of the behaviour of reinforced concrete slabs under concentrated loads, adopting finite element modeling and comparison with experimental results. These elements are characterized by shear behaviour and crisi, whose modeling is, from a computational point of view, a hard challeng, due to the brittle behavior combined with three-dimensional effects. The numerical modeling of the failure is studied through Sequentially Linear Analysis (SLA), an alternative Finite Element method, with respect to traditional incremental and iterative approaches. The comparison between the two different numerical techniques represents one of the first works and comparisons in a three-dimensional environment. It's carried out adopting one of the experimental test executed on reinforced concrete slabs as well. The advantage of the SLA is to avoid the well known problems of convergence of typical non-linear analysis, by directly specifying a damage increment, in terms of reduction of stiffness and resistance in particular finite element, instead of load or displacement increasing on the whole structure . For the first time, particular attention has been paid to specific aspects of the slabs, like an accurate constraints modeling and sensitivity of the solution with respect to the mesh density. This detailed analysis with respect to the main parameters proofed a strong influence of the tensile fracture energy, mesh density and chosen model on the solution in terms of force-displacement diagram, distribution of the crack patterns and shear failure mode. The SLA showed a great potential, but it requires a further developments for what regards two aspects of modeling: load conditions (constant and proportional loads) and softening behaviour of brittle materials (like concrete) in the three-dimensional field, in order to widen its horizons in these new contexts of study.

"Offshore tower or platform foundations: numerical analysis of a laterally loaded single pile or pile group in soft clay and analysis of actions on a jacket structure"

Laterally loaded piles are a typical situation for a large number of cases in which deep foundations are used. Dissertation herein reported, is a focus upon the numerical simulation of laterally loaded piles. In the first chapter the best model settings are largely discussed, so a clear idea about the effects of interface adoption, model dimension, refinement cluster and mesh coarseness is reached. At a second stage, there are three distinct parametric analyses, in which the model response sensibility is studied for variation of interface reduction factor, Eps50 and tensile cut-off. In addition, the adoption of an advanced soil model is analysed (NGI-ADP). This was done in order to use the complex behaviour (different undrained shear strengths are involved) that governs the resisting process of clay under short time static loads. Once set a definitive model, a series of analyses has been carried out with the objective of defining the resistance-deflection (P-y) curves for Plaxis3D (2013) data. Major results of a large number of comparisons made with curves from API (America Petroleum Institute) recommendation are that the empirical curves have almost the same ultimate resistance but a bigger initial stiffness. In the second part of the thesis a simplified structural preliminary design of a jacket structure has been carried out to evaluate the environmental forces that act on it and on its piles foundation. Finally, pile lateral response is studied using the empirical curves.

Analysis and numerical solution of the Peterlin viscoelastic model

Liquids and gasses form a vital part of nature. Many of these are complex fluids with non-Newtonian behaviour. We introduce a mathematical model describing the unsteady motion of an incompressible polymeric fluid. Each polymer molecule is treated as two beads connected by a spring. For the nonlinear spring force it is not possible to obtain a closed system of equations, unless we approximate the force law. The Peterlin approximation replaces the length of the spring by the length of the average spring. Consequently, the macroscopic dumbbell-based model for dilute polymer solutions is obtained. The model consists of the conservation of mass and momentum and time evolution of the symmetric positive definite conformation tensor, where the diffusive effects are taken into account. In two space dimensions we prove global in time existence of weak solutions. Assuming more regular data we show higher regularity and consequently uniqueness of the weak solution. For the Oseen-type Peterlin model we propose a linear pressure-stabilized characteristics finite element scheme. We derive the corresponding error estimates and we prove, for linear finite elements, the optimal first order accuracy. Theoretical error of the pressure-stabilized characteristic finite element scheme is confirmed by a series of numerical experiments.