Award Winner in the Young Investigator Category, 2014 Society for Biomaterials Annual Meeting and Exposition, Denver, Colorado, April 16-19, 2014: Periodically perforated core-shell collagen biomaterials balance cell infiltration, bioactivity, and mechanical properties.

Orthopedic tissue engineering requires biomaterials with robust mechanics as well as adequate porosity and permeability to support cell motility, proliferation, and new extracellular matrix (ECM) synthesis. While collagen-glycosaminoglycan (CG) scaffolds have been developed for a range of tissue engineering applications, they exhibit poor mechanical properties. Building on previous work in our lab that described composite CG biomaterials containing a porous scaffold core and nonporous CG membrane shell inspired by mechanically efficient core-shell composites in nature, this study explores an approach to improve cellular infiltration and metabolic health within these core-shell composites. We use indentation analyses to demonstrate that CG membranes, while less permeable than porous CG scaffolds, show similar permeability to dense materials such as small intestine submucosa (SIS). We also describe a simple method to fabricate CG membranes with organized arrays of microscale perforations. We demonstrate that perforated membranes support improved tenocyte migration into CG scaffolds, and that migration is enhanced by platelet-derived growth factor BB-mediated chemotaxis. CG core-shell composites fabricated with perforated membranes display scaffold-membrane integration with significantly improved tensile properties compared to scaffolds without membrane shells. Finally, we show that perforated membrane-scaffold composites support sustained tenocyte metabolic activity as well as improved cell infiltration and reduced expression of hypoxia-inducible factor 1α compared to composites with nonperforated membranes. These results will guide the design of improved biomaterials for tendon repair that are mechanically competent while also supporting infiltration of exogenous cells and other extrinsic mediators of wound healing.

Computational engineering design for micro-scale combustion devices: A thermally improved configuration

A multi-objective design optimisation study has been carried out with the objectives to improve the overall efficiency of the device and to reduce the fuel consumption for the proposed micro-scale combustor design configuration. In a previous study we identified the topology of the combustion chamber that produced improved behaviour of the device in terms of the above design criteria. We now extend our design approach, and we propose a new configuration by the addition of a micro-cooling channel that will improve the thermal behaviour of the design as previously suggested in literature. Our initial numerical results revealed an improvement of 2.6% in the combustion efficiency when we applied the micro-cooling channel to an optimum design configuration we identified from our earlier multi-objective optimisation study, and under the same operating conditions. The computational modelling of the combustion process is implemented in the commercial computational fluid dynamics package ANSYS-CFX using Finite Rate Chemistry and a single step hydrogen-air reaction. With this model we try to balance good accuracy of the combustion solution and at the same time practicality within the context of an optimisation process. The whole design system comprises also the ANSYS-ICEM CFD package for the automatic geometry and mesh generation and the Multi-Objective Tabu Search algorithm for the design space exploration. We model the design problem with 5 geometrical parameters and 3 operational parameters subject to 5 design constraints that secure practicality and feasibility of the new optimum design configurations. The final results demonstrate the reliability and efficiency of the developed computational design system and most importantly we assess the practicality and manufacturability of the revealed optimum design configurations of micro-combustor devices. Copyright © 2013 by ASME.

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.

Carbon source/sink function of a subtropical, eutrophic lake determined from an overall mass balance and a gas exchange and carbon burial balance

Although studies on carbon burial in lake sediments have shown that lakes are disproportionately important carbon sinks, many studies on gaseous carbon exchange across the water-air interface have demonstrated that lakes are supersaturated with CO2 and CH4 causing a net release of CO2 and CH4 to the atmosphere. In order to more accurately estimate the net carbon source/sink function of lake ecosystems, a more comprehensive carbon budget is needed, especially for gaseous carbon exchange across the water-air interface. Using two methods, overall mass balance and gas exchange and carbon burial balance, we assessed the carbon source/sink function of Lake Donghu, a subtropical, eutrophic take, from April 2003 to March 2004. With the overall mass balance calculations, total carbon input was 14 905 t, total carbon output was 4950 1, and net carbon budget was +9955 t, suggesting that Lake Donghu was a great carbon sink. For the gas exchange and carbon burial balance, gaseous carbon (CO2 and CH4) emission across the water-air interface totaled 752 t while carbon burial in the lake sediment was 9477 t. The ratio of carbon emission into the atmosphere to carbon burial into the sediment was only 0.08. This low ratio indicates that Lake Donghu is a great carbon sink. Results showed good agreement between the two methods with both showing Lake Donghu to be a great carbon sink. This results from the high primary production of Lake Donghu, substantive allochthonous carbon inputs and intensive anthropogenic activity. Gaseous carbon emission accounted for about 15% of the total carbon output, indicating that the total output would be underestimated without including gaseous carbon exchange. (C) 2007 Elsevier Ltd. All rights reserved.

Novel fluid grid and voidage calculation techniques for a discrete element model of a 3D cylindrical fluidized bed

A discrete element model (DEM) combined with computational fluid dynamics (CFD) was developed to model particle and fluid behaviour in 3D cylindrical fluidized beds. Novel techniques were developed to (1) keep fluid cells, defined in cylindrical coordinates, at a constant volume in order to ensure the conditions for validity of the volume-averaged fluid equations were satisfied and (2) smoothly and accurately measure voidage in arbitrarily shaped fluid cells. The new technique for calculating voidage was more stable than traditional techniques, also examined in the paper, whilst remaining computationally-effective. The model was validated by quantitative comparison with experimental results from the magnetic resonance imaging of a fluidised bed analysed to give time-averaged particle velocities. Comparisons were also made between theoretical determinations of slug rise velocity in a tall bed. It was concluded that the DEM-CFD model is able to investigate aspects of the underlying physics of fluidisation not readily investigated by experiment. © 2014 The Authors.

Comparison of moment-curvature models for fiber-reinforced polymer plate-end debonding studies using global energy balance approach

This paper compares a number of different moment-curvature models for cracked concrete sections that contain both steel and external fiber-reinforced polymer (FRP) reinforcement. The question of whether to use a whole-section analysis or one that considers the FRP separately is discussed. Five existing and three new models are compared with test data for moment-curvature or load deflection behavior, and five models are compared with test results for plate-end debonding using a global energy balance approach (GEBA). A proposal is made for the use of one of the simplified models. The availability of a simplified model opens the way to the production of design aids so that the GEBA can be made available to practicing engineers through design guides and parametric studies. Copyright © 2014, American Concrete Institute.

Scaling arguments for the fluxes in turbulent miscible fountains

For established axisymmetric turbulent miscible Boussinesq fountains in quiescent uniform environments, expressions are developed for the fluxes of volume, momentum and buoyancy at the outflow from the fountain: the outflow referring to the counterflow at the horizontal plane of the source. The fluxes are expressed in terms of the fountain source conditions and two dimensionless functions of the source Froude number, Fr0: a radial function (relating a horizontal scale of the outflow to the source radius) and a volume flux function (relating the outflow and source volume fluxes). The forms taken by these two functions at low Fr0 and high Fr0 are deduced, thereby providing the outflow fluxes and outflow Froude number solely in terms of the source conditions. For high Fr0, the outflow Froude number, Frout, is shown to be invariant, indicating (by analogy with plumes for which the 'far-field' Froude number is invariant with source Froude number) that the outflow may be regarded as 'far-field' since the fluxes within the fountain have adjusted to attain a balance which is independent of the source conditions. Based on Frout, the fluxes in the plume that forms beyond the fountain outflow are deduced. Finally, from the results of previously published studies, we show that the scalings deduced for fountains are valid for 0.0025 ≲ Fr0 ≲ 1.0 for low Fr0 and Fr0≳ 3.0 for high Fr0. © 2014 Cambridge University Press.

Wetness loss prediction for a low pressure steam turbine using computational fluid dynamics

Two-phase computational fluid dynamics modelling is used to investigate the magnitude of different contributions to the wet steam losses in a three-stage model low pressure steam turbine. The thermodynamic losses (due to irreversible heat transfer across a finite temperature difference) and the kinematic relaxation losses (due to the frictional drag of the drops) are evaluated directly from the computational fluid dynamics simulation using a concept based on entropy production rates. The braking losses (due to the impact of large drops on the rotor) are investigated by a separate numerical prediction. The simulations show that in the present case, the dominant effect is the thermodynamic loss that accounts for over 90% of the wetness losses and that both the thermodynamic and the kinematic relaxation losses depend on the droplet diameter. The numerical results are brought into context with the well-known Baumann correlation, and a comparison with available measurement data in the literature is given. The ability of the numerical approach to predict the main wetness losses is confirmed, which permits the use of computational fluid dynamics for further studies on wetness loss correlations. © IMechE 2013 Reprints and permissions: sagepub.co.uk/journalsPermissions.nav.

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.

simulation and interaction of fluid dynamics

In the fluid simulation, the fluids and their surroundings may greatly change properties such as shape and temperature simultaneously, and different surroundings would characterize different interactions, which would change the shape and motion of the fluids in different ways. On the other hand, interactions among fluid mixtures of different kinds would generate more comprehensive behavior. To investigate the interaction behavior in physically based simulation of fluids, it is of importance to build physically correct models to represent the varying interactions between fluids and the environments, as well as interactions among the mixtures. In this paper, we will make a simple review of the interactions, and focus on those most interesting to us, and model them with various physical solutions. In particular, more detail will be given on the simulation of miscible and immiscible binary mixtures. In some of the methods, it is advantageous to be taken with the graphics processing unit (GPU) to achieve real-time computation for middle-scale simulation.

Gas-bearing fluid influx sub-system for gas hydrate geological system in Shenhu Area, Northern South China Sea

Based on the comprehensive interpretation and study of the Neogene fracture system and diapiric structure, it can be concluded that the diapiric structures, high-angle fractures and vertical fissure system are the main gas-bearing fluid influx sub-system for gas hydrate geological system in Shenhu Area, northern South China Sea. The Neogene fractures widely developed in the study area may be classed into two groups: NW (NNW)-trending and NE (NNE)-trending. The first group was active in the Late Miocene, while the second one was active since the Pliocene. The NE (NNE)-trending fractures were characterized by lower activity strength and larger scale, and cut through the sediment layers deposited since the Pliocene. Within the top sediment layers, the high-angle fracture and vertical fissure system was developed. The diapiric structures display various types such as a turtle-back-like arch, weak piercing, gas chimney, and fracture (or crack, fissure). On the seismic profile, some diapiric structures show the vertical chimney pathway whose top is narrow and the bottom is wide, where some ones extend horizontally into pocket or flower-shaped structures and formed the seismic reflection chaotic zones. Within the overlying sediment layers of the diapiric structure, the tree branch, flower-shaped high-angle fractures and vertical fissures were developed and became the pathway and migration system of the gas-bearing fluid influx. In the study area, the diapiric structures indicate a high temperature/over pressure system ever developed. Closely associated and abundant bright-spots show the methane-bearing fluid influx migrated vertically or horizontally through the diapiric structures, high-angle fractures and vertical fissures. In the place where the temperature and pressure conditions were favor for the formation of gas hydrate, the hydrate reservoir deposition sub-system was developed.