Effect of mesh distortion on the accuracy of high order vorticity-velocity CFD approaches

The numerical solution of the incompressible Navier-Stokes equations offers an alternative to experimental analysis of fluid-structure interaction (FSI). We would save a lot of time and effort and help cut back on costs, if we are able to accurately model systems by these numerical solutions. These advantages are even more obvious when considering huge structures like bridges, high rise buildings or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the Kinematic Laplacian Equation (KLE) to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ordinary differential equations (ODE) time integration schemes, allowing us to tackle each problem as a separate module. The current algortihm for the KLE uses an unstructured quadrilateral mesh, formed by dividing each triangle of an unstructured triangular mesh into three quadrilaterals for spatial discretization. This research deals with determining a suitable measure of mesh quality based on the physics of the problems being tackled. This is followed by exploring methods to improve the quality of quadrilateral elements obtained from the triangles and thereby improving the overall mesh quality. A series of numerical experiments were designed and conducted for this purpose and the results obtained were tested on different geometries with varying degrees of mesh density.

Evaluating stress patterns for Plio-Quaternary brittle deformation along the Calama-Olacapato-El Toro fault zone, Central Andes

Understanding the geometry and kinematics of the major structures of an orogen is important to elucidate its style of deformation, as well as its tectonic evolution. We describe the temporal and spatial changes in the state of stress of the trans-orogen area of the Calama-Olacapato-El Toro (COT) Fault Zone in the Central Andes, at about 24°S within the northern portion of the Puna Plateau between the Argentina-Chile border. The importance of the COT derives principally from the Quaternary-Holocene activity recognized on some segments, which may shed new light on its possible control on Quaternary volcanism and on the seismic hazard evaluation of the area. Field geological surveys along with kinematic analysis and numerical inversion of ∼140 new fault-slip measurements have revealed that this portion of the COT zone, previously considered a continuous, long-lived lineament, in reality has been subjected to three different kinematic regimes: 1) a Miocene transpressional phase with the maximum principal stress (σ1) chiefly trending NNE-SSW; 2) an extensional phase that started by 9 Ma, with a horizontal NW-SE-striking minimum principal stress (σ3) – permutations between σ2 and σ3 axes have been recognized at two sites – and 3) a left-lateral strike-slip phase with a horizontal ∼E-W &sigma1 and ∼N-S σ3 dating to the Late Pliocene-Quaternary. Spatially, in the Quaternary, the left-lateral component decreases toward the westernmost tip of the COT, where it transitions to extension; this produced to a N-S horst and graben structure. Hence, even if transcurrence is still active in the eastern portion of the COT, as focal mechanisms of crustal earthquakes indicate, our study demonstrates that extension is becoming the predominant structural style of deformation, at least in the western region. These major temporal and spatial changes in the tectonic regimes are attributed in part to changes in the magnitude of the boundary forces due to subduction processes. The overall orogen-perpendicular extension might be the result of vertical stress larger than both the horizontal stresses induced by gravitational effect of a thickened crust.

Swimming as a model of task-specific locomotor retraining after spinal cord injury in the rat

BACKGROUND: The authors have shown that rats can be retrained to swim after a moderately severe thoracic spinal cord contusion. They also found that improvements in body position and hindlimb activity occurred rapidly over the first 2 weeks of training, reaching a plateau by week 4. Overground walking was not influenced by swim training, suggesting that swimming may be a task-specific model of locomotor retraining. OBJECTIVE: To provide a quantitative description of hindlimb movements of uninjured adult rats during swimming, and then after injury and retraining. METHODS: The authors used a novel and streamlined kinematic assessment of swimming in which each limb is described in 2 dimensions, as 3 segments and 2 angles. RESULTS: The kinematics of uninjured rats do not change over 4 weeks of daily swimming, suggesting that acclimatization does not involve refinements in hindlimb movement. After spinal cord injury, retraining involved increases in hindlimb excursion and improved limb position, but the velocity of the movements remained slow. CONCLUSION: These data suggest that the activity pattern of swimming is hardwired in the rat spinal cord. After spinal cord injury, repetition is sufficient to bring about significant improvements in the pattern of hindlimb movement but does not improve the forces generated, leaving the animals with persistent deficits. These data support the concept that force (load) and pattern generation (recruitment) are independent and may have to be managed together with respect to postinjury rehabilitation.

Der Einfluss kinematischer Parameter auf den Energiebedarf eines Regalbediengerätes

Bisher wurde noch nicht erforscht wie sich die Variation kinematischer Parameter konkret auf den Energiebedarf von Regalbediengeräten auswirkt. Dieser Beitrag untersucht daher die energetischen Folgen bei unterschiedlichen Geschwindigkeiten und Beschleunigungen der Antriebe eines Regalbediengerätes. Veränderte Geschwindigkeiten und Beschleunigungen beeinflussen jedoch auch die Fahr- und Hubzeiten des Regalbediengerätes. Daher wird ebenfalls eine Berechnungsmethode vorgestellt, mit der der optimale Startpunkt des Hubwerkes und die optimale Fahrgeschwindigkeit auf Basis der eingestellten kinematischen Parameter berechnet werden können.

Simulationsgestützte Untersuchung der Spurtreue von Routenzügen

Für den innerbetrieblichen Transport und besonders für die Materialversorgung in der Produktion werden zunehmend Routenzüge eingesetzt, die aus einem Schleppfahrzeug und bis zu fünf Anhängern bestehen. Um gute Manövrierbarkeit und einen geringen Platzbedarf zu gewährleisten, sollten Routenzuganhänger möglichst spurtreu sein. In diesem Beitrag wird die Spurtreue exemplarisch für zwei am Markt verfügbare Anhänger für Routenzüge mit zwei ungelenkten und vier gelenkten Rädern untersucht. Hierfür wird zunächst ein Gütekriterium definiert, das die maximale Abweichung von der Spurtreue quantitativ erfasst. Außerdem werden Testszenarien vorgeschlagen, um die verschiedenen Fahrwerks- und Lenkkonzepte vergleichen zu können. Mit Hilfe eines entwickelten analytischen Modells werden die Abweichungen in der Spurtreue bei stationärer Kreisfahrt für die zwei gewählten Konzepte berechnet und dargestellt. Zusätzlich wird eine Mehrkörper-Simulation, die eine tiefere physikalische Modellierung und die Untersuchung komplexerer Fahrmanöver erlaubt, durchgeführt und mit den analytischen Ergebnissen verglichen. Es kann gezeigt werden, dass neben der Lenkkinematik weitere Parameter wie das Schräglaufverhalten der Reifen Einfluss auf die Spurtreue haben.

GPS-only gravity field determination from GOCE data

The Gravity field and steady-state Ocean Circulation Explorer (GOCE) is now in orbit for more than four years. This is longer than the originally planned lifetime of the satellite and after three years on the same altitude the satellite has been lowered to 235 km in several steps. In the frame of the GOCE High-level Processing Facility the Astronomical Institute of the University of Bern (AIUB) is responsible for the determination of the official Precise Science Orbit (PSO) product. Kinematic GOCE orbits are part of this product and are used by several institutions in- and outside the HPF for determining the low degrees of the Earth’s gravity field. AIUB GOCE GPS-only gravity field solutions using the Celestial Mechanics Approach and covering the Release 4 period as well as a more recent time interval at the lower orbit altitude are shown and discussed. Special attention is paid to the impact of systematic deficiencies in the kinematic orbits on the resulting gravity fields, e.g., related to the geomagnetic equator, and on possibilities to get rid of them.

Morphology–elasticity relationships using decreasing fabric information of human trabecular bone from three major anatomical locations

With improving clinical CT scanning technology, the accuracy of CT-based finite element (FE) models of the human skeleton may be ameliorated by an enhanced description of apparent level bone mechanical properties. Micro-finite element (μFE) modeling can be used to study the apparent elastic behavior of human cancellous bone. In this study, samples from the femur, radius and vertebral body were investigated to evaluate the predictive power of morphology–elasticity relationships and to compare them across different anatomical regions. μFE models of 701 trabecular bone cubes with a side length of 5.3 mm were analyzed using kinematic boundary conditions. Based on the FE results, four morphology–elasticity models using bone volume fraction as well as full, limited or no fabric information were calibrated for each anatomical region. The 5 parameter Zysset–Curnier model using full fabric information showed excellent predictive power with coefficients of determination ( r2adj ) of 0.98, 0.95 and 0.94 of the femur, radius and vertebra data, respectively, with mean total norm errors between 14 and 20%. A constant orthotropy model and a constant transverse isotropy model, where the elastic anisotropy is defined by the model parameters, yielded coefficients of determination between 0.90 and 0.98 with total norm errors between 16 and 25%. Neglecting fabric information and using an isotropic model led to r2adj between 0.73 and 0.92 with total norm errors between 38 and 49%. A comparison of the model regressions revealed minor but significant (p<0.01) differences for the fabric–elasticity model parameters calibrated for the different anatomical regions. The proposed models and identified parameters can be used in future studies to compute the apparent elastic properties of human cancellous bone for homogenized FE models.

Estimating the patient's contribution during robot-assisted therapy

Robot-assisted therapy has become increasingly common in neurorehabilitation. Sophisticated controllers have been developed for robots to assist and cooperate with the patient. It is difficult for the patient to judge to what extent the robot contributes to the execution of a movement. Therefore, methods to comprehensively quantify the patient's contribution and provide feedback are of key importance. We developed a method comprehensively to estimate the patient's contribution by combining kinematic measures and the motor assistance applied. Inverse dynamic models of the robot and the passive human arm calculate the required torques to move the robot and the arm and build, together with the recorded motor torque, a metric (in percentage) that represents the patient's contribution to the movement. To evaluate the developed metric, 12 nondisabled subjects and 7 patients with neurological problems simulated instructed movement contributions. The results are compared with a common performance metric. The estimation shows very satisfying results for both groups, even though the arm model used was strongly simplified. Displaying this metric to patients during therapy can potentially motivate them to actively participate in the training.

Quantitative 3D strain analysis in analogue experiments: Integration of X-ray computed tomography and digital volume correlation techniques

The combination of scaled analogue experiments, material mechanics, X-ray computed tomography (XRCT) and Digital Volume Correlation techniques (DVC) is a powerful new tool not only to examine the 3 dimensional structure and kinematic evolution of complex deformation structures in scaled analogue experiments, but also to fully quantify their spatial strain distribution and complete strain history. Digital image correlation (DIC) is an important advance in quantitative physical modelling and helps to understand non-linear deformation processes. Optical non-intrusive (DIC) techniques enable the quantification of localised and distributed deformation in analogue experiments based either on images taken through transparent sidewalls (2D DIC) or on surface views (3D DIC). X-ray computed tomography (XRCT) analysis permits the non-destructive visualisation of the internal structure and kinematic evolution of scaled analogue experiments simulating tectonic evolution of complex geological structures. The combination of XRCT sectional image data of analogue experiments with 2D DIC only allows quantification of 2D displacement and strain components in section direction. This completely omits the potential of CT experiments for full 3D strain analysis of complex, non-cylindrical deformation structures. In this study, we apply digital volume correlation (DVC) techniques on XRCT scan data of “solid” analogue experiments to fully quantify the internal displacement and strain in 3 dimensions over time. Our first results indicate that the application of DVC techniques on XRCT volume data can successfully be used to quantify the 3D spatial and temporal strain patterns inside analogue experiments. We demonstrate the potential of combining DVC techniques and XRCT volume imaging for 3D strain analysis of a contractional experiment simulating the development of a non-cylindrical pop-up structure. Furthermore, we discuss various options for optimisation of granular materials, pattern generation, and data acquisition for increased resolution and accuracy of the strain results. Three-dimensional strain analysis of analogue models is of particular interest for geological and seismic interpretations of complex, non-cylindrical geological structures. The volume strain data enable the analysis of the large-scale and small-scale strain history of geological structures.

Jet energy measurement with the ATLAS detector in proton-proton collisions at sqrt(s) = 7 TeV

The jet energy scale (JES) and its systematic uncertainty are determined for jets measured with the ATLAS detector at the LHC in proton-proton collision data at a centre-of-mass energy of sqrt(s) = 7 TeV corresponding to an integrated luminosity of 38 inverse pb. Jets are reconstructed with the anti-kt algorithm with distance parameters R=0.4 or R=0.6. Jet energy and angle corrections are determined from Monte Carlo simulations to calibrate jets with transverse momenta pt > 20 GeV and pseudorapidities eta<4.5. The JES systematic uncertainty is estimated using the single isolated hadron response measured in situ and in test-beams. The JES uncertainty is less than 2.5% in the central calorimeter region (eta<0.8) for jets with 60 < pt < 800 GeV, and is maximally 14% for pt < 30 GeV in the most forward region 3.2 50 GeV after a dedicated correction for this effect. The JES is validated for jet transverse momenta up to 1 TeV to the level of a few percent using several in situ techniques by comparing a well-known reference such as the recoiling photon pt, the sum of the transverse momenta of tracks associated to the jet, or a system of low-pt jets recoiling against a high-pt jet. More sophisticated jet calibration schemes are presented based on calorimeter cell energy density weighting or hadronic properties of jets, providing an improved jet energy resolution and a reduced flavour dependence of the jet response. The JES systematic uncertainty determined from a combination of in situ techniques are consistent with the one derived from single hadron response measurements over a wide kinematic range. The nominal corrections and uncertainties are derived for isolated jets in an inclusive sample of high-pt jets.

Search for new phenomena in events with three charged leptons at a center-of-mass energy of 7 TeV with the ATLAS detector

A generic search for anomalous production of events with at least three charged leptons is presented. The search uses a pp-collision data sample at a center-of-mass energy of root s = 7 TeV corresponding to 4.6 fb(-1) of integrated luminosity collected in 2011 by the ATLAS detector at the CERN Large Hadron Collider. Events are required to contain at least two electrons or muons, while the third lepton may either be an additional electron or muon, or a hadronically decaying tau lepton. Events are categorized by the presence or absence of a reconstructed tau-lepton or Z-boson candidate decaying to leptons. No significant excess above backgrounds expected from Standard Model processes is observed. Results are presented as upper limits on event yields from non-Standard-Model processes producing at least three prompt, isolated leptons, given as functions of lower bounds on several kinematic variables. Fiducial efficiencies for model testing are also provided. The use of the results is illustrated by setting upper limits on the production of doubly charged Higgs bosons decaying to same-sign lepton pairs.

Measurement of ℤℤ production in pp collisions at sqrts=7 TeV and limits on anomalous ℤℤℤ ℤℤgamma couplings with the ATLAS detector

A measurement of the ZZ production cross section in proton-proton collisions at root s = 7 TeV using data recorded by the ATLAS experiment at the Large Hadron Collider is presented. In a data sample corresponding to an integrated luminosity of 4.6 fb(-1) collected in 2011, events are selected that are consistent either with two Z bosons decaying to electrons or muons or with one Z boson decaying to electrons or muons and a second Z boson decaying to neutrinos. The ZZ((*)) -> l(+)l(-)l'(+)l'(-) and ZZ -> l(+)l(-) nu(nu) over bar cross sections are measured in restricted phase-space regions. These results are then used to derive the total cross section for ZZ events produced with both Z bosons in the mass range 66 to 116 GeV, sigma(tot)(ZZ) = 6.7 +/- 0.7 (stat.) (+0.4)(-0.3) (syst.) +/- 0.3 (lumi.) pb, which is consistent with the Standard Model prediction of 5.89(-0.18)(+0.22) pb calculated at next-to-leading order in QCD. The normalized differential cross sections in bins of various kinematic variables are presented. Finally, the differential event yield as a function of the transverse momentum of the leading Z boson is used to set limits on anomalous neutral triple gauge boson couplings in ZZ production.

Measurement of the flavour composition of dijet events in pp collisions at sqrts=7 TeV with the AT detector

This paper describes a measurement of the flavour composition of dijet events produced in pp collisions at root s = 7 TeV using the ATLAS detector. The measurement uses the full 2010 data sample, corresponding to an integrated luminosity of 39 pb(-1). Six possible combinations of light, charm and bottom jets are identified in the dijet events, where the jet flavour is defined by the presence of bottom, charm or solely light flavour hadrons in the jet. Kinematic variables, based on the properties of displaced decay vertices and optimised for jet flavour identification, are used in a multidimensional template fit to measure the fractions of these dijet flavour states as functions of the leading jet transverse momentum in the range 40 GeV to 500 GeV and jet rapidity vertical bar y vertical bar < 2.1. The fit results agree with the predictions of leading-and next-to-leading-order calculations, with the exception of the dijet fraction composed of bottom and light flavour jets, which is underestimated by all models at large transverse jet momenta. The ability to identify jets containing two b-hadrons, originating from e. g. gluon splitting, is demonstrated. The difference between bottom jet production rates in leading and subleading jets is consistent with the next-to-leading-order predictions.

Measurement of the cross-section for W boson production in association with b-jets in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

This paper reports a measurement of the W+b-jets (W+b+X and W+b (b) over bar +X) production cross-section in proton-proton collisions at a centre-of-mass energy of 7 TeV at the LHC. These results are based on data corresponding to an integrated luminosity of 4.6 fb(-1), collected with the ATLAS detector. Cross-sections are presented as a function of jet multiplicity and of the transverse momentum of the leading b-jet for both the muon and electron decay modes of the W boson. The W+b-jets cross-section, corrected for all known detector effects, is quoted in a limited kinematic range. Combining the muon and electron channels, the fiducial cross-section for W+b-jets is measured to be 7.1 +/- 0.5 (stat) +/- 1.4 (syst) pb, consistent with the next-to-leading order QCD prediction, corrected for non-perturbative and double-parton interactions (DPI) contributions, of 4.70 +/- 0.09 (stat) (+0.60)(-0.49) (scale) +/- 0.06 (PDF) +/- 0.16 (non-pert) (+0.52)(-0.38) (DPI) pb.

Kinetics and kinematics of the tölt: effects of rider interaction and shoeing manipulations

Introduction . Compared to most equine horse breeds which are able to walk, trot and canter /gallop, the gait repertoire of the Icelandic horses additionally includes the lateral gait tölt and frequently also the pace. With respect to the tölt gait, special shoeing, saddling and riding techniques have been developed for Icelandic horses in order to enhance its expressiveness and regularity. Toes are left unnaturally long and heavy shoes and paddings, as well as weighted boots are used to enforce the individual gait predisposition. For the same reason, the rider is placed more caudally to the horse's centre of mass as compared to other riding techniques. The biomechanical impact of these methods on the health of the locomotor system has so far never been subject of systematic research. Objectives . The aims of the presented study are (1) to describe the kinetic and kinematic characteristics of the tölt performed on a treadmill, (2) to understand the mechanical consequences of shoeing manipulation (long hooves, weighted boots) on the loading and protraction movement of the limbs, as well as (3) to study the pressure distribution and effects on the gait pattern of 3 different saddle types used for riding Icelandic horses. Materials and methods . Gait analysis was carried out in 13 Icelandic horses at walk and at slow and medium tölting and trotting speeds on a high-speed treadmill instrumented for measuring vertical ground reaction forces as well as temporal and spatial gait variables. Kinematic data of horse, rider and saddle were measured simultaneously. Gait analysis was first carried out with high, long hooves (SH) without and in combination with weighted boots (ad aim (2)). Afterwards, horses were re-shod according to current horseshoeing standards (SN) and gait analysis was repeated (ad aims (1) and (2)). In a second trial, horses were additionally equipped with a pressure sensitive saddle mat and were ridden with a dressage-like saddle (SDres), an Icelandic saddle (Slcel) and a saddle cushion (SCush) in the standard saddle position (ad aim 3). Results and conclusions . Compared to trot at the same speed, tölting horses had a higher stride rate and lower stride impulses. At the tölt loading of the forelimbs was increased in form of higher peak vertical forces (Fzpeak) due to shorter relative stance durations (StDrel). Conversely, in the hindlimbs, longer StDrel resulted in lower Fzpeak. Despite the higher head-neck position at tölt, there was no measurable shift in weight to the hindlimbs. Footfall rhythm was in most horses laterally coupled at the tölt and frequently had a slight fourbeat and a very short suspension phase at trot; underlining the fact that performance of correct gaits in Icelandic horses needs special training. Gait performance as it is currently judged in competition could be improved using a shoeing with SH, resulting in a 21 ± 5 mm longer dorsal hoof wall, but also a weight gain of 273 ± 50 g at the distal limb due to heavier shoeing material. Compared to SN, SH led to a lower stride rate, a longer stride length and a higher, but not wider, forelimb protraction arc, which were also positively associated with speed. At the tölt, the footfall rhythm showed less tendency to lateral couplets and at the trot, the suspension phase was longer. However, on the long term, SH may have negative implications for the health of the palmar structures of the distal foot by increased limb impulses, higher torques at breakover (up to 20%); as well as peak vertical forces at faster speeds. Compared to the shoeing style, the saddle type had less influence on limb forces or movements. The slight weight shift to the rear with SCush and Slcel may be explained by the more caudal position of the rider relative to the horse's back. With SCush, pressure was highest under the cranial part of the saddle, whereas the saddles with trees had more pressure under the caudal area.