An objective anisotropic elastic plastic model and Algorithm applicable to bone mechanics

The implicit projection algorithm of isotropic plasticity is extended to an objective anisotropic elastic perfectly plastic model. The recursion formula developed to project the trial stress on the yield surface, is applicable to any non linear elastic law and any plastic yield function.A curvilinear transverse isotropic model based on a quadratic elastic potential and on Hill's quadratic yield criterion is then developed and implemented in a computer program for bone mechanics perspectives. The paper concludes with a numerical study of a schematic bone-prosthesis system to illustrate the potential of the model.

Permeability estimates from poro-elastic dispersion analyses of sonic logs

Modern sonic logging tools designed for shallow environmental and engineering applications allow for P-wave phase velocity measurements over a wide frequency band. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and source frequencies ranging from approximately 1 to 30 kHz, the observable poro-elastic P-wave velocity dispersion is sufficiently pronounced to allow for reliable first-order estimations of the underlying permeability structure. These predictions have been tested on and verified for a surficial alluvial aquifer. Our results indicate that, even without any further calibration, the thus obtained permeability estimates as well as their variabilities within the pertinent lithological units are remarkably close to those expected based on the corresponding granulometric characteristics.

Mendacastrum n. gen. and Domuzdagia n. gen., two Jurassic spherical Spumellaria (Radiolaria) with hagiastrid medullary shell

Two new spumellarian radiolarian genera, Mendacastrum and Domuzdagia, are described from the lower Tithonian and lower Pliensbachian respectively. Both have a spherical cortical shell of actinommid type and a spherical or subspherical double medullary shell with the inner medullary shell of hagiastrid s.l. type. The inner medullary shell of Mendacastrum is of dactyliosphaerid or higumastrid s. sit. type, whereas that of Domuzdagia is of angulobracchiid type. Since they cannot be assigned to any described Mesozoic pyloniacean families, they are considered as type genera of two new families: Mendacastridae and Domuzdagiidae respectively.

Transhepatic ipsilateral right portal vein embolization extended to segment IV: improving hypertrophy and resection outcomes with spherical particles and coils.

PURPOSE: To analyze outcomes after right portal vein embolization extended to segment IV (right PVE + IV) before extended right hepatectomy, including liver hypertrophy, resection rates, and complications after embolization and resection, and to assess differences in outcomes with two different particulate embolic agents. MATERIALS AND METHODS: Between 1998 and 2004, transhepatic ipsilateral right PVE + IV with particles and coils was performed in 44 patients with malignant hepatobiliary disease, including metastases (n = 24), biliary cancer (n = 14), and hepatocellular carcinoma (n = 6). Right PVE + IV was considered if the future liver remnant (FLR; segments II/III with or without I) was less than 25% of the total estimated liver volume (TELV). Tris-acryl microspheres (100-700 microm; n = 21) or polyvinyl alcohol (PVA) particles (355-1,000 microm; n = 23) were administered in a stepwise fashion. Smaller particles were used to occlude distal branches, followed by larger particles to occlude proximal branches until near-complete stasis. Coils were then placed in secondary portal branches. Computed tomographic volumetry was performed before and 3-4 weeks after right PVE + IV to assess FLR hypertrophy. Liver volumes and postembolization and postoperative outcomes were measured. RESULTS: After right PVE + IV with PVA particles, FLR volume increased 45.5% +/- 40.9% and FLR/TELV ratio increased 6.9% +/- 5.6%. After right PVE + IV with tris-acryl microspheres, FLR volume increased 69.0% +/- 30.7% and FLR/TELV ratio increased 9.7% +/- 3.3%. Differences in FLR volume (P = .0011), FLR/TELV ratio (P = .027), and resection rates (P = .02) were statistically significant. Seventy-one percent of patients underwent extended right hepatectomy (86% after receiving tris-acryl microspheres, 57% after receiving PVA). Thirteen patients (29%) did not undergo resection (extrahepatic spread [n = 9], inadequate hypertrophy [n = 3], other reasons [n = 1]). No patient developed postembolization syndrome or progressive liver insufficiency after embolization or resection. One death after resection occurred as a result of sepsis and hemorrhage. Median hospital stays were 1 day after right PVE + IV and 7 days after resection. CONCLUSION: Transhepatic ipsilateral right PVE + IV with use of particles and coils is a safe, effective method for inducing contralateral hypertrophy before extended right hepatectomy. Embolization with small spherical particles provides improved hypertrophy and resection rates compared with larger, nonspherical particles.

Traitement des fractures de l'enfant par enclouage centro-médullaire élastique stable (ECMES) [Treatment of pediatric fractures by intramedullary stable elastic pinning]

During the past decade several new techniques for the treatment of children's fractures respecting the specificity of the growing bone have been described. The goal of all these techniques was to mechanically stabilise the fracture however to preserve a certain instability of the fracture gap itself inducing early callus formation and subsequent consolidation. The dynamic external fixation as well as the elastic stable intramedullary pinning have become accepted means in the treatment of long bone fractures in the paediatric age group. We report our experience of the last seven years with the intramedullary pinning of 105 fractures. Eighty-four were fractures of the femur, 9 of the humerus, 8 of the forearm, and a further 4 of the tibial shaft. The intramedullary elastic pinning represents a simple technique which supports or even enhances the natural process of fracture healing of the growing bone. The method is not very invasive, is cost effective, and allows short hospitalisation. Early physical activity is guaranteed due to early consolidation of the fracture. Complications are rare and the final orthopedic and cosmetic outcome is excellent.

A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid-solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently bench-marked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

Quantitative comparison of simulations of seismic wave propagation in heterogeneous poro-elastic media involving fluid-solid interfaces and in equivalent visco-elastic solids

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes the predominant attenuation mechanism at seismic frequencies. As a consequence, centimeter-scale perturbations of the subsurface physical properties should be taken into account for seismic modeling whenever detailed and accurate responses of the target structures are desired. This is, however, computationally prohibitive since extremely small grid spacings would be necessary. A convenient way to circumvent this problem is to use an upscaling procedure to replace the heterogeneous porous media by equivalent visco-elastic solids. In this work, we solve Biot's equations of motion to perform numerical simulations of seismic wave propagation through porous media containing mesoscopic heterogeneities. We then use an upscaling procedure to replace the heterogeneous poro-elastic regions by homogeneous equivalent visco-elastic solids and repeat the simulations using visco-elastic equations of motion. We find that, despite the equivalent attenuation behavior of the heterogeneous poro-elastic medium and the equivalent visco-elastic solid, the seismograms may differ due to diverging boundary conditions at fluid-solid interfaces, where there exist additional options for the poro-elastic case. In particular, we observe that the seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an interesting result, which has potentially important implications for wave-equation-based algorithms in exploration geophysics involving fluid-solid interfaces, such as, for example, wave field decomposition.

Chemical labelling of oyster shells used for time-calibrated high-resolution Mg/Ca ratios: A tool for estimation of past seasonal temperature variations

The geochemical compositions of biogenic carbonates are increasingly used for palaeoenvironmental reconstructions. The skeletal delta O-18 temperature relationship is dependent on water salinity, so many recent studies have focused on the Mg/Ca and Sr/Ca ratios because those ratios in water do not change significantly on short time scales. Thus, those elemental ratios are considered to be good palaeotemperature proxies in many biominerals, although their use remains ambiguous in bivalve shells. Here, we present the high-resolution Mg/Ca ratios of two modern species of juvenile and adult oyster shells, Crassostrea gigas and Ostrea edulis. These specimens were grown in controlled conditions for over one year in two different locations. In situ monthly Mn-marking of the shells has been used for day calibration. The daily Mg/Ca.ratios in the shell have been measured with an electron microprobe. The high frequency Mg/Ca variation of all specimens displays good synchronism with lunar cycles, suggesting that tides strongly influence the incorporation of Mg/Ca into the shells. Highly significant correlation coefficients (0.70<R<0.83, p<0.0001) between the Mg/Ca ratios and the seawater temperature are obtained only for juvenile C. gigas samples, while metabolic control of Mg/Ca incorporation and lower shell growth rates preclude the use of the Mg/Ca ratio in adult shells as a palaeothermometer. Data from three juvenile C. gigas shells from the two study sites are selected to establish a relationship: T = 3.77Mg/Ca + 1.88, where T is in degrees C and Mg/Ca in mmol/mol. (c) 2012 Elsevier B.V. All rights reserved.

A pseudo-spectral method for simulating of poro-elastic seismic wave propagation in complex borehole environments

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in cylindrical coordinates. An important application of this method is the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh consisting of three concentric domains representing the borehole fluid in the center, the borehole casing and the surrounding porous formation. The spatial discretization is based on a Chebyshev expansion in the radial direction, Fourier expansions in the other directions, and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method based on the method of characteristics is used to match the boundary conditions at the fluid/porous-solid and porous-solid/porous-solid interfaces. The viability and accuracy of the proposed method has been tested and verified in 2D polar coordinates through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently benchmarked solution for 2D Cartesian coordinates. The proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is handled adequately.

Constraints on the permeability structure of alluvial aquifers from the poro-elastic inversion of multi-frequency P-wave velocity logs

We have explored the possibility of obtaining first-order permeability estimates for saturated alluvial sediments based on the poro-elastic interpretation of the P-wave velocity dispersion inferred from sonic logs. Modern sonic logging tools designed for environmental and engineering applications allow one for P-wave velocity measurements at multiple emitter frequencies over a bandwidth covering 5 to 10 octaves. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and typical emitter frequencies ranging from approximately 1 to 30 kHz, the observable velocity dispersion should be sufficiently pronounced to allow one for reliable first-order estimations of the permeability structure. The corresponding predictions have been tested on and verified for a borehole penetrating a typical surficial alluvial aquifer. In addition to multifrequency sonic logs, a comprehensive suite of nuclear and electrical logs, an S-wave log, a litholog, and a limited number laboratory measurements of the permeability from retrieved core material were also available. This complementary information was found to be essential for parameterizing the poro-elastic inversion procedure and for assessing the uncertainty and internal consistency of corresponding permeability estimates. Our results indicate that the thus obtained permeability estimates are largely consistent with those expected based on the corresponding granulometric characteristics, as well as with the available evidence form laboratory measurements. These findings are also consistent with evidence from ocean acoustics, which indicate that, over a frequency range of several orders-of-magnitude, the classical theory of poro-elasticity is generally capable of explaining the observed P-wave velocity dispersion in medium- to fine-grained seabed sediments

Quantitative comparison between simulations of seismic wave propagation in heterogeneous poro-elastic media and equivalent visco-elastic solids for marine-type environments

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes an important seismic attenuation mechanism in porous rocks. As a consequence, centimetre-scale perturbations of the rock physical properties should be taken into account for seismic modelling whenever detailed and accurate responses of specific target structures are desired, which is, however, computationally prohibitive. A convenient way to circumvent this problem is to use an upscaling procedure to replace each of the heterogeneous porous media composing the geological model by corresponding equivalent visco-elastic solids and to solve the visco-elastic equations of motion for the inferred equivalent model. While the overall qualitative validity of this procedure is well established, there are as of yet no quantitative analyses regarding the equivalence of the seismograms resulting from the original poro-elastic and the corresponding upscaled visco-elastic models. To address this issue, we compare poro-elastic and visco-elastic solutions for a range of marine-type models of increasing complexity. We found that despite the identical dispersion and attenuation behaviour of the heterogeneous poro-elastic and the equivalent visco-elastic media, the seismograms may differ substantially due to diverging boundary conditions, where there exist additional options for the poro-elastic case. In particular, we observe that at the fluid/porous-solid interface, the poro- and visco-elastic seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an important result which has potentially far-reaching implications for wave-equation-based algorithms in exploration geophysics involving fluid/porous-solid interfaces, such as, for example, wavefield decomposition.