The spatial distribution of water in the inner coma of Comet 9P/Tempel 1: Comparison between models and observations

The near nucleus coma of Comet 9P/Tempel 1 has been simulated with the 3D Direct Simulation Monte Carlo (DSMC) code PDSC++ (Su, C.-C. [2013]. Parallel Direct Simulation Monte Carlo (DSMC) Methods for Modeling Rarefied Gas Dynamics. PhD Thesis, National Chiao Tung University, Taiwan) and the derived column densities have been compared to observations of the water vapour distribution found by using infrared imaging spectrometer on the Deep Impact spacecraft (Feaga, L.M., A’Hearn, M.F., Sunshine, J.M., Groussin, O., Farnham, T.L. [2007]. Icarus 191(2), 134–145. http://dx.doi.org/10.1016/j.icarus.2007.04.038). Modelled total production rates are also compared to various observations made at the time of the Deep Impact encounter. Three different models were tested. For all models, the shape model constructed from the Deep Impact observations by Thomas et al. (Thomas, P.C., Veverka, J., Belton, M.J.S., Hidy, A., A’Hearn, M.F., Farnham, T.L., et al. [2007]. Icarus, 187(1), 4–15. http://dx.doi.org/10.1016/j.icarus.2006.12.013) was used. Outgassing depending only on the cosine of the solar insolation angle on each shape model facet is shown to provide an unsatisfactory model. Models constructed on the basis of active areas suggested by Kossacki and Szutowicz (Kossacki, K., Szutowicz, S. [2008]. Icarus, 195(2), 705–724. http://dx.doi.org/10.1016/j.icarus.2007.12.014) are shown to be superior. The Kossacki and Szutowicz model, however, also shows deficits which we have sought to improve upon. For the best model we investigate the properties of the outflow.

Epithelial transformations in the establishment of the blood-gas barrier in the developing chick embryo lung

The tall epithelium of the developing chick embryo lung is converted to a squamous one, which participates in formation of the thin blood-gas barrier. We show that this conversion occurred through processes resembling exocrine secretion. Initially, cells formed intraluminal protrusions (aposomes), and then transcellular double membranes were established. Gaps between the membranes opened, thus, severing the aposome from the cell. Alternatively, aposomes were squeezed out by adjacent cells or were spontaneously constricted and extruded. As a third mechanism, formation and fusion of severed vesicles or vacuoles below the aposome and their fusion with the apicolateral plasma membrane resulted in severing of the aposome. The atria started to form by progressive epithelial attenuation and subsequent invasion of the surrounding mesenchyme at regions delineated by subepithelial alpha-smooth muscle actin-positive cells. Further epithelial attenuation was achieved by vacuolation; rupture of such vacuoles with resultant numerous microfolds and microvilli, which were abscised to accomplish a smooth squamous epithelium just before hatching.

Influence of inspired oxygen on glucose-lactate dynamics after subdural hematoma in the rat

The mechanisms causing brain damage after acute subdural hematoma (SDH) are poorly understood. A decrease in cerebral blood flow develops immediately after the hematoma forms, thus reducing cerebral oxygenation. This in turn may activate mitochondrial failure and tissue damage leading to ionic imbalance and possibly to cellular breakdown. The purpose of this study was to test whether a simple therapeutic measure, namely increased fraction of inspired oxygen (FiO2 100), and hence increased arterial and brain tissue oxygen tension, can influence brain glucose and lactate dynamics acutely after subdural hematoma in the rat. Twenty-five male Sprague-Dawley anesthetized rats were studied before, during and after induction of the SDH in two separate groups. The Oxygen group (n = 10) was ventilated with 100% oxygen immediately after induction of the SDH. The Air group (n = 10) was ventilated during the entire study with 21% oxygen. Brain microdialysate samples were analyzed for glucose and lactate. All rats were monitored with femoral arterial blood pressure catheters, arterial blood gas analysis, arterial glucose, lactate and end tidal CO2 (EtCO2). Five male Sprague-Dawley rats were sham operated to measure the effect of oxygen challenge on glucose-lactate dynamics without injury. Arterial oxygen tension in the Oxygen group was 371 +/- 30 mmHg and was associated with significantly greater increase in dialysate lactate in the first 30 min after induction of SDH. Dialysate glucose initially dropped in both groups, after SDH, but then reverted significantly faster to values above baseline in the Oxygen group. Changes in ventilatory parameters had no significant effect on dialysate glucose and lactate parameters in the sham group. Extracellular dialysate lactate and glucose are influenced by administration of 100% O2 after SDH. Dialysate glucose normalizes significantly quicker upon 100% oxygen ventilation. We hypothesize that increased neural tissue oxygen tension, in presence of reduced regional CBF, and possibly compromised mitochondrial function, after acute SDH results in upregulation of rate-limiting enzyme systems responsible for both glycolytic and aerobic metabolism. Similar changes have been seen in severe human head injury, and suggest that a simple therapeutic measure, such as early ventilation with 100% O2, may improve cerebral energy metabolism, early after SDH. Further studies to measure the generation of adenosine triphosphate (ATP) are needed to validate the hypothesis.

A universal description of ultraslow glass dynamics

The dynamics of glass is of importance in materials science but its nature has not yet been fully understood. Here we report that a verification of the temperature dependencies of the primary relaxation time or viscosity in the ultraslowing/ultraviscous domain of glass-forming systems can be carried out via the analysis of the inverse of the Dyre-Olsen temperature index. The subsequent analysis of experimental data indicates the possibility of the self-consistent description of glass-forming low-molecular-weight liquids, polymers, liquid crystals, orientationally disordered crystals and Ising spin-glass-like systems, as well as the prevalence of equations associated with the 'finite temperature divergence'. All these lead to a new formula for the configurational entropy in glass-forming systems. Furthermore, a link to the dominated local symmetry for a given glass former is identified here. Results obtained show a new relationship between the glass transition and critical phenomena.

Glacial–interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ15N measurements

Correct estimation of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice core studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: outputs of a firn densification model, and measurements of δ15N of N2 in air trapped in ice core, assuming that δ15N is only affected by gravitational fractionation in the firn column. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi-coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available ice core air-δ15N measurements from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rates and temperature conditions. Our δ15N profiles reveal a heterogeneous response of the firn structure to glacial–interglacial climatic changes. While firn densification simulations correctly predict TALDICE δ15N variations, they systematically fail to capture the large millennial-scale δ15N variations measured at BI and the δ15N glacial levels measured at JRI and EDML – a mismatch previously reported for central East Antarctic ice cores. New constraints of the EDML gas–ice depth offset during the Laschamp event (~41 ka) and the last deglaciation do not favour the hypothesis of a large convective zone within the firn as the explanation of the glacial firn model–δ15N data mismatch for this site. While we could not conduct an in-depth study of the influence of impurities in snow for firnification from the existing datasets, our detailed comparison between the δ15N profiles and firn model simulations under different temperature and accumulation rate scenarios suggests that the role of accumulation rate may have been underestimated in the current description of firnification models.

Transient simulations of the carbon and nitrogen dynamics in northern peatlands: from the Last Glacial Maximum to the 21st century

The development of northern high-latitude peatlands played an important role in the carbon (C) balance of the land biosphere since the Last Glacial Maximum (LGM). At present, carbon storage in northern peatlands is substantial and estimated to be 500 ± 100 Pg C (1 Pg C = 1015 g C). Here, we develop and apply a peatland module embedded in a dynamic global vegetation and land surface process model (LPX-Bern 1.0). The peatland module features a dynamic nitrogen cycle, a dynamic C transfer between peatland acrotelm (upper oxic layer) and catotelm (deep anoxic layer), hydrology- and temperature-dependent respiration rates, and peatland specific plant functional types. Nitrogen limitation down-regulates average modern net primary productivity over peatlands by about half. Decadal acrotelm-to-catotelm C fluxes vary between −20 and +50 g C m−2 yr−1 over the Holocene. Key model parameters are calibrated with reconstructed peat accumulation rates from peat-core data. The model reproduces the major features of the peat core data and of the observation-based modern circumpolar soil carbon distribution. Results from a set of simulations for possible evolutions of northern peat development and areal extent show that soil C stocks in modern peatlands increased by 365–550 Pg C since the LGM, of which 175–272 Pg C accumulated between 11 and 5 kyr BP. Furthermore, our simulations suggest a persistent C sequestration rate of 35–50 Pg C per 1000 yr in present-day peatlands under current climate conditions, and that this C sink could either sustain or turn towards a source by 2100 AD depending on climate trajectories as projected for different representative greenhouse gas concentration pathways.

Three-dimensional Vortex Structures on Heated Micro-lattice in a Gas

This paper addresses the microscale heat transfer problem from heated lattice to the gas. A micro-device for enhanced heat transfer is presented and numerically investigated. Thermal creep induces 3-D vortex structures in the vicinity of the lattice. The gas flow is in the slip flow regime (Knudsen number Kn⩽0.1Kn⩽0.1). The simulations are performed using slip flow Navier–Stokes equations with boundary condition formulations proposed by Maxwell and Smoluchowski. In this study the wire thicknesses and distances of the heated lattice are varied. The surface geometrical properties alter significantly heat flux through the surface.

Terahertz macrospin dynamics in insulating ferrimagnets

We investigate numerically the excitation of nonlinear magnetic interactions in a ferrite material by an energetic pump pulse of terahertz (THz) radiation. The calculations are performed by solving the coupled Maxwell and Landau-Lifshitz-Gilbert differential equations. In a time-resolved THz pump/THz probe scheme, it is demonstrated that Faraday rotation of a delayed THz probe pulse can be used to map these interactions. Our study is motivated by the ability of soft x-ray free electron lasers to perform time-resolved imaging of the magnetization process at the submicrometer and subpicosecond length and time scales.

ANALYTICAL MODELS OF EXOPLANETARY ATMOSPHERES. I. ATMOSPHERIC DYNAMICS VIA THE SHALLOW WATER SYSTEM

Within the context of exoplanetary atmospheres, we present a comprehensive linear analysis of forced, damped, magnetized shallow water systems, exploring the effects of dimensionality, geometry (Cartesian, pseudo-spherical, and spherical), rotation, magnetic tension, and hydrodynamic and magnetic sources of friction. Across a broad range of conditions, we find that the key governing equation for atmospheres and quantum harmonic oscillators are identical, even when forcing (stellar irradiation), sources of friction (molecular viscosity, Rayleigh drag, and magnetic drag), and magnetic tension are included. The global atmospheric structure is largely controlled by a single key parameter that involves the Rossby and Prandtl numbers. This near-universality breaks down when either molecular viscosity or magnetic drag acts non-uniformly across latitude or a poloidal magnetic field is present, suggesting that these effects will introduce qualitative changes to the familiar chevron-shaped feature witnessed in simulations of atmospheric circulation. We also find that hydrodynamic and magnetic sources of friction have dissimilar phase signatures and affect the flow in fundamentally different ways, implying that using Rayleigh drag to mimic magnetic drag is inaccurate. We exhaustively lay down the theoretical formalism (dispersion relations, governing equations, and time-dependent wave solutions) for a broad suite of models. In all situations, we derive the steady state of an atmosphere, which is relevant to interpreting infrared phase and eclipse maps of exoplanetary atmospheres. We elucidate a pinching effect that confines the atmospheric structure to be near the equator. Our suite of analytical models may be used to develop decisively physical intuition and as a reference point for three-dimensional magnetohydrodynamic simulations of atmospheric circulation.

A numerical model for inert gas transport in the lung based on fractal airway morphology

Patients suffering from cystic fibrosis (CF) show thick secretions, mucus plugging and bronchiectasis in bronchial and alveolar ducts. This results in substantial structural changes of the airway morphology and heterogeneous ventilation. Disease progression and treatment effects are monitored by so-called gas washout tests, where the change in concentration of an inert gas is measured over a single or multiple breaths. The result of the tests based on the profile of the measured concentration is a marker for the severity of the ventilation inhomogeneity strongly affected by the airway morphology. However, it is hard to localize underlying obstructions to specific parts of the airways, especially if occurring in the lung periphery. In order to support the analysis of lung function tests (e.g. multi-breath washout), we developed a numerical model of the entire airway tree, coupling a lumped parameter model for the lung ventilation with a 4th-order accurate finite difference model of a 1D advection-diffusion equation for the transport of an inert gas. The boundary conditions for the flow problem comprise the pressure and flow profile at the mouth, which is typically known from clinical washout tests. The natural asymmetry of the lung morphology is approximated by a generic, fractal, asymmetric branching scheme which we applied for the conducting airways. A conducting airway ends when its dimension falls below a predefined limit. A model acinus is then connected to each terminal airway. The morphology of an acinus unit comprises a network of expandable cells. A regional, linear constitutive law describes the pressure-volume relation between the pleural gap and the acinus. The cyclic expansion (breathing) of each acinus unit depends on the resistance of the feeding airway and on the flow resistance and stiffness of the cells themselves. Special care was taken in the development of a conservative numerical scheme for the gas transport across bifurcations, handling spatially and temporally varying advective and diffusive fluxes over a wide range of scales. Implicit time integration was applied to account for the numerical stiffness resulting from the discretized transport equation. Local or regional modification of the airway dimension, resistance or tissue stiffness are introduced to mimic pathological airway restrictions typical for CF. This leads to a more heterogeneous ventilation of the model lung. As a result the concentration in some distal parts of the lung model remains increased for a longer duration. The inert gas concentration at the mouth towards the end of the expirations is composed of gas from regions with very different washout efficiency. This results in a steeper slope of the corresponding part of the washout profile.

Morphology and dynamics of the jets of comet 67P/Churyumov-Gerasimenko: Early-phase development

Aims. The OSIRIS camera onboard the Rosetta spacecraft obtained close-up views of the dust coma of comet 67P. The jet structures can be used to trace their source regions and to examine the possible effect of gas-surface interaction. Methods. We analyzed the wide-angle images obtained in the special dust observation sequences between August and September 2014. The jet features detected in different images were compared to study their time variability. The locations of the potential source regions of some of the jets are identified by ray tracing. We used a ring-masking technique to calculate the brightness distribution of dust jets along the projected distance. Results. The jets detected between August and September 2014 mostly originated in the Hapi region. Morphological changes appeared over a timescale of several days in September. The brightness slope of the dust jets is much steeper than the background coma. This might be related to the sublimation or fragmentation of the emitted dust grains. Interaction of the expanding gas flow with the cliff walls on both sides of Hapi could lead to erosion and material down-fall to the nucleus surface.