Fabric-mechanical property relationships of trabecular bone allografts are altered by supercritical CO2 treatment and gamma sterilization

Tissue grafts are implanted in orthopedic surgery every day. In order to minimize infection risk, bone allografts are often delipidated with supercritical CO2 and sterilized prior to implantation. This treatment may, however, impair the mechanical behavior of the bone graft tissue. The goal of this study was to determine clinically relevant mechanical properties of treated/sterilized human trabecular bone grafts, e.g. the apparent modulus, strength, and the ability to absorb energy during compaction. They were compared with results of identical experiments performed previously on untreated/fresh frozen human trabecular bone from the same anatomical site (Charlebois, 2008). We tested the hypothesis that the morphology–mechanical property relationships of treated cancellous allografts are similar to those of fresh untreated bone. The morphology of the allografts was determined by μCT. Subsequently, cylindrical samples were tested in unconfined and confined compression. To account for various morphologies, the experimental data was fitted to phenomenological mechanical models for elasticity, strength, and dissipated energy density based on bone volume fraction (BV/TV) and the fabric tensor determined by MIL. The treatment/sterilization process does not appear to influence bone graft stiffness. However, strength and energy dissipation of the bone grafts were found to be significantly reduced by 36% to 47% and 66% to 81%, respectively, for a broad range of volume fraction (0.14 < BV/TV < 0.39) and degree of anisotropy (1.24 < DA < 2.18). Since the latter properties are strongly dominated by BV/TV, the clinical consequences of this reduction can be compensated by using grafts with lower porosity. The data of this study suggests that an increase of 5–10% in BV/TV is sufficient to compensate for the reduced post-yield mechanical properties of treated/sterilized bone in monotonic compression. In applications where graft stiffness needs to be matched and strength is not a concern, treated allograft with the same BV/TV as an appropriate fresh bone graft may be used.

Prediction of water–rock interaction and porosity evolution in a granitoid-hosted enhanced geothermal system, using constraints from the 5 km Basel-1 well

Numerical simulations based on plans for a deep geothermal system in Basel, Switzerland are used here to understand chemical processes that occur in an initially dry granitoid reservoir during hydraulic stimulation and long-term water circulation to extract heat. An important question regarding the sustainability of such enhanced geothermal systems (EGS), is whether water–rock reactions will eventually lead to clogging of flow paths in the reservoir and thereby reduce or even completely block fluid throughput. A reactive transport model allows the main chemical reactions to be predicted and the resulting evolution of porosity to be tracked over the expected 30-year operational lifetime of the system. The simulations show that injection of surface water to stimulate fracture permeability in the monzogranite reservoir at 190 °C and 5000 m depth induces redox reactions between the oxidised surface water and the reduced wall rock. Although new calcite, chlorite, hematite and other minerals precipitate near the injection well, their volumes are low and more than compensated by those of the dissolving wall-rock minerals. Thus, during stimulation, reduction of injectivity by mineral precipitation is unlikely. During the simulated long-term operation of the system, the main mineral reactions are the hydration and albitization of plagioclase, the alteration of hornblende to an assemblage of smectites and chlorites and of primary K-feldspar to muscovite and microcline. Within a closed-system doublet, the composition of the circulated fluid changes only slightly during its repeated passage through the reservoir, as the wall rock essentially undergoes isochemical recrystallization. Even after 30 years of circulation, the calculations show that porosity is reduced by only ∼0.2%, well below the expected fracture porosity induced by stimulation. This result suggests that permeability reduction owing to water–rock interaction is unlikely to jeopardize the long-term operation of deep, granitoid-hosted EGS systems. A peculiarity at Basel is the presence of anhydrite as fracture coatings at ∼5000 m depth. Simulated exposure of the circulating fluid to anhydrite induces a stronger redox disequilibrium in the reservoir, driving dissolution of ferrous minerals and precipitation of ferric smectites, hematite and pyrite. However, even in this scenario the porosity reduction is at most 0.5%, a value which is unproblematic for sustainable fluid circulation through the reservoir.

Characterization of granite matrix porosity and pore-space geometry by in situ and laboratory methods

Most available studies of interconnected matrix porosity of crystalline rocks are based on laboratory investigations; that is, work on samples that have undergone stress relaxation and were affected by drilling and sample preparation. The extrapolation of the results to in situ conditions is therefore associated with considerable uncertainty, and this was the motivation to conduct the ‘in situ Connected Porosity’ experiment at the Grimsel Test Site (Central Swiss Alps). An acrylic resin doped with fluorescent agents was used to impregnate the microporous granitic matrix in situ around an injection borehole, and samples were obtained by overcoring. The 3-D structure of the porespace, represented by microcracks, was studied by U-stage fluorescence microscopy. Petrophysical methods, including the determination of porosity, permeability and P -wave velocity, were also applied. Investigations were conducted both on samples that were impregnated in situ and on non-impregnated samples, so that natural features could be distinguished from artefacts. The investigated deformed granites display complex microcrack populations representing a polyphase deformation at varying conditions. The crack population is dominated by open cleavage cracks in mica and grain boundary cracks. The porosity of non-impregnated samples lies slightly above 1 per cent, which is 2–2.5 times higher than the in situ porosity obtained for impregnated samples. Measurements of seismic velocities (Vp ) on spherical rock samples as a function of confining pressure, spatial direction and water saturation for both non-impregnated and impregnated samples provide further constraints on the distinction between natural and induced crack types. The main conclusions are that (1) an interconnected network of microcracks exists in the whole granitic matrix, irrespective of the distance to ductile and brittle shear zones, and (2) conventional laboratory methods overestimate the matrix porosity. Calculations of contaminant transport through fractured media often rely on matrix diffusion as a retardation mechanism.

Mixing-cell boundary conditions and apparent mass balance errors for advective-dispersive solute transport

In many field or laboratory situations, well-mixed reservoirs like, for instance, injection or detection wells and gas distribution or sampling chambers define boundaries of transport domains. Exchange of solutes or gases across such boundaries can occur through advective or diffusive processes. First we analyzed situations, where the inlet region consists of a well-mixed reservoir, in a systematic way by interpreting them in terms of injection type. Second, we discussed the mass balance errors that seem to appear in case of resident injections. Mixing cells (MC) can be coupled mathematically in different ways to a domain where advective-dispersive transport occurs: by assuming a continuous solute flux at the interface (flux injection, MC-FI), or by assuming a continuous resident concentration (resident injection). In the latter case, the flux leaving the mixing cell can be defined in two ways: either as the value when the interface is approached from the mixing-cell side (MC-RT -), or as the value when it is approached from the column side (MC-RT +). Solutions of these injection types with constant or-in one case-distance-dependent transport parameters were compared to each other as well as to a solution of a two-layer system, where the first layer was characterized by a large dispersion coefficient. These solutions differ mainly at small Peclet numbers. For most real situations, the model for resident injection MC-RI + is considered to be relevant. This type of injection was modeled with a constant or with an exponentially varying dispersion coefficient within the porous medium. A constant dispersion coefficient will be appropriate for gases because of the Eulerian nature of the usually dominating gaseous diffusion coefficient, whereas the asymptotically growing dispersion coefficient will be more appropriate for solutes due to the Lagrangian nature of mechanical dispersion, which evolves only with the fluid flow. Assuming a continuous resident concentration at the interface between a mixing cell and a column, as in case of the MC-RI + model, entails a flux discontinuity. This flux discontinuity arises inherently from the definition of a mixing cell: the mixing process is included in the balance equation, but does not appear in the description of the flux through the mixing cell. There, only convection appears because of the homogeneous concentration within the mixing cell. Thus, the solute flux through a mixing cell in close contact with a transport domain is generally underestimated. This leads to (apparent) mass balance errors, which are often reported for similar situations and erroneously used to judge the validity of such models. Finally, the mixing cell model MC-RI + defines a universal basis regarding the type of solute injection at a boundary. Depending on the mixing cell parameters, it represents, in its limits, flux as well as resident injections. (C) 1998 Elsevier Science B.V. All rights reserved.

Investigation of young water inflow in karst aquifers using SF6–CFC–3H/He–85Kr–39Ar and stable isotope components

Karst aquifers are known for their wide distribution of water transfer velocities. From this observation, a multiple geochemical tracer approach seems to be particularly well suited to provide a significant assessment of groundwater flows, but the choice of adapted tracers is essential. In this study, several common tracers in karst aquifers such as physicochemical parameters, major ions, stable isotopes, and d13C to more specific tracers such as dating tracers – 14C, 3H, 3H–3He, CFC-12, SF6 and 85Kr, and 39Ar – were used, in a fractured karstic carbonated aquifer located in Burgundy (France). The information carried by each tracer and the best sampling strategy are compared on the basis of geochemical monitoring done during several recharge events and over longer time periods (months to years). This study’s results demonstrate that at the seasonal and recharge event time scale, the variability of concentrations is low for most tracers due to the broad spectrum of groundwater mixings. The tracers used traditionally for the study of karst aquifers, i.e., physicochemical parameters and major ions, efficiently describe hydrological processes such as the direct and differed recharge, but require being monitored at short time steps during recharge events to be maximized. From stable isotopes, tritium, and Cl� contents, the proportion of the fast direct recharge by the largest porosity was estimated using a binary mixing model. The use of tracers such as CFC-12, SF6, and 85Kr in karst aquifers provides additional information, notably an estimation of apparent age, but they require good preliminary knowledge of the karst system to interpret the results suitably. The CFC-12 and SF6 methods efficiently determine the apparent age of baseflow, but it is preferable to sample the groundwater during the recharge event. Furthermore, these methods are based on different assumptions such as regional enrichment in atmospheric SF6, excess air, and flow models among others. 85Kr and 39Ar concentrations can potentially provide a more direct estimation of groundwater residence time. Conversely, the 3H–3He method is inefficient in the karst aquifer for dating due to 3He degassing.

Diffusion-weighted MR imaging of upper abdominal organs: field strength and intervendor variability of apparent diffusion coefficients

PURPOSE To determine the variability of apparent diffusion coefficient (ADC) values in various anatomic regions in the upper abdomen measured with magnetic resonance (MR) systems from different vendors and with different field strengths. MATERIALS AND METHODS Ten healthy men (mean age, 36.6 years ± 7.7 [standard deviation]) gave written informed consent to participate in this prospective ethics committee-approved study. Diffusion-weighted (DW) MR imaging was performed in each subject with 1.5- and 3.0-T MR systems from each of three vendors at two institutions. Two readers independently measured ADC values in seven upper abdominal regions (left and right liver lobe, gallbladder, pancreas, spleen, and renal cortex and medulla). ADC values were tested for interobserver differences, as well as for differences related to field strength and vendor, with repeated-measures analysis of variance; coefficients of variation (CVs) and variance components were calculated. RESULTS Interreader agreement was excellent (intraclass coefficient, 0.876). ADC values were (77.5-88.8) ×10(-5) mm(2)/sec in the spleen and (250.6-278.5) ×10(-5) mm(2)/sec in the gallbladder. There were no significant differences between ADC values measured at 1.5 T and those measured at 3.0 T in any anatomic region (P >.10 for all). In two of seven regions at 1.5 T (left and right liver lobes, P < .023) and in four of seven regions at 3.0 T (left liver lobe, pancreas, and renal cortex and medulla, P < .008), intervendor differences were significant. CVs ranged from 7.0% to 27.1% depending on the anatomic location. CONCLUSION Despite significant intervendor differences in ADC values of various anatomic regions of the upper abdomen, ADC values of the gallbladder, pancreas, spleen, and kidney may be comparable between MR systems from different vendors and between different field strengths.

Dissecting abdominal aortic aneurysm in Ang II-infused mice: suprarenal branch ruptures and apparent luminal dilatation.

AIMS In this work, we provide novel insight into the morphology of dissecting abdominal aortic aneurysms in angiotensin II-infused mice. We demonstrate why they exhibit a large variation in shape and, unlike their human counterparts, are located suprarenally rather than infrarenally. METHODS AND RESULTS We combined synchrotron-based, ultra-high resolution ex vivo imaging (phase contrast X-Ray tomographic microscopy) with in vivo imaging (high-frequency ultrasound and contrast-enhanced micro-CT) and image-guided histology. In all mice, we observed a tear in the tunica media of the abdominal aorta near the ostium of the celiac artery. Independently we found that, unlike the gradual luminal expansion typical for human aneurysms, the outer diameter increase of angiotensin II-induced dissecting aneurysms in mice was related to one or several intramural haematomas. These were caused by ruptures of the tunica media near the ostium of small suprarenal side branches, which had never been detected by the established small animal imaging techniques. The tear near the celiac artery led to apparent luminal dilatation, while the intramural haematoma led to a dissection of the tunica adventitia on the left suprarenal side of the aorta. The number of ruptured branches was higher in those aneurysms that extended into the thoracic aorta, which explained the observed variability in aneurysm shape. CONCLUSION Our results are the first to describe apparent luminal dilatation, suprarenal branch ruptures, and intramural haematoma formation in dissecting abdominal aortic aneurysms in mice. Moreover, we validate and demonstrate the vast potential of phase contrast X-ray tomographic microscopy in cardiovascular small animal applications.

Implementation and evaluation of permeability-porosity and tortuosity-porosity relationships linked to mineral dissolution-precipitation

Changes of porosity, permeability, and tortuosity due to physical and geochemical processes are of vital importance for a variety of hydrogeological systems, including passive treatment facilities for contaminated groundwater, engineered barrier systems (EBS), and host rocks for high-level nuclear waste (HLW) repositories. Due to the nonlinear nature and chemical complexity of the problem, in most cases, it is impossible to verify reactive transport codes analytically, and code intercomparisons are the most suitable method to assess code capabilities and model performance. This paper summarizes model intercomparisons for six hypothetical scenarios with generally increasing geochemical or physical complexity using the reactive transport codes CrunchFlow, HP1, MIN3P, PFlotran, and TOUGHREACT. Benchmark problems include the enhancement of porosity and permeability through mineral dissolution, as well as near complete clogging due to localized mineral precipitation, leading to reduction of permeability and tortuosity. Processes considered in the benchmark simulations are advective-dispersive transport in saturated media, kinetically controlled mineral dissolution-precipitation, and aqueous complexation. Porosity changes are induced by mineral dissolution-precipitation reactions, and the Carman-Kozeny relationship is used to describe changes in permeability as a function of porosity. Archie’s law is used to update the tortuosity and the pore diffusion coefficient as a function of porosity. Results demonstrate that, generally, good agreement is reached amongst the computer models despite significant differences in model formulations. Some differences are observed, in particular for the more complex scenarios involving clogging; however, these differences do not affect the interpretation of system behavior and evolution.

Reconstruction of in-situ porosity and porewater compositions of low-permeability crystalline rocks: Magnitude of artefacts induced by drilling and sample recovery

Geological site characterisation programmes typically rely on drill cores for direct information on subsurface rocks. However, porosity, transport properties and porewater composition measured on drill cores can deviate from in-situ values due to two main artefacts caused by drilling and sample recovery: (1) mechanical disruption that increases porosity and (2) contamination of the porewater by drilling fluid. We investigated the effect and magnitude of these perturbations on large drill core samples (12–20 cm long, 5 cmdiameter) of high-grade, granitic gneisses obtained from 350 to 600 m depth in a borehole on Olkiluoto Island (SW Finland). The drilling fluid was traced with sodium–iodide. By combining out-diffusion experiments, gravimetry, UV-microscopy and iodide mass balance calculations, we successfully quantified the magnitudes of the artefacts: 2–6% increase in porosity relative to the bulk connected porosity and 0.9 to 8.9 vol.% contamination by drilling fluid. The spatial distribution of the drilling-induced perturbations was revealed by numerical simulations of 2D diffusion matched to the experimental data. This showed that the rims of the samples have a mechanically disrupted zone 0.04 to 0.22 cm wide, characterised by faster transport properties compared to the undisturbed centre (1.8 to 7.7 times higher pore diffusion coefficient). Chemical contamination was shown to affect an even wider zone in all samples, ranging from 0.15 to 0.60 cm, inwhich iodide enrichmentwas up to 180 mg/kgwater, compared to 0.5 mg/kgwater in the uncontaminated centre. For all samples in the present case study, it turned out that the magnitude of the artefacts caused by drilling and sample recovery is so small that no correction is required for their effects. Therefore, the standard laboratory measurements of porosity, transport properties and porewater composition can be taken as valid in-situ estimates. However, it is clear that the magnitudes strongly depend on site- and drilling-specific factors and therefore our results cannot be transferred simply to other locations. We recommend the approach presented in this study as a route to obtain reliable values in future drilling campaigns aimed at characterising in-situ bedrock properties.