Three-dimensional models of mantle convection : influence of phase transitions and temperature-dependent viscosity

Two of the most important questions in mantle dynamics are investigated in three separate studies: the influence of phase transitions (studies 1 and 2), and the influence of temperature-dependent viscosity (study 3).

(1) Numerical modeling of mantle convection in a three-dimensional spherical shell incorporating the two major mantle phase transitions reveals an inherently three-dimensional flow pattern characterized by accumulation of cold downwellings above the 670 km discontinuity, and cylindrical 'avalanches' of upper mantle material into the lower mantle. The exothermic phase transition at 400 km depth reduces the degree of layering. A region of strongly-depressed temperature occurs at the base of the mantle. The temperature field is strongly modulated by this partial layering, both locally and in globally-averaged diagnostics. Flow penetration is strongly wavelength-dependent, with easy penetration at long wavelengths but strong inhibition at short wavelengths. The amplitude of the geoid is not significantly affected.

(2) Using a simple criterion for the deflection of an upwelling or downwelling by an endothermic phase transition, the scaling of the critical phase buoyancy parameter with the important lengthscales is obtained. The derived trends match those observed in numerical simulations, i.e., deflection is enhanced by (a) shorter wavelengths, (b) narrower up/downwellings (c) internal heating and (d) narrower phase loops.

(3) A systematic investigation into the effects of temperature-dependent viscosity on mantle convection has been performed in three-dimensional Cartesian geometry, with a factor of 1000-2500 viscosity variation, and Rayleigh numbers of 10^5-10^7. Enormous differences in model behavior are found, depending on the details of rheology, heating mode, compressibility and boundary conditions. Stress-free boundaries, compressibility, and temperature-dependent viscosity all favor long-wavelength flows, even in internally heated cases. However, small cells are obtained with some parameter combinations. Downwelling plumes and upwelling sheets are possible when viscosity is dependent solely on temperature. Viscous dissipation becomes important with temperature-dependent viscosity.

The sensitivity of mantle flow and structure to these various complexities illustrates the importance of performing mantle convection calculations with rheological and thermodynamic properties matching as closely as possible those of the Earth.

Comparison of polar ozone loss rates simulated by one-dimensional and three-dimensional models with Match observations in recent Antarctic and Arctic winters

Simulations of ozone loss rates using a three-dimensional chemical transport model and a box model during recent Antarctic and Arctic winters are compared with experimental loss rates. The study focused on the Antarctic winter 2003, during which the first Antarctic Match campaign was organized, and on Arctic winters 1999/2000, 2002/2003. The maximum ozone loss rates retrieved by the Match technique for the winters and levels studied reached 6 ppbv/sunlit hour and both types of simulations could generally reproduce the observations at 2-sigma error bar level. In some cases, for example, for the Arctic winter 2002/2003 at 475 K level, an excellent agreement within 1-sigma standard deviation level was obtained. An overestimation was also found with the box model simulation at some isentropic levels for the Antarctic winter and the Arctic winter 1999/2000, indicating an overestimation of chlorine activation in the model. Loss rates in the Antarctic show signs of saturation in September, which have to be considered in the comparison. Sensitivity tests were performed with the box model in order to assess the impact of kinetic parameters of the ClO-Cl2O2 catalytic cycle and total bromine content on the ozone loss rate. These tests resulted in a maximum change in ozone loss rates of 1.2 ppbv/sunlit hour, generally in high solar zenith angle conditions. In some cases, a better agreement was achieved with fastest photolysis of Cl2O2 and additional source of total inorganic bromine but at the expense of overestimation of smaller ozone loss rates derived later in the winter.

Automated detection of the osseous acetabular rim using three-dimensional models of the pelvis

An automated algorithm for detection of the acetabular rim was developed. Accuracy of the algorithm was validated in a sawbone study and compared against manually conducted digitization attempts, which were established as the ground truth. The latter proved to be reliable and reproducible, demonstrated by almost perfect intra- and interobserver reliability. Validation of the automated algorithm showed no significant difference compared to the manually acquired data in terms of detected version and inclination. Automated detection of the acetabular rim contour and the spatial orientation of the acetabular opening plane can be accurately achieved with this algorithm.

Three dimensional models of water flow in the Persian Gulf

Persian Gulf region is globally of great importance due to its economical and political reasons. The importance lies in oil sources and sea exports. Geophysical phenomena dominated in the water circulation affected this region is called Monsoon it stretches from African coasts to the half way of Red Seal affected all coasts of Persian Gulf and goes toward east to the Indian ocean. Other essential factors in the water circulation in this region are net evaporation (several meters in per year), high density and high salinity. In this article the effects of wind stress and evaporation in the water circulation in the region will be considered and model equations for wind forces, density, pressure, gradient, and bottom friction for Persian Gulf will be discussed.

Three-dimensional printing of models for preoperative planning and simulation of transcatheter valve replacement

In this study, we show the use of three-dimensional printing models for preoperative planning of transcatheter valve replacement in a patient with an extreme porcelain aorta. A 70-year-old man with severe aortic stenosis and a porcelain aorta was referred to our center for transcatheter aortic valve replacement. Unfortunately, the patient died after the procedure because of a potential ischemic event. Therefore, we decided to fabricate three-dimensional models to evaluate the potential effects of these constructs for previous surgical planning and simulation of the transcatheter valve replacement.

Automated fit quantification of tibial nail designs during the insertion using computer three-dimensional modelling

Intramedullary nailing is the standard fixation method for displaced diaphyseal fractures of the tibia. An optimal nail design should both facilitate insertion and anatomically fit the bone geometry at its final position in order to reduce the risk of stress fractures and malalignments. Due to the nonexistence of suitable commercial software, we developed a software tool for the automated fit assessment of nail designs. Furthermore, we demonstrated that an optimised nail, which fits better at the final position, is also easier to insert. Three-dimensional models of two nail designs and 20 tibiae were used. The fitting was quantified in terms of surface area, maximum distance, sum of surface areas and sum of maximum distances by which the nail was protruding into the cortex. The software was programmed to insert the nail into the bone model and to quantify the fit at defined increment levels. On average, the misfit during the insertion in terms of the four fitting parameters was smaller for the Expert Tibial Nail Proximal bend (476.3 mm2, 1.5 mm, 2029.8 mm2, 6.5 mm) than the Expert Tibial Nail (736.7 mm2, 2.2 mm, 2491.4 mm2, 8.0 mm). The differences were statistically significant (p ≤ 0.05). The software could be used by nail implant manufacturers for the purpose of implant design validation.

CFD-DEM simulation of three-dimensional aeolian sand movement

A three-dimensional CFD-DEM model is proposed to investigate the aeolian sand movement. The results show that the mean particle horizontal velocity can be expressed by a power function of heights. The probability distribution of the impact and lift-off velocities of particles can be described by a log-normal function, and that of the impact and lift-off angles can be expressed by an exponential function. The probability distribution of particle horizontal velocity at different heights can be described as a lognormal function, while the probability distribution of longitudinal and vertical velocity can be described as a normal function. The comparison with previous two-dimensional calculations shows that the variations of mean particle horizontal velocity along the heights in two-dimensional and three-dimensional models are similar. However, the mean particle density of the two-dimensional model is larger than that in reality, which will result in the overestimation of sand transportation rate in the two-dimensional calculation. The study also shows that the predicted probability distributions of particle velocities are in good agreement with the experimental results.

MAGNETOHYDRODYNAMIC SIMULATIONS OF RECONNECTION AND PARTICLE ACCELERATION: THREE-DIMENSIONAL EFFECTS

Magnetic fields can change their topology through a process known as magnetic reconnection. This process in not only important for understanding the origin and evolution of the large-scale magnetic field, but is seen as a possibly efficient particle accelerator producing cosmic rays mainly through the first-order Fermi process. In this work we study the properties of particle acceleration inserted in reconnection zones and show that the velocity component parallel to the magnetic field of test particles inserted in magnetohydrodynamic (MHD) domains of reconnection without including kinetic effects, such as pressure anisotropy, the Hall term, or anomalous effects, increases exponentially. Also, the acceleration of the perpendicular component is always possible in such models. We find that within contracting magnetic islands or current sheets the particles accelerate predominantly through the first-order Fermi process, as previously described, while outside the current sheets and islands the particles experience mostly drift acceleration due to magnetic field gradients. Considering two-dimensional MHD models without a guide field, we find that the parallel acceleration stops at some level. This saturation effect is, however, removed in the presence of an out-of-plane guide field or in three-dimensional models. Therefore, we stress the importance of the guide field and fully three-dimensional studies for a complete understanding of the process of particle acceleration in astrophysical reconnection environments.

Influence of Cusp Inclination on Stress Distribution in Implant-Supported Prostheses. A Three-Dimensional Finite Element Analysis

Purpose: The aim of this study was to assess the influence of cusp inclination on stress distribution in implant-supported prostheses by 3D finite element method.Materials and Methods: Three-dimensional models were created to simulate a mandibular bone section with an implant (3.75 mm diameter x 10 mm length) and crown by means of a 3D scanner and 3D CAD software. A screw-retained single crown was simulated using three cusp inclinations (10 degrees, 20 degrees, 30 degrees). The 3D models (model 10d, model 20d, and model 30d) were transferred to the finite element program NeiNastran 9.0 to generate a mesh and perform the stress analysis. An oblique load of 200 N was applied on the internal vestibular face of the metal ceramic crown.Results: The results were visualized by means of von Mises stress maps. Maximum stress concentration was located at the point of application. The implant showed higher stress values in model 30d (160.68 MPa). Cortical bone showed higher stress values in model 10d (28.23 MPa).Conclusion: Stresses on the implant and implant/abutment interface increased with increasing cusp inclination, and stresses on the cortical bone decreased with increasing cusp inclination.

Passivity Versus Unilateral Angular Misfit: Evaluation of Stress Distribution on Implant-Supported Single Crowns: Three-Dimensional Finite Element Analysis

The aim of this study was to evaluate the effect of unilateral angular misfit of 100 Km on stress distribution of implant-supported single crowns with ceramic veneering and gold framework by three-dimensional finite element analysis. Two three-dimensional models representing a maxillary section of premolar region were constructed: group 1 (control)-crown completely adapted to the implant and group 2-crown with unilateral angular misfit of 100 Km. A vertical force of 100 N was applied on 2 centric points of the crown. The von Mises stress was used as an analysis criterion. The stress values and distribution in the main maps (204.4 MPa for group 1 and 205.0 MPa for group 2) and in the other structures (aesthetic veneering, framework, retention screw, implant, and bone tissue) were similar for both groups. The highest stress values were observed between the first and second threads of the retention screw. Considering the bone tissue, the highest stress values were exhibited in the peri-implant cortical bone. The unilateral angular misfit of 100 Km did not influence the stress distribution on the implant-supported prosthesis under static loading.

Three-Dimensional Finite Element Analysis of Vertical and Angular Misfit in Implant-Supported Fixed Prostheses

Purpose: Three-dimensional finite element analysis was used to evaluate the effect of vertical and angular misfit in three-piece implant-supported screw-retained fixed prostheses on the biomechanical response in the peri-implant bone, implants, and prosthetic components. Materials and Methods: Four three-dimensional models were fabricated to represent a right posterior mandibular section with one implant in the region of the second premolar (2PM) and another in the region of the second molar (2M). The implants were splinted by a three-piece implant-supported metal-ceramic prosthesis and differed according to the type of misfit, as represented by four different models: Control = prosthesis with complete fit to the implants; UAM (unilateral angular misfit) = prosthesis presenting unilateral angular misfit of 100 pm in the mesial region of the 2M; UVM (unilateral vertical misfit) = prosthesis presenting unilateral vertical misfit of 100 pm in the mesial region of the 2M; and TVM (total vertical misfit) = prosthesis presenting total vertical misfit of 100 pm in the platform of the framework in the 2M. A vertical load of 400 N was distributed and applied on 12 centric points by the software Ansys, ie, a vertical load of 150 N was applied to each molar in the prosthesis and a vertical load of 100 N was applied at the 2PM. Results: The stress values and distribution in peri-implant bone tissue were similar for all groups. The models with misfit exhibited different distribution patterns and increased stress magnitude in comparison to the control. The highest stress values in group UAM were observed in the implant body and retention screw. The groups UVM and TVM exhibited high stress values in the platform of the framework and the implant hexagon, respectively. Conclusions: The three types of misfit influenced the magnitude and distribution of stresses. The influence of misfit on peri-implant bone tissue was modest. Each type of misfit increased the stress values in different regions of the system. INT J ORAL MAXILLOFAC IMPLANTS 2011;26:788-796

Three-dimensional printing for perioperative planning of complex aortic arch surgery

Purpose In this study, we show the use of three-dimensional printing models for preoperative planning of surgery for patients with complex aortic arch anomalies. Description A 70-year-old man with an extensively arteriosclerotic aneurysm reaching from the ascending aorta to the descending aorta was referred to our center for complete aortic arch replacement. We visualized and reconstructed computed tomography data of the patient and fabricated a flexible three-dimensional model of the aortic arch including the aneurysm. Evaluation This model was very helpful for the preoperative decision making and planning of the frozen elephant trunk procedure owing to the exact and lifelike illustration of the native aortic arch. Conclusions Three-dimensional models are helpful in preoperative planning and postoperative evaluation of frozen elephant trunk procedures in patients with complex aortic anatomy.

Simulations of pollutant dispersion at toll plazas using three-dimensional CFD models

Toll plazas are particularly susceptible to build-ups of vehicle-emitted pollutants because vehicles pass through in low gear. To look at this, three-dimensional computational fluid dynamics simulations of pollutant dispersion are used on the standard k e turbulence model. The effects of wind speed, wind direction and topography on pollutant dispersion were discussed. The Wuzhuang toll plaza on the Hefei-Nanjing expressway is considered, and the effects of the retaining walls along both sides of the plaza on pollutant dispersion is analysed. There are greater pollutant concentrations near the tollbooths as the angle between the direction of the wind and traffic increases implying that retaining walls impede dispersion. The slope of the walls has little influence on the variations in pollutant concentration.

Three dimensional in-vitro models for studying cancer angiogenesis

Introduction Hydrogels prepared from star-shaped poly(ethylene glycol) (PEG) and maleimide-functionalized heparin provide a potential matrix for use in developing three dimensional (3D) models. We have previously demonstrated that these hydrogels support the cultivation of human umbilical vein endothelial cells (HUVECs). We extend this body of work to study the ability to create an extracellular matrix (ECM)-like model to study breast and prostate cancer cell growth in 3D. Also, we investigate the ability to produce a tri-culture mimicking tumour angiogenesis with cancer spheroids, HUVECs and mesenchymal stem cells (MSCs). Materials and Methods The breast cancer cell lines, MCF-7 and MDA-MB-231, and prostate cancer cell lines, LNCaP and PC3, were seeded into starPEG-heparin hydrogels and grown for 14 Days to analyze the effects of varying hydrogel stiffness on spheroid development. Resulting hydrogel constructs were analyzed via proliferation assays, light microscopy, and immunostaining. Cancer cell lines were then seeded into starPEG-heparin hydrogels functionalized with growth factors as spheroids with HUVECs and MSCs and grown as a tri-culture. Cultures were analyzed via immunostaining and observed using confocal microscopy. Results Cultures prepared in MMP-cleavable starPEG-heparin hydrogels display spheroid formation in contrast to adherent growth on tissue culture plastic. Small differences were visualized in cancer spheroid growth between different gel stiffness across the range of cell lines. Cancer cell lines were able to be co-cultivated with HUVECs and MSC. Interaction was visualized between tumours and HUVECs via confocal microscopy. Further studies intend to further optimize and mimic the ECM environment of in-situ tumour angiogenesis. Discussion Our results confirm the suitability of hydrogels constructed from starPEG-heparin for HUVEC and MSC co-cultivation with cancer cell lines to study cell-cell and cell-matrix interactions in a 3D environment. This represents a step forward in the development of 3D culture models to study the pathomechanisms of breast and prostate cancer.