A new Lagrangian-Eulerian shell-fluid coupling algorithm based on level sets

We propose a computational method for the coupled simulation of a compressible flow interacting with a thin-shell structure undergoing large deformations. An Eulerian finite volume formulation is adopted for the fluid and a Lagrangian formulation based on subdivision finite elements is adopted for the shell response. The coupling between the fluid and the solid response is achieved via a novel approach based on level sets. The basic approach furnishes a general algorithm for coupling Lagrangian shell solvers with Cartesian grid based Eulerian fluid solvers. The efficiency and robustness of the proposed approach is demonstrated with a airbag deployment simulation. It bears emphasis that in the proposed approach the solid and the fluid components as well as their coupled interaction are considered in full detail and modeled with an equivalent level of fidelity without any oversimplifying assumptions or bias towards a particular physical aspect of the problem.

Techniques to improve the shear strength of impacted bone graft: the effect of particle size and washing of the graft.

BACKGROUND: When fresh morselized graft is compacted, as in impaction bone-grafting for revision hip surgery, fat and marrow fluid is either exuded or trapped in the voids between particles. We hypothesized that the presence of incompressible fluid damps and resists compressive forces during impaction and prevents the graft particles from moving into a closer formation, thus reducing the graft strength. In addition, viscous fluid such as fat may act as an interparticle lubricant, thus reducing the interlocking of the particles. METHODS: We performed mechanical shear testing in the laboratory with use of fresh-frozen human femoral-head allografts that had been passed through different orthopaedic bone mills to produce graft of differing particle-size distributions (grading). RESULTS: After compaction of fresh graft, fat and marrow fluid continued to escape on application of normal loads. Washed graft, however, had little lubricating fluid and better contact between the particles, increasing the shear resistance. On mechanical testing, washed graft was significantly (p < 0.001) more resistant to shearing forces than fresh graft was. This feature was consistent for different bone mills that produced graft of different particle-size distributions and shear strengths. CONCLUSIONS: Removal of fat and marrow fluid from milled human allograft by washing the graft allows the production of stronger compacted graft that is more resistant to shear, which is the usual mode of failure. Further research into the optimum grading of particle sizes from bone mills is required.

Techniques to improve the shear strength of impacted bone graft: The effect of particle size and washing of the graft

Background: When fresh morselized graft is compacted, as in impaction bone-grafting for revision hip surgery, fat and marrow fluid is either exuded or trapped in the voids between particles. We hypothesized that the presence of incompressible fluid damps and resists compressive forces during impaction and prevents the graft particles from moving into a closer formation, thus reducing the graft strength. In addition, viscous fluid such as fat may act as an interparticle lubricant, thus reducing the interlocking of the particles. Methods: We performed mechanical shear testing in the laboratory with use of fresh-frozen human femoral-head allografts that had been passed through different orthopaedic bone mills to produce graft of differing particle-size distributions (grading). Results: After compaction of fresh graft, fat and marrow fluid continued to escape on application of normal loads. Washed graft, however, had little lubricating fluid and better contact between the particles, increasing the shear resistance. On mechanical testing, washed graft was significantly (p < 0.001) more resistant to shearing forces than fresh graft was. This feature was consistent for different bone mills that produced graft of different particle-size distributions and shear strengths. Conclusions: Removal of fat and marrow fluid from milled human allograft by washing the graft allows the production of stronger compacted graft that is more resistant to shear, which is the usual mode of failure. Further research into the optimum grading of particle sizes from bone mills is required. Clinical Relevance: Understanding the mechanical properties of milled human allograft is important when impaction grafting is used for mechanical support. A simple means of improving the mechanical strength of graft produced by currently available bone mills, including an intraoperative washing technique, is described.

A fixed-grid b-spline finite element technique for fluid-structure interaction

We present a fixed-grid finite element technique for fluid-structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b-spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision-stabilisation technique is used to ensure inf-sup stability. The beam equations are discretised with b-splines and the shell equations with subdivision basis functions, both leading to a rotation-free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet-Robin partitioning scheme, and the fluid equations are solved with a pressure-correction method. Auxiliary techniques employed for improving numerical robustness include the level-set based implicit representation of the structure interface on the fluid grid, a cut-cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. © 2013 John Wiley & Sons, Ltd.

Finite element solution of Navier Stokes equations adapted to a priori error estimates

The paper is based on qualitative properties of the solution of the Navier-Stokes equations for incompressible fluid, and on properties of their finite element solution. In problems with corner-like singularities (e.g. on the well-known L-shaped domain) usually some adaptive strategy is used. In this paper we present an alternative approach. For flow problems on domains with corner singularities we use the a priori error estimates and asymptotic expansion of the solution to derive an algorithm for refining the mesh near the corner. It gives very precise solution in a cheap way. We present some numerical results.

ANALYTICAL DEVICE FOR TEST FLUIDS, COMPRISES A CONDUCTIVE POLYMER COMPOSITION

A sensor for chemical species or biological species or radiation presenting to test fluid a polymer composition comprises polymer and conductive filler metal, alloy or reduced metal oxide and having a first level of electrical conductance when quiescent and being convertible to a second level of conductance by change of stress applied by stretching or compression or electric field, in which the polymer composition is characterised by at least one of the features in the form of particles at least 90% w/w held on a 100 mesh sieve; and/or comprising a permeable body extending across a channel of fluid flow; and/or affording in-and-out diffusion of test fluid and/or mechanically coupled to a workpiece of polymer swellable by a constituent of test fluid.

Velocity profiles in a cylindrical liquid jet by reconstructed velocimetry

An experimental setup and a simple reconstruction method are presented to measure velocity fields inside slightly tapering cylindrical liquid jets traveling through still air. Particle image velocimetry algorithms are used to calculate velocity fields from high speed images of jets of transparent liquid containing seed particles. An inner central plane is illuminated by a laser sheet pointed at the center of the jet and visualized through the jet by a high speed camera. Optical distortions produced by the shape of the jet and the difference between the refractive index of the fluid and the surrounding air are corrected by using a ray tracing method. The effect of the jet speed on the velocity fields is investigated at four jet speeds. The relaxation rate for the velocity profile downstream of the nozzle exit is reasonably consistent with theoretical expectations for the low Reynolds numbers and the fluid used, although the velocity profiles are considerably flatter than expected. © 2012 American Society of Mechanical Engineers.

A comparison of the value of viscosity for several water models using Poiseuille flow in a nano-channel.

The viscosity-temperature relation is determined for the water models SPC/E, TIP4P, TIP4P/Ew, and TIP4P/2005 by considering Poiseuille flow inside a nano-channel using molecular dynamics. The viscosity is determined by fitting the resulting velocity profile (away from the walls) to the continuum solution for a Newtonian fluid and then compared to experimental values. The results show that the TIP4P/2005 model gives the best prediction of the viscosity for the complete range of temperatures for liquid water, and thus it is the preferred water model of these considered here for simulations where the magnitude of viscosity is crucial. On the other hand, with the TIP4P model, the viscosity is severely underpredicted, and overall the model performed worst, whereas the SPC/E and TIP4P/Ew models perform moderately.

Thylakoid terminal oxidases are essential for the cyanobacterium Synechocystis sp. PCC 6803 to survive rapidly changing light intensities.

Cyanobacteria perform photosynthesis and respiration in the thylakoid membrane, suggesting that the two processes are interlinked. However, the role of the respiratory electron transfer chain under natural environmental conditions has not been established. Through targeted gene disruption, mutants of Synechocystis sp. PCC 6803 were generated that lacked combinations of the three terminal oxidases: the thylakoid membrane-localized cytochrome c oxidase (COX) and quinol oxidase (Cyd) and the cytoplasmic membrane-localized alternative respiratory terminal oxidase. All strains demonstrated similar growth under continuous moderate or high light or 12-h moderate-light/dark square-wave cycles. However, under 12-h high-light/dark square-wave cycles, the COX/Cyd mutant displayed impaired growth and was completely photobleached after approximately 2 d. In contrast, use of sinusoidal light/dark cycles to simulate natural diurnal conditions resulted in little photobleaching, although growth was slower. Under high-light/dark square-wave cycles, the COX/Cyd mutant suffered a significant loss of photosynthetic efficiency during dark periods, a greater level of oxidative stress, and reduced glycogen degradation compared with the wild type. The mutant was susceptible to photoinhibition under pulsing but not constant light. These findings confirm a role for thylakoid-localized terminal oxidases in efficient dark respiration, reduction of oxidative stress, and accommodation of sudden light changes, demonstrating the strong selective pressure to maintain linked photosynthetic and respiratory electron chains within the thylakoid membrane. To our knowledge, this study is the first to report a phenotypic difference in growth between terminal oxidase mutants and wild-type cells and highlights the need to examine mutant phenotypes under a range of conditions.

Winglets for Improved Aerothermal Performance of High Pressure Turbines

This paper investigates the design of winglet tips for unshrouded high pressure turbine rotors, considering aerodynamic and thermal performance simultaneously. A novel parameterization method has been developed to alter the tip geometry of a rotor blade. A design survey of un-cooled, flat-tipped winglets is performed using RANS calculations for a single rotor at engine representative operating conditions. Compared to a plain tip, large efficiency gains can be realized by employing an overhang around the full perimeter of the blade, but the overall heat load rises significantly. By employing an overhang on only the early suction surface, significant efficiency improvements can be obtained without increasing the overall heat transfer to the blade. The flow physics are explored in detail to explain the results. For a plain tip, the leakage and passage vortices interact to create a three-dimensional impingement onto the blade suction surface, causing high heat transfer. The addition of an overhang on the early suction surface displaces the tip leakage vortex away from the blade, weakening the impingement effect and reducing the heat transfer on the blade. The winglets reduce the aerodynamic losses by unloading the tip section, reducing the leakage flow rate, turning the leakage flow in a more streamwise direction and reducing the interaction between the leakage fluid and endwall flows. Generally these effects are most effective close to the leading edge of the tip, where the leakage flow is subsonic.

Initiation of immersed granular avalanches.

By means of coupled molecular dynamics-computational fluid dynamics simulations, we analyze the initiation of avalanches in a granular bed of spherical particles immersed in a viscous fluid and inclined above its angle of repose. In quantitative agreement with experiments, we find that the bed is unstable for a packing fraction below 0.59 but is stabilized above this packing fraction by negative excess pore pressure induced by the effect of dilatancy. From detailed numerical data, we explore the time evolution of shear strain, packing fraction, excess pore pressures, and granular microstructure in this creeplike pressure redistribution regime, and we show that they scale excellently with a characteristic time extracted from a model based on the balance of granular stresses in the presence of a negative excess pressure and its interplay with dilatancy. The cumulative shear strain at failure is found to be ≃ 0.2, in close agreement with the experiments, irrespective of the initial packing fraction and inclination angle. Remarkably, the avalanche is triggered when dilatancy vanishes instantly as a result of fluctuations while the average dilatancy is still positive (expanding bed) with a packing fraction that declines with the initial packing fraction. Another nontrivial feature of this creeplike regime is that, in contrast to dry granular materials, the internal friction angle of the bed at failure is independent of dilatancy but depends on the inclination angle, leading therefore to a nonlinear dependence of the excess pore pressure on the inclination angle. We show that this behavior may be described in terms of the contact network anisotropy, which increases with a nearly constant connectivity and levels off at a value (critical state) that increases with the inclination angle. These features suggest that the behavior of immersed granular materials is controlled not only directly by hydrodynamic forces acting on the particles but also by the influence of the fluid on the granular microstructure.