Full scale test for explosion water barriers in small cross-section galleries

Underground coal mines explosions generally arise from the inflammation of a methane/air mixture. This explosion can also generate a subsequent coal dust explosion. Traditionally such explosions have being fought eliminating one or several of the factors needed by the explosion to take place. Although several preventive measures are taken to prevent explosions, other measures should be considered to reduce the effects or even to extinguish the flame front. Unlike other protection methods that remove one or two of the explosion triangle elements, namely; the ignition source, the oxidizing agent and the fuel, explosion barriers removes all of them: reduces the quantity of coal in suspension, cools the flame front and the steam generated by vaporization removes the oxygen present in the flame. Passive water barriers are autonomous protection systems against explosions that reduce to a satisfactory safety level the effects of methane and/or flammable dust explosions. The barriers are activated by the pressure wave provoked in the explosion destroying the barrier troughs and producing a uniform dispersion of the extinguishing agent throughout the gallery section in quantity enough to extinguish the explosion flame. Full scale tests have been carried out in Polish Barbara experimental mine at GIG Central Mining Institute in order to determine the requirements and the optimal installation conditions of these devices for small sections galleries which are very frequent in the Spanish coal mines. Full scale tests results have been analyzed to understand the explosion timing and development, in order to assess on the use of water barriers in the typical small crosssection Spanish galleries. Several arrangements of water barriers have been designed and tested to verify the effectiveness of the explosion suppression in each case. The results obtained demonstrate the efficiency of the water barriers in stopping the flame front even with smaller amounts of water than those established by the European standard. According to the tests realized, water barriers activation times are between 0.52 s and 0.78 s and the flame propagation speed are between 75 m/s and 80 m/s. The maximum pressures (Pmax) obtained in the full scale tests have varied between 0.2 bar and 1.8 bar. Passive barriers protect effectively against the spread of the flame but cannot be used as a safeguard of the gallery between the ignition source and the first row of water troughs or bags, or even after them, as the pressure could remain high after them even if the flame front has been extinguished.

Application of protection means against explosions in underground mines

Underground coal mines explosions generally arise from the inflammation of a methane/air mixture. This explosion can also generate a subsequent coal dust explosion. Traditionally such explosions have being fought eliminating one or several of the factors needed by the explosion to take place. Although several preventive measures are taken to prevent explosions, other measures should be considered to reduce the effects or even to extinguish the flame front. Unlike other protection methods that remove one or two of the explosion triangle elements, namely; the ignition source, the oxidizing agent and the fuel, explosion barriers removes all of them: reduces the quantity of coal in suspension, cools the flame front and the steam generated by vaporization removes the oxygen present in the flame. The present paper is essentially based on the comprehensive state-of–the-art of Protective Systems in underground coal mines, and particularly on the application of Explosion Barriers to improve safety level in Spanish coal mining industry. After an exhaustive study of series EN 14591 standards covering explosion prevention and protection in underground mines, authors have proven explosion barriers effectiveness in underground galleries by Full Scale Tests performed in Polish Barbara experimental mine, showing that the barriers can reduce the effects of methane and/or flammable coal dust explosions to a satisfactory safety level.

Experimental and numerical study of cw green laser crystallization of a-Si:H thin films

Crystallization and grain growth technique of thin film silicon are among the most promising methods for improving efficiency and lowering cost of solar cells. A major advantage of laser crystallization and annealing over conventional heating methods is its ability to limit rapid heating and cooling to thin surface layers. Laser energy is used to heat the amorphous silicon thin film, melting it and changing the microstructure to polycrystalline silicon (poly-Si) as it cools. Depending on the laser density, the vaporization temperature can be reached at the center of the irradiated area. In these cases ablation effects are expected and the annealing process becomes ineffective. The heating process in the a-Si thin film is governed by the general heat transfer equation. The two dimensional non-linear heat transfer equation with a moving heat source is solve numerically using the finite element method (FEM), particularly COMSOL Multiphysics. The numerical model help to establish the density and the process speed range needed to assure the melting and crystallization without damage or ablation of the silicon surface. The samples of a-Si obtained by physical vapour deposition were irradiated with a cw-green laser source (Millennia Prime from Newport-Spectra) that delivers up to 15 W of average power. The morphology of the irradiated area was characterized by confocal laser scanning microscopy (Leica DCM3D) and Scanning Electron Microscopy (SEM Hitachi 3000N). The structural properties were studied by micro-Raman spectroscopy (Renishaw, inVia Raman microscope).

Dynamics of thermal ignition of spray flames in mixing layers

Conditions are identified under which analyses of laminar mixing layers can shed light on aspects of turbulent spray combustion. With this in mind, laminar spray-combustion models are formulated for both non-premixed and partially premixed systems. The laminar mixing layer separating a hot-air stream from a monodisperse spray carried by either an inert gas or air is investigated numerically and analytically in an effort to increase understanding of the ignition process leading to stabilization of high-speed spray combustion. The problem is formulated in an Eulerian framework, with the conservation equations written in the boundary-layer approximation and with a one-step Arrhenius model adopted for the chemistry description. The numerical integrations unveil two different types of ignition behaviour depending on the fuel availability in the reaction kernel, which in turn depends on the rates of droplet vaporization and fuel-vapour diffusion. When sufficient fuel is available near the hot boundary, as occurs when the thermochemical properties of heptane are employed for the fuel in the integrations, combustion is established through a precipitous temperature increase at a well-defined thermal-runaway location, a phenomenon that is amenable to a theoretical analysis based on activation-energy asymptotics, presented here, following earlier ideas developed in describing unsteady gaseous ignition in mixing layers. By way of contrast, when the amount of fuel vapour reaching the hot boundary is small, as is observed in the computations employing the thermochemical properties of methanol, the incipient chemical reaction gives rise to a slowly developing lean deflagration that consumes the available fuel as it propagates across the mixing layer towards the spray. The flame structure that develops downstream from the ignition point depends on the fuel considered and also on the spray carrier gas, with fuel sprays carried by air displaying either a lean deflagration bounding a region of distributed reaction or a distinct double-flame structure with a rich premixed flame on the spray side and a diffusion flame on the air side. Results are calculated for the distributions of mixture fraction and scalar dissipation rate across the mixing layer that reveal complexities that serve to identify differences between spray-flamelet and gaseous-flamelet problems.

Flamelet structures in spray ignition

In typical liquid-fueled burners the fuel is injected as a high-velocity liquid jet that breaks up to form the spray. The initial heating and vaporization of the liquid fuel rely on the relatively large temperatures of the sourrounding gas, which may include hot combustion products and preheated air. The heat exchange between the liquid and the gas phases is enhanced by droplet dispersion arising from the turbulent motion. Chemical reaction takes place once molecular mixing between the fuel vapor and the oxidizer has occurred in mixing layers separating the spray flow from the hot air stream. Since in most applications the injection velocities are much larger than the premixed-flame propagation velocity, combustion stabilization relies on autoignition of the fuel-oxygen mixture, with the combustion stand-off distance being controlled by the interaction of turbulent transport, droplet heating and vaporization, and gas-phase chemical reactions. In this study, conditions are identified under which analyses of laminar flamelets canshed light on aspects of turbulent spray ignition. This study extends earlier fundamental work by Liñan & Crespo (1976) on ignition in gaseous mixing layers to ignition of sprays. Studies of laminar mixing layers have been found to be instrumental in developing un-derstanding of turbulent combustion (Peters 2000), including the ignition of turbulent gaseous diffusion flames (Mastorakos 2009). For the spray problem at hand, the configuration selected, shown in Figure 1, involves a coflow mixing layer formed between a stream of hot air moving at velocity UA and a monodisperse spray moving at velocity USUA. The boundary-layer approximation will be used below to describe the resulting sl ender flow, which exhibits different igniting behaviors depending on the characteristics of t he fuel. In this approximation, consideration of the case U A = U S enables laminar ignition distances to be related to ignition times of unstrained spray flamelets, thereby pro viding quantitative information of direct applicability in regions of low scala r dissipation-rate in turbulent reactive flows (see the discussion in pp. 181–186 of Peters (2000)) . This report is organized as follows. Effects of droplet dispersion dynamics on ignition of sprays in turbulent mixing layers are discussed in Section 2. The formulation f or ignition in laminar mixing layers is outlined in Sections 3 and 4. The results are presented in Section 5. In Section 6, the mixture-fraction field and associated scalar dissipat ion rates for spray ignition are discussed. Finally, some brief conclusions are drawn in Section 7.

Metodología de simulación numérica de la transferencia de calor por ebullición en película en torno a superficies esféricas refrigeradas por líquido saturado en convección libre: aplicación a reactores nucleares de lecho de bolas (PBR)

En la presente Tesis se realizó un análisis numérico, usando el código comercial Ansys-Fluent, de la refrigeración de una bola de combustible de un reactor de lecho de bolas (PBR, por sus siglas en inglés), ante un escenario de emergencia en el cual el núcleo sea desensamblado y las bolas se dejen caer en una piscina de agua, donde el mecanismo de transferencia de calor inicialmente sería la ebullición en película, implicando la convección y la radiación al fluido. Previamente se realizaron pruebas de validación, comparando los resultados numéricos con datos experimentales disponibles en la literatura para tres geometrías diferentes, lo cual permitió seleccionar los esquemas y modelos numéricos con mejor precisión y menor costo computacional. Una vez identificada la metodología numérica, todas las pruebas de validación fueron ampliamente satisfactorias, encontrándose muy buena concordancia en el flujo de calor promedio con los datos experimentales. Durante estas pruebas de validación se lograron caracterizar numéricamente algunos parámetros importantes en la ebullición en película con los cuales existen ciertos niveles de incertidumbre, como son el factor de acoplamiento entre convección y radiación, y el factor de corrección del calor latente de vaporización. El análisis térmico de la refrigeración de la bola del reactor por ebullición en película mostró que la misma se enfría, a pesar del calor de decaimiento, con una temperatura superficial de la bola que desciende de forma oscilatoria, debido al comportamiento inestable de la película de vapor. Sin embargo, la temperatura de esta superficie tiene una buena uniformidad, notándose que las áreas mejor y peor refrigeradas están localizadas en la parte superior de la bola. Se observó la formación de múltiples domos de vapor en diferentes posiciones circunferenciales, lo cual causa que el área más caliente de la superficie se localice donde se forman los domos más grandes. La separación entre los domos de vapor fue consistente con la teoría hidrodinámica, con la adición de que la separación entre domos se reduce a medida que evolucionan y crecen, debido a la curvatura de la superficie. ABSTRACT A numerical cooling analysis of a PBR fuel pebble, after an emergency scenario in which the nucleus disassembly is made and the pebbles are dropped into a water pool, transmitting heat by film boiling, involving convection and radiation to the fluid, is carried out in this Thesis. First, were performed validation tests comparing the numerical results with experimental works available for three different geometries, which allowed the selection of numerical models and schemes with better precision and lower computational cost. Once identified the numerical methodology, all validation tests were widely satisfactory, finding very good agreement with experimental works in average heat flux. During these validation tests were achieved numerically characterize some important parameters in film boiling with which there are certain levels of uncertainty, such as the coupling factor between convection and radiation, and the correction factor of the latent heat of vaporization. The thermal analysis of pebble cooling by film boiling shows that despite its decay heat, cooling occurs, with pebble surface temperature descending from an oscillatory manner, due to the instability of the vapor film. However, the temperature of this surface has a good uniformity, noting that the best and worst refrigerated area is located at the top of the pebble. The formation of multiple vapor domes at different circumferential positions is observed, which cause that the hottest area of the surface was located where biggest vapor domes were formed. The separation between vapor domes was consistent with the hydrodynamic theory, with the addition that the separation is reduced as the vapor dome evolves and grows, due to the surface curvature.

(Table 2) Concentrations of dissolved aluminum in interstitial water from sediments of the Indian Ocean

Results are presented of application of laser stepwise photoionization of atoms in combination with thermal atomization of matter in vacuum for direct determination of aluminum dissolved in sea and interstitial waters. Dry residue from evaporation of 40 ?l sea water was atomized in a crucible at 1800°C, and aluminum atoms in the beam thus formed were energized into Rydberg state in two steps by two tunable dye laser beams; the atoms were then ionized by an electric pulse and resulting ions were recorded by secondary emission electron multiplier (ion detector). Ionic signal dependence on sample vaporization time was studied. The procedure is suggested for separating out a selective signal in a single measurement. Dissolved aluminum concentrations in interstitial waters of the Indian Ocean and in waters of the river-sea zone were determined using preliminarily plotted calibration characteristics for aluminum solutions in deionized and sea waters. The minimum detectable Al concentration in seawater was 1 ?g/l that corresponds to 40 pg of Al in a sample.

Geochemistry of microtektites from Victoria Land Transantarctic Mountains

We extended the petrographic and geochemical dataset for the recently discovered Transantarctic Mountain microtektites in order to check our previous claim that they are related to the Australasian strewn field. Based on color and composition, the 465 microtektites so far identified include two groups of transparent glass spheres less than ca. 800 µm in diameter: the most abundant pale-yellow, or normal, microtektites, and the rare pale-green, or high-Mg, microtektites. The major element composition of normal microtektites determined through electron microprobe analysis is characterized by high contents of silica (SiO2 = 71.5 ± 3.6 (1 sigma) wt%) and alumina (Al2O3 = 15.5 ± 2.2 (1 sigma) wt%), low total alkali element contents (0.50-1.85 wt%), and MgO abundances <6 wt%. The high-Mg microtektites have a distinctly higher MgO content >10 wt%. Transantarctic Mountain microtektites contain rare silica-rich (up to 93 wt% SiO2) glassy inclusions similar to those found in two Australasian microtektites analyzed here for comparison. These inclusions are interpreted as partially digested, lechatelierite-like inclusions typically found in tektites and microtektites. The major and trace element (by laser ablation - inductively coupled plasma - mass spectrometry) abundance pattern of the Transantarctic Mountain microtektites matches the average upper continental crust composition for most elements. Major deviations include a strong to moderate depletion in volatile elements including Pb, Zn, Na, K, Rb, Sr and Cs, as a likely result of severe volatile loss during the high temperature melting and vaporization of crustal target rocks. The normal and high-Mg Transantarctic Mountain microtektites have compositions similar to the most volatile-poor normal and high-Mg Australasian microtektites reported in the literature. Their very low H2O and B contents (by secondary ion mass spectrometry) of 85 ± 58 (1 sigma) ?g/g and 0.53 ± 0.21 ?g/g, respectively, evidence the extreme volatile loss characteristically observed in tektites. The Sr and Nd isotopic compositions of multigrain samples of Transantarctic Mountain microtektites are 87Sr/86Sr ~ 0.71629 and 143Nd/144Nd ~ 0.51209, and fall into the Australasian tektite compositional field. The Nd model age calculated with respect to the chondritic uniform reservoir (CHUR) is TNdCHUR ~ 1.1 Ga, indicating a Meso-Proterozoic crustal source rock, as was derived for Australasian tektites as well. Coupled with the Quaternary age from the literature, the extended dataset presented in this work strengthens our previous conclusion that Transantarctic Mountain microtektites represent a major southward extension of the Australasian tektite/microtektite strewn field. Furthermore, the significant depletion in volatile elements (i.e., Pb, B, Na, K, Zn, Rb, Sr and Cs) of both normal and high-Mg Transantarctic Mountain microtektites relative to the Australasian ones provide us with further confirmation of a possible relationship between high temperature-time regimes in the microtektite-forming process and ejection distance.

Combining Sr and Na additions in hypoeutectic Al-Si foundry alloys

From recent published data, it is still unclear whether combining additions of Na and Sr have synergistic effects or deleterious interactions, This paper clarifies the interactions between these two modifiers and investigates the effects of such interactions on alloy solidification and castability. It was found that combined additions of Sr and Na do not appear to cause improvement of the modification of the eutectic microstructure even after only a short period after addition. Na addition may promote Sr vaporization and/or oxidation kinetically. leading to a quicker loss of both modifiers, which is blamed for the rapid loss of the modification effect during melt holding. Quenching trials during the eutectic arrest indicate that addition of Sr into Na-modified melts does not alter the eutectic solidification behaviour The effect of Na on eutectic solidification dominates, and the eutectic is observed to evolve with a significant dependency on the thermal gradient, Combining Sr and Na additions produced no beneficial effects on porosity and casting defects. (c) 2005 Elsevier B.V. All rights reserved.

Multi-objective optimization of the molecular structure of refrigerants

The aim of this work was to develop a new methodology, which can be used to design new refrigerants that are better than the currently used refrigerants. The methodology draws some parallels with the general approach of computer aided molecular design. However, the mathematical way of representing the molecular structure of an organic compound and the use of meta models during the optimization process make it different. In essence, this approach aimed to generate molecules that conform to various property requirements that are known and specified a priori. A modified way of mathematically representing the molecular structure of an organic compound having up to four carbon atoms, along with atoms of other elements such as hydrogen, oxygen, fluorine, chlorine and bromine, was developed. The normal boiling temperature, enthalpy of vaporization, vapor pressure, tropospheric lifetime and biodegradability of 295 different organic compounds, were collected from open literature and data bases or estimated. Surrogate models linking the previously mentioned quantities with the molecular structure were developed. Constraints ensuring the generation of structurally feasible molecules were formulated and used in commercially available optimization algorithms to generate molecular structures of promising new refrigerants. This study was intended to serve as a proof-of-concept of designing refrigerants using the newly developed methodology.

(Table 3, page 8) Chemical composition and physical cahracteristics of manganese nodules from the Pacific Ocean

This paper reviews the state of the art in processing and extraction of ocean floor manganese nodules. It briefly reviews the mining sites where the abundant rich nodules occur and also discusses the metal distribution in nodules in view of economical processing and extraction of these metal values. The paper discloses in a detailed manner the physical and chemical characteristics of nodules, including porosity, surface area, water content and the effect of temperature on crystal structure of major constituents of nodules. In the extraction aspect of nodules, the paper reviews two different extraction schemes revealed in the literature, namely hydrometallurgical treatment and pyrometallurgical treatment. The hydrometallurgical treatments include acid leaching, ammonia leaching, leaching with reducing agents and leaching after high temperature pre-treatments such as in sulfating rousting, while the pyrometallurgical processes include smelting, chlorination-vaporization and segregation. The paper also covers metal recovery processes from leach liquor. An economic survey of processing nodules has been made in terms of problems associated with metal-marketing, and impact of metal production from nodules on mineral industries.

Estudo da gasificação do glicerol em reator de leito fixo

As metas da União Europeia para 2020 em termos de biocombustíveis e biolíquidos traduziram-se, na última década, num destaque da indústria de biodiesel em Portugal. Inerente ao processo de produção biodiesel está um subproduto, o glicerol bruto, cujo estudo tem vindo a ser alvo de interesse na comunidade científica. O objetivo principal deste trabalho consistiu no estudo da gasificação do glicerol técnico e do glicerol bruto, usando vapor como agente oxidante. Pretendeu-se avaliar a composição do gás de produção obtido e os parâmetros de gasificação, como a percentagem de conversão de carbono e de hidrogénio, o rendimento de gás seco, a eficiência de gás frio e o poder calorífico do gás produzido. No estudo da gasificação do glicerol técnico avaliou-se o efeito da temperatura na performance do processo, entre 750 – 1000 ºC, e estudou-se ainda o efeito do caudal de alimentação ao reator (3,8 mL/min, 6,5 mL/min e 10,0 mL/min). Para o caudal mais baixo, estudou-se o efeito da razão de mistura glicerol/água (25/75, 40/60, 60/40 e 75/25) e para a razão de mistura 60/40 foi avaliada a influência da adição de ar como agente gasificante. O estudo da gasificação do glicerol bruto foi feito realizando ensaios de gasificação numa gama de temperaturas de 750 ºC a 1000 ºC, para uma razão de mistura glicerol/água (60/40) com o caudal de 3,8 mL/min e usando apenas vapor de água como agente de gasificação. Os ensaios foram realizados num reator de leito fixo de 500 mm de comprimento e 90 mm de diâmetro interno, composto por um leito de alumina com partículas de 5 mm de diâmetro. O aquecimento foi realizado com um forno elétrico de 4 kW. A amostra de gás de produção recolhida foi analisada por cromatografia gasosa com detector de termocondutividade. Os resultados obtidos na gasificação do glicerol técnico, revelaram que a temperatura é uma variável preponderante no desempenho do processo de gasificação. À exceção do poder calorífico superior, para o qual se obteve uma ligeira diminuição de valores com o aumento da temperatura, os valores mais elevados dos parâmetros de gasificação foram obtidos para temperaturas superiores a 900 ºC. Esta temperatura parece ser determinante no modelo cinético de gasificação do glicerol, condicionando a composição do gás de produção obtido. Concluiu-se ainda que, na gama de caudais testada, o caudal de alimentação ao reator não teve influência no processo de gasificação. Os ensaios realizados para avaliar o efeito da razão de mistura permitiram verificar que, o aumento da adição de água à alimentação se traduz na redução do teor de CO e de CH4 e no aumento do teor de H2 e CO2, no gás de produção. Para a razão de mistura 25/75 foram obtidos valores de 1,3 para o rácio H2/CO para temperaturas superiores a 900 ºC. A influência da adição de água tornou-se mais evidente nos ensaios de gasificação realizados a temperaturas superiores a 900 ºC. Verificou-se um aumento da conversão de carbono, do rendimento de gás seco e da eficiência do gás frio e uma ligeira diminuição do poder calorífico e da potência disponível, no gás de produção. Para as razões de misturas 60/40 e 40/60 obtiveram-se resultados, para os parâmetros de gasificação, da mesma ordem de grandeza e com valores intermédios entre os obtidos para as razões de mistura 25/75 e 75/25. Porém, quanto maior o teor de água alimentado maior o consumo de energia associado à vaporização da água. Assim, o aumento do teor de água na mistura só apresentará interesse industrial se o objetivo passar pela produção de hidrogénio. Quanto ao efeito da adição de ar como agente de gasificação, os resultados obtidos dão indicação que se poderão potenciar algumas reações exotérmicas que contribuirão para a redução do consumo energético global do processo. Por outro lado, o gás de produção apresentou um rácio H2/CO interessante do ponto de vista da sua aplicação industrial, superior em 35 % ao verificado para a gasificação efetuada apenas na presença de vapor. À exceção do decréscimo no valor do poder calorífico superior do gás de produção, os restantes parâmetros estudados apresentaram a mesma ordem de grandeza, dos obtidos para o estudo da mesma razão de mistura na ausência de ar. Relativamente ao estudo da gasificação do glicerol bruto, obtiveram-se valores de rácio H2/CO e eficiência de gás frio mais elevados que os valores obtidos para a mesma razão de mistura usando glicerol técnico. Os demais parâmetros de gasificação avaliados mostraram-se semelhantes entre as duas matérias-primas, verificando-se apenas uma ligeira diminuição no valor do poder calorífico superior do gás produzido com glicerol bruto. Os resultados obtidos demonstram a possibilidade de valorização energética do glicerol bruto resultante da produção de biodiesel.

Diseño del proceso de una torre de vacío. Ventajas de la simulación

Los hidrocarburos pesados son el mayor recurso del petróleo en el mundo, sin embargo en el pasado se habían dejado de lado como recurso energético debido a las dificultades y costos asociados de su producción [1]. La industria financia estas investigaciones por la importancia del tema en producción y caracterización. Al trabajar con una torre de vacio los datos necesarios para los cálculos son las temperaturas ASTM (10mmHg) y la densidad del crudo con la cual se obtiene la curva TBP760 (True Boiling Point), también se necesita las especificaciones de los productos y los rendimientos respecto de la alimentación. Para poder correlacionar los distintos puntos de ebullición con los porcentajes de vaporizado para cada cambio de presión de los distintos productos, se construye un diagrama de fases con las temperaturas EFV760 (Equilibrium Flash Vaporization) y EFV10. El simulador a través de cálculos internos resuelve automáticamente el diagrama de fases, en comparación con la dificultad que representan los cálculos manuales del mismo, tal como son explicitados precedentemente. En este trabajo se desarrolla la simulación de una torre de vacío mediante el simulador Aspen HYSYS V8.3, empleando como alimentación un crudo pesado. Lo antes expuesto constituye una importante ventaja el uso del simulador frente al cálculo convencional, considerando los tiempos de resolución de los diseños de procesos.

Vapor pressure and liquid density of fluorinated alcohols:Experimental, simulation and GC-SAFT-VR predictions

The vapor pressure of four liquid 1H,1H-perfluoroalcohols (CF3(CF2)n(CH2)OH, n ¼ 1, 2, 3, 4), often called odd-fluorotelomer alcohols, was measured as a function of temperature between 278 K and 328 K. Liquid densities were also measured for a temperature range between 278 K and 353 K. Molar enthalpies of vaporization were calculated from the experimental data. The results are compared with data from the literature for other perfluoroalcohols as well as with the equivalent hydrogenated alcohols. The results were modeled and interpreted using molecular dynamics simulations and the GC-SAFT-VR equation of state.