Comparison of friction models applied to a control valve

Eight different models to represent the effect of friction in control valves are presented: four models based on physical principles and four empirical ones. The physical models, both static and dynamic, have the same structure. The models are implemented in Simulink/Matlab (R) and compared, using different friction coefficients and input signals. Three of the models were able to reproduce the stick-slip phenomenon and passed all the tests, which were applied following ISA standards. (C) 2008 Elsevier Ltd. All rights reserved.

Computational fluid dynamics investigation of a centrifugal blood pump

In the development of a ventricular assist device, computational fluid dynamics (CFD) analysis is an efficient tool to obtain the best design before making the final prototype. In this study, different designs of a centrifugal blood pump were developed to investigate flow characteristics and performance. This study assumed the blood flow as being an incompressible homogeneous Newtonian fluid. A constant velocity was applied at the inlet; no slip boundary conditions were applied at device wall; and pressure boundary conditions were applied at the outlet. The CFD code used in this work was based on the finite volume method. In the future, the results of CFD analysis can be compared with flow visualization and hemolysis tests.

A demonstration of artificial dissipation effects in some finite volume solution methods for the Navier-Stokes equations

The artificial dissipation effects in some solutions obtained with a Navier-Stokes flow solver are demonstrated. The solvers were used to calculate the flow of an artificially dissipative fluid, which is a fluid having dissipative properties which arise entirely from the solution method itself. This was done by setting the viscosity and heat conduction coefficients in the Navier-Stokes solvers to zero everywhere inside the flow, while at the same time applying the usual no-slip and thermal conducting boundary conditions at solid boundaries. An artificially dissipative flow solution is found where the dissipation depends entirely on the solver itself. If the difference between the solutions obtained with the viscosity and thermal conductivity set to zero and their correct values is small, it is clear that the artificial dissipation is dominating and the solutions are unreliable.

A relative gradient theory for layered materials

It is possible to remedy certain difficulties with the description of short wave length phenomena and interfacial slip in standard models of a laminated material by considering the bending stiffness of the layers. If the couple or moment stresses are assumed to be proportional to the relative deformation gradient, then the bending effect disappears for vanishing interface slip, and the model correctly reduces to an isotropic standard continuum. In earlier Cosserat-type models this was not the case. Laminated materials of the kind considered here occur naturally as layered rock, or at a different scale, in synthetic layered materials and composites. Similarities to the situation in regular dislocation structures with couple stresses, also make these ideas relevant to single slip in crystalline materials. Application of the theory to a one-dimensional model for layered beams demonstrates agreement with exact results at the extremes of zero and infinite interface stiffness. Moreover, comparison with finite element calculations confirm the accuracy of the prediction for intermediate interfacial stiffness.

Modelling the large deformations in stratified media - the Cosserat continuum approach

Methods employing continuum approximation in describing the deformation of layered materials possess a clear advantage over explicit models, However, the conventional implicit models based on the theory of anisotropic continua suffers from certain difficulties associated with interface slip and internal instabilities. These difficulties can be remedied by considering the bending stiffness of the layers. This implies the introduction of moment (couple) stresses and internal rotations, which leads to a Cosserat-type theory. In the present model, the behaviour of the layered material is assumed to be linearly elastic; the interfaces are assumed to be elastic perfectly plastic. Conditions of slip or no slip at the interfaces are detected by a Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformation analysis. The model is incorporated into the finite element program AFENA and validated against analytical solutions of elementary buckling problems in layered medium. A problem associated with buckling of the roof and the floor of a rectangular excavation in jointed rock mass under high horizontal in situ stresses is considered as the main application of the theory. Copyright (C) 1999 John Wiley & Sons, Ltd.

A numerical study of flexural buckling of foliated rock slopes

The occurrence of foliated rock masses is common in mining environment. Methods employing continuum approximation in describing the deformation of such rock masses possess a clear advantage over methods where each rock layer and each inter-layer interface (joint) is explicitly modelled. In devising such a continuum model it is imperative that moment (couple) stresses and internal rotations associated with the bending of the rock layers be properly incorporated in the model formulation. Such an approach will lead to a Cosserat-type theory. In the present model, the behaviour of the intact rock layer is assumed to be linearly elastic and the joints are assumed to be elastic-perfectly plastic. Condition of slip at the interfaces are determined by a Mohr-Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformations. The model is incorporated into the finite element program AFENA and validated against an analytical solution of elementary buckling problems of a layered medium under gravity loading. A design chart suitable for assessing the stability of slopes in foliated rock masses against flexural buckling failure has been developed. The design chart is easy to use and provides a quick estimate of critical loading factors for slopes in foliated rock masses. It is shown that the model based on Euler's buckling theory as proposed by Cavers (Rock Mechanics and Rock Engineering 1981; 14:87-104) substantially overestimates the critical heights for a vertical slope and underestimates the same for sub-vertical slopes. Copyright (C) 2001 John Wiley & Sons, Ltd.

Relative motions of Africa, Iberia and Europe during Alpine orogeny

A revised kinematic model for the motions of Africa and Iberia relative to Europe since the Middle Jurassic is presented in order to provide boundary conditions for Alpine-Mediterranean reconstructions. These motions were calculated using up-to-date kinematic data predominantly based on magnetic isochrons in the Atlantic Ocean and published by various authors during the last 15 years. It is shown that convergence of Africa with respect to Europe commenced during the Cretaceous Normal Superchron (CNS), between chrons MO and 34 (120-83 Ma). This motion was subjected to fluctuations in convergence rates characterised by two periods of relatively rapid convergence (during Late Cretaceous and Eocene-Oligocene times) that alternated with periods of slower convergence (during the Paleocene and since the Early Miocene). Distinct changes in plate kinematics are recognised in the motion of Iberia with respect to Europe, indicated by: (1) a Late Jurassic-Early Cretaceous left-lateral strike-slip motion; (2) Late Cretaceous convergence; (3) Paleocene quiescence; (4) a short period of right-lateral strike-slip motion; and (5) final Eocene-Oligocene convergence. Based on these results, it is speculated that a collisional episode in the Alpine orogeny at ca. 65 Ma resulted in a dramatic decrease in the relative plate motions and that a slower motion since the Early Miocene promoted extension in the Mediterranean back-arc basins. (C) 2002 Elsevier Science B.V. All rights reserved.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.

Neogene and Quaternary rollback evolution of the Tyrrhenian Sea, the Apennines, and the Sicilian Maghrebides

Reconstruction of the evolution of the Tyrrhenian Sea shows that the major stage of rifting associated with the opening of this basin began at similar to10 Ma. It involved two episodes of back arc extension, which were induced by the rollback of a west dipping subducting slab. The first period of extension (10-6 Ma) was prominent in the northern Tyrrhenian Sea and in the western part of the southern Tyrrhenian Sea. The second period of extension, mainly affected the southern Tyrrhenian Sea, began in the latest Messinian (6-5 Ma) and has been accompanied by subduction rollback at rates of 60-100 km Myr(-1). Slab reconstruction, combined with paleomagnetic and paleogeographic constraints, indicates that in the central Apennines, the latest Messinian (6-5 Ma) arrival of a carbonate platform at the subduction zone impeded subduction and initiated a slab tear and major strike-slip faults. These processes resulted in the formation of a narrow subducting slab beneath the Ionian Sea that has undergone faster subduction rollback and induced extreme rates of back arc extension.

The ambivalent practices of reflexivity

Reflexivity involves turning one's reflexive gaze oil discourse-turning language back on itself to see the Work it does in constituting the world. The subject/researcher sees simultaneously the object of her or his gaze and the means by which the object (which may include oneself as subject) is being constituted. The consciousness of self that reflexive writing sometimes entails may be seen to slip inadvertently into constituting the very (real) self that seems to contradict a focus on the constitutive power of discourse. This article explores this site of slippage and of ambivalence. In a collective biography oil the topic of reflexivity, the authors tell and write stories about reflexivity and in a doubled reflexive arc, examine themselves at work during the workshop. Examining their own memories and reflexive practices, they explore this place of slippage and provide theoretical and practical insight into what is going on in reflexive research and writing.

Satellite observation and climate system model simulation of the St. Lawrence Island polynya

The St. Lawrence Island polynya (SLIP) is a commonly occurring winter phenomenon in the Bering Sea, in which dense saline water produced during new ice formation is thought to flow northward through the Bering Strait to help maintain the Arctic Ocean halocline. Winter darkness and inclement weather conditions have made continuous in situ and remote observation of this polynya difficult. However, imagery acquired from the European Space Agency ERS-1 Synthetic Aperture Radar (SAR) has allowed observation of the St. Lawrence Island polynya using both the imagery and derived ice displacement products. With the development of ARCSyM, a high resolution regional model of the Arctic atmosphere/sea ice system, simulation of the SLIP in a climate model is now possible. Intercomparisons between remotely sensed products and simulations can lead to additional insight into the SLIP formation process. Low resolution SAR, SSM/I and AVHRR infrared imagery for the St. Lawrence Island region are compared with the results of a model simulation for the period of 24-27 February 1992. The imagery illustrates a polynya event (polynya opening). With the northerly winds strong and consistent over several days, the coupled model captures the SLIP event with moderate accuracy. However, the introduction of a stability dependent atmosphere-ice drag coefficient, which allows feedbacks between atmospheric stability, open water, and air-ice drag, produces a more accurate simulation of the SLIP in comparison to satellite imagery. Model experiments show that the polynya event is forced primarily by changes in atmospheric circulation followed by persistent favorable conditions: ocean surface currents are found to have a small but positive impact on the simulation which is enhanced when wind forcing is weak or variable.

Complex distal insertions of the tibialis posterior tendon: detailed anatomic and MR imaging investigation in cadavers

Purpose The purpose of this report was to demonstrate the normal complex insertional anatomy of the tibialis posterior tendon (TPT) in cadavers using magnetic resonance (MR) imaging with anatomic and histologic correlation. Material and methods Ten cadaveric ankles were used according to institutional guidelines. MR T1-weighted spin echo imaging was performed to demonstrate aspects of the complex anatomic distal insertions of the TPT in cadaveric specimens. Findings on MR imaging were correlated with those derived from anatomic and histologic study. Reults Generally, the TPT revealed a low signal in all MR images, except near the level of the medial malleolus, where the TPT suddenly changed direction and ""magic angle"" artifact could be observed. In five out of ten specimens (50%), a type I accessory navicular bone was found in the TPT. In all cases with a type I accessory navicular bone, the TPT had an altered signal in this area. Axial and coronal planes on MR imaging were the best in identifying the distal insertions of the TPT. A normal division of the TPT was observed just proximal to the insertion into the navicular bone in five specimens (100%) occurring at a maximum proximal distance from its attachment to the navicular bone of approximately 1.5 to 2 cm. In the other five specimens, in which a type I accessory navicular bone was present, the TPT directly inserted into the accessory bone and a slip less than 1.5 mm in thickness could be observed attaching to the medial aspect of the navicular bone (100%). Anatomic inspection confirmed the sites of the distal insertions of the components of the TPT. Conclusion MR imaging enabled detailed analysis of the complex distal insertions of the TPT as well as a better understanding of those features of its insertion that can simulate a lesion.

Effect of processing induced particle alignment on the fracture toughness and fracture behavior of multiphase dental ceramics

Objective. To investigate the processing induced particle alignment on fracture behavior of four multiphase dental ceramics (one porcelain, two glass-ceramics and a glass-infiltrated-alumina composite). Methods. Disks (empty set12mm x 1.1 mm-thick) and bars (3 mm x 4 mm x 20 mm) of each material were processed according to manufacturer instructions, machined and polished. Fracture toughness (K(IC)) was determined by the indentation strength method using 3-point bending and biaxial flexure fixtures for the fracture of bars and disks, respectively. Microstructural and fractographic analyses were performed with scanning electron microscopy, energy dispersive spectroscopy and X-ray diffraction. Results. The isotropic microstructure of the porcelain and the leucite-based glass-ceramic resulted in similar fracture toughness values regardless of the specimen geometry. On the other hand, materials containing second-phase particles with high aspect ratio (lithium disilicate glass-ceramic and glass-infiltrated-alumina composite) showed lower fracture toughness for disk specimens compared to bars. For the lithium disilicate glass-ceramic disks, it was demonstrated that the occurrence of particle alignment during the heat-pressing procedure resulted in an unfavorable pattern that created weak microstructural paths during the biaxial test. For the glass-infiltrated-alumina composite, the microstructural analysis showed that the large alumina platelets tended to align their large surfaces perpendicularly to the direction of particle deposition during slip casting of green preforms. Significance. The fracture toughness of dental ceramics with anisotropic microstructure should be determined by means of biaxial testing, since it results in lower values. (C) 2009 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

The effect of channel geometry and wall boundary conditions on the formation of extrusion surface instabilities for LLDPE

It is believed that surface instabilities can occur during the extrusion of linear low density polyethylene due to high extensional stresses at the exit of the die. Local crack development can occur at a critical stress level when melt rupture is reached. This high extensional stress results from the rearrangement of the flow at the boundary transition between the wall exit and the free surface. The stress is highest at the extrudate surface and decreases into the bulk of the material. The location of the region where the critical level is reached can determine the amplitude of the extrudate surface distortion, This paper studies the effect of wall slip on the numerically simulated extensional stress level at the die exit and correlates this to the experimentally determined amplitude of the surface instability. The effect of die exit radius and die wall roughness on extrusion surface instabilities is also correlated to the exit stress level in the same way. Whereas full slip may completely suppress the surface instability, a reduction in the exit stress level and instability amplitude is also shown for a rounded die exit and a slight increase in instability is shown to result from a rough die wall. A surface instability map demonstrates how the shear rate for onset of extrusion surface instabilities can be predicted on the basis of melt strength measurements and simulated stress peaks at the exit of the die. (C) 2001 Elsevier Science B.V. All rights reserved.