Correlating flow induced shear stress and chondrocytes activity in micro-porous scaffold using computational fluid dynamics and rapid prototyping

Flow induced shear stress plays an important role in regulating cell growth and distribution in scaffolds. This study sought to correlate wall shear stress and chondrocytes activity for engineering design of micro-porous osteochondral grafts based on the hypothesis that it is possible to capture and discriminate between the transmitted force and cell response at the inner irregularities. Unlike common tissue engineering therapies with perfusion bioreactors in which flow-mediated stress is the controlling parameter, this work assigned the associated stress as a function of porosity to influence in vitro proliferation of chondrocytes. D-optimality criterion was used to accommodate three pore characteristics for appraisal in a mixed level fractional design of experiment (DOE); namely, pore size (4 levels), distribution pattern (2 levels) and density (3 levels). Micro-porous scaffolds (n=12) were fabricated according to the DOE using rapid prototyping of an acrylic-based bio-photopolymer. Computational fluid dynamics (CFD) models were created correspondingly and used on an idealized boundary condition with a Newtonian fluid domain to simulate the dynamic microenvironment inside the pores. In vitro condition was reproduced for the 3D printed constructs seeded by high pellet densities of human chondrocytes and cultured for 72 hours. The results showed that cell proliferation was significantly different in the constructs (p<0.05). Inlet fluid velocity of 3×10-2mms-1 and average shear stress of 5.65×10-2 Pa corresponded with increased cell proliferation for scaffolds with smaller pores in hexagonal pattern and lower densities. Although the analytical solution of a Poiseuille flow inside the pores was found insufficient for the description of the flow profile probably due to the outside flow induced turbulence, it showed that the shear stress would increase with cell growth and decrease with pore size. This correlation demonstrated the basis for determining the relation between the induced stress and chondrocyte activity to optimize microfabrication of engineered cartilaginous constructs.

Effects of scaffold architecture on mechanical characteristics and osteoblast response to static and perfusion bioreactor cultures

Tissue engineering focuses on the repair and regeneration of tissues through the use of biodegradable scaffold systems that structurally support regions of injury whilst recruiting and/or stimulating cell populations to rebuild the target tissue. Within bone tissue engineering, the effects of scaffold architecture on cellular response have not been conclusively characterized in a controlled-density environment. We present a theoretical and practical assessment of the effects of polycaprolactone (PCL) scaffold architectural modifications on mechanical and flow characteristics as well as MC3T3-E1 preosteoblast cellular response in an in vitro static plate and custom-designed perfusion bioreactor model. Four scaffold architectures were contrasted, which varied in inter-layer lay-down angle and offset between layers, whilst maintaining a structural porosity of 60 ± 5%. We established that as layer angle was decreased (90° vs. 60°) and offset was introduced (0 vs. 0.5 between layers), structural stiffness, yield stress, strength, pore size and permeability decreased, whilst computational fluid dynamics-modeled wall shear stress was increased. Most significant effects were noted with layer offset. Seeding efficiencies in static culture were also dramatically increased due to offset (~45% to ~86%), with static culture exhibiting a much higher seeding efficiency than perfusion culture. Scaffold architecture had minimal effect on cell response in static culture. However, architecture influenced osteogenic differentiation in perfusion culture, likely by modifying the microfluidic environment.

Strong and bioactive tri-calcium phosphate scaffolds with tube-like macropores

Calcium phosphate ceramic scaffolds have been widely investigated for bone tissue engineering due to their excellent biocompatibility and biodegradation. Unfortunately, they have the shortcoming of low mechanical properties. In order to provide strong, bioactive, and biodegradable scaffolds, a new approach of infiltrating the macro-tube ABS (acrylontrile butadiene styrene) templates with a hydroxyapatite/bioactive glass mixed slurry was developed to fabricate porous Si-doped TCP (tri-calcium phosphate) scaffolds. The porous Si-doped TCP ceramics with a high porosity (~65%) and with interconnected macrotubes (~0.8mm in diameter) and micropores (5-100 m) had a high compressive strength (up to 14.68+0.2MPa), which was comparable to that of a trabecular bone and was much higher than those of pure TCP scaffolds. Additional cell attachment study and MTT cytotoxicity assay proved the bioactivity and biocompatibility of the new scaffolds. Thus a potential bioceramic material and a new approach to make the potential scaffolds were developed for bone tissue engineering.

Quantitative estimation of bioturbation based on digital image analysis

Quantitative determination of modification of primary sediment features, by the activity of organisms (i.e., bioturbation) is essential in geosciences. Some methods proposed since the 1960s are mainly based on visual or subjective determinations. The first semiquantitative evaluations of the Bioturbation Index, Ichnofabric Index, or the amount of bioturbation were attempted, in the best cases using a series of flashcards designed in different situations. Recently, more effective methods involve the use of analytical and computational methods such as X-rays, magnetic resonance imaging or computed tomography; these methods are complex and often expensive. This paper presents a compilation of different methods, using Adobe® Photoshop® software CS6, for digital estimation that are a part of the IDIAP (Ichnological Digital Analysis Images Package), which is an inexpensive alternative to recently proposed methods, easy to use, and especially recommended for core samples. The different methods — “Similar Pixel Selection Method (SPSM)”, “Magic Wand Method (MWM)” and the “Color Range Selection Method (CRSM)” — entail advantages and disadvantages depending on the sediment (e.g., composition, color, texture, porosity, etc.) and ichnological features (size of traces, infilling material, burrow wall, etc.). The IDIAP provides an estimation of the amount of trace fossils produced by a particular ichnotaxon, by a whole ichnocoenosis or even for a complete ichnofabric. We recommend the application of the complete IDIAP to a given case study, followed by selection of the most appropriate method. The IDIAP was applied to core material recovered from the IODP Expedition 339, enabling us, for the first time, to arrive at a quantitative estimation of the discrete trace fossil assemblage in core samples.

Gas sensing of ruthenium implanted tungsten oxide thin films

Different amounts of Ru were implanted into thermally evaporated WO3 thin films by ion implantation. The films were subsequently annealed at 600oC for 2 hours in air to remove defects generated during the ion implantation. The Ru concentrations of four samples have been quantified by Rutherford Backscattering Spectrometry as 0.8, 5.5, 9 and 11.5 at%. The un-implanted WO3 films were highly porous but the porosity decreased significantly after ion implantation as observed by Transmission Electron Microscopy and Scanning Electron Microscopy. The thickness of the films also decreased with increasing Ru-ion dose, which is mainly due to densification of the porous films during ion implantation. From Raman spectroscopy two peaks at 408 and 451 cm-1 (in addition to the typical vibrational peaks of the monoclinic WO3 phase) associated with Ru were observed. Their intensity increased with increasing Ru concentration. X-Ray Photoelectron Spectroscopy showed a metallic state of Ru with binding energy of Ru 3d5/2 at 280.1 eV. This peak position remained almost unchanged with increasing Ru concentration. The resistances of the Ru-implanted films were found to increase in the presence of NO2 and NO with higher sensor response to NO2. The effect of Ru concentration on the sensing performance of the films was not explicitly observed due to reduced film thickness and porosity with increasing Ru concentration. However, the results indicate that the implantation of Ru into WO3 films with sufficient film porosity and film thickness can be beneficial for NO2 sensing at temperatures in the range of 250°C to 350°C.

Pore size distribution and accessible pore size distribution in bituminous coals

The porosity and pore size distribution of coals determine many of their properties, from gas release to their behavior on carbonization, and yet most methods of determining pore size distribution can only examine a restricted size range. Even then, only accessible pores can be investigated with these methods. Small-angle neutron scattering (SANS) and ultra small-angle neutron scattering (USANS) are increasingly used to characterize the size distribution of all of the pores non-destructively. Here we have used USANS/SANS to examine 24 well-characterized bituminous and subbituminous coals: three from the eastern US, two from Poland, one from New Zealand and the rest from the Sydney and Bowen Basins in Eastern Australia, and determined the relationships of the scattering intensity corresponding to different pore sizes with other coal properties. The range of pore radii examinable with these techniques is 2.5nm to 7μm. We confirm that there is a wide range of pore sizes in coal. The pore size distribution was found to be strongly affected by both rank and type (expressed as either hydrogen or vitrinite content) in the size range 250nm to 7μm and 5 to 10nm, but weakly in intermediate regions. The results suggest that different mechanisms control coal porosity on different scales. Contrast-matching USANS and SANS were also used to determine the size distribution of the fraction of the pores in these coals that are inaccessible to deuterated methane, CD4, at ambient temperature. In some coals most of the small (~10nm) pores were found to be inaccessible to CD4 on the time scale of the measurement (~30min–16h). This inaccessibility suggests that in these coals a considerable fraction of inherent methane may be trapped for extended periods of time, thus reducing the effectiveness of methane release from (or sorption by) these coals. Although the number of small pores was less in higher rank coals, the fraction of total pores that was inaccessible was not rank dependent. In the Australian coals, at the 10nm to 50nm size scales the pores in inertinites appeared to be completely accessible to CD4, whereas the pores in the vitrinite were about 75% inaccessible. Unlike the results for total porosity that showed no regional effects on relationships between porosity and coal properties, clear regional differences in the relationships between fraction of closed porosity and coal properties were found. The 10 to 50nm-sized pores of inertinites of the US and Polish coals examined appeared less accessible to methane than those of the inertinites of Australian coals. This difference in pore accessibility in inertinites may explain why empirical relationships between fluidity and coking properties developed using Carboniferous coals do not apply to Australian coals.

Pore accessibility by methane and carbon dioxide in coal as determined by neutron scattering

Contrast-matching ultrasmall-angle neutron scattering (USANS) and small-angle neutron scattering (SANS) techniques were used for the first time to determine both the total pore volume and the fraction of the pore volume that is inaccessible to deuterated methane, CD4, in four bituminous coals in the range of pore sizes between ∼10 Å and ∼5 μm. Two samples originated from the Illinois Basin in the U.S.A., and the other two samples were commercial Australian bituminous coals from the Bowen Basin. The total and inaccessible porosity were determined in each coal using both Porod invariant and the polydisperse spherical particle (PDSP) model analysis of the scattering data acquired from coals both in vacuum and at the pressure of CD4, at which the scattering length density of the pore-saturating fluid is equal to that of the solid coal matrix (zero average contrast pressure). The total porosity of the coals studied ranged from 7 to 13%, and the volume of pores inaccessible to CD4 varied from ∼13 to ∼36% of the total pore volume. The volume fraction of inaccessible pores shows no correlation with the maceral composition; however, it increases with a decreasing total pore volume. In situ measurements of the structure of one coal saturated with CO2 and CD4 were conducted as a function of the pressure in the range of 1−400 bar. The neutron scattering intensity from small pores with radii less than 35 Å in this coal increased sharply immediately after the fluid injection for both gases, which demonstrates strong condensation and densification of the invading subcritical CO2 and supercritical methane in small pores.

Accessibility of pores in coal to methane and carbon dioxide

Fluid–solid interactions in natural and engineered porous solids underlie a variety of technological processes, including geological storage of anthropogenic greenhouse gases, enhanced coal bed methane recovery, membrane separation, and heterogeneous catalysis. The size, distribution and interconnectivity of pores, the chemical and physical properties of the solid and fluid phases collectively dictate how fluid molecules migrate into and through the micro- and meso-porous media, adsorb and ultimately react with the solid surfaces. Due to the high penetration power and relatively short wavelength of neutrons, smallangle neutron scattering (SANS) as well as ultra small-angle scattering (USANS) techniques are ideally suited for assessing the phase behavior of confined fluids under pressure as well as for evaluating the total porosity in engineered and natural porous systems including coal. Here we demonstrate that SANS and USANS can be also used for determining the fraction of the pore volume that is actually accessible to fluids as a function of pore sizes and study the fraction of inaccessible pores as a function of pore size in three coals from the Illinois Basin (USA) and Bowen Basin (Australia). Experiments were performed at CO2 and methane pressures up to 780 bar, including pressures corresponding to zero average contrast condition (ZAC), which is the pressure where no scattering from the accessible pores occurs. Scattering curves at the ZAC were compared with the scattering from same coals under vacuum and analysed using a newly developed approach that shows that the volume fraction of accessible pores in these coals varies between �90% in the macropore region to �30% in the mesopore region and the variation is distinctive for each of the examined coals. The developed methodology may be also applied for assessing the volume of accessible pores in other natural underground formations of interest for CO2 sequestration, such as saline aquifers as well as for estimating closed porosity in engineered porous solids of technological importance.

Two-scale computational modelling of unsaturated flow in soils exhibiting small scale heterogeneities

Unsaturated water flow in soil is commonly modelled using Richards’ equation, which requires the hydraulic properties of the soil (e.g., porosity, hydraulic conductivity, etc.) to be characterised. Naturally occurring soils, however, are heterogeneous in nature, that is, they are composed of a number of interwoven homogeneous soils each with their own set of hydraulic properties. When the length scale of these soil heterogeneities is small, numerical solution of Richards’ equation is computationally impractical due to the immense effort and refinement required to mesh the actual heterogeneous geometry. A classic way forward is to use a macroscopic model, where the heterogeneous medium is replaced with a fictitious homogeneous medium, which attempts to give the average flow behaviour at the macroscopic scale (i.e., at a scale much larger than the scale of the heterogeneities). Using the homogenisation theory, a macroscopic equation can be derived that takes the form of Richards’ equation with effective parameters. A disadvantage of the macroscopic approach, however, is that it fails in cases when the assumption of local equilibrium does not hold. This limitation has seen the introduction of two-scale models that include at each point in the macroscopic domain an additional flow equation at the scale of the heterogeneities (microscopic scale). This report outlines a well-known two-scale model and contributes to the literature a number of important advances in its numerical implementation. These include the use of an unstructured control volume finite element method and image-based meshing techniques, that allow for irregular micro-scale geometries to be treated, and the use of an exponential time integration scheme that permits both scales to be resolved simultaneously in a completely coupled manner. Numerical comparisons against a classical macroscopic model confirm that only the two-scale model correctly captures the important features of the flow for a range of parameter values.

A novel route in bone tissue engineering : magnetic biomimetic scaffolds

In recent years, interest in tissue engineering and its solutions has increased considerably. In particular, scaffolds have become fundamental tools in bone graft substitution and are used in combination with a variety of bio-agents. However, a long-standing problem in the use of these conventional scaffolds lies in the impossibility of re-loading the scaffold with the bio-agents after implantation. This work introduces the magnetic scaffold as a conceptually new solution. The magnetic scaffold is able, via magnetic driving, to attract and take up in vivo growth factors, stem cells or other bio-agents bound to magnetic particles. The authors succeeded in developing a simple and inexpensive technique able to transform standard commercial scaffolds made of hydroxyapatite and collagen in magnetic scaffolds. This innovative process involves dip-coating of the scaffolds in aqueous ferrofluids containing iron oxide nanoparticles coated with various biopolymers. After dip-coating, the nanoparticles are integrated into the structure of the scaffolds, providing the latter with magnetization values as high as 15 emu g�1 at 10 kOe. These values are suitable for generating magnetic gradients, enabling magnetic guiding in the vicinity and inside the scaffold. The magnetic scaffolds do not suffer from any structural damage during the process, maintaining their specific porosity and shape. Moreover, they do not release magnetic particles under a constant flow of simulated body fluids over a period of 8 days. Finally, preliminary studies indicate the ability of the magnetic scaffolds to support adhesion and proliferation of human bone marrow stem cells in vitro. Hence, this new type of scaffold is a valuable candidate for tissue engineering applications, featuring a novel magnetic guiding option.

Optimization approaches for the design of additively manufactured scaffolds

Scaffolds play a pivotal role in tissue engineering, promoting the synthesis of neo extra-cellular matrix (ECM), and providing temporary mechanical support for the cells during tissue regeneration. Advances introduced by additive manufacturing techniques have significantly improved the ability to regulate scaffold architecture, enhancing the control over scaffold shape and porosity. Thus, considerable research efforts have been devoted to the fabrication of 3D porous scaffolds with optimized micro-architectural features. This chapter gives an overview of the methods for the design of additively manufactured scaffolds and their applicability in tissue engineering (TE). Along with a survey of the state of the art, the Authors will also present a recently developed method, called Load-Adaptive Scaffold Architecturing (LASA), which returns scaffold architectures optimized for given applied mechanical loads systems, once the specific stress distribution is evaluated through Finite Element Analysis (FEA).

A meshfree particle based model for microscale shrinkage mechanisms of food materials in drying conditions

Cells are the fundamental building block of plant based food materials and many of the food processing born structural changes can fundamentally be derived as a function of the deformations of the cellular structure. In food dehydration the bulk level changes in porosity, density and shrinkage can be better explained using cellular level deformations initiated by the moisture removal from the cellular fluid. A novel approach is used in this research to model the cell fluid with Smoothed Particle Hydrodynamics (SPH) and cell walls with Discrete Element Methods (DEM), that are fundamentally known to be robust in treating complex fluid and solid mechanics. High Performance Computing (HPC) is used for the computations due to its computing advantages. Comparing with the deficiencies of the state of the art drying models, the current model is found to be robust in replicating drying mechanics of plant based food materials in microscale.

Effect of ion current density on the properties of vacuum arc-deposited TiN coatings

The influence of ion current density on the thickness of coatings deposited in a vacuum arc setup has been investigated to optimize the coating porosity. A planar probe was used to measure the ion current density distribution across plasma flux. A current density from 20 to 50 A/m2 was obtained, depending on the probe position relative to the substrate center. TiN coatings were deposited onto the cutting inserts placed at different locations on the substrate, and SEM was used to characterize the surfaces of the coatings. It was found that lowdensity coatings were formed at the decreased ion current density. A quantitative dependence of the coating thickness on the ion current density in the range of 20-50 A/m2 were obtained for the films deposited at substrate bias of 200 V and nitrogen pressure 0.1 Pa, and the coating porosity was calculated. The coated cutting inserts were tested by lathe machining of the martensitic stainless steel AISI 431. The results may be useful for controlling ion flux distribution over large industrial-scale substrates.

Durability performances of concrete with nano-silica

This study investigated the durability properties of concrete containing nano-silica at dosages of 0.3% and 0.9%, respectively. Due to the nano-filler effect and the pozzolanic reaction, the microstructure became more homogeneous and less porous, especially at the interfacial transition zone (ITZ), which led to reduced permeability. Tests on the durability properties verified the beneficial effects of nano-silica. The channels for harmful agents through the cement composites were partially filled and blocked. The pore size distribution also indicated that the large capillary pores were refined by the nano-silica, due to the combined contribution of the nano-filler effect and the pozzolanic reaction.

Theoretical modelling of momentum transfer function of bi-disperse porous media

The aim of this paper is to obtain the momentum transfer coefficient between the two phases, denoted by f and p, occupying a bi-disperse porous medium by mapping the available experimental data to the theoretical model proposed by Nield and Kuznetsov. Data pertinent to plate-fin heat exchangers, as bi-disperse porous media, were used. The measured pressure drops for such heat exchangers are then used to give the overall permeability which is linked to the porosity and permeability of each phase as well as the interfacial momentum transfer coefficient between the two phases. Accordingly, numerical values are obtained for the momentum transfer coefficient for three different fin spacing values considered in the heat exchanger experiments.