Formulação teórica dos fundamentos da otimização global topográfica com análise de desempenho e aplicações à estabilidade de fases de misturas termodinâmicas

Métodos de otimização que utilizam condições de otimalidade de primeira e/ou segunda ordem são conhecidos por serem eficientes. Comumente, esses métodos iterativos são desenvolvidos e analisados à luz da análise matemática do espaço euclidiano n-dimensional, cuja natureza é de caráter local. Consequentemente, esses métodos levam a algoritmos iterativos que executam apenas as buscas locais. Assim, a aplicação de tais algoritmos para o cálculo de minimizadores globais de uma função não linear,especialmente não-convexas e multimodais, depende fortemente da localização dos pontos de partida. O método de Otimização Global Topográfico é um algoritmo de agrupamento, que utiliza uma abordagem baseada em conceitos elementares da teoria dos grafos, a fim de gerar bons pontos de partida para os métodos de busca local, a partir de pontos distribuídos de modo uniforme no interior da região viável. Este trabalho tem dois objetivos. O primeiro é realizar uma nova abordagem sobre método de Otimização Global Topográfica, onde, pela primeira vez, seus fundamentos são formalmente descritos e suas propriedades básicas são matematicamente comprovadas. Neste contexto, propõe-se uma fórmula semi-empírica para calcular o parâmetro chave deste algoritmo de agrupamento, e, usando um método robusto e eficiente de direções viáveis por pontos-interiores, estendemos o uso do método de Otimização Global Topográfica a problemas com restrições de desigualdade. O segundo objetivo é a aplicação deste método para a análise de estabilidade de fase em misturas termodinâmicas,o qual consiste em determinar se uma dada mistura se apresenta em uma ou mais fases. A solução deste problema de otimização global é necessária para o cálculo do equilíbrio de fases, que é um problema de grande importância em processos da engenharia, como, por exemplo, na separação por destilação, em processos de extração e simulação da recuperação terciária de petróleo, entre outros. Além disso, afim de ter uma avaliação inicial do potencial dessa técnica, primeiro vamos resolver 70 problemas testes, e então comparar o desempenho do método proposto aqui com o solver MIDACO, um poderoso software recentemente introduzido no campo da otimização global.

Determinação experimental e modelagem termodinâmica do equilíbrio de fases em sistemas com CO2, n-hexano, n-hexadecano e tetralina

A descoberta de petróleo na camada de Pré-Sal possibilita a geração de ganhos em relação à dependência energética do país, mas também grandes desafios econômicos e tecnológicos. Os custos de extração são maiores devido a vários fatores como a exigência de equipamentos de exploração que suportem elevadas pressões, altas temperaturas e grandes concentrações de gases ácidos, tais quais, dióxido de carbono (CO2) e sulfeto de hidrogênio (H2S). Uma das principais preocupações com o CO2 é evitar liberá-lo para a atmosfera durante a produção. Com a modelagem termodinâmica de dados de equilíbrio de sistemas envolvendo CO2 supercrítico e hidrocarbonetos é possível projetar equipamentos utilizados em processos de separação. A principal motivação do trabalho é o levantamento de dados de equilíbrio de fases de sistemas compostos de CO2 e hidrocarbonetos, possibilitando assim prever o comportamento dessas misturas. Os objetivos específicos são a avaliação do procedimento experimental, a estimação e predição dos parâmetros de interação binários para assim prever o comportamento de fases dos sistemas ternários envolvendo CO2 e hidrocarbonetos. Duas metodologias foram utilizadas para obtenção dos dados de equilíbrio: método estático sintético (visual) e método dinâmico analítico (recirculação das fases). Os sistemas avaliados foram: CO2 + n-hexano, CO2 + tetralina, CO2 + n-hexadecano, CO2 + n-hexano + tetralina e CO2 + tetralina + n-hexadecano à alta pressão; tetralina + n-hexadecano à baixa pressão. Para o tratamento dos dados foi utilizada equação de estado cúbica de Peng-Robinson e a regra de mistura clássica

Effect of growth conditions on the GaN thin film by sputtering deposition

The phase transition between thermodynamically stable hexagonal wurtzite (h-WZ) gallium nitride (GaN) and metastable cubic zinc-blende (c-ZB) GaN during growth by radio-frequency planar magnetron sputtering is studied. GaN films grown on substrates with lower mismatches tend to have a h-WZ structure, but when grown on substrates with higher mismatches, a c-ZB structure is preferred. GaN films grown under high nitrogen pressure also tend to have a h-WZ structure, whereas a c-ZB structure is preferred when grown under low nitrogen pressure. In addition, low target-power growth not only helps to improve hexagonal GaN (h-GaN) crystalline quality at high nitrogen pressure on low-mismatch substrates, but also enhances cubic GaN (c-GaN) quality at low nitrogen pressure on high-mismatch substrates. (c) 2007 Elsevier B.V. All rights reserved.

Investigation of Interfacial Tensions for Carbon Dioxide Aqueous Solutions by Perturbed-Chain Statistical Associating Fluid Theory Combined with Density-Gradient Theory

The perturbed-chain statistical associating fluid theory and density-gradient theory are used to construct an equation of state (EOS) applicable for the phase behaviors of carbon dioxide aqueous solutions. With the molecular parameters and influence parameters respectively regressed from bulk properties and surface tensions of pure fluids as input, both the bulk and interfacial properties of carbon dioxide aqueous solutions are satisfactorily correlated by adjusting the binary interaction parameter (k(ij)). Our results show that the constructed EOS is able to describe the interfacial properties of carbon dioxide aqueous solutions in a wide temperature range, and illustrate the influences of temperature, pressure, and densities in each phase on the interfacial properties.

Hydrate dissociation conditions for gas mixtures containing carbon dioxide, hydrogen, hydrogen sulfide, nitrogen, and hydrocarbons using SAFT

A new method, a molecular thermodynamic model based on statistical mechanics, is employed to predict the hydrate dissociation conditions for binary gas mixtures with carbon dioxide, hydrogen, hydrogen sulfide, nitrogen, and hydrocarbons in the presence of aqueous solutions. The statistical associating fluid theory (SAFT) equation of state is employed to characterize the vapor and liquid phases and the statistical model of van der Waals and Platteeuw for the hydrate phase. The predictions of the proposed model were found to be in satisfactory to excellent agreement with the experimental data.

Separation of ethyl acetate-ethanol azeotropic mixture using hydrophilic ionic liquids

The separation of ethyl acetate and ethanol (EtOH) is important but difficult due to their close boiling points and formation of an azeotropic mixture. The separation of the azeotropic mixture of ethyl acetate and EtOH using the hydrophilic ionic liquids (ILs) 1-alkyl-3-methylimidazolium chloride (alkyl = butyl, hexyl, and octyl) ([C(n)mim]Cl, n = 4, 6, 8) and 1-allyl-3-methylimidazolium chloride and bromide ([Amim]Cl and [Amim]Br) has been investigated. Triangle phase diagrams of five ILs with ethyl acetate and EtOH were constructed, and the biphasic regions were found as follows: [Amim]Cl > [Amim]Br > [C(4)mim]Cl > [C(6)mim]Cl > [C(8)mim]Cl. The mechanisms of the ILs including cation, anion, and polarity effect were discussed.

Calculating the equation of state parameters and predicting the spinodal curve of isotactic polypropylene/poly(ethylene-co-octene) blend by molecular dynamics simulations combined with Sanchez-Lacombe lattice fluid theory

Molecular dynamics simulations are adopted to calculate the equation of state characteristic parameters P*, rho*, and T* of isotactic polypropylene (iPP) and poly(ethylene-co-octene) (PEOC), which can be further used in the Sanchez-Lacombe lattice fluid theory (SLLFT) to describe the respective physical properties. The calculated T* is a function of the temperature, which was also found in the literature. To solve this problem, we propose a Boltzmann fitting of the data and obtain T* at the high-temperature limit. With these characteristic parameters, the pressure-volume-temperature (PVT) data of iPP and PEOC are predicted by the SLLFT equation of state. To justify the correctness of our results, we also obtain the PVT data for iPP and PEOC by experiments. Good agreement is found between the two sets of data. By integrating the Euler-Lagrange equation and the Cahn-Hilliard relation, we predict the density profiles and the surface tensions for iPP and PEOC, respectively. Furthermore, a recursive method is proposed to obtain the characteristic interaction energy parameter between iPP and PEOC. This method, which does not require fitting to the experimental phase equilibrium data, suggests an alternative way to predict the phase diagrams that are not easily obtained in experiments.

Compositional dependence of the microphase separation transition in styrene-butadiene-styrene triblock copolymer/polystyrene blends

The change in the microphase separation transition (MST) temperature of a styrene-butadiene-styrene (SBS) triblock copolymer induced by the addition of polystyrene (PS) was investigated by small-angle X-ray scattering. It was found that the transition temperature was determined from the molecular weight (M(H)) Of the added PS in relation to that of the corresponding blocks (M(A)) in the copolymer. The MST temperature decreased with added PS if M(H)/M(A) < 1/4, while it increased with added PS when M(H)/M(A) > 1/4 Analysis of the theoretical expression based on the random phase approximation showed exactly the same tendency of change in the transition temperatures as that observed experimentally. The interaction parameter, chi(SB), obtained by nonlinear fitting of the scattering profiles of SBS/PS blends in the disordered state, was found to be a function of temperature and composition. Composition fluctuations were found to exist in SBS/PS blends, increasing with increasing addition of PS but diminishing with increasing molecular weight of the added PS.

COMPATIBILITY AND CRYSTALLIZATION OF TETRAHYDROFURAN-METHYL METHACRYLATE DIBLOCK COPOLYMER/POLYTETRAHYDROFURAN BLENDS

The compatibility and crystallization of tetrahydrofuran-methyl methacrylate diblock copolymer (PTHF-b-PMMA)/tetrahydrofuran homopolymer (PTHF) blends were studied. Our results showed that the crystallization and morphology of compatible PTHF-b-PMMA/PTHF

水-盐-气天然流体的物理化学性质模拟与流体包裹体应用研究

Natural fluids with water-salt-gas are often found in every sphere of the Earth, whose physicochemical properties and geochemical behaviors are complicated. To study these properties and behaviors turns out to be one of the challenging issues in geosciences. Traditional approaches mainly depend on experiments and observations. However, it is impossible to obtain a large number of data covering a large T-P space of the Earth by experimental methods in the near future, which will hinder the advance of the theoretical study. Therefore, it is important to model natural fluids by advanced theoretical methods, by which limited experimental data can be extended to a large temperature-pressure-composition space. Physicochemical models developed in this dissertation are not only more accurate, but also extend the applied T-P-m region of the experimental data of the multi-fluid systems by about two times. These models provide the new and accurate theoretical tools for the geochemical research, especially for the water-rock interactions and the study of the fluid inclusions. The main achievements can be summarized as follows: (1) A solubility model on components of natural gases is presented. The solubility model on the systems of CH4-H2O-NaCl, C2H6-H2O-NaCl or N2-H2O-NaCl takes advantage of modern physicochemical theory and methods, and is an improvement over previous models whose prediction and precision are relatively poor. The model can predict not only the gas solubility in liquid phase but also water content in the gas phase. In addition, it can predict gases (methane or nitrogen) solubility in seawater and brine. Isochores can be determined, which are very important in the interpretation of fluid inclusions. (2) A density model on common aqueous salt solutions is developed. The density models with high precision for common aqueous salt solutions (H2O-NaCl, H2O-LiCl, H2O-KCl, H2O-MgCl2, H2O-CaCl2, H2O-SrCl2 or H2O-BaCl2) are absent in the past. Previous density models are limited to the relatively small range of experimental data, and cannot meet the requirement of the study of natural fluids. So a general density model of the above systems is presented by us based on the international standard density model of the water. The model exceeds the other models in both precision and prediction. (3) A viscosity model on common aqueous alkali-chloride solutions is proposed. Dynamic viscosity of water-salt systems, an important physics variable, is widely used in three-dimension simulation of the fluids. But in most cases, due to the lack of viscosity models with a wide T-P range, the viscosity of aqueous salt solutions is replaced by that of the water, giving rise to a relatively large uncertainty. A viscosity model with good prediction for the systems (H2O-NaCl, H2O-LiCl or H2O-KCl) is presented on the base of the international standard viscosity model of water and the density model developed before. (4) Equation of State applied in fluid inclusions. The best Equations of State in the world developed by others or us recently are applied in the study of the fluid inclusions. Phase equilibria and isochores of unitary system (e.g. H2O, CO2, CH4, O2, N2, C2H6 or H2S), binary H2O-NaCl system and ternary H2O-CH4-NaCl system are finished. From these programs and thermodynamic equations of coexisting ores, the physicochemical conditions before or after the deposits form can be determined. To some extent, it is a better tool.

分子水平上研究地质流体的物理化学性质

Geological fluids exist in every geosphere of the Earth and play important roles in many processes of material transformations, energetic interchanges and geochemical interactions. To study the physicochemical properties and geochemical behaviors of geological fluids turn Girt to be one of the challenging issues in geosciences. Compared with conventional approaches of experiments and semi-theoretical modeling, computer simulation on molecular level shows its advantages on quantitative predictions of the physicochemical properties of geological fluids under extreme conditions and emerges as a promising approach to find the characteristics of geological fluids and their interactions in different geospheres of the Earth interior.This dissertation systematically discusses the physicochemical properties of typical geological fluids with state-of-the-art computer simulation techniques. The main results can be summarized as follows: (1) The experimental phase behaviors of the systems CH4-C2H6 and. CO2 have been successfully reproduced with Monte Carlo simulations. (2) Through comprehensive isothermal-isobaric molecular dynamics simulations, the PVT data of water hia^e been extended beyond experimental range to about 2000 K and 20 GPa and an improved equation of state for water has been established. (3) Based on extensive computer simulations, am optimized molecular potential for carbon dioxide have been proposed, this model is expected to predict different properties of carbon dioxide (volumetric properties, phase equilibria, heat of vaporization, structural and dynamical properties) with improved accuracies. (4) On the basis of the above researches of the end-members, a set of parameters for unlike interactions has been proposed by non-linear fitting to the ab initio potential surface of CO2-H2O and is superior to the common used mixing rule and the results of prior workers vs/Ith remarkable accuracies, then a number of simulations of the mixture have been carried out to generate data under high temperatures and pressures as an important complement to the limited experiments. (5) With molecular dynamics simulations, various structural, dynamical and thermodynamical properties of ionic solvations and associations have been oomprehensively analyzed, these results not only agree well with experimental data and first principle calculation results, but also reveal some new insights into the microscopic ionic solvation and association processes.

Speeds of sound and isentropic compressibilities of binary mixtures containing trialkylamines with alkanes and mono-alkylamines at 303.15 and 313.15 K

Isentropic compressibilities ?S, and excess isentropic compressibilities ?SE have been determined from measurements of speeds of sound u and densities ? of 14 binary mixtures of triethylamine (TEA) and tri-n-butylamine (TBA) with n-hexane, n-octane, iso-octane, n-propylamine, n-butylamine, n-hexylamine and n-octylamine. The relative magnitude and sign of ?SE have been interpreted in terms of molecular interactions and interstitial accommodation. The values of ?SE for TEA + alkane are positive while for TBA + alkane are negative. The values of ?SE for TEA + primary amine become progressively less positive and eventually to negative with the increase in chain length of alkylamine. In case of TBA + primary amine, the values of ?SE increase from n-propylamine to n-butylamine, and then decrease with chain length of primary amine. The experimental speeds of sound u have been analyzed in terms of collision factor theory, free length theory and Prigogine–Flory–Patterson statistical theory of solutions.

Acoustic, volumetric, compressibility and refractivity properties and reduction parameters for the ERAS and Flory models of some homologous series of amines from 298.15 to 328.15 K

The speeds of sound u, densities ? and refractive indices nD of homologous series of mono-, di-, and tri-alkylamines were measured in the temperature range from 298.15 to 328.15 K. Isentropic and isothermal compressibilities ?S and ?T, molar refraction Rm, Eykman’s constant Cm, Rao’s molar sound function R, thermal expansion coefficient a, thermal pressure coefficient ?, and reduction parameters P*, V*, and T* in frameworks of the ERAS model for associated amines and Flory model for tertiary amines have been calculated from the measured experimental data. Applicability of the Rao theory and the ERAS and Flory models have been examined and discussed for the alkyl amines.