An unstructured CVFEM and moving interface algorithm for non-Newtonian Hele-Shaw flows in injection molding

Purpose - The purpose of this paper is to develop a novel unstructured simulation approach for injection molding processes described by the Hele-Shaw model. Design/methodology/approach - The scheme involves dual dynamic meshes with active and inactive cells determined from an initial background pointset. The quasi-static pressure solution in each timestep for this evolving unstructured mesh system is approximated using a control volume finite element method formulation coupled to a corresponding modified volume of fluid method. The flow is considered to be isothermal and non-Newtonian. Findings - Supporting numerical tests and performance studies for polystyrene described by Carreau, Cross, Ellis and Power-law fluid models are conducted. Results for the present method are shown to be comparable to those from other methods for both Newtonian fluid and polystyrene fluid injected in different mold geometries. Research limitations/implications - With respect to the methodology, the background pointset infers a mesh that is dynamically reconstructed here, and there are a number of efficiency issues and improvements that would be relevant to industrial applications. For instance, one can use the pointset to construct special bases and invoke a so-called ""meshless"" scheme using the basis. This would require some interesting strategies to deal with the dynamic point enrichment of the moving front that could benefit from the present front treatment strategy. There are also issues related to mass conservation and fill-time errors that might be addressed by introducing suitable projections. The general question of ""rate of convergence"" of these schemes requires analysis. Numerical results here suggest first-order accuracy and are consistent with the approximations made, but theoretical results are not available yet for these methods. Originality/value - This novel unstructured simulation approach involves dual meshes with active and inactive cells determined from an initial background pointset: local active dual patches are constructed ""on-the-fly"" for each ""active point"" to form a dynamic virtual mesh of active elements that evolves with the moving interface.

Biocompatible in-tube solid phase microextraction coupled with liquid chromatography-fluorescence detection for determination of interferon alpha in plasma samples

The present work demonstrates the successful application of automated biocompatible in-tube solid-phase microextraction coupled with liquid chromatography (in-tube SPME/LC) for determination of interferon alpha(2a) (IFN alpha(2a)) in plasma samples for therapeutic drug monitoring. A restricted access material (RAM, protein-coated silica) was employed for preparation of a lab-made biocompatible in-tube SPME capillary that enables the direct injection of biological fluids as well as the simultaneous exclusion of macromolecules by chemical diffusion barrier and drug pre-concentration. The in-tube SPME variables, such as sample volume, draw/eject volume, number of draw-eject cycles, and desorption mode were optimized, to improve the sensitivity of the proposed method. The IFN alpha(2a) analyses in plasma sample were carried out within 25 min (sample preparation and LC analyses). The response of the proposed method was linear over a dynamic range, from 0.06 to 3.0 MIU mL(-1), with correlation coefficient equal to 0.998. The interday precision of the method presented coefficient of variation lower than 8%. The proposed automated method has adequate analytical sensitivity and selectivity for determination of IFN alpha(2a) in plasma samples for therapeutic drug monitoring. (C) 2010 Elsevier B.V. All rights reserved.

Análise de modelos geológicos utilizando percolação dinâmica

In the present study we elaborated algorithms by using concepts from percolation theory which analyze the connectivity conditions in geological models of petroleum reservoirs. From the petrophysical parameters such as permeability, porosity, transmittivity and others, which may be generated by any statistical process, it is possible to determine the portion of the model with more connected cells, what the interconnected wells are, and the critical path between injector and source wells. This allows to classify the reservoir according to the modeled petrophysical parameters. This also make it possible to determine the percentage of the reservoir to which each well is connected. Generally, the connected regions and the respective minima and/or maxima in the occurrence of the petrophysical parameters studied constitute a good manner to characterize a reservoir volumetrically. Therefore, the algorithms allow to optimize the positioning of wells, offering a preview of the general conditions of the given model s connectivity. The intent is not to evaluate geological models, but to show how to interpret the deposits, how their petrophysical characteristics are spatially distributed, and how the connections between the several parts of the system are resolved, showing their critical paths and backbones. The execution of these algorithms allows us to know the properties of the model s connectivity before the work on reservoir flux simulation is started

Estudo comparativo da utilização do nitrogênio como fluido alternativo após a injeção de vapor

Currently, due to part of world is focalized to petroleum, many researches with this theme have been advanced to make possible the production into reservoirs which were classified as unviable. Because of geological and operational challenges presented to oil recovery, more and more efficient methods which are economically successful have been searched. In this background, steam flood is in evidence mainly when it is combined with other procedures to purpose low costs and high recovery factors. This work utilized nitrogen as an alternative fluid after steam flood to adjust the best combination of alternation between these fluids in terms of time and rate injection. To describe the simplified economic profile, many analysis based on liquid cumulative production were performed. The completion interval and injection fluid rates were fixed and the oil viscosity was ranged at 300 cP, 1.000 cP and 3.000 cP. The results defined, for each viscosity, one specific model indicating the best period to stop the introduction of steam and insertion of nitrogen, when the first injected fluid reached its economic limit. Simulations in physics model defined from one-eighth nine-spot inverted were realized using the commercial simulator Steam, Thermal and Advanced Processes Reservoir Simulator STARS of Computer Modelling Group CMG

Análise paramétrica do método de injeção alternada de água e CO2(WAG) em reservatórios de petróleo

After the decline of production from natural energy of the reservoir, the methods of enhanced oil recovery, which methods result from the application of special processes such as chemical injection, miscible gases, thermal and others can be applied. The advanced recovery method with alternating - CO2 injection WAG uses the injection of water and gas, normally miscible that will come in contact with the stock oil. In Brazil with the discovery of pre-salt layer that gas gained prominence. The amount of CO2 present in the oil produced in the pre-salt layer, as well as some reservoirs is one of the challenges to be overcome in relation to sustainable production once this gas needs to be processed in some way. Many targets for CO2 are proposed by researchers to describe some alternatives to the use of CO2 gas produced such as enhanced recovery, storage depleted fields, salt caverns storage and marketing of CO2 even in plants. The largest oil discoveries in Brazil have recently been made by Petrobras in the pre -salt layer located between the states of Santa Catarina and Espírito Santo, where he met large volumes of light oil with a density of approximately 28 ° API, low acidity and low sulfur content. This oil that has a large amount of dissolved CO2 and thus a pioneering solution for the fate of this gas comes with an advanced recovery. The objective of this research is to analyze which parameters had the greatest influence on the enhanced recovery process. The simulations were performed using the "GEM" module of the Computer Modelling Group, with the aim of studying the advanced recovery method in question. For this work, semi - synthetic models were used with reservoir and fluid data that can be extrapolated to practical situations in the Brazilian Northeast. The results showed the influence of the alternating injection of water and gas on the recovery factor and flow rate of oil production process, when compared to primary recovery and continuous water injection or continuous gas injection

Simulador computacional do comportamento em regime de poços de petróleo equipados com gas-lift contínuo

Amongst the results of the AutPoc Project - Automation of Wells, established between UFRN and Petrobras with the support of the CNPq, FINEP, CTPETRO, FUNPEC, was developed a simulator for equipped wells of oil with the method of rise for continuous gas-lift. The gas-lift is a method of rise sufficiently used in production offshore (sea production), and its basic concept is to inject gas in the deep one of the producing well of oil transform it less dense in order to facilitate its displacement since the reservoir until the surface. Based in the use of tables and equations that condense the biggest number of information on characteristics of the reservoir, the well and the valves of gas injection, it is allowed, through successive interpolations, to simulate representative curves of the physical behavior of the existing characteristic variable. With a simulator that approaches a computer of real the physical conditions of an oil well is possible to analyze peculiar behaviors with very bigger speeds, since the constants of time of the system in question well are raised e, moreover, to optimize costs with assays in field. The simulator presents great versatility, with prominance the analysis of the influence of parameters, as the static pressure, relation gas-liquid, pressure in the head of the well, BSW (Relation Basic Sediments and Water) in curves of request in deep of the well and the attainment of the curve of performance of the well where it can be simulated rules of control and otimization. In moving the rules of control, the simulator allows the use in two ways of simulation: the application of the control saw software simulated enclosed in the proper simulator, as well as the use of external controllers. This implies that the simulator can be used as tool of validation of control algorithms. Through the potentialities above cited, of course one another powerful application for the simulator appears: the didactic use of the tool. It will be possible to use it in formation courses and recycling of engineers

Avaliação do impacto do reservatório de Três Irmãos sobre a superfície potenciométrica do aquífero livre na cidade de Pereira Barreto (SP): uma abordagem numérica e geoestatística

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Estimativas da condutividade térmica dos minerais e rochas e influência de parâmetros térmicos e petrofísicos na resistividade aparente da formação

O presente estudo realiza estimativas da condutividade térmica dos principais minerais formadores de rochas, bem como estimativas da condutividade média da fase sólida de cinco litologias básicas (arenitos, calcários, dolomitos, anidritas e litologias argilosas). Alguns modelos térmicos foram comparados entre si, possibilitando a verificação daquele mais apropriado para representar o agregado de minerais e fluidos que compõem as rochas. Os resultados obtidos podem ser aplicados a modelamentos térmicos os mais variados. A metodologia empregada baseia-se em um algoritmo de regressão não-linear denominado de Busca Aleatória Controlada. O comportamento do algoritmo é avaliado para dados sintéticos antes de ser usado em dados reais. O modelo usado na regressão para obter a condutividade térmica dos minerais é o modelo geométrico médio. O método de regressão, usado em cada subconjunto litológico, forneceu os seguintes valores para a condutividade térmica média da fase sólida: arenitos 5,9 ± 1,33 W/mK, calcários 3.1 ± 0.12 W/mK, dolomitos 4.7 ± 0.56 W/mK, anidritas 6.3 ± 0.27 W/mK e para litologias argilosas 3.4 ± 0.48 W/mK. Na sequência, são fornecidas as bases para o estudo da difusão do calor em coordenadas cilíndricas, considerando o efeito de invasão do filtrado da lama na formação, através de uma adaptação da simulação de injeção de poços proveniente das teorias relativas à engenharia de reservatório. Com isto, estimam-se os erros relativos sobre a resistividade aparente assumindo como referência a temperatura original da formação. Nesta etapa do trabalho, faz-se uso do método de diferenças finitas para avaliar a distribuição de temperatura poço-formação. A simulação da invasão é realizada, em coordenadas cilíndricas, através da adaptação da equação de Buckley-Leverett em coordenadas cartesianas. Efeitos como o aparecimento do reboco de lama na parede do poço, gravidade e pressão capilar não são levados em consideração. A partir das distribuições de saturação e temperatura, obtém-se a distribuição radial de resistividade, a qual é convolvida com a resposta radial da ferramenta de indução (transmissor-receptor) resultando na resistividade aparente da formação. Admitindo como referência a temperatura original da formação, são obtidos os erros relativos da resistividade aparente. Através da variação de alguns parâmetros, verifica-se que a porosidade e a saturação original da formação podem ser responsáveis por enormes erros na obtenção da resistividade, principalmente se tais "leituras" forem realizadas logo após a perfuração (MWD). A diferença de temperatura entre poço e formação é a principal causadora de tais erros, indicando que em situações onde esta diferença de temperatura seja grande, perfilagens com ferramentas de indução devam ser realizadas de um a dois dias após a perfuração do poço.

Estudo Comparativo da Taxa de Injeção em Injetor do Tipo Venturi com e Sem Válvula de Retenção

Para se realizar a técnica da quimigação é necessário que o sistema de irrigação possua um injetor para incorporar os produtos na água de irrigação. Dentre os métodos existentes, o que tem se destacado nos últimos anos é o injetor do tipo Venturi. Essa preferência tem sido motivada pela sua simplicidade, por não possuir peças móveis e ser de baixo custo. Apesar dessas vantagens, o uso do injetor Venturi, assim como os demais injetores, necessita de dispositivos de segurança que minimizem o risco de contaminação ambiental. Neste trabalho, objetivou-se analisar a influência da válvula de retenção na taxa de injeção de água de um injetor do tipo Venturi instalado em pivô central e seu rendimento. O dispositivo de injeção era constituído de um injetor Venturi de 1,5 polegadas associado a uma motobomba centrífuga de 2,2 kW. O equipamento foi instalado na base do pivô central. Foram tomados tempos de esvaziamento de um reservatório de volume de 20 litros. O delineamento experimental utilizado foi o inteiramente casualizado, sendo realizadas dez repetições em cada situação avaliada, ou seja, injetor equipado com válvula de retenção e sem a presença da mesma. Aplicou-se a análise de variância, seguida do teste de Tukey. Concluiu-se que a presença de válvula de retenção em injetores do tipo Venturi não influenciou de forma significativa as taxas de injeção e o rendimento, tendo este último apresentado valores da ordem de 10%.

Cross linked Poly (4-Vinylpyridine-Ethylene Glycol Dimethacrylate) Used for Preconcentration of Cd(II) and its Determination by Flow Injection Flame Atomic Absorption Spectrometry

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Hydroxocobalamin quantification in human plasma by high-performance liquid chromatography coupled with electrospray tandem mass spectrometry in a pharmacokinetic study

A rapid, sensitive and specific method for quantifying hydroxocobalamin in human plasma using paracetamol as the internal standard (IS) is described. The analyte and the IS were extracted from plasma by liquid-liquid extraction using an organic solvent (ethanol 100%; -20°C). The extracts were analyzed by high performance liquid chromatography coupled with electrospray tandem mass spectrometry (HPLC-MS-MS). Chromatography was performed on Prevail C8 3 μm, analytical column (2.1×100 mm i.d.). The method had a chromatographic run time of 3.4 min and a linear calibration curve over the range 5-400 ng.mL-1 (r>0.9983). The limit of quantification was 5 ng.mL-1. The method was also validated without the use of the internal standard. The precision in the intra-batch validation with IS was 9.6%, 8.9%, 1.0% and 2.8% whereas without IS was 9.2%, 8.2%, 1.8% and 1.5% for 5, 15, 80 and 320 ng/mL, respectively. The accuracy in intra-batch validation with IS was 108.9%, 99.9%, 98.9% and 99.0% whereas without IS was 101.1%, 99.3%, 97.5% and 92.5% for 5, 15, 80 and 320 ng/mL, respectively. The precision in the inter-batch validation with IS was 9.4%, 6.9%, 4.6% and 5.5% whereas without IS was 10.9%, 6.4%, 5.0% and 6.2% for 5, 15, 80 and 320 ng/mL, respectively. The accuracy in inter-batch validation with IS was 101.9%, 104.1%, 103.2% and 99.7% whereas without IS was 94.4%, 101.2%, 101.6% and 96.0% for 5, 15, 80 and 320 ng/mL, respectively. This HPLC-MS-MS procedure was used to assess the pharmacokinetics of Hydroxo cobalamin following intramuscular injection 5000 μg in healthy volunteers of both sexes (10 males and 10 females). The volunteers had the following clinical characteristics (according to gender and expressed as mean ± SD [range]): males: age: 32.40 ± 8.00 y [23.00-46.00], height: 1.73 ± 0.07 m [1.62-1.85], body weight: 72.48 ± 10.22 Kg [60.20- 88.00]; females: age: 28.60 ± 9.54 y [18.00-44.00], height: 1.60 ± 0.05 m [1.54-1.70], body weight: 58.64 ± 6.09 Kg [51.70- 66.70]. The following pharmacokinetic parameters were obtained from the hydroxocobalamin plasma concentration vs. time curves: AUClast, T1/2, Tmax, Vd, Cl, Cmax and Clast. The pharmacokinetic parameters were 120 (± 25) ng/mL for Cmax, 2044 (± 641) ng.h/mL for AUClast, 8 (± 3.2) ng.mL-1 for Clast, 38 (± 15.8) hr for T1/2 and 2.5 (range 1-6) hr for Tmax. Female volunteers presented significant (p=0.0136) lower AUC (1706 ± 704) ng.h/mL) and larger (p=0.0205) clearance (2.91 ± 1.41 L/hr), as compared to male 2383 ± 343 ng.h/mL and 1.76 ± 0.23 L/hr, respectively. These pharmacokinetic differences could explain the higher prevalence of vitamin B12 deficiency in female patients. The method described validated well without the use of the internal standard and this approach should be investigated in other HPLC-MS-MS methods.

Effect of resonant tunneling on exciton dynamics in coupled dot-well nanostructures

Excitonic dynamics in a hybrid dot-well system composed of InAs quantum dots (QDs) and an InGaAs quantum well (QW) is studied by means of femtosecond pump-probe reflection and continuous wave (cw) photoluminescence (PL) spectroscopy. The system is engineered to bring the QW ground exciton state into resonance with the third QD excited state. The resonant tunneling rate is varied by changing the effective barrier thickness between the QD and QW layers. This strongly affects the exciton dynamics in these hybrid structures as compared to isolated QW or QD systems. Optically measured decay times of the coupled system demonstrate dramatically different response to temperature change depending on the strength of the resonant tunneling or coupling strength. This reflects a competition between purely quantum mechanical and thermodynamical processes.

Nonlinear/electromagnetic co-simulation of microwaves and millimeter-waves links

The Ph.D. thesis describes the simulations of different microwave links from the transmitter to the receiver intermediate-frequency ports, by means of a rigorous circuit-level nonlinear analysis approach coupled with the electromagnetic characterization of the transmitter and receiver front ends. This includes a full electromagnetic computation of the radiated far field which is used to establish the connection between transmitter and receiver. Digitally modulated radio-frequency drive is treated by a modulation-oriented harmonic-balance method based on Krylov-subspace model-order reduction to allow the handling of large-size front ends. Different examples of links have been presented: an End-to-End link simulated by making use of an artificial neural network model; the latter allows a fast computation of the link itself when driven by long sequences of the order of millions of samples. In this way a meaningful evaluation of such link performance aspects as the bit error rate becomes possible at the circuit level. Subsequently, a work focused on the co-simulation an entire link including a realistic simulation of the radio channel has been presented. The channel has been characterized by means of a deterministic approach, such as Ray Tracing technique. Then, a 2x2 multiple-input multiple-output antenna link has been simulated; in this work near-field and far-field coupling between radiating elements, as well as the environment factors, has been rigorously taken into account. Finally, within the scope to simulate an entire ultra-wideband link, the transmitting side of an ultrawideband link has been designed, and an interesting Front-End co-design technique application has been setup.

Novel simulation methods for Coulomb and hydrodynamic interactions

This thesis presents new methods to simulate systems with hydrodynamic and electrostatic interactions. Part 1 is devoted to computer simulations of Brownian particles with hydrodynamic interactions. The main influence of the solvent on the dynamics of Brownian particles is that it mediates hydrodynamic interactions. In the method, this is simulated by numerical solution of the Navier--Stokes equation on a lattice. To this end, the Lattice--Boltzmann method is used, namely its D3Q19 version. This model is capable to simulate compressible flow. It gives us the advantage to treat dense systems, in particular away from thermal equilibrium. The Lattice--Boltzmann equation is coupled to the particles via a friction force. In addition to this force, acting on {it point} particles, we construct another coupling force, which comes from the pressure tensor. The coupling is purely local, i.~e. the algorithm scales linearly with the total number of particles. In order to be able to map the physical properties of the Lattice--Boltzmann fluid onto a Molecular Dynamics (MD) fluid, the case of an almost incompressible flow is considered. The Fluctuation--Dissipation theorem for the hybrid coupling is analyzed, and a geometric interpretation of the friction coefficient in terms of a Stokes radius is given. Part 2 is devoted to the simulation of charged particles. We present a novel method for obtaining Coulomb interactions as the potential of mean force between charges which are dynamically coupled to a local electromagnetic field. This algorithm scales linearly, too. We focus on the Molecular Dynamics version of the method and show that it is intimately related to the Car--Parrinello approach, while being equivalent to solving Maxwell's equations with freely adjustable speed of light. The Lagrangian formulation of the coupled particles--fields system is derived. The quasi--Hamiltonian dynamics of the system is studied in great detail. For implementation on the computer, the equations of motion are discretized with respect to both space and time. The discretization of the electromagnetic fields on a lattice, as well as the interpolation of the particle charges on the lattice is given. The algorithm is as local as possible: Only nearest neighbors sites of the lattice are interacting with a charged particle. Unphysical self--energies arise as a result of the lattice interpolation of charges, and are corrected by a subtraction scheme based on the exact lattice Green's function. The method allows easy parallelization using standard domain decomposition. Some benchmarking results of the algorithm are presented and discussed.