Implementação paralelizada de métodos de resolução de sistemas algébricos na simulação de reservatórios de gás

Desde a década de 1960, devido à pertinência para a indústria petrolífera, a simulação numérica de reservatórios de petróleo tornou-se uma ferramenta usual e uma intensa área de pesquisa. O principal objetivo da modelagem computacional e do uso de métodos numéricos, para a simulação de reservatórios de petróleo, é o de possibilitar um melhor gerenciamento do campo produtor, de maneira que haja uma maximização na recuperação de hidrocarbonetos. Este trabalho tem como objetivo principal paralelizar, empregando a interface de programação de aplicativo OpenMP (Open Multi-Processing), o método numérico utilizado na resolução do sistema algébrico resultante da discretização da equação que descreve o escoamento monofásico em um reservatório de gás, em termos da variável pressão. O conjunto de equações governantes é formado pela equação da continuidade, por uma expressão para o balanço da quantidade de movimento e por uma equação de estado. A Equação da Difusividade Hidráulica (EDH), para a variável pressão, é obtida a partir deste conjunto de equações fundamentais, sendo então discretizada pela utilização do Método de Diferenças Finitas, com a escolha por uma formulação implícita. Diferentes testes numéricos são realizados a fim de estudar a eficiência computacional das versões paralelizadas dos métodos iterativos de Jacobi, Gauss-Seidel, Sobre-relaxação Sucessiva, Gradientes Conjugados (CG), Gradiente Biconjugado (BiCG) e Gradiente Biconjugado Estabilizado (BiCGStab), visando a uma futura aplicação dos mesmos na simulação de reservatórios de gás. Ressalta-se que a presença de heterogeneidades na rocha reservatório e/ou às não-linearidades presentes na EDH para o escoamento de gás aumentam a necessidade de métodos eficientes do ponto de vista de custo computacional, como é o caso de estratégias usando OpenMP.

Application of fast Pade approximation in simulating photonic crystal nanocavities by FDTD technology

The exact calculation of mode quality factor Q is a key problem in the design of high-Q photonic crystal nanocavity. On the basis of further investigation on conventional Pade approximation, FDM and DFT, Pade approximation with Baker's algorithm is enhanced through introducing multiple frequency search and parabola interpolation. Though Pade approximation is a nonlinear signal processing method and only short time sequence is needed, we find the different length of sequence requirements for 2D and 3D FDTD, which is very important to obtain convergent and accurate results. By using the modified Pade approximation method and 3D FDTD, the 2D slab photonic crystal nanocavity is analyzed and high-Q multimode can be solved quickly instead of large range high-resolution scanning. Monitor position has also been investigated. These results are very helpful to the design of photonic crystal nanocavity devices. (C) 2008 Elsevier B.V. All rights reserved.

Two-dimensional simulation of high-order laterally-coupled GaAs-AlGaAs DFB laser diodes

Laterally-coupled distributed feedback (LC-DFB) laser diodes made without an epitaxial re-growth process have the advantage of a simple fabrication process. In this paper, two-dimensional optical field distribution of the fundamental quasi TE (transverse electric) mode is calculated by means of a semivectorial finite-difference method (SV-FDM). The dependence of the effective coupling coefficient (kappa(eff)) on the dutycycle of first-, second- and third-order LC-DFB LDs is investigated using modified coupled wave equations.

地震波谱元法数值模拟和偏移研究

With the development of both seismic theory and computer technology, numerical modeling technology of seismic wave has achieved great advancement during the past half century. The current methods under development include finite differentiation method (FDM), finite element method (FEM), pseudospectral method (PSM), integral equation method (IEM) and spectral element method (SEM). They exert their very important roles in every corner of seismology and seismic prospecting. Large quantity of researches towards spectral element method in the end of last century bring this method to a new era, which results in perfect solution of many difficult problems. However, parts of posterior works such as seismic migration and inversion which base on spectral element method have never been studied widely at least up to the present whereas are of importance to seismic imaging and seismic wave propagation. Based on previous work, this paper uses spectral element method to investigate the characteristics and laws of the seismic wave propagation in isotropic and anisotropic media. By thoroughly studying this high-accuracy method, we implement a kind of reverse-time pre- and post-stack migration based on SEM. In order to verify the validity of the SEM method, we have simulated the propagation of seismic wave in several different models. The simulation results show that: (1) spectral element method can be used to model any complex models and the computational results are comparable with the expected results and the analytic results; (2) the optimum accuracy can be achieved when the rank is between 4 and 9. When it is below 4, the dispersion may occur; and when it is above 9, the time step-length will be changed accordingly with the reducing space step-length in order to keep the computation stability. This will exponentially increase the computation time and at the same time the memory even if simulating the same media. This paper also applies explosive reflection surface imaging technology, time constancy principle of wave-filed extrapolation and least travetime raytracing technology of surface source to SEM pre- and post-stack migration of isotropic and anisotropic media. All imaging results derived by the above methods agree well with the real geological models and the position of interface and inflexions can also return to their right location well. This indicates that the method proposed in this paper is a kind of technology with high accuracy and robust stability. It can serve as an alternative method in real seismic data processing. All these work can boost the development of high-accuracy seismic imaging, and therefore have significant inference value.

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

Efectes de la competici?? sobre la for??a din??mica m??xima en el jugador de tennis d'elit : estudi d'un cas

Contiene tablas y gr??ficos de resultados de ejercicios y de entrenamientos. Resumen tomado del autor

Revisión sistemática : uso de imágenes por difusión en RMN, para predecir sobrevida en pacientes adultos con diagnostico de glioblastoma multiforme

El Glioblastoma multiforme (GBM), es el tumor cerebral más frecuente, con pronóstico grave y baja sensibilidad al tratamiento inicial. El propósito de este estudio fue evaluar si la Difusión en RM (IDRM), es un biomarcador temprano de respuesta tumoral, útil para tomar decisiones tempranas de tratamiento y para obtener información pronostica. Metodología La búsqueda se realizo en las bases de datos EMBASE, CENTRAL, MEDLINE; las bibliografías también fueron revisadas. Los artículos seleccionados fueron estudios observacionales (casos y controles, cohortes, corte transversal), no se encontró ningún ensayo clínico; todos los participante tenían diagnostico histopatológico de GBM, sometidos a resección quirúrgica y/o radio-quimioterapia y seguimiento de respuesta al tratamiento con IDRM por al menos 6 meses. Los datos extraídos de forma independiente fueron tipo de estudio, participantes, intervenciones, seguimiento, desenlaces (sobrevida, progresión/estabilización de la enfermedad, muerte) Resultados Quince estudios cumplieron los criterios de inclusión. Entre las técnicas empleadas de IDRM para evaluar respuesta radiológica al tratamiento, fueron histogramas del coeficiente aparente de difusion ADC (compararon valores inferiores a la media y el percentil 10 de ADC, con los valores superiores); encontrando en términos generales que un ADC bajo es un fuerte predictor de sobrevida y/o progresión del tumor. (Esto fue significativo en 5 estudios); mapas funcionales de difusion (FDM) (midieron el porcentaje de cambio de ADC basal vs pos tratamiento) que mostro ser un fuerte predictor de sobrevida en pacientes con progresión tumoral. DISCUSION Desafortunadamente la calidad de los estudios fue intermedia-baja lo que hace que la aplicabilidad de los estudios sea limitada.

Effects of fat on the properties of halloumi cheese

Halloumi cheese was produced from 11 bovine milks with fat contents of 1.61-4.04%, giving a range of 32-53% fat in dry matter (FDM) in the cheeses. Starter culture and/or microparticulated whey protein (Simplesse((R)) 100(E)) was also added to selected batches of milk. Hardness decreased with increasing FDM, with increase in moisture and with lower pH. On sensory evaluation, there was an increase in preference score with FDM (R-2 = 0.8). Inclusion of microparticulated whey protein may have had a fat mimetic effect, as preference scores otherwise decreased with increasing protein levels (R-2 = 0.75).

Os discursos sobre língua e ensino no Brasil da 1a. e da 2a. república : o duplo lugar da determinação e da contradição

A presente pesquisa, filiada à Análise de Discurso de linha francesa, trata dos processos de determinação na regulamentação das questões de ensino de língua portuguesa no Brasil no período histórico compreendido entre a 1ª e a 2ª Repúblicas (1889-1945). Por estar filiado à referida teoria, este trabalho não busca esgotar o recorte cronológico e nem tampouco abarcar a totalidade das enunciações que foram feitas, nessa época, em torno da questão selecionada para estudo. Os gestos de análise ora desenhados recortaram as enunciações em torno das questões de ensino e de língua, tecidas a partir de quatro formações discursivas diferentes, desdobradas em suas contradições, através: a) da posição-sujeito dos operários-anarquistas, identificada à Formação Discursiva Libertária (FDL); b)da posição-sujeito das enunciações das propostas de Reforma do Ensino na 1ª República, identificada à Formação Discursiva Estatal (FDE); c) de enunciações da posição sujeito situada no âmago das formulações jurídicas elaboradas no Governo Vargas, identificada à Formação Discursiva Ministerial (FDM); d) e por último as enunciações constantes em gramáticas, identificadas à a Formação Discursiva Pedagógica (FDP). As análises deste trabalho de investigação estão articuladas de modo que os saberes das FDs em questão são olhados de formas diferentes, no sentido de que se situando em lugares de antagonismo e de contradição, em muitos momentos faz-se possível observar essas relações que, no interior de uma dada FD, remetem aos saberes do discursooutro. No recorte específico dos saberes sobre a língua privilegiamos a constituição de um imaginário de língua derivado do imaginário de nação, buscando investigar em que medida e de que formas eles (re) produziram posturas xenofobistas que ressoaram nos modos como foi tipificado o sujeito imigrante operário, a partir de sua língua e de sua prática política. Assim, recorremos ao interdiscurso, lugar onde os enunciados se articularam, descrevendo os diferentes modos como esses foram linearizados e, assim, produziram sentidos, no embate tenso entre os jogos de aliança e refutação. Para proceder à fundamentação teórica que deu suporte às análises, tecemos diálogos com outras regiões do conhecimento, desde o campo da filosofia, nas discussões envolvendo as questões de designação e determinação, passando pelo olhar das gramáticas (cartesiana, filosófica, histórica e normativa) sobre a língua, para então chegarmos à ordem do discurso, lugar onde conceitos foram ressignificados, para podermos tratar da língua, enquanto objeto nunca imune às condições de produção dos discursos por ela materializados.

Spatial variability of forage yield and soil physical attributes of a Brachiaria decumbens pasture in the Brazilian Cerrado

The objective of this study was to analyze variability, linear and spatial correlations of forage dry mass yield (FDM) and dry matter percentage (DM%) of Brachiaria decumbens with the bulk density (BD), gravimetric (GM) and volumetric (VM) moisture, mechanical resistance to penetration (RP) and organic matter content (OM), at depths 1 (0-0.10 m) and 2 (0.10-0.20 m), in a Red Latosol (Oxisol), in order to select an indicator of soil physical quality and identify possible causes of pasture degradation. The geostatistical grid was installed to collect soil and plant data, with 121 sampling points, over an area of 2.56 ha. The linear correlation between FDM × DM% and FDM × BD2 was low, but highly significant. Spatial correlations varied inversely and positively, respectively. Except for DM% and BD, at both depths, the other attributes showed average to high variability, indicating a heterogeneous environment. Thus, geostatistics emerges as an important tool in understanding the interactions in pasture ecosystems, in order to minimize possible causes of degradation and indicate better alternatives for soil-plant-animal management. The decrease in FDM and increased BD1 are indicators of physical degradation (compaction) of Red Latosol (Oxisol), particularly in the places with the highest concentration of animals and excessive trampling, in Cerrado conditions, in the municipality of Selvíria, Mato Grosso do Sul State, Brazil.

Spatial and linear correlations between soil and corn

O conjunto de tecnologias aplicadas ao sistema de produção agrícola tem como características principais a verticalização da produtividade, diminuição de custos, melhoria nas características físicas, químicas e biológicas do solo para proporcionarem o crescimento sustentável do meio de produção. Desta forma, o trabalho teve como objetivo determinar a variabilidade e as correlações lineares e especiais entre os atributos da planta e do solo, visando selecionar um indicador da qualidade física do solo de boa representatividade para produção de forragem. No ano agrícola de 2006, na Fazenda Bonança, município de Pereira Barreto (SP), foram analisados a produtividade de forragem do milho outonal (FDM) em sistema plantio direto irrigado e os atributos físicos do solo, num Latossolo Vermelho distrófico. O objetivo foi estudar a variabilidade e as correlações lineares e espaciais entre os atributos da planta e os do solo, visando selecionar um indicador da qualidade física do solo de boa representatividade para a produtividade da forragem. Foi instalada a malha geoestatística, para coleta de dados do solo e planta, contendo 125 pontos amostrais, numa área de 2.500 m². Os atributos estudados, além de não terem variado aleatoriamente, apresentaram variabilidade dos dados entre baixa e muito alta e seguiram padrões espaciais bem definidos, com alcances entre 7,8 e 38,0 m. Por outro lado, a correlação linear entre os atributos da planta com o do solo foi baixa e extremamente significativa. Os pares Massa Seca de forragem versus Microporosidade e Diâmetro do colmo versus Densidade do Solo foram melhor correlacionados na camada de 0-0.10m, enquanto os outros pares - Massa Seca de Forragem versus Macroporosidade - e Porosidade Total - apresentaram correlação inversa para a mesma camada. Entretanto, do ponto de vista espacial, houve uma alta correlação inversa entre Massa Seca de Forragem com Microporosidade, de modo que a microporosidade na camada de 0-0.10m pode ser considerada um bom indicador de qualidade física do solo, tendo em vista a produção de forragem de milho.

Adaptive control charts: A Markovian approach for processes subject to independent disturbances

Recent studies have shown that adaptive X control charts are quicker than traditional X charts in detecting small to moderate shifts in a process. In this article, we propose a joint statistical design of adaptive X and R charts having all design parameters varying adaptively. The process is subjected to two independent assignable causes. One cause changes the process mean and the other changes the process variance. However, the occurrence of one kind of assignable cause does not preclude the occurrence of the other. It is assumed that the quality characteristic is normally distributed and the time that the process remains in control has exponential distribution. Performance measures of these adaptive control charts are obtained through a Markov chain approach. (c) 2005 Elsevier B.V. All rights reserved.

Carboxylic acids of marine evaporitic oils from Sergipe-Alagoas Basin, Brazil

A. suite of 10 different marine evaporitic oil samples from Sergipe-Alagoas Basin, Brazil was studied for its biomarker content, in particular its acidic constituents. The oils showed different molecular distributions and relative abundances of n-alkanoic, isoprenoid and hopanoic acids. The observed differences were assigned to the incorporation of immature organic matter in the oils and fractionation along the migration pathway. The diagenetic precursor functionality (alcohol/ether or acid) was proposed based on the comparison of the relative abundances of the neutral and acidic biomarkers (hopanoids, isoprenoids, alkyl-steranes, monoaromatic alkyl-steroids). In the acidic fraction, 3 series of steroid-alkanoic acids and monoaromatic steroid-alkanoic acids (steroid-methanoic, ethanoic and propanoic acids and monoaromatic steroid-methanoic, ethanoic and propanoic acids) were detected, while in the neutral fraction only 2 series of each corresponding class could be observed (methyl and ethyl-steranes and monoaromatic methyl and ethyl-steroids). These carbon shifts suggest that decarboxylation is an important process in the formation of the alkyrsteranes and monoaromatic alkyl-steroids, and we infer that carboxylic acids are the diagenetic precursors of these classes of compounds. When alcohol or ether are the diagenetic precursors (isoprenoids and hopanoids), no significant differences in the molecular distributions between neutral and acidic fractions were observed. (C) 2000 Elsevier B.V. Ltd. All rights reserved.

Predicting breeding values for milk yield of Guzera (Bos indicus) cows using random regression models

Given the importance of Guzera breeding programs for milk production in the tropics, the objective of this study was to compare alternative random regression models for estimation of genetic parameters and prediction of breeding values. Test-day milk yields records (TDR) were collected monthly, in a maximum of 10 measurements. The database included 20,524 records of first lactation from 2816 Guzera cows. TDR data were analyzed by random regression models (RRM) considering additive genetic, permanent environmental and residual effects as random and the effects of contemporary group (CG), calving age as a covariate (linear and quadratic effects) and mean lactation curve as fixed. The genetic additive and permanent environmental effects were modeled by RRM using Wilmink, All and Schaeffer and cubic B-spline functions as well as Legendre polynomials. Residual variances were considered as heterogeneous classes, grouped differently according to the model used. Multi-trait analysis using finite-dimensional models (FDM) for testday milk records (TDR) and a single-trait model for 305-days milk yields (default) using the restricted maximum likelihood method were also carried out as further comparisons. Through the statistical criteria adopted, the best RRM was the one that used the cubic B-spline function with five random regression coefficients for the genetic additive and permanent environmental effects. However, the models using the Ali and Schaeffer function or Legendre polynomials with second and fifth order for, respectively, the additive genetic and permanent environmental effects can be adopted, as little variation was observed in the genetic parameter estimates compared to those estimated by models using the B-spline function. Therefore, due to the lower complexity in the (co)variance estimations, the model using Legendre polynomials represented the best option for the genetic evaluation of the Guzera lactation records. An increase of 3.6% in the accuracy of the estimated breeding values was verified when using RRM. The ranks of animals were very close whatever the RRM for the data set used to predict breeding values. Considering P305, results indicated only small to medium difference in the animals' ranking based on breeding values predicted by the conventional model or by RRM. Therefore, the sum of all the RRM-predicted breeding values along the lactation period (RRM305) can be used as a selection criterion for 305-day milk production. (c) 2014 Elsevier B.V. All rights reserved.

Deconvolution of the EPR spectra of vanadium oxide nanotubes

In this work we report results of continuous wave (CW) electron paramagnetic resonance (EPR) spectroscopy of vanadium oxide nanotubes. The observed EPR spectra are composed of a weak well-resolved spectrum of isolated V4+ ions on top of an intense and broad structure-less line shape, attributed to spin-spin exchanged V4+ clusters. With the purpose to deconvolute the structured weak spectrum from the composed broad line, a new approach based on the Krylov basis diagonalization method (KBDM) is introduced. It is based on the discrimination between broad and sharp components with respect to a selectable threshold and can be executed with few adjustable parameters, without the need of a priori information on the shape and structure of the lines. This makes the method advantageous with respect to other procedures and suitable for fast and routine spectral analysis, which, in conjunction with simulation techniques based on the spin Hamiltonian parameters, can provide a full characterization of the EPR spectrum. Results demonstrate and characterize the coexistence of two V4+ species in the nanotubes and show good progress toward the goal of obtaining high fidelity deconvoluted spectra from complex signals with overlapping broader line shapes. (C) 2012 Elsevier Inc. All rights reserved.