Conselho Nacional de Desenvolvimento Científico e Tecnológico / Neste trabalho é apresentado um novo procedimento numérico para o upscaling de permeabilidade utilizando condições de contorno periódicas. Este procedimento combina decomposição de domínio com elementos finitos mistos na discretização do problema local de pressão-velocidade necessário para se encontrar as permeabilidades equivalentes. / A new numerical method is proposed for the permeabilities upscaling take into consideration periodic boundary conditions. This method combines domain decomposition with mixed finite elements in discretization of the local problem of pressure-velocity necessary to meet the equivalent permeabilities.

Conselho Nacional de Desenvolvimento Científico e Tecnológico / Neste trabalho é apresentado um novo procedimento numérico para o upscaling de permeabilidade utilizando condições de contorno periódicas. Este procedimento combina decomposição de domínio com elementos finitos mistos na discretização do problema local de pressão-velocidade necessário para se encontrar as permeabilidades equivalentes. / A new numerical method is proposed for the permeabilities upscaling take into consideration periodic boundary conditions. This method combines domain decomposition with mixed finite elements in discretization of the local problem of pressure-velocity necessary to meet the equivalent permeabilities.

Cette thèse a pour but de développer des approches déterministes et heuristiques pour résoudre certaines classes des problèmes d'optimisation en télécommunication et la mobilité d'un réseau de transport : problèmes de routage, problèmes de covoiturage, problèmes de contrôle de l'alimentation dans un réseau sans fil, problèmes d'équilibrage du spectre dans les réseaux DSL. Il s'agit des problèmes d'optimisation non convexe de très grande taille. Nos approches sont basées sur la programmation DC&DCA, méthode de décomposition proximale et la méthode d'étiquetage des graphes. Grâce aux techniques de formulation/reformulation et de pénalité exacte, nous avons établi des programmes DC équivalents en vue de leur résolution par DCA. Selon la structure de ces problèmes, on peut fournir des décompositions DC appropriées ou de bons points initiaux de DCA. Nos méthodes ont été programmées sous MATLAB, C/C++. Ils montrent la performance de nos algorithmes par rapport à des méthodes existantes. / In this thesis, we focus on developing deterministic and heuristic approaches for solving some classes of optimization problems in Telecommunication and Mobility & Transport domain: Routing problems, Car pooloing problems, Power control problems in wireless network, Optimal spectrum balancing problems in DSL networks. They are large-scale nonconvex optimization problems. Our methodologies are focus on DC programming and DCA, Proximal decomposition method and Labeling method in graph theory. They are well-known as powerful tools in optimization. The considered problems were reformulated using the DC formulation/reformulation and exact penalty techniques and the DCA was used to obtain the solution. Also, depending on the structure of considered problems, we can provide appropriate DE decompositions or good initial points for DCA. All these proposed methods have been implemented with MATLAB, C/C++ to confirm the practical aspects and enhance our research works.

The thesis deals with qualitative properties of solutions of singular initial value problems for ordinary differential and integrodifferential equations which occur in the theory of linear and nonlinear electrical circuits and the theory of therminionic currents. The research is concentrated especially on questions of existence and uniqueness of solutions, asymptotic estimates of solutions and modications of Adomian decomposition method for singular initial problems. Solution algoritms are derived for scalar differential equations of Lane-Emden type using Taylor series and modication of the Adomian decomposition method. For certain classes of nonlinear of integrodifferential equations asymptotic expansions of solutions are constructed in a neighbourhood of a singular point. By means of the combination of Wazewski's topological method and Schauder xed-point theorem there are proved asymptotic estimates of solutions in a region which is homeomorphic to a cone having vertex coinciding with the initial point. Using Banach xed-point theorem the uniqueness of a solution of the singular initial value problem is proved for systems of integrodifferential equations of Volterra and Fredholm type including implicit systems. Moreover, conditions of continuous dependence of a solution on a parameter are determined. Obtained results are presented in illustrative examples.

Le contexte d’étude est celui de la Compatibilité ÉlectroMagnétique (CEM). L’objectif de la CEM est, comme son nom l’indique, d’assurer la compatibilité entre une source de perturbation électromagnétique et un système électronique victime. Or, la prédiction de ces niveaux de perturbation ne peut pas s’effectuer à l’aide d’un simple calcul analytique, en raison de la géométrie qui est généralement complexe pour le système que l’on étudie, tel que le champ à l’intérieur d’un cockpit d’avion par exemple. En conséquence, nous sommes contraints d’employer des méthodes numériques, dans le but de prédire ce niveau de couplage entre les sources et les victimes. Parmi les nombreuses méthodes numériques existantes à ce jour, les méthodes Multi-Domaines (MD) sont très prisées. En effet, elles offrent la liberté aux utilisateurs de choisir la méthode numérique la plus adaptée, en fonction de la zone géométrique à calculer. Au sein de ces méthodes MD, la « Domain Decomposition Method » (DDM) présente l’avantage supplémentaire de découpler chacun de ces domaines. En conséquence, la DDM est particulièrement intéressante, vis-à-vis des méthodes concurrentes, en particulier sur l’aspect du coût numérique. Pour preuve, l’ONERA continue de développer cette méthode qui ne cesse de montrer son efficacité depuis plusieurs années, notamment pour le domaine des Surfaces Équivalentes Radar (SER) et des antennes. L’objectif de l’étude est de tirer profit des avantages de cette méthode pour des problématiques de CEM. Jusqu’à maintenant, de nombreuses applications de CEM, traitées par le code DDM, fournissaient des résultats fortement bruités. Même pour des problématiques électromagnétiques très simples, des problèmes subsistaient, sans explication convaincante. Ceci justifie cette étude. Le but de cette thèse est de pouvoir appliquer ce formalisme DDM à des problématiques de CEM. Dans cette optique, nous avons été amenés à redéfinir un certain nombre de conventions, qui interviennent au sein de la DDM. Par ailleurs, nous avons développé un modèle spécifique pour les ouvertures, qui sont des voies de couplage privilégiées par les ondes, à l’intérieur des cavités que représentent les blindages. Comme les ouvertures sont, en pratique, de petites dimensions devant la longueur d’onde, on s’est intéressé à un modèle quasi-statique. Nous proposons alors un modèle, qui a été implémenté, puis validé. Suite à ce modèle, nous avons développé une méthode originale, basée sur un calcul en deux étapes, permettant de ne plus discrétiser le support des ouvertures dans les calculs 3D. / The context of the study is the ElectroMagnetic Compatibility (EMC). Principal aim of the EMC is to ensure the compatibility between an electromagnetic perturbance source and an electronic device victim. Unfortunately, the perturbation levels prediction can not be made using an analytic formula, because the geometry which is generally complex for the interesting system, for example the field inside an aircraft’s cockpit. Therefore, we are contrained to use numerical methods, to be able to evaluate this coupling level between sources and victims. Among several existing numerical methods, Multi-Domains (MD) methods are very interesting. They offer to users the freedom to choose the most powerfull numerical method, in terms of the geometrical zone evaluated. With the MD methods, « Domain Decomposition Method » (DDM) has the avantage of decouplingeach of theses areas. Therefore, DDM is very interesting, compared to other methods, in particular on the numerical cost. ONERA keeps on developing this method, which has not stop showing his efficiency since several years, in particular in Radar Cross Section (RCS) and antennas. The objective of this study is to take the benefits of this method for EMC problems. Up to now, several EMC applications treated by the DDM code provided results strongly noisy. Even for with very simple electromagnetic cases, some problems remained without convincing explanations. This justifies this study. The aim of this thesis is to can be able to apply DDM formalism to EMC problems. Then, we have been induced to redefine a number of conventions which are involved in the DDM. Otherwise, we have developed a specific model for the apertures which are privilegied tracts of the coupling by the penetration of waves inside cavities (shieldings). As the apertures have in practice smaller dimensions compared to the wavelength, we have been interested to a quasistatic model which was developped, implemented and validated. Following this model, we have developed an original method, based on a two step calculation, able to do not discretize the apertures support in 3D computations.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / Considera-se neste trabalho um modelo matemático para escoamentos com duas e três fases em reservatórios petrolíferos e a modelagem computacional do sistema de equações governantes para a sua solução numérica. Os fluidos são imiscíveis e incompressíveis e as heterogeneidades da rocha reservatório são modeladas estocasticamente. Além disso, é modelado o fenômeno de histerese para a fase óleo via funções de permeabilidades relativas. No caso de escoamentos trifásicos água-óleo-gás a escolha de expressões gerais para as funções de permeabilidades relativas pode levar à perda de hiperbolicidade estrita e, desta maneira, à existência de uma região elíptica ou de pontos umbílicos para o sistema não linear de leis de conservação hiperbólicas que descreve o transporte convectivo das fases fluidas. Como conseqüência, a perda de hiperbolicidade estrita pode levar à existência de choques não clássicos (também chamados de choques transicionais ou choques subcompressivos) nas soluções de escoamentos trifásicos, de difícil simulação numérica. Indica-se um método numérico com passo de tempo fracionário, baseado em uma técnica de decomposição de operadores, para a solução numérica do sistema governante de equações diferenciais parciais que modela o escoamento bifásico água-óleo imiscível em reservatórios de petróleo heterogêneos. Um simulador numérico bifásico água-óleo eficiente desenvolvido pelo grupo de pesquisa no qual o autor está inserido foi modificado com sucesso para incorporar a histerese sob as hipóteses consideradas. Os resultados numéricos obtidos para este caso indicam fortes evidências que o método proposto pode ser estendido para o caso trifásico água-óleo-gás. A técnica de decomposição de operadores em dois níveis permite o uso de passos de tempo distintos para os quatro problemas definidos pelo procedimento de decomposição: convecção, difusão, pressão-velocidade e relaxação para histerese. O problema de transporte convectivo (hiperbólico) das fases fluidas é aproximado por um esquema central de diferenças finitas explícito, conservativo, não oscilatório e de segunda ordem. Este esquema é combinado com elementos finitos mistos, localmente conservativos, para a aproximação dos problemas de transporte difusivo (parabólico) e de pressão-velocidade (elíptico). O operador temporal associado ao problema parabólico de difusão é resolvido fazendo-se uso de uma estratégia implícita de solução (Backward Euler). Uma equação diferencial ordinária é resolvida (analiticamente) para a relaxação relacionada à histerese. Resultados numéricos para o problema bifásico água-óleo em uma dimensão espacial em concordância com resultados semi-analíticos disponíveis na literatura foram reproduzidos e novos resultados em meios heterogêneos, em duas dimensões espaciais, são apresentados e a extensão desta técnica para o caso de problemas trifásicos água-óleo-gás é proposta. / We consider in this work a mathematical model for two- and three-phase flow problems in petroleum reservoirs and the computational modeling of the governing equations for its numerical solution. We consider two- (water-oil) and three-phase (water-gas-oil) incompressible, immiscible flow problems and the reservoir rock is considered to be heterogeneous. In our model, we also take into account the hysteresis effects in the oil relative permeability functions. In the case of three-phase flow, the choice of general expressions for the relative permeability functions may lead to the loss of strict hyperbolicity and, therefore, to the existence of an elliptic region or umbilic points for the system of nonlinear hyperbolic conservation laws describing the convective transport of the fluid phases. As a consequence, the loss of hyperbolicity may lead to the existence of nonclassical shocks (also called transitional shocks or undercompressive shocks) in three-phase flow solutions. We present a new, accurate fractional time-step method based on an operator splitting technique for the numerical solution of a system of partial differential equations modeling two-phase, immiscible water-oil flow problems in heterogeneous petroleum reservoirs. An efficient two-phase water-oil numerical simulator developed by our research group was sucessfuly extended to take into account hysteresis effects under the hypotesis previously annouced. The numerical results obtained by the procedure proposed indicate numerical evidence the method at hand can be extended for the case of related three-phase water-gas-oil flow problems. A two-level operator splitting technique allows for the use of distinct time steps for the four problems defined by the splitting procedure: convection, diffusion, pressure-velocity and relaxation for hysteresis. The convective transport (hyperbolic) of the fluid phases is approximated by a high resolution, nonoscillatory, second-order, conservative central difference scheme in the convection step. This scheme is combined with locally conservative mixed finite elements for the numerical solution of the diffusive transport (parabolic) and the pressure-velocity (elliptic) problems. The time discretization of the parabolic problem is performed by means of the implicit Backward Euler method. An ordinary diferential equation is solved (analytically) for the relaxation related to hysteresis. Two-phase water-oil numerical results in one space dimensional, in which are in a very good agreement with semi-analitycal results available in the literature, were computationaly reproduced and new numerical results in two dimensional heterogeneous media are also presented and the extension of this technique to the case of three-phase water-oil-gas flows problems is proposed.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / Considera-se neste trabalho um modelo matemático para escoamentos com duas e três fases em reservatórios petrolíferos e a modelagem computacional do sistema de equações governantes para a sua solução numérica. Os fluidos são imiscíveis e incompressíveis e as heterogeneidades da rocha reservatório são modeladas estocasticamente. Além disso, é modelado o fenômeno de histerese para a fase óleo via funções de permeabilidades relativas. No caso de escoamentos trifásicos água-óleo-gás a escolha de expressões gerais para as funções de permeabilidades relativas pode levar à perda de hiperbolicidade estrita e, desta maneira, à existência de uma região elíptica ou de pontos umbílicos para o sistema não linear de leis de conservação hiperbólicas que descreve o transporte convectivo das fases fluidas. Como conseqüência, a perda de hiperbolicidade estrita pode levar à existência de choques não clássicos (também chamados de choques transicionais ou choques subcompressivos) nas soluções de escoamentos trifásicos, de difícil simulação numérica. Indica-se um método numérico com passo de tempo fracionário, baseado em uma técnica de decomposição de operadores, para a solução numérica do sistema governante de equações diferenciais parciais que modela o escoamento bifásico água-óleo imiscível em reservatórios de petróleo heterogêneos. Um simulador numérico bifásico água-óleo eficiente desenvolvido pelo grupo de pesquisa no qual o autor está inserido foi modificado com sucesso para incorporar a histerese sob as hipóteses consideradas. Os resultados numéricos obtidos para este caso indicam fortes evidências que o método proposto pode ser estendido para o caso trifásico água-óleo-gás. A técnica de decomposição de operadores em dois níveis permite o uso de passos de tempo distintos para os quatro problemas definidos pelo procedimento de decomposição: convecção, difusão, pressão-velocidade e relaxação para histerese. O problema de transporte convectivo (hiperbólico) das fases fluidas é aproximado por um esquema central de diferenças finitas explícito, conservativo, não oscilatório e de segunda ordem. Este esquema é combinado com elementos finitos mistos, localmente conservativos, para a aproximação dos problemas de transporte difusivo (parabólico) e de pressão-velocidade (elíptico). O operador temporal associado ao problema parabólico de difusão é resolvido fazendo-se uso de uma estratégia implícita de solução (Backward Euler). Uma equação diferencial ordinária é resolvida (analiticamente) para a relaxação relacionada à histerese. Resultados numéricos para o problema bifásico água-óleo em uma dimensão espacial em concordância com resultados semi-analíticos disponíveis na literatura foram reproduzidos e novos resultados em meios heterogêneos, em duas dimensões espaciais, são apresentados e a extensão desta técnica para o caso de problemas trifásicos água-óleo-gás é proposta. / We consider in this work a mathematical model for two- and three-phase flow problems in petroleum reservoirs and the computational modeling of the governing equations for its numerical solution. We consider two- (water-oil) and three-phase (water-gas-oil) incompressible, immiscible flow problems and the reservoir rock is considered to be heterogeneous. In our model, we also take into account the hysteresis effects in the oil relative permeability functions. In the case of three-phase flow, the choice of general expressions for the relative permeability functions may lead to the loss of strict hyperbolicity and, therefore, to the existence of an elliptic region or umbilic points for the system of nonlinear hyperbolic conservation laws describing the convective transport of the fluid phases. As a consequence, the loss of hyperbolicity may lead to the existence of nonclassical shocks (also called transitional shocks or undercompressive shocks) in three-phase flow solutions. We present a new, accurate fractional time-step method based on an operator splitting technique for the numerical solution of a system of partial differential equations modeling two-phase, immiscible water-oil flow problems in heterogeneous petroleum reservoirs. An efficient two-phase water-oil numerical simulator developed by our research group was sucessfuly extended to take into account hysteresis effects under the hypotesis previously annouced. The numerical results obtained by the procedure proposed indicate numerical evidence the method at hand can be extended for the case of related three-phase water-gas-oil flow problems. A two-level operator splitting technique allows for the use of distinct time steps for the four problems defined by the splitting procedure: convection, diffusion, pressure-velocity and relaxation for hysteresis. The convective transport (hyperbolic) of the fluid phases is approximated by a high resolution, nonoscillatory, second-order, conservative central difference scheme in the convection step. This scheme is combined with locally conservative mixed finite elements for the numerical solution of the diffusive transport (parabolic) and the pressure-velocity (elliptic) problems. The time discretization of the parabolic problem is performed by means of the implicit Backward Euler method. An ordinary diferential equation is solved (analytically) for the relaxation related to hysteresis. Two-phase water-oil numerical results in one space dimensional, in which are in a very good agreement with semi-analitycal results available in the literature, were computationaly reproduced and new numerical results in two dimensional heterogeneous media are also presented and the extension of this technique to the case of three-phase water-oil-gas flows problems is proposed.

No description available.