Ovine bone and marrow derived progenitor cells : isolation, characterization, and osteogenic potential

Recently, research has focused on bone marrow derived multipotent mesenchymal precursor cells (MPC) for their potential clinical use in bone engineering. Prior to clinical application, MPC-based treatment concepts need to be evaluated in preclinical, immunocompetent, large animal models. Sheep in particular are considered a valid model for orthopaedic and trauma related research. However, ovine MPC and their osteogenic potential remain poorly characterized. In the present study, ex vivo expanded MPC isolated from ovine bone marrow proliferated at a higher rate than osteoblasts (OB) derived from tibial compact bone as assessed using standard 2D culture. MPC expressed the respective phenotypic profile typical for different mesenchymal cell populations (CD14-/CD31-/CD45- /CD29+/CD44+/CD166+) and showed a multilineage differentiation potential. When compared to OB, MPC had a higher mineralization potential under standard osteogenic culture conditions and expressed typical markers such as osteocalcin, osteonectin and type I collagen at the mRNA and protein level. After 4 weeks in 3D culture, MPC constructs demonstrated higher cell density and mineralization, whilst cell viability on the scaffolds was assessed >90%. Cells displayed a spindle-like morphology and formed an interconnected network. Implanted subcutaneously into NOD/SCID mice on type I collagen coated polycaprolactone-tricalciumphosphate (mPCL-TCP) scaffolds, MPC presented a higher developmental potential than osteoblasts. In summary, this study provides a detailed in vitro characterisation of ovine MPC from a bone engineering perspective and suggests that MPC provide promising means for future bone disease related treatment applications.

Numerical evaluation of pollutant dispersion at a toll plaza based on system dynamics and Computational Fluid Dynamics models

The health of tollbooth workers is seriously threatened by long-term exposure to polluted air from vehicle exhausts. Using traffic data collected at a toll plaza, vehicle movements were simulated by a system dynamics model with different traffic volumes and toll collection procedures. This allowed the average travel time of vehicles to be calculated. A three-dimension Computational Fluid Dynamics (CFD) model was used with a k–ε turbulence model to simulate pollutant dispersion at the toll plaza for different traffic volumes and toll collection procedures. It was shown that pollutant concentration around tollbooths increases as traffic volume increases. Whether traffic volume is low or high (1500 vehicles/h or 2500 vehicles/h), pollutant concentration decreases if electronic toll collection (ETC) is adopted. In addition, pollutant concentration around tollbooths decreases as the proportion of ETC-equipped vehicles increases. However, if the proportion of ETC-equipped vehicles is very low and the traffic volume is not heavy, then pollutant concentration increases as the number of ETC lanes increases.

Executive function effects and numerical development in children : behavioural and ERP evidence from a numerical Stroop paradigm

Most research on numerical development in children is behavioural, focusing on accuracy and response time in different problem formats. However, Temple and Posner (1998) used ERPs and the numerical distance task with 5-year-olds to show that the development of numerical representations is difficult to disentangle from the development of the executive components of response organization and execution. Here we use the numerical Stroop paradigm (NSP) and ERPs to study possible executive interference in numerical processing tasks in 6–8-year-old children. In the NSP, the numerical magnitude of the digits is task-relevant and the physical size of the digits is task-irrelevant. We show that younger children are highly susceptible to interference from irrelevant physical information such as digit size, but that access to the numerical representation is almost as fast in young children as in adults. We argue that the developmental trajectories for executive function and numerical processing may act together to determine numerical development in young children.

Atomistic numerical investigation of single-crystal copper nanowire with surface defects

Based on the embedded atom method (EAM), a molecular dynamics (MD) simulation is performed to study the single-crystal copper nanowire with surface defects through tension. The tension simulations for nanowire without defect are first carried out under different temperatures, strain rates and time steps and then surface defect effects for nanowire are investigated. The stress-strain curves obtained by the MD simulations of various strain rates show a rate below 1 x 10(9) s-1 will exert less effect on the yield strength and yield point, and the Young's modulus is independent of strain rate. a time step below 5 fs is recommend for the atomic model during the MD simulation. It is observed that high temperature leads to low Young's modulus, as well as the yield strength. The surface defects on nanowires are systematically studied in considering different defect orientations. It is found that the surface defect serves as a dislocation source, and the yield strength shows 34.20% decresse with 45 degree surface defect. Both yield strength and yield point are significantly influenced by the surface defects, except the Young's modulus.

Bioactive SrO–SiO2 glass with well-ordered mesopores : characterization, physiochemistry and biological properties

For a biomaterial to be considered suitable for bone repair it should ideally be both bioactive and have a capacity for controllable drug delivery; as such, mesoporous SiO2 glass has been proposed as a new class of bone regeneration material by virtue of its high drug-loading ability and generally good biocompatibility. It does, however, have less than optimum bioactivity and controllable drug delivery properties. In this study, we incorporated strontium (Sr) into mesoporous SiO2 in an effort to develop a bioactive mesoporous SrO–SiO2 (Sr–Si) glass with the capacity to deliver Sr2+ ions, as well as a drug, at a controlled rate, thereby producing a material better suited for bone repair. The effects of Sr2+ on the structure, physiochemistry, drug delivery and biological properties of mesoporous Sr–Si glass were investigated. The prepared mesoporous Sr–Si glass was found to have an excellent release profile of bioactive Sr2+ ions and dexamethasone, and the incorporation of Sr2+ improved structural properties, such as mesopore size, pore volume and specific surface area, as well as rate of dissolution and protein adsorption. The mesoporous Sr–Si glass had no cytotoxic effects and its release of Sr2+ and SiO44− ions enhanced alkaline phosphatase activity – a marker of osteogenic cell differentiation – in human bone mesenchymal stem cells. Mesoporous Sr–Si glasses can be prepared to porous scaffolds which show a more sustained drug release. This study suggests that incorporating Sr2+ into mesoporous SiO2 glass produces a material with a more optimal drug delivery profile coupled with improved bioactivity, making it an excellent material for bone repair applications. Keywords: Mesoporous Sr–Si glass; Drug delivery; Bioactivity; Bone repair; Scaffolds

Numerical solution to the Saffman-Taylor finger problem with kinetic undercooling regularisation

The Saffman-Taylor finger problem is to predict the shape and,in particular, width of a finger of fluid travelling in a Hele-Shaw cell filled with a different, more viscous fluid. In experiments the width is dependent on the speed of propagation of the finger, tending to half the total cell width as the speed increases. To predict this result mathematically, nonlinear effects on the fluid interface must be considered; usually surface tension is included for this purpose. This makes the mathematical problem suffciently diffcult that asymptotic or numerical methods must be used. In this paper we adapt numerical methods used to solve the Saffman-Taylor finger problem with surface tension to instead include the effect of kinetic undercooling, a regularisation effect important in Stefan melting-freezing problems, for which Hele-Shaw flow serves as a leading order approximation when the specific heat of a substance is much smaller than its latent heat. We find the existence of a solution branch where the finger width tends to zero as the propagation speed increases, disagreeing with some aspects of the asymptotic analysis of the same problem. We also find a second solution branch, supporting the idea of a countably infinite number of branches as with the surface tension problem.

Shrinking neighborhood evolution : a novel stochastic algorithm for numerical optimization

Focusing on the conditions that an optimization problem may comply with, the so-called convergence conditions have been proposed and sequentially a stochastic optimization algorithm named as DSZ algorithm is presented in order to deal with both unconstrained and constrained optimizations. The principle is discussed in the theoretical model of DSZ algorithm, from which we present the practical model of DSZ algorithm. Practical model efficiency is demonstrated by the comparison with the similar algorithms such as Enhanced simulated annealing (ESA), Monte Carlo simulated annealing (MCS), Sniffer Global Optimization (SGO), Directed Tabu Search (DTS), and Genetic Algorithm (GA), using a set of well-known unconstrained and constrained optimization test cases. Meanwhile, further attention goes to the strategies how to optimize the high-dimensional unconstrained problem using DSZ algorithm.

Numerical simulation of flow in a photocatalytic reactor containing baffles

In this study a new immobilized flat plate photocatalytic reactor for wastewater treatment has been investigated using computational fluid dynamics (CFD). The reactor consists of a reactor inlet, a reactive section where the catalyst is coated, and outlet parts. For simulation, the reactive section of the reactor was modelled with an array of baffles. In order to optimize the fluid mixing and reactor design, this study attempts to investigate the influence of baffles with differing heights on the flow field of the flat plate reactor. The results obtained from the simulation of a baffled flat plate reactor hydrodynamics for differing baffle heights for certain positions are presented. Under the conditions simulated, the qualitative flow features, such as the distribution of local stream lines, velocity contours, and high shear region, boundary layers separation, vortex formation, and the underlying mechanism are examined. At low and high Re numbers, the influence of baffle heights on the distribution of species mass fraction of a model pollutant are also highlighted. The simulation of qualitative and quantitative properties of fluid dynamics in a baffled reactor provides valuable insight to fully understand the effect of baffles and their role on the flow pattern, behaviour, and features of wastewater treatment using a photocatalytic reactor.

The effect of friction on indenter force and pile-up in numerical simulations of bone nanoindentation

Nanoindentation is a useful technique for probing the mechanical properties of bone, and finite element (FE) modeling of the indentation allows inverse determination of elasto-plastic constitutive properties. However, all but one FE study to date have assumed frictionless contact between indenter and bone. The aim of this study was to explore the effect of friction in simulations of bone nanoindentation. Two dimensional axisymmetric FE simulations were performed using a spheroconical indenter of tip radius 0.6 m and angle 90°. The coefficient of friction between indenter and bone was varied between 0.0 (frictionless) and 0.3. Isotropic linear elasticity was used in all simulations, with bone elastic modulus E=13.56GPa and Poisson‟s ratio f 0.3. Plasticity was incorporated using both Drucker-Prager and von Mises yield surfaces. Friction had a modest effect on the predicted force-indentation curve for both von Mises and Drucker-Prager plasticity, reducing maximum indenter displacement by 10% and 20% respectively as friction coefficient was increased from zero to 0.3 (at a maximum indenter force of 5mN). However, friction has a much greater effect on predicted pile-up after indentation, reducing predicted pile-up from 0.27 to 0.11 m with a von Mises model, and from 0.09 to 0.02 m with Drucker-Prager plasticity. We conclude that it is potentially important to include friction in nanoindentation simulations of bone if pile-up is used to compare simulation results with experiment.

Numerical study of two ill-posed one-phase Stefan problems

We treat two related moving boundary problems. The first is the ill-posed Stefan problem for melting a superheated solid in one Cartesian coordinate. Mathematically, this is the same problem as that for freezing a supercooled liquid, with applications to crystal growth. By applying a front-fixing technique with finite differences, we reproduce existing numerical results in the literature, concentrating on solutions that break down in finite time. This sort of finite-time blow-up is characterised by the speed of the moving boundary becoming unbounded in the blow-up limit. The second problem, which is an extension of the first, is proposed to simulate aspects of a particular two-phase Stefan problem with surface tension. We study this novel moving boundary problem numerically, and provide results that support the hypothesis that it exhibits a similar type of finite-time blow-up as the more complicated two-phase problem. The results are unusual in the sense that it appears the addition of surface tension transforms a well-posed problem into an ill-posed one.

Numerical simulation for the variable-order Galilei invariant advection diffusion equation with a nonlinear source term

In this paper, we consider the variable-order Galilei advection diffusion equation with a nonlinear source term. A numerical scheme with first order temporal accuracy and second order spatial accuracy is developed to simulate the equation. The stability and convergence of the numerical scheme are analyzed. Besides, another numerical scheme for improving temporal accuracy is also developed. Finally, some numerical examples are given and the results demonstrate the effectiveness of theoretical analysis. Keywords: The variable-order Galilei invariant advection diffusion equation with a nonlinear source term; The variable-order Riemann–Liouville fractional partial derivative; Stability; Convergence; Numerical scheme improving temporal accuracy