Comparative study of the different models for the calculation of BLEVEs overpressure effects

The BLEVE, acronym for Boiling Liquid Expanding Vapour Explosion, is one of the most dangerous accidents that can occur in pressure vessels. It can be defined as an explosion resulting from the failure of a vessel containing a pressure liquefied gas stored at a temperature significantly above its boiling point at atmospheric pressure. This phenomenon frequently appears when a vessel is engulfed by a fire: the heat causes the internal pressure to raise and the mechanical proprieties of the wall to decrease, with the consequent rupture of the tank and the instantaneous release of its whole content. After the breakage, the vapour outflows and expands and the liquid phase starts boiling due to the pressure drop. The formation and propagation of a distructive schock wave may occur, together with the ejection of fragments, the generation of a fireball if the stored fluid is flammable and immediately ignited or the atmospheric dispersion of a toxic cloud if the fluid contained inside the vessel is toxic. Despite the presence of many studies on the BLEVE mechanism, the exact causes and conditions of its occurrence are still elusive. In order to better understand this phenomenon, in the present study first of all the concept and definition of BLEVE are investigated. A historical analysis of the major events that have occurred over the past 60 years is described. A research of the principal causes of this event, including the analysis of the substances most frequently involved, is presented too. Afterwards a description of the main effects of BLEVEs is reported, focusing especially on the overpressure. Though the major aim of the present thesis is to contribute, with a comparative analysis, to the validation of the main models present in the literature for the calculation and prediction of the overpressure caused by BLEVEs. In line with this purpose, after a short overview of the available approaches, their ability to reproduce the trend of the overpressure is investigated. The overpressure calculated with the different models is compared with values deriving from events happened in the past and ad-hoc experiments, focusing the attention especially on medium and large scale phenomena. The ability of the models to consider different filling levels of the reservoir and different substances is analyzed too. The results of these calculations are extensively discussed. Finally some conclusive remarks are reported.

An Investigation On The Use Of Microalgae For Biogas Enrichment

There is a need for biomethane capture and carbon dioxide sequestration to mitigate evident global climate change. This research work investigated the potential for microalgae to remove CO2 from biogas as a biotechnical method for upgrading the thermal value for subsequent compression, liquification, or introduction to natural gas pipelines. Because biogas is largely methane, the effect of high methane environments on mixed microalgae was explored and found that specific carbon utilization rates were not statistically different when microalgae were exposed to biogas environments (70% v/v CH4) , relative to high CO2 environment. The uses of conventional bubbled column photobioreactors (PBR) were assessed for CO2 removal and subsequent CH4 enrichment. A continuously-bubbled biogas PBR (cB-PBR5) and intermittently-bubbled biogas PBR (iB-PBR) experienced CO2 loading rates of about 1664 and 832 mg C/L*day and showed 30.0 and 60.1 % carbon removal, respectively. However, a lack of biogas enrichment and issues associated growth inhibition due to high CO2 environments as well as stripping the dissolved gases, namely oxygen and nitrogen, from the bulk liquid and introduction to the outlet gas prompted the consideration for gas/liquid separation using nonporous hollow-fiber (HF) membranes for CO2 transfer. The potential for two non-porous HF membrane materials [polydimethylsiloxane (PDMS) and composite polyurethane (PU)] were modeled along fiber length using a mechanistic model based on polymeric material transport properties (Gilmore et al., 2009). Based on a high CO2:CH4 permeability selectivity for PU of 76.2 the model predicted gas enrichment along an 8.5 cm fiber length. Because PDMS permeability selectivity is low (3.5), evident gas transfer was not predicated along a 34.3 cm length. Both of these HF materials were implemented in hollow-fiber membrane-carbonated biofilm (HFMcB) PBRs for microalgal-mediated biogas enrichment. Phototrophic biofilm colonization occurred on the membrane, where CO2 concentration was greatest. The presence of a biofilm demonstrated greater resiliency to high CO2 environments, compared to the conventional PBRs. However, as the PDMS model predicted, the PDMS HFMcBs did not demonstrate gas enrichment. These reactors received CO2 loading rates of 200 mg C/L*day based on PDMS permeability flux and showed approximately 65% removal of the total C transferred across the membrane. Thus, the HFMcBs demonstrated controlled carbonation of the bulk liquid via a nonporous HF membrane. Likewise, the experimental PU HFMcB did not show gas enrichment yet this result should be further explored due to the high permeability selectivity of the polymeric material. Chemical stratifications, namely pH and dissolved O2, present in a PDMS membrane-carbonated biofilm were analyzed using electrochemical microsensors. Results indicated that high DO (20 mg L-1) exists at surface of the biofilm where light availability is greatest and low pH microenvironments (pH=5.40) exist deep in the biofilm where the diffusive flux of CO2 drives transfer through the biofilm. The presence of a 400-600 ¿m liquid phase boundary layer was evident from microsensor profiles. Cryosectioning of the biofilm samples showed the biofilm to be approximately 1.17 ± 0.07 mm thick, suggesting that the high localized concentration of biomass associated with the phototrophic biofilm aided in overcoming inhibition in a microenvironment dominated by CO2(aq). Challenges of biofilm detachment and PBR fouling as well as microalgal growth inhibition in the presence of high CO2 content remain for applications of microalgae for biogas enrichment.

Diesel spray mixing limited vaporization with non-ideal and multi-component fuel thermophysical property effects

Internal combustion engines are, and will continue to be, a primary mode of power generation for ground transportation. Challenges exist in meeting fuel consumption regulations and emission standards while upholding performance, as fuel prices rise, and resource depletion and environmental impacts are of increasing concern. Diesel engines are advantageous due to their inherent efficiency advantage over spark ignition engines; however, their NOx and soot emissions can be difficult to control and reduce due to an inherent tradeoff. Diesel combustion is spray and mixing controlled providing an intrinsic link between spray and emissions, motivating detailed, fundamental studies on spray, vaporization, mixing, and combustion characteristics under engine relevant conditions. An optical combustion vessel facility has been developed at Michigan Technological University for these studies, with detailed tests and analysis being conducted. In this combustion vessel facility a preburn procedure for thermodynamic state generation is used, and validated using chemical kinetics modeling both for the MTU vessel, and institutions comprising the Engine Combustion Network international collaborative research initiative. It is shown that minor species produced are representative of modern diesel engines running exhaust gas recirculation and do not impact the autoignition of n-heptane. Diesel spray testing of a high-pressure (2000 bar) multi-hole injector is undertaken including non-vaporizing, vaporizing, and combusting tests, with sprays characterized using Mie back scatter imaging diagnostics. Liquid phase spray parameter trends agree with literature. Fluctuations in liquid length about a quasi-steady value are quantified, along with plume to plume variations. Hypotheses are developed for their causes including fuel pressure fluctuations, nozzle cavitation, internal injector flow and geometry, chamber temperature gradients, and turbulence. These are explored using a mixing limited vaporization model with an equation of state approach for thermopyhysical properties. This model is also applied to single and multi-component surrogates. Results include the development of the combustion research facility and validated thermodynamic state generation procedure. The developed equation of state approach provides application for improving surrogate fuels, both single and multi-component, in terms of diesel spray liquid length, with knowledge of only critical fuel properties. Experimental studies are coupled with modeling incorporating improved thermodynamic non-ideal gas and fuel

COMPUTATIONAL STUDY OF MICROSTRUCTURE-PROPERTYMECHANISM RELATIONS IN FERROIC COMPOSITES

Ferroic materials, as notable members of smart materials, have been widely used in applications that perform sensing, actuation and control. The macroscopic property change of ferroic materials may become remarkably large during ferroic phase transition, leading to the fact that the macroscopic properties can be tuned by carefully applying a suitable external field (electric, magnetic, stress). To obtain an enhancement in physical and/or mechanical properties, different kinds of ferroic composites have been fabricated. The properties of a ferroic composite are determined not only by the properties and relative amounts of the constituent phases, but also by the microstructure of individual phase such as the phase connectivity, phase size, shape and spatial arrangement. This dissertation mainly focuses on the computational study of microstructure – property – mechanism relations in two representative ferroic composites, i.e., two-phase particulate magnetoelectric (ME) composite and polymer matrix ferroelectric composite. The former is a great example of ferroic composite exhibiting a new property and functionality that neither of the constituent phases possesses individually. The latter well represents the kind of ferroic composites having property combinations that are better than the existing materials. Phase field modeling was employed as the computing tool, and the required models for ferroic composites were developed based on existing models for monolithic materials. Extensive computational simulations were performed to investigate the microstructure-property relations and the underlying mechanism in ferroic composites. In particulate, it is found that for ME composite 0-3 connectivity (isolated magnetostrictive phase) is necessary to exhibit ME effect, and small but finite electrical conductivity of isolated magnetic phase can beneficially enhance ME effect. It is revealed that longitudinal and transverse ME coefficients of isotropic 0-3 particulate composites can be effectively tailored by controlling magnetic domain structures without resort to anisotropic two-phase microstructures. Simulations also show that the macroscopic properties of the ferroelectricpolymer composites critically depend on the ferroelectric phase connectivity while are not sensitive to the sizes and internal grain structures of the ceramic particles. Texturing is found critical to exploit the paraelectric«ferroelectric phase transition and nonlinear polarization behavior in paraelectric polycrystal and its polymer matrix composite. Additionally, a Diffuse Interface Field model was developed to simulate packing and motion in liquid phase which is promising for studying the fabrication of particulatepolymer composites.

MAGNETO-PHOTONIC CRYSTALS FOR OPTICAL SENSING APPLICATIONS

Among the optical structures investigated for optical sensing purpose, a significant amount of research has been conducted on photonic crystal based sensors. A particular advantage of photonic crystal based sensors is that they show superior sensitivity for ultra-small volume sensing. In this study we investigate polarization changes in response to the changes in the cover index of magneto-optic active photonic band gap structures. One-dimensional photonic-band gap structures fabricated on iron garnet materials yield large polarization rotations at the band gap edges. The enhanced polarization effects serve as an excellent tool for chemical sensing showing high degree of sensitivity for photonic crystal cover refractive index changes. The one dimensional waveguide photonic crystals are fabricated on single-layer bismuth-substituted rare earth iron garnet films ((Bi, Y, Lu)3(Fe, Ga)5O12 ) grown by liquid phase epitaxy on gadolinium gallium garnet substrates. Band gaps have been observed where Bragg scattering conditions links forward-going fundamental waveguide modes to backscattered high-order waveguide modes. Large near-band-edge polarization rotations which increase progressively with backscattered-mode order have been experimentally demonstrated for multiple samples with different composition, film thickness and fabrication parameters. Experimental findings are supported by theoretical analysis of Bloch modes polarization states showing that large near stop-band edge rotations are induced by the magneto-photonic crystal. Theoretical and experimental analysis conducted on polarization rotation sensitivity to waveguide photonic crystal cover refractive index changes shows a monotonic enhancement of the rotation with cover index. The sensor is further developed for selective chemical sensing by employing Polypyrrole as the photonic crystal cover layer. Polypyrrole is one of the extensively studied conducting polymers for selective analyte detection. Successful detection of aqueous ammonia and methanol has been achieved with Polypyrrole deposited magneto-photonic crystals.

High-resolution computer simulations of EKC

The electrophoresis simulation software, GENTRANS, has been modified to include the interaction of analytes with an electrolyte additive to allow the simulation of liquid-phase EKC separations. The modifications account for interaction of weak and strong acid and base analytes with a single weak or strong acid or base background electrolyte additive and can be used to simulate a range of EKC separations with both charged and neutral additives. Simulations of separations of alkylphenyl ketones under real experimental conditions were performed using mobility and interaction constant data obtained from the literature and agreed well with experimental separations. Migration times in fused-silica capillaries and linear polyacrylamide-coated capillaries were within 7% of the experimental values, while peak widths were always narrower than the experimental values, but were still within 50% of those obtained by experiment. Simulations of sweeping were also performed; although migration time agreement was not as good as for simple EKC separations, peak widths were in good agreement, being within 1-50% of the experimental values. All simulations for comparison with experimental data were performed under real experimental conditions using a 47 cm capillary and a voltage of 20 kV and represent the first quantitative attempt at simulating EKC separations with and without sweeping.

MODELING REACTIVE GAS UPTAKE, TRANSPORT, AND TRANSFORMATION IN AGGREGATED SOILS

Gas diffusion research in soils covers, to a large extent, the transport behavior of practically insoluble gases. We extend the mathematical description of gas transport to include reactive gaseous components that hydrolyze in water such as SO2 and CO2. The path between the free atmosphere and the microporous niches is modeled by assuming penetration through gas-filled macropores, air-water phase transfer, and diffusion and speciation in the liquid phase. For hydrolyzable gases, the rate of mass transfer into and the total absorption capacity of the soil solution may be high. Both the capacity and the transfer rate are influenced by the soil-solution pH; for high pH, they become extremely high for SO2. The soil absorption of such gases is also influenced by soil structure. Well-aerated, near-neutral soils are a potentially important sink for SO2.

Water–montmorillonite systems: Neutron scattering and tracer throughdiffusion studies

Designs for deep geological respositories of nuclear waste include bentonite as a hydraulic and chemisorption buffer material to protect the biosphere from leakage of radionuclides. Bentonite is chosen because it is a cheap, naturally occurring material with the required properties. It consists essentially of montmorillonite, a swelling clay mineral. Upon contact with groundwater such clays can seal the repository by incorporating water in the interlayers of their crystalline structure. The intercalated water exhibits significantly different properties to bulk water in the surrounding interparticle pores, such as lower diffusion coefficients (González Sánchez et. al. 2008). This doctoral thesis presents water distribution and diffusion behavior on various time and space scales in montmorillonite. Experimental results are presented for Na- and Cs-montmorillonite samples with a range of bulk dry densities (0.8 to 1.7 g/cm3). The experimental methods employed were neutron scattering (backscattering, diffraction, time-of-flight), adsorption measurements (water, nitrogen) and tracer-through diffusion. For the tracer experiments the samples were fully saturated via the liquid phase under volume-constrained conditions. In contrast, for the neutron scattering experiments, the samples were hydrated via the vapor phase and subsequently compacted, leaving a significant fraction of interparticle pores unfilled with water. Owing to these differences in saturation, the water contents of the samples for neutron scattering were characterized by gravimetry whereas those for the tracer experiments were obtained from the bulk dry density. The amount of surface water in interlayer pores could be successfully discriminated from the amount of bulk-like water in interparticle pores in Na- and Csmontmorillonite using neutron spectroscopy. For the first time in the literature, the distribution of water between these two pore environments was deciphered as a function of gravimetric water content. The amount was compared to a geometrical estimation of the amount of interlayer and interparticle water determined by neutron diffraction and adsorption measurements. The relative abundances of the 1 to 4 molecular water layers in the interlayer were determined from the area ratios of the (001)-diffraction peaks. Depending on the characterization method, different fractions of surface water and interlayer water were obtained. Only surface and interlayer water exists in amontmorillonite with water contents up to 0.18 g/g according to spectroscopic measurements and up to 0.32 g/g according to geometrical estimations, respectively. At higher water contents, bulk-like and interparticle water also exists. The amounts increase monotonically, but not linearly, from zero to 0.33 g/g for bulk-like water and to 0.43 g/g for interparticle water. It was found that water most likely redistributes between the surface and interlayer sites during the spectroscopic measurements and therefore the reported fraction is relevant only below about -10 ºC (Anderson, 1967). The redistribution effect can explain the discrepancy in fractions between the methods. In a novel approach the fractions of water in different pore environments were treated as a fixed parameter to derive local diffusion coefficients for water from quasielastic neutron scattering data, in particular for samples with high water contents. Local diffusion coefficients were obtained for the 1 to 4 molecular water layers in the interlayer of 0.5·10–9, 0.9·10–9, 1.5·10–9 and 1.4·10–9 m²/s, respectively, taking account of the different water fractions (molecular water layer, bulk-like water). The diffusive transport of 22Na and HTO through Na-montmorillonite was measured on the laboratory experimental scale (i.e. cm, days) by tracer through-diffusion experiments. We confirmed that diffusion of HTO is independent of the ionic strength of the external solution in contact with the clay sample but dependent on the bulk dry density. In contrast, the diffusion of 22Na was found to depend on both the ionic strength of the pore solution and on the bulk dry density. The ratio of the pore and surface diffusion could be experimentally determined for 22Na from the dependence of the diffusion coefficient on the ionic strength. Activation energies were derived from the temperaturedependent diffusion coefficients via the Arrhenius relation. In samples with high bulk dry density the activation energies are slightly higher than those of bulk water whereas in low density samples they are lower. The activation energies as a function of ionic strengths of the pore solutions are similar for 22Na and HTO. The facts that (i) the slope of the logarithmic effective diffusion coefficients as a function of the logarithmic ionic strength is less than unity for low bulk dry densities and (ii) two water populations can be observed for high gravimetric water contents (low bulk dry densities) support the interlayer and interparticle porosity model proposed by Glaus et al. (2007), Bourg et al. (2006, 2007) and Gimmi and Kosakowski (2011).

The isotopic composition of pore-water from Tournemire argillite (France): An inter-comparison study.

The determination of stable isotope contents of pore-water from consolidated argillaceous rocks remains a critical issue. In order to understand the processes involved in techniques developed for acquiring stable isotope compositions of pore-water, a comparative study between different methods was based on core samples of the Tournemire argillite. It concerns two water extraction techniques based on vacuum distillation and two pore-water equilibration techniques (radial diffusion in liquid phase and diffusive exchange in vapor phase). The water-content values obtained from vacuum distillation at 50 °C are always the lowest, on average 8% lower than the values obtained by heating at 105 °C and 17% lower than the values obtained by heating at 150 °C. The amounts of pore-water estimated from vacuum distillation at 105 °C and 150 °C and from radial diffusion method are in good agreement with those determined by heating. On the contrary, the vapor exchange method provides the highest values of water contents. Concerning stable isotope data, a good agreement was found between those obtained by equilibration techniques and those of fracture water, especially for 2H. Vacuum distillation at high temperature (particularly at 150 °C) also provided results consistent with data of fracture fluids. On the other hand, distillation at 50 °C provides a systematic depletion in heavy isotopes (about –20‰ for 2H and –2.7‰ for 18O) that can be modelled by an incomplete Rayleigh-type distillation process.

Minimally invasive photopolymerization in intervertebral disc tissue cavities

Photopolymerized hydrogels are commonly used for a broad range of biomedical applications. As long as the polymer volume is accessible, gels can easily be hardened using light illumination. However, in clinics, especially for minimally invasive surgery, it becomes highly challenging to control photopolymerization. The ratios between polymerization- volume and radiating-surface-area are several orders of magnitude higher than for ex-vivo settings. Also tissue scattering occurs and influences the reaction. We developed a Monte Carlo model for photopolymerization, which takes into account the solid/liquid phase changes, moving solid/liquid-boundaries and refraction on these boundaries as well as tissue scattering in arbitrarily designable tissue cavities. The model provides a tool to tailor both the light probe and the scattering/absorption properties of the photopolymer for applications such as medical implants or tissue replacements. Based on the simulations, we have previously shown that by adding scattering additives to the liquid monomer, the photopolymerized volume was considerably increased. In this study, we have used bovine intervertebral disc cavities, as a model for spinal degeneration, to study photopolymerization in-vitro. The cavity is created by enzyme digestion. Using a custom designed probe, hydrogels were injected and photopolymerized. Magnetic resonance imaging (MRI) and visual inspection tools were employed to investigate the successful photopolymerization outcomes. The results provide insights for the development of novel endoscopic light-scattering polymerization probes paving the way for a new generation of implantable hydrogels.

A review of air–ice chemical and physical interactions (AICI): liquids, quasi-liquids, and solids in snow

Snow in the environment acts as a host to rich chemistry and provides a matrix for physical exchange of contaminants within the ecosystem. The goal of this review is to summarise the current state of knowledge of physical processes and chemical reactivity in surface snow with relevance to polar regions. It focuses on a description of impurities in distinct compartments present in surface snow, such as snow crystals, grain boundaries, crystal surfaces, and liquid parts. It emphasises the microscopic description of the ice surface and its link with the environment. Distinct differences between the disordered air–ice interface, often termed quasi-liquid layer, and a liquid phase are highlighted. The reactivity in these different compartments of surface snow is discussed using many experimental studies, simulations, and selected snow models from the molecular to the macro-scale. Although new experimental techniques have extended our knowledge of the surface properties of ice and their impact on some single reactions and processes, others occurring on, at or within snow grains remain unquantified. The presence of liquid or liquid-like compartments either due to the formation of brine or disorder at surfaces of snow crystals below the freezing point may strongly modify reaction rates. Therefore, future experiments should include a detailed characterisation of the surface properties of the ice matrices. A further point that remains largely unresolved is the distribution of impurities between the different domains of the condensed phase inside the snowpack, i.e. in the bulk solid, in liquid at the surface or trapped in confined pockets within or between grains, or at the surface. While surface-sensitive laboratory techniques may in the future help to resolve this point for equilibrium conditions, additional uncertainty for the environmental snowpack may be caused by the highly dynamic nature of the snowpack due to the fast metamorphism occurring under certain environmental conditions. Due to these gaps in knowledge the first snow chemistry models have attempted to reproduce certain processes like the long-term incorporation of volatile compounds in snow and firn or the release of reactive species from the snowpack. Although so far none of the models offers a coupled approach of physical and chemical processes or a detailed representation of the different compartments, they have successfully been used to reproduce some field experiments. A fully coupled snow chemistry and physics model remains to be developed.

Dissolution-Precipitation in Porous Media: Experiments and Modelling

The evolution of porosity due to dissolution/precipitation processes of minerals and the associated change of transport parameters are of major interest for natural geological environments and engineered underground structures. We designed a reproducible and fast to conduct 2D experiment, which is flexible enough to investigate several process couplings implemented in the numerical code OpenGeosys-GEM (OGS-GEM). We investigated advective-diffusive transport of solutes, effect of liquid phase density on advective transport, and kinetically controlled dissolution/precipitation reactions causing porosity changes. In addition, the system allowed to investigate the influence of microscopic (pore scale) processes on macroscopic (continuum scale) transport. A Plexiglas tank of dimension 10 × 10 cm was filled with a 1 cm thick reactive layer consisting of a bimodal grain size distribution of celestite (SrSO4) crystals, sandwiched between two layers of sand. A barium chloride solution was injected into the tank causing an asymmetric flow field to develop. As the barium chloride reached the celestite region, dissolution of celestite was initiated and barite precipitated. Due to the higher molar volume of barite, its precipitation caused a porosity decrease and thus also a decrease in the permeability of the porous medium. The change of flow in space and time was observed via injection of conservative tracers and analysis of effluents. In addition, an extensive post-mortem analysis of the reacted medium was conducted. We could successfully model the flow (with and without fluid density effects) and the transport of conservative tracers with a (continuum scale) reactive transport model. The prediction of the reactive experiments initially failed. Only the inclusion of information from post-mortem analysis gave a satisfactory match for the case where the flow field changed due to dissolution/precipitation reactions. We concentrated on the refinement of post-mortem analysis and the investigation of the dissolution/precipitation mechanisms at the pore scale. Our analytical techniques combined scanning electron microscopy (SEM) and synchrotron X-ray micro-diffraction/micro-fluorescence performed at the XAS beamline (Swiss Light Source). The newly formed phases include an epitaxial growth of barite micro-crystals on large celestite crystals (epitaxial growth) and a nano-crystalline barite phase (resulting from the dissolution of small celestite crystals) with residues of celestite crystals in the pore interstices. Classical nucleation theory, using well-established and estimated parameters describing barite precipitation, was applied to explain the mineralogical changes occurring in our system. Our pore scale investigation showed limits of the continuum scale reactive transport model. Although kinetic effects were implemented by fixing two distinct rates for the dissolution of large and small celestite crystals, instantaneous precipitation of barite was assumed as soon as oversaturation occurred. Precipitation kinetics, passivation of large celestite crystals and metastability of supersaturated solutions, i.e. the conditions under which nucleation cannot occur despite high supersaturation, were neglected. These results will be used to develop a pore scale model that describes precipitation and dissolution of crystals at the pore scale for various transport and chemical conditions. Pore scale modelling can be used to parameterize constitutive equations to introduce pore-scale corrections into macroscopic (continuum) reactive transport models. Microscopic understanding of the system is fundamental for modelling from the pore to the continuum scale.

Concentration of nutrients and primary production in sea ice of the high-latitudinal Arctic

Dynamics of phosphates and silicates in sea ice of the high-latitudinal Arctic are considered for period from November 2005 to May 2006. It is shown that, during ice formation, silicates are included into it in the same ratio to salinity that is characteristic of under-ice water. Further dynamics of silicates are determined by their bioassimilation with beginning of the polar day and by biogenic silicon accumulation at bottom meltwater pools with subsequent leaching. Phosphates are included into ice in a ratio higher than that occurring in the under-ice water. This is caused by the fact that liquid phase of sea ice represents composition of the surface microlayer at the ice-water interface, which is enriched in organic matter and in products of its destruction (particularly in phosphates). With onset of the polar day, content of phosphates first markedly increases (due to photo oxidation of biogenic organic matter) and then decreases because of bioassimilation. At the beginning of the polar day, primary production of diatoms was estimated to be ~0.3 mg C/m**2/day.

Li & Pb MD Sinulations

Pb17Li is today a reference breeder material in diverse fusion R&D programs worldwide. Extracting dynamic and structural properties of liquid LiPb mixtures via molecular dynamics simulations, represent a crucial step for multiscale modeling efforts in order to understand the suitability of this compound for future Nuclear Fusion technologies. At present a Li-Pb cross potential is not available in the literature. Here we present our first results on the validation of two semi-empirical potentials for Li and Pb in liquid phase. Our results represent the establishment of a solid base as a previous but crucial step to implement a LiPb cross potential. Structural and thermodynamical analyses confirm that the implemented potentials for Li and Pb are realistic to simulate both elements in the liquid phase.

Mechanical properties of tungsten alloys with Y2O3 and titanium additions

In this research the mechanical behaviour of pure tungsten (W) and its alloys (2 wt.% Ti–0.47 wt.% Y2O3 and 4 wt.% Ti–0.5 wt.% Y2O3) is compared. These tungsten alloys, have been obtained by powder metallurgy. The yield strength, fracture toughness and elastic modulus have been studied in the temperature interval of 25 °C to 1000 °C. The results have shown that the addition of Ti substantially improves the bending strength and toughness of W, but it also dramatically increases the DBTT. On the other hand, the addition of 0.5% Y2O3, is enough to improve noticeably the oxidation behaviour at the higher temperatures. The grain size, fractography and microstructure are studied in these materials. Titanium is a good grain growth inhibitor and effective precursor of liquid phase in HIP. The simultaneous presence of Y2O3 and Ti permits to obtain materials with low pores presence