Mechanobiology of the meniscus.

The meniscus plays a critical biomechanical role in the knee, providing load support, joint stability, and congruity. Importantly, growing evidence indicates that the mechanobiologic response of meniscal cells plays a critical role in the physiologic, pathologic, and repair responses of the meniscus. Here we review experimental and theoretical studies that have begun to directly measure the biomechanical effects of joint loading on the meniscus under physiologic and pathologic conditions, showing that the menisci are exposed to high contact stresses, resulting in a complex and nonuniform stress-strain environment within the tissue. By combining microscale measurements of the mechanical properties of meniscal cells and their pericellular and extracellular matrix regions, theoretical and experimental models indicate that the cells in the meniscus are exposed to a complex and inhomogeneous environment of stress, strain, fluid pressure, fluid flow, and a variety of physicochemical factors. Studies across a range of culture systems from isolated cells to tissues have revealed that the biological response of meniscal cells is directly influenced by physical factors, such as tension, compression, and hydrostatic pressure. In addition, these studies have provided new insights into the mechanotransduction mechanisms by which physical signals are converted into metabolic or pro/anti-inflammatory responses. Taken together, these in vivo and in vitro studies show that mechanical factors play an important role in the health, degeneration, and regeneration of the meniscus. A more thorough understanding of the mechanobiologic responses of the meniscus will hopefully lead to therapeutic approaches to prevent degeneration and enhance repair of the meniscus.

In vitro fluid dynamics of the Ahmed glaucoma valve modified with expanded polytetrafluoroethylene.

PURPOSE: Long-term intraocular pressure reduction by glaucoma drainage devices (GDDs) is often limited by the fibrotic capsule that forms around them. Prior work demonstrates that modifying a GDD with a porous membrane promotes a vascularized and more permeable capsule. This work examines the in vitro fluid dynamics of the Ahmed valve after enclosing the outflow tract with a porous membrane of expanded polytetrafluoroethylene (ePTFE). MATERIALS AND METHODS: The control and modified Ahmed implants (termed porous retrofitted implant with modified enclosure or PRIME-Ahmed) were submerged in saline and gelatin and perfused in a system that monitored flow (Q) and pressure (P). Flow rates of 1-50 μl/min were applied and steady state pressure recorded. Resistance was calculated by dividing pressure by flow. RESULTS: Modifying the Ahmed valve implant outflow with expanded ePTFE increased pressure and resistance. Pressure at a flow of 2 μl/min was increased in the PRIME-Ahmed (11.6 ± 1.5 mm Hg) relative to the control implant (6.5 ± 1.2 mm Hg). Resistance at a flow of 2 μl/min was increased in the PRIME-Ahmed (5.8 ± 0.8 mm Hg/μl/min) when compared to the control implant (3.2 ± 0.6 mm Hg/μl/min). CONCLUSIONS: Modifying the outflow tract of the Ahmed valve with a porous membrane adds resistance that decreases with increasing flow. The Ahmed valve implant behaves as a variable resistor. It is partially open at low pressures and provides reduced resistance at physiologic flow rates.

EAE cerebrospinal fluid augments in vitro phagocytosis and metabolism of CNS myelin by macrophages.

Previous studies from this laboratory have shown that CNS myelin is phagocytized and metabolized by cultured rat macrophages to a much larger extent when myelin is pretreated with serum containing antibodies to myelin constituents than when it is left untreated or pretreated with non-specific serum. In this study the effect of cerebrospinal fluid (CSF) from rabbits with experimental allergic encephalomyelitis (EAE) in promoting myelin phagocytosis was examined. Fourteen rabbits were immunized with purified myelin in Freund's complete adjuvant, seven of which developed clinical EAE symptoms. Serum and CSF were collected from EAE and control rabbits, and the CSF was centrifuged to remove cells. Sera and CSF from these rabbits and from Freund's adjuvant-immunized controls and untreated controls were measured for IgG content by radial diffusion assay, their myelin antibody characteristics were analyzed by immunoblots, and the ability of these serum and CSF samples to promote myelin phagocytosis when used for myelin opsonization was examined. The ability of a CSF sample to enhance radioactive myelin uptake and phagocytosis by cultured macrophages as measured by the appearance of radioactive cholesterol ester was linearly proportional to its total IgG titer, and correlated approximately both with clinical symptoms of the animal and the presence of antibody against the myelin constituents myelin basic protein, proteolipid protein, and galactocerebroside. The cholesterol esterification activities of EAE sera correlated to a lesser extent with IgG levels and clinical symptoms.(ABSTRACT TRUNCATED AT 250 WORDS)

Optical-CT 3D Dosimetry Using Fresnel Lenses with Minimal Refractive-Index Matching Fluid.

Telecentric optical computed tomography (optical-CT) is a state-of-the-art method for visualizing and quantifying 3-dimensional dose distributions in radiochromic dosimeters. In this work a prototype telecentric system (DFOS-Duke Fresnel Optical-CT Scanner) is evaluated which incorporates two substantial design changes: the use of Fresnel lenses (reducing lens costs from $10-30K t0 $1-3K) and the use of a 'solid tank' (which reduces noise, and the volume of refractively matched fluid from 1 ltr to 10 cc). The efficacy of DFOS was evaluated by direct comparison against commissioned scanners in our lab. Measured dose distributions from all systems were compared against the predicted dose distributions from a commissioned treatment planning system (TPS). Three treatment plans were investigated including a simple four-field box treatment, a multiple small field delivery, and a complex IMRT treatment. Dosimeters were imaged within 2 h post irradiation, using consistent scanning techniques (360 projections acquired at 1 degree intervals, reconstruction at 2mm). DFOS efficacy was evaluated through inspection of dose line-profiles, and 2D and 3D dose and gamma maps. DFOS/TPS gamma pass rates with 3%/3mm dose difference/distance-to-agreement criteria ranged from 89.3% to 92.2%, compared to from 95.6% to 99.0% obtained with the commissioned system. The 3D gamma pass rate between the commissioned system and DFOS was 98.2%. The typical noise rates in DFOS reconstructions were up to 3%, compared to under 2% for the commissioned system. In conclusion, while the introduction of a solid tank proved advantageous with regards to cost and convenience, further work is required to improve the image quality and dose reconstruction accuracy of the new DFOS optical-CT system.

Modelling meso-scale diffusion processes in stochastic fluid bio-membranes

The space–time dynamics of rigid inhomogeneities (inclusions) free to move in a randomly fluctuating fluid bio-membrane is derived and numerically simulated as a function of the membrane shape changes. Both vertically placed (embedded) inclusions and horizontally placed (surface) inclusions are considered. The energetics of the membrane, as a two-dimensional (2D) meso-scale continuum sheet, is described by the Canham–Helfrich Hamiltonian, with the membrane height function treated as a stochastic process. The diffusion parameter of this process acts as the link coupling the membrane shape fluctuations to the kinematics of the inclusions. The latter is described via Ito stochastic differential equation. In addition to stochastic forces, the inclusions also experience membrane-induced deterministic forces. Our aim is to simulate the diffusion-driven aggregation of inclusions and show how the external inclusions arrive at the sites of the embedded inclusions. The model has potential use in such emerging fields as designing a targeted drug delivery system.

A non-Newtonian computational fluid dynamics study of the stencil printing process

This paper describes the application of computational fluid dynamics (CFD) to simulate the macroscopic bulk motion of solder paste ahead of a moving squeegee blade in the stencil printing process during the manufacture of electronic components. The successful outcome of the stencil printing process is dependent on the interaction of numerous process parameters. A better understanding of these parameters is required to determine their relation to print quality and improve guidelines for process optimization. Various modelling techniques have arisen to analyse the flow behaviour of solder paste, including macroscopic studies of the whole mass of paste as well as microstructural analyses of the motion of individual solder particles suspended in the carrier fluid. This work builds on the knowledge gained to date from earlier analytical models and CFD investigations by considering the important non-Newtonian rheological properties of solder pastes which have been neglected in previous macroscopic studies. Pressure and velocity distributions are obtained from both Newtonian and non-Newtonian CFD simulations and evaluated against each other as well as existing established analytical models. Significant differences between the results are observed, which demonstrate the importance of modelling non-Newtonian properties for realistic representation of the flow behaviour of solder paste.

SMARTFIRE: an integrated computational fluid dynamics code and expert system for fire field modelling

SMARTFIRE, an open architecture integrated CFD code and knowledge based system attempts to make fire field modeling accessible to non-experts in Computational Fluid Dynamics (CFD) such as fire fighters, architects and fire safety engineers. This is achieved by embedding expert knowledge into CFD software. This enables the 'black-art' associated with the CFD analysis such as selection of solvers, relaxation parameters, convergence criteria, time steps, grid and boundary condition specification to be guided by expert advice from the software. The user is however given the option of overriding these decisions, thus retaining ultimate control. SMARTFIRE also makes use of recent developments in CFD technology such as unstructured meshes and group solvers in order to make the CFD analysis more efficient. This paper describes the incorporation within SMARTFIRE of the expert fire modeling knowledge required for automatic problem setup and mesh generation as well as the concept and use of group solvers for automatic and manual dynamic control of the CFD code.

Magnetic levitation fluid dynamics

This work is concerned with the accurate computation of flow in a rapidly deforming liquid metal droplet, suspended in an AC magnetic field. Intense flow motion due to the induced electromagnetic force distorts dynamically the droplet envelope, which is initially spherical. The relative positional change between the liquid metal surface and the surrounding coil means that fluid flow and magnetic field computations need to be closely coupled. A spectral technique is used to solve this problem, which is assumed axisymmetric. The computed results are compared against a physical experiment and "ideal sphere" analytic solutions. A comparison between the "magnetic pressure" approximation and the full electromagnetic force solutions, shows fundamental differences; the full electromagnetic force solution is necessary for accurate results in most practical applications of this technique. The physical reason for the fundamental discrepancy is the difference in the electromagnetic force representation: only the gradient part of the full force is accounted for in the "magnetic pressure" approximation. Figs 9, Refs 13.

Utilising computational fluid dynamics (CFD) for the modelling of granular material in large-scale engineering processes

In this paper, the framework is described for the modelling of granular material by employing Computational Fluid Dynamics (CFD). This is achieved through the use and implementation in the continuum theory of constitutive relations, which are derived in a granular dynamics framework and parametrise particle interactions that occur at the micro-scale level. The simulation of a process often met in bulk solids handling industrial plants involving granular matter, (i.e. filling of a flat-bottomed bin with a binary material mixture through pneumatic conveying-emptying of the bin in core flow mode-pneumatic conveying of the material coming out of a the bin) is presented. The results of the presented simulation demonstrate the capability of the numerical model to represent successfully key granular processes (i.e. segregation/degradation), the prediction of which is of great importance in the process engineering industry.

Modelling the fluid dynamics and coupled phenomena in arc weld pools

A 3D model of melt pool created by a moving arc type heat sources has been developed. The model solves the equations of turbulent fluid flow, heat transfer and electromagnetic field to demonstrate the flow behaviour phase-change in the pool. The coupled effects of buoyancy, capillary (Marangoni) and electromagnetic (Lorentz) forces are included within an unstructured finite volume mesh environment. The movement of the welding arc along the workpiece is accomplished via a moving co-ordinator system. Additionally a method enabling movement of the weld pool surface by fluid convection is presented whereby the mesh in the liquid region is allowed to move through a free surface. The surface grid lines move to restore equilibrium at the end of each computational time step and interior grid points then adjust following the solution of a Laplace equation.

A volume of fluid numerical model for wave impacts at coastal structures

The first stages in the development of a new design tool, to be used by coastal engineers to improve the efficiency, analysis, design, management and operation of a wide range of coastal and harbour structures, are described. The tool is based on a two-dimensional numerical model, NEWMOTICS-2D, using the volume of fluid (VOF) method, which permits the rapid calculation of wave hydrodynamics at impermeable natural and man-made structures. The critical hydrodynamic flow processes and forces are identified together with the equations that describe these key processes. The different possible numerical approaches for the solution of these equations, and the types of numerical models currently available, are examined and assessed. Preliminary tests of the model, using comparisons with results from a series of hydraulic model test cases, are described. The results of these tests demonstrate that the VOF approach is particularly appropriate for the simulation of the dynamics of waves at coastal structures because of its flexibility in representing the complex free surfaces encountered during wave impact and breaking. The further programme of work, required to develop the existing model into a tool for use in routine engineering design, is outlined.

Dynamics of magnetically suspended fluid

Magnetic suspension is a technique for processing pure or reactive materials without contact to walls. This work is concerned with the flow in the rapidly deforming liquid volume, suspended in an AC magnetic field. Intense flow motion due to the induced electromagnetic force distorts dynamically the droplet envelope. The relative positional change between the liquid surface and the surrounding coil means that fluid flow and magnetic field computations need to be closely coupled. The computed results are compared against a physical experiment and nearly spherical analytic solutions. A comparison between the "magetic pressure" approximation and the full electromagnetic force solutions shows fundamental differences; the full electromagnetic force is necessary for accurate results in most practical applications of this technique. The physical reason for the fundamental discrepancy is the difference in the electromagnetic force representation: only the gradient part of the full force is accounted for in the "magnetic pressure" approximation.

Details of an integrated approach to three-dimensional dynamic fluid structure interaction

Fluid structure interaction, as applied to flexible structures, has wide application in diverse areas such as flutter in aircraft, wind response of buildings, flows in elastic pipes and blood vessels. Numerical modelling of dynamic fluid-structure interaction (DFSI) involves the coupling of fluid flow and structural mechanics, two fields that are conventionally modelled using two dissimilar methods, thus a single comprehensive computational model of both phenomena is a considerable challenge and until recently work in this area focused on one phenomenon and represented the behaviour of the other more simply. A single, finite volume unstructured mesh (FV-UM) spatial discretisation method has been employed on a single mesh for the entire domain. The Navier Stokes equations for fluid flow are solved using a SIMPLE type procedure and the Newmark b algorithm is employed for solving the dynamic equilibrium equations for linear elastic solid mechanics and mesh movement is achieved using a spring based mesh procedure for dynamic mesh movement. In the paper we describe a number of additional computation issues for the efficient and accurate modelling of three-dimensional, dynamic fluid-structure interaction problems.

A finite volume unstructured mesh approach to dynamic fluid-structure interaction: an assessment of the challenge of flutter analysis

Computational modelling of dynamic fluid-structure interaction (DFSI) is problematical since conventionally computational fluid dynamics (CFD) is solved using finite volume (FV) methods and computational structural mechanics (CSM) is based entirely on finite element (FE) methods. Hence, progress in modelling the emerging multi-physics problem of dynamic fluid-structure interaction in a consistent manner is frustrated and significant problems in computation convergence may be encountered in transferring and filtering data from one mesh and solution procedure to another, unless the fluid-structure coupling is either one way, very weak or both. This paper sets out the solution procedure for modelling the multi-physics dynamic fluid-structure interaction problem within a single software framework PHYSICA, using finite volume, unstructured mesh (FV-UM) procedures and will focus upon some of the problems and issues that have to be resolved for time accurate closely coupled dynamic fluid-structure flutter analysis.

Dynamic fluid-structure interactions using finite volume unstructured mesh procedures

A three-dimensional finite volume, unstructured mesh (FV-UM) method for dynamic fluid–structure interaction (DFSI) is described. Fluid structure interaction, as applied to flexible structures, has wide application in diverse areas such as flutter in aircraft, wind response of buildings, flows in elastic pipes and blood vessels. It involves the coupling of fluid flow and structural mechanics, two fields that are conventionally modelled using two dissimilar methods, thus a single comprehensive computational model of both phenomena is a considerable challenge. Until recently work in this area focused on one phenomenon and represented the behaviour of the other more simply. More recently, strategies for solving the full coupling between the fluid and solid mechanics behaviour have been developed. A key contribution has been made by Farhat et al. [Int. J. Numer. Meth. Fluids 21 (1995) 807] employing FV-UM methods for solving the Euler flow equations and a conventional finite element method for the elastic solid mechanics and the spring based mesh procedure of Batina [AIAA paper 0115, 1989] for mesh movement. In this paper, we describe an approach which broadly exploits the three field strategy described by Farhat for fluid flow, structural dynamics and mesh movement but, in the context of DFSI, contains a number of novel features: a single mesh covering the entire domain, a Navier–Stokes flow, a single FV-UM discretisation approach for both the flow and solid mechanics procedures, an implicit predictor–corrector version of the Newmark algorithm, a single code embedding the whole strategy.