Computer simulation methods to study interfacial tensions : from the Ising model to colloidal crystals

In condensed matter systems, the interfacial tension plays a central role for a multitude of phenomena. It is the driving force for nucleation processes, determines the shape and structure of crystalline structures and is important for industrial applications. Despite its importance, the interfacial tension is hard to determine in experiments and also in computer simulations. While for liquid-vapor interfacial tensions there exist sophisticated simulation methods to compute the interfacial tension, current methods for solid-liquid interfaces produce unsatisfactory results.rnrnAs a first approach to this topic, the influence of the interfacial tension on nuclei is studied within the three-dimensional Ising model. This model is well suited because despite its simplicity, one can learn much about nucleation of crystalline nuclei. Below the so-called roughening temperature, nuclei in the Ising model are not spherical anymore but become cubic because of the anisotropy of the interfacial tension. This is similar to crystalline nuclei, which are in general not spherical but more like a convex polyhedron with flat facets on the surface. In this context, the problem of distinguishing between the two bulk phases in the vicinity of the diffuse droplet surface is addressed. A new definition is found which correctly determines the volume of a droplet in a given configuration if compared to the volume predicted by simple macroscopic assumptions.rnrnTo compute the interfacial tension of solid-liquid interfaces, a new Monte Carlo method called ensemble switch method'' is presented which allows to compute the interfacial tension of liquid-vapor interfaces as well as solid-liquid interfaces with great accuracy. In the past, the dependence of the interfacial tension on the finite size and shape of the simulation box has often been neglected although there is a nontrivial dependence on the box dimensions. As a consequence, one needs to systematically increase the box size and extrapolate to infinite volume in order to accurately predict the interfacial tension. Therefore, a thorough finite-size scaling analysis is established in this thesis. Logarithmic corrections to the finite-size scaling are motivated and identified, which are of leading order and therefore must not be neglected. The astounding feature of these logarithmic corrections is that they do not depend at all on the model under consideration. Using the ensemble switch method, the validity of a finite-size scaling ansatz containing the aforementioned logarithmic corrections is carefully tested and confirmed. Combining the finite-size scaling theory with the ensemble switch method, the interfacial tension of several model systems, ranging from the Ising model to colloidal systems, is computed with great accuracy.

The analysis of the collapse of a precast r.c. industrial building during the 29th may 2012 emilia romagna earthquake

Analysis of the collapse of a precast r.c. industrial building during the 2012 Emilia earthquake, focus on the failure mechanisms in particular on the flexure-shear interactions. Analysis performed by a time history analysis using a FEM model with the software SAP2000. Finally a reconstruction of the collapse on the basis of the numerical data coming from the strength capacity of the elements failed, using formulation for lightly reinforced columns with high shear and bending moment.

Analytical model for the prediction of wrinkling of rear pressure bulkheads

In Airbus GmbH (Hamburg) has been developed a new design of Rear Pressure Bulkhead (RPB) for the A320-family. The new model has been formed with vacuum forming technology. During this process the wrinkling phenomenon occurs. In this thesis is described an analytical model for prediction of wrinkling based on the energetic method of Timoshenko. Large deflection theory has been used for analyze two cases of study: a simply supported circular thin plate stamped by a spherical punch and a simply supported circular thin plate formed with vacuum forming technique. If the edges are free to displace radially, thin plates will develop radial wrinkles near the edge at a central deflection approximately equal to four plate thicknesses w0/ℎ≈4 if they’re stamped by a spherical punch and w0/ℎ≈3 if they’re formed with vacuum forming technique. Initially, there are four symmetrical wrinkles, but the number increases if the central deflection is increased. By using experimental results, the “Snaptrhough” phenomenon is described.

An Urn Model for Cascading Failures on a Lattice

A cascading failure is a failure in a system of interconnected parts, in which the breakdown of one element can lead to the subsequent collapse of the others. The aim of this paper is to introduce a simple combinatorial model for the study of cascading failures. In particular, having in mind particle systems and Markov random fields, we take into consideration a network of interacting urns displaced over a lattice. Every urn is Pólya-like and its reinforcement matrix is not only a function of time (time contagion) but also of the behavior of the neighboring urns (spatial contagion), and of a random component, which can represent either simple fate or the impact of exogenous factors. In this way a non-trivial dependence structure among the urns is built, and it is used to study default avalanches over the lattice. Thanks to its flexibility and its interesting probabilistic properties, the given construction may be used to model different phenomena characterized by cascading failures such as power grids and financial networks.

Functional and structural characterization of the first prokaryotic member of the L-amino acid transporter (LAT) family: a model for APC transporters

We have identified YkbA from Bacillus subtilis as a novel member of the L-amino acid transporter (LAT) family of amino acid transporters. The protein is approximately 30% identical in amino acid sequence to the light subunits of human heteromeric amino acid transporters. Purified His-tagged YkbA from Escherichia coli membranes reconstituted in proteoliposomes exhibited sodium-independent, obligatory exchange activity for L-serine and L-threonine and also for aromatic amino acids, albeit with less activity. Thus, we propose that YkbA be renamed SteT (Ser/Thr exchanger transporter). Kinetic analysis supports a sequential mechanism of exchange for SteT. Freeze-fracture analysis of purified, functionally active SteT in proteoliposomes, together with blue native polyacrylamide gel electrophoresis and transmission electron microscopy of detergent-solubilized purified SteT, suggest that the transporter exists in a monomeric form. Freeze-fracture analysis showed spherical particles with a diameter of 7.4 nm. Transmission electron microscopy revealed elliptical particles (diameters 6 x 7 nm) with a distinct central depression. To our knowledge, this is the first functional characterization of a prokaryotic member of the LAT family and the first structural data on an APC (amino acids, polyamines, and choline for organocations) transporter. SteT represents an excellent model to study the molecular architecture of the light subunits of heteromeric amino acid transporters and other APC transporters.

Effect of chest compressions only during experimental basic life support on alveolar collapse and recruitment

AIM: The importance of ventilatory support during cardiac arrest and basic life support is controversial. This experimental study used dynamic computed tomography (CT) to assess the effects of chest compressions only during cardiopulmonary resuscitation (CCO-CPR) on alveolar recruitment and haemodynamic parameters in porcine model of ventricular fibrillation. MATERIALS AND METHODS: Twelve anaesthetized pigs (26+/-1kg) were randomly assigned to one of the following groups: (1) intermittent positive pressure ventilation (IPPV) both during basic life support and advanced cardiac life support, or (2) CCO during basic life support and IPPV during advanced cardiac life support. Measurements were acquired at baseline prior to cardiac arrest, during basic life support, during advanced life support, and after return of spontaneous circulation (ROSC), as follows: dynamic CT series, arterial and central venous pressures, blood gases, and regional organ blood flow. The ventilated and atelectatic lung area was quantified from dynamic CT images. Differences between groups were analyzed using the Kruskal-Wallis test, and a p<0.05 was considered statistically significant. RESULTS: IPPV was associated with cyclic alveolar recruitment and de-recruitment. Compared with controls, the CCO-CPR group had a significantly larger mean fractional area of atelectasis (p=0.009), and significantly lower PaO(2) (p=0.002) and mean arterial pressure (p=0.023). The increase in mean atelectatic lung area observed during basic life support in the CCO-CPR group remained clinically relevant throughout the subsequent advanced cardiac life support period and following ROSC, and was associated with prolonged impaired haemodynamics. No inter-group differences in myocardial and cerebral blood flow were observed. CONCLUSION: A lack of ventilation during basic life support is associated with excessive atelectasis, arterial hypoxaemia and compromised CPR haemodynamics. Moreover, these detrimental effects remain evident even after restoration of IPPV.

Visualization of alveolar recruitment in a porcine model of unilateral lung lavage using 3He-MRI

BACKGROUND: In the acute respiratory distress syndrome potentially recruitable lung volume is currently discussed. (3)He-magnetic resonance imaging ((3)He-MRI) offers the possibility to visualize alveolar recruitment directly. METHODS: With the approval of the state animal care committee, unilateral lung damage was induced in seven anesthetized pigs by saline lavage of the right lungs. The left lung served as an intraindividual control (healthy lung). Unilateral lung damage was confirmed by conventional proton MRI and spiral-CT scanning. The total aerated lung volume was determined both at a positive end-expiratory pressure (PEEP) of 0 and 10 mbar from three-dimensionally reconstructed (3)He images, both for healthy and damaged lungs. The fractional increase of aerated volume in damaged and healthy lungs, followed by a PEEP increase from 0 to 10 mbar, was compared. RESULTS: Aerated gas space was visualized with a high spatial resolution in the three-dimensionally reconstructed (3)He-MR images, and aeration defects in the lavaged lung matched the regional distribution of atelectasis in proton MRI. After recruitment and PEEP increase, the aerated volume increased significantly both in healthy lungs from 415 ml [270-445] (median [min-max]) to 481 ml [347-523] and in lavaged lungs from 264 ml [71-424] to 424 ml [129-520]. The fractional increase in lavaged lungs was significantly larger than that in healthy lungs (healthy: 17% [11-38] vs. lavage: 42% [14-90] (P=0.031). CONCLUSION: The (3)He-MRI signal might offer an experimental approach to discriminate atelectatic vs. poor aerated lung areas in a lung damage animal model. Our results confirm the presence of potential recruitable lung volume by either alveolar collapse or alveolar flooding, in accordance with previous reports by computed tomography.

Experimental Validation of Finite Element Analysis of Human Vertebral Collapse Under Large Compressive Strains

Osteoporosis-related vertebral fractures represent a major health problem in elderly populations. Such fractures can often only be diagnosed after a substantial deformation history of the vertebral body. Therefore, it remains a challenge for clinicians to distinguish between stable and progressive potentially harmful fractures. Accordingly, novel criteria for selection of the appropriate conservative or surgical treatment are urgently needed. Computer tomography-based finite element analysis is an increasingly accepted method to predict the quasi-static vertebral strength and to follow up this small strain property longitudinally in time. A recent development in constitutive modeling allows us to simulate strain localization and densification in trabecular bone under large compressive strains without mesh dependence. The aim of this work was to validate this recently developed constitutive model of trabecular bone for the prediction of strain localization and densification in the human vertebral body subjected to large compressive deformation. A custom-made stepwise loading device mounted in a high resolution peripheral computer tomography system was used to describe the progressive collapse of 13 human vertebrae under axial compression. Continuum finite element analyses of the 13 compression tests were realized and the zones of high volumetric strain were compared with the experiments. A fair qualitative correspondence of the strain localization zone between the experiment and finite element analysis was achieved in 9 out of 13 tests and significant correlations of the volumetric strains were obtained throughout the range of applied axial compression. Interestingly, the stepwise propagating localization zones in trabecular bone converged to the buckling locations in the cortical shell. While the adopted continuum finite element approach still suffers from several limitations, these encouraging preliminary results towardsthe prediction of extended vertebral collapse may help in assessing fracture stability in future work.

Alveolar derecruitment and collapse induration as crucial mechanisms in lung injury and fibrosis.

Idiopathic pulmonary fibrosis (IPF) and bleomycin-induced pulmonary fibrosis are associated with surfactant system dysfunction, alveolar collapse (derecruitment), and collapse induration (irreversible collapse). These events play undefined roles in the loss of lung function. The purpose of this study was to quantify how surfactant inactivation, alveolar collapse, and collapse induration lead to degradation of lung function. Design-based stereology and invasive pulmonary function tests were performed 1, 3, 7, and 14 days after intratracheal bleomycin-instillation in rats. The number and size of open alveoli was correlated to mechanical properties. Active surfactant subtypes declined by Day 1, associated with a progressive alveolar derecruitment and a decrease in compliance. Alveolar epithelial damage was more pronounced in closed alveoli compared with ventilated alveoli. Collapse induration occurred on Day 7 and Day 14 as indicated by collapsed alveoli overgrown by a hyperplastic alveolar epithelium. This pathophysiology was also observed for the first time in human IPF lung explants. Before the onset of collapse induration, distal airspaces were easily recruited, and lung elastance could be kept low after recruitment by positive end-expiratory pressure (PEEP). At later time points, the recruitable fraction of the lung was reduced by collapse induration, causing elastance to be elevated at high levels of PEEP. Surfactant inactivation leading to alveolar collapse and subsequent collapse induration might be the primary pathway for the loss of alveoli in this animal model. Loss of alveoli is highly correlated with the degradation of lung function. Our ultrastructural observations suggest that collapse induration is important in human IPF.

Shape model, reference system definition, and cartographic mapping standards for comet 67P/Churyumov-Gerasimenko Stereo-photogrammetric analysis of Rosetta/OSIRIS image data

We analyzed more than 200 OSIRIS NAC images with a pixel scale of 0.9-2.4 m/pixel of comet 67P/Churyumov-Gerasimenko (67P) that have been acquired from onboard the Rosetta spacecraft in August and September 2014 using stereo-photogrammetric methods (SPG). We derived improved spacecraft position and pointing data for the OSIRIS images and a high-resolution shape model that consists of about 16 million facets (2 m horizontal sampling) and a typical vertical accuracy at the decimeter scale. From this model, we derive a volume for the northern hemisphere of 9.35 km(3) +/- 0.1 km(3). With the assumption of a homogeneous density distribution and taking into account the current uncertainty of the position of the comet's center-of-mass, we extrapolated this value to an overall volume of 18.7 km(3) +/- 1.2 km(3), and, with a current best estimate of 1.0 X 10(13) kg for the mass, we derive a bulk density of 535 kg/m(3) +/- 35 kg/m(3). Furthermore, we used SPG methods to analyze the rotational elements of 67P. The rotational period for August and September 2014 was determined to be 12.4041 +/- 0.0004 h. For the orientation of the rotational axis (z-axis of the body-fixed reference frame) we derived a precession model with a half-cone angle of 0.14 degrees, a cone center position at 69.54 degrees/64.11 degrees (RA/Dec J2000 equatorial coordinates), and a precession period of 10.7 days. For the definition of zero longitude (x-axis orientation), we finally selected the boulder-like Cheops feature on the big lobe of 67P and fixed its spherical coordinates to 142.35 degrees right-hand-rule eastern longitude and -0.28 degrees latitude. This completes the definition of the new Cheops reference frame for 67P. Finally, we defined cartographic mapping standards for common use and combined analyses of scientific results that have been obtained not only within the OSIRIS team, but also within other groups of the Rosetta mission.

Gravity potential spherical harmonic series

Time variable gravity fields, reflecting variations of mass distribution in the system Earth is one of the key parameters to understand the changing Earth. Mass variations are caused either by redistribution of mass in, on or above the Earth's surface or by geophysical processes in the Earth's interior. The first set of observations of monthly variations of the Earth gravity field was provided by the US/German GRACE satellite mission beginning in 2002. This mission is still providing valuable information to the science community. However, as GRACE has outlived its expected lifetime, the geoscience community is currently seeking successor missions in order to maintain the long time series of climate change that was begun by GRACE. Several studies on science requirements and technical feasibility have been conducted in the recent years. These studies required a realistic model of the time variable gravity field in order to perform simulation studies on sensitivity of satellites and their instrumentation. This was the primary reason for the European Space Agency (ESA) to initiate a study on ''Monitoring and Modelling individual Sources of Mass Distribution and Transport in the Earth System by Means of Satellites''. The goal of this interdisciplinary study was to create as realistic as possible simulated time variable gravity fields based on coupled geophysical models, which could be used in the simulation processes in a controlled environment. For this purpose global atmosphere, ocean, continental hydrology and ice models were used. The coupling was performed by using consistent forcing throughout the models and by including water flow between the different domains of the Earth system. In addition gravity field changes due to solid Earth processes like continuous glacial isostatic adjustment (GIA) and a sudden earthquake with co-seismic and post-seismic signals were modelled. All individual model results were combined and converted to gravity field spherical harmonic series, which is the quantity commonly used to describe the Earth's global gravity field. The result of this study is a twelve-year time-series of 6-hourly time variable gravity field spherical harmonics up to degree and order 180 corresponding to a global spatial resolution of 1 degree in latitude and longitude. In this paper, we outline the input data sets and the process of combining these data sets into a coherent model of temporal gravity field changes. The resulting time series was used in some follow-on studies and is available to anybody interested.

Late glacial and deglacial age model and biomarker analyses of sediment core MSM05/5_712-2 from the western continental slope of Svalbard

The transition from last glacial to deglacial and subsequently to modern interglacial climate conditions was accompanied by abrupt shifts in the palaeoceanographic setting in the subpolar North Atlantic. Knowledge about the role that sea ice coverage played during these rapid climate reversals is limited since most marine sediment cores from the higher latitudes provide only a coarse temporal resolution and often poorly preserved microfossils. Here we present a highly resolved reconstruction of sea ice conditions that characterised the eastern Fram Strait - a key area for water mass exchange between the Arctic Ocean and the North Atlantic - for the past 30 ka BP. This reconstruction is based on the distribution of the sea ice biomarker IP25 and phytoplankton derived biomarkers in a sediment core from the continental slope of western Svalbard. During the late glacial (30 ka to 19 ka BP), recurrent advances and retreats of sea ice characterised the study area and point to a hitherto less considered oceanic (and/or atmospheric) variability. A long-lasting perennial sea ice coverage in eastern Fram Strait persisted only at the very end of the Last Glacial Maximum (i.e. from 19.2 to 17.6 ka BP) and was abruptly reduced at the onset of Heinrich Event 1 - coincident with or possibly even inducing the collapse of the Atlantic Meridional Overturning Circulation (AMOC). Further maximum sea ice conditions prevailed during the Younger Dryas cooling event and support the assumption of an AMOC reduction due to increased formation and export of Arctic sea ice through Fram Strait. A significant retreat of sea ice and sea surface warming are observed for the Early Holocene.