Nesta contribuição apresentamos uma solução para a equação de transferência radiativa condutiva em geometria cilíndrica. Esta solução é aplicada para simular a radiação e campo de temperatura juntamente com o transporte de energia radiativa e condutiva proveniente do escape liberado em lançamentos de foguetes. Para este fim, discutimos uma abordagem semianalítica reduzindo a equação original, que é contínua nas variáveis angulares, numa equação semelhante ao problema SN da transferência radiativa condutiva. A solução é construída usando um método de composição por transformada de Laplace e o método da decomposição de Adomian. O esquema recursivo ´e apresentado para o sistema de equações de ordenadas duplamente discretas juntamente com as dependências dos parâmetros e suas influências sobre a convergência heurística da solução. A solução obtida, em seguida, permite construir o campo próximo relevante para caracterizar o termo fonte para problemas de dispersão ao ajustar os parâmetros do modelo, tais como, emissividade, refletividade, albedo e outros, em comparação com a observação, que são relevantes para os processos de dispersão de campo distante e podem ser manipulados de forma independente do presente problema. Além do método de solução, também relatamos sobre algumas soluções e simulações numéricas. / In this contribution we present a solution for the radiative conductive transfer equation in cylinder geometry. This solution is applied to simulate the radiation and temperature field together with conductive and radiative energy transport originated from the exhaust released in rocket launches. To this end we discuss a semi-analytical approach reducing the original equation, which is continuous in the angular variables, into an equation similar to the SN radiative conductive transfer problem. The solution is constructed using a composite method by Laplace transform and Adomian decomposition method. The recursive scheme is presented for the doubly discrete ordinate equations system together with parameter dependencies and their influence on heuristic convergence of the solution. The obtained solution allows then to construct the relevant near field to characterize the source term for dispersion problems when adjusting the model parameters such as emissivity, reflectivity, albedo and others in comparison to the observation, that are relevant for far field dispersion processes and may be handled independently from the present problem. In addition to the solution method we also report some solutions and numerical simulations.

Made available in DSpace on 2016-06-02T19:51:37Z (GMT). No. of bitstreams: 1 2228.pdf: 979050 bytes, checksum: ffa6f96c8eada124b6f1e6ba3ebe02da (MD5) Previous issue date: 2008-12-11 / Financiadora de Estudos e Projetos / Cutting and packing problems are important in the production planning of various industrial segments involving goals such as minimizing the negative efects generated by waste of materials or idle spaces. The loss of material due to an inadequate programming of the cutting or packing patterns, can be substantial, and, in general, parts of these losses can be avoided only with a more eficient production planning, not resulting in additional investments in production processes. This study aimed at evaluating the performance of the Benders decomposition method, applied to the manufacturer and distributor pallet loading models. The manufacturer pallet loading model involves packing equal boxes on a pallet, so as to optimize its use. The distributor pallet loading model involves packing boxes of diferent sizes on a pallet, also a way to optimize its use. The approach based on Benders decomposition, defines a relaxation algorithm that partitions the original problem in two other problems easier to be solved. To check the effectiveness of the approach, computational tests were carried out by comparing the results with those obtained by a computational package composed of a modeling language (GAMS) and a last generation optimization solver (CPLEX ). / Os problemas de corte e empacotamento são importantes no planejamento da produção de vários segmentos industriais envolvendo objetivos como, por exemplo, minimizar os efeitos negativos gerados por desperdício de materiais ou espaços ociosos. As perdas de material, devido a uma programação pouco adequada dos padrões de corte ou empacotamento, podem ser substanciais, sendo que, em geral, parte destas perdas pode ser evitada apenas com uma programação da produção mais eficiente, não implicando em investimentos adicionais nos processos de produção. O objetivo deste estudo é verificar o desempenho do método de decomposição de Benders aplicado a modelos de carregamento de paletes do produtor e do distribuidor. O problema de carregamento de paletes do produtor envolve empacotar caixas iguais sobre um palete, de maneira a otimizar o aproveitamento deste. O problema de carregamento de paletes do distribuidor envolve empacotar caixas de tamanhos diferentes sobre um palete, também de maneira a otimizar o aproveitamento deste. A abordagem baseada na reformulação de Benders define um algoritmo de relaxação que particiona o problema original em dois outros problemas mais simples de serem resolvidos. Para verificar a eficiência da abordagem, realizaram-se testes computacionais, comparando os resultados obtidos com os obtidos pelo pacote computacional composto de uma linguagem de modelagem (GAMS) e um software de otimização de última geração (CPLEX).

Nesta contribuição apresentamos uma solução para a equação de transferência radiativa condutiva em geometria cilíndrica. Esta solução é aplicada para simular a radiação e campo de temperatura juntamente com o transporte de energia radiativa e condutiva proveniente do escape liberado em lançamentos de foguetes. Para este fim, discutimos uma abordagem semianalítica reduzindo a equação original, que é contínua nas variáveis angulares, numa equação semelhante ao problema SN da transferência radiativa condutiva. A solução é construída usando um método de composição por transformada de Laplace e o método da decomposição de Adomian. O esquema recursivo ´e apresentado para o sistema de equações de ordenadas duplamente discretas juntamente com as dependências dos parâmetros e suas influências sobre a convergência heurística da solução. A solução obtida, em seguida, permite construir o campo próximo relevante para caracterizar o termo fonte para problemas de dispersão ao ajustar os parâmetros do modelo, tais como, emissividade, refletividade, albedo e outros, em comparação com a observação, que são relevantes para os processos de dispersão de campo distante e podem ser manipulados de forma independente do presente problema. Além do método de solução, também relatamos sobre algumas soluções e simulações numéricas. / In this contribution we present a solution for the radiative conductive transfer equation in cylinder geometry. This solution is applied to simulate the radiation and temperature field together with conductive and radiative energy transport originated from the exhaust released in rocket launches. To this end we discuss a semi-analytical approach reducing the original equation, which is continuous in the angular variables, into an equation similar to the SN radiative conductive transfer problem. The solution is constructed using a composite method by Laplace transform and Adomian decomposition method. The recursive scheme is presented for the doubly discrete ordinate equations system together with parameter dependencies and their influence on heuristic convergence of the solution. The obtained solution allows then to construct the relevant near field to characterize the source term for dispersion problems when adjusting the model parameters such as emissivity, reflectivity, albedo and others in comparison to the observation, that are relevant for far field dispersion processes and may be handled independently from the present problem. In addition to the solution method we also report some solutions and numerical simulations.

Neste trabalho desenvolveu-se e implementou-se um modelo computacional paralelo multifísica para a simulação do transporte de substâncias e do escoamento hidrodinâmico, bidimensional (2D) e tridimensional (3D), em corpos de água. Sua motivação está centrada no fato de que as margens e zonas costeiras de rios, lagos, estuários, mares e oceanos são locais de aglomerações de seres humanos, dada a sua importância para as atividades econômica, de transporte e de lazer, causando desequilíbrios a esses ecossistemas. Esse fato impulsiona o desenvolvimento de pesquisas relativas a esta temática. Portanto, o objetivo deste trabalho é o de construir um modelo computacional com alta qualidade numérica, que possibilite simular os comportamentos da hidrodinâmica e do transporte escalar de substâncias em corpos de água com complexa configuração geométrica, visando a contribuir para seu manejo racional. Visto que a ênfase nessa tese são os aspectos numéricos e computacionais dos algoritmos, analisaram-se as características e propriedades numérico-computacionais que as soluções devem contemplar, tais como a estabilidade, a monotonicidade, a positividade e a conservação da massa. As estratégias de soluções enfocam os termos advectivos e difusivos, horizontais e verticais, da equação do transporte. Desse modo, a advecção horizontal é resolvida empregando o método da limitação dos fluxos de Sweby, e o transporte vertical (advecção e difusão) é resolvido com os métodos beta de Gross e de Crank-Nicolson. São empregadas malhas com distintas resoluções para a solução do problema multifísica. O esquema numérico resultante é semi-implícito, computacionalmente eficiente, estável e fornece acurácia espacial e temporal de segunda ordem. Os sistemas de equações resultantes da discretização, em diferenças finitas, das equações do escoamento e do transporte 3D, são de grande porte, lineares, esparsos e simétricos definidos-positivos (SDP). No caso 2D os sistemas são lineares, mas os sistemas de equações para a equação do transporte não são simétricos. Assim, para a solução de sistemas de equações SDP e dos sistemas não simétricos empregam-se, respectivamente, os métodos do subespaço de Krylov do gradiente conjugado e do resíduo mínimo generalizado. No caso da solução dos sistemas 3-diagonal, utiliza-se o algoritmo de Thomas e o algoritmo de Cholesky. A solução paralela foi obtida sob duas abordagens. A decomposição ou particionamento de dados, onde as operações e os dados são distribuídos entre os processos disponíveis e são resolvidos em paralelo. E, a decomposição de domínio, onde obtém-se a solução do problema global combinando as soluções de subproblemas locais. Em particular, emprega-se neste trabalho, o método de decomposição de domínio aditivo de Schwarz, como método de solução, e como pré-condicionador. Para maximizar a relação computação/comunicação, visto que a eficiência computacional da solução paralela depende diretamente do balanceamento de carga e da minimização da comunicação entre os processos, empregou-se algoritmos de particionamento de grafos para obter localmente os subproblemas, ou as partes dos dados. O modelo computacional paralelo resultante mostrou-se computacionalmente eficiente e com alta qualidade numérica. / A multi-physics parallel computational model was developed and implemented for the simulation of substance transport and for the two-dimensional (2D) and threedimensional (3D) hydrodynamic flow in water bodies. The motivation for this work is focused in the fact that the margins and coastal zones of rivers, lakes, estuaries, seas and oceans are places of human agglomeration, because of their importance for economic, transport, and leisure activities causing ecosystem disequilibrium. This fact stimulates the researches related to this topic. Therefore, the goal of this work is to build a computational model of high numerical quality, that allows the simulation of hydrodynamics and of scalar transport of substances behavior in water bodies of complex configuration, aiming at their rational management. Since the focuses of this thesis are the numerical and computational aspects of the algorithms, the main numerical-computational characteristics and properties that the solutions need to fulfill were analyzed. That is: stability, monotonicity, positivity and mass conservation. Solution strategies focus on advective and diffusive terms, horizontal and vertical terms of the transport equation. In this way, horizontal advection is solved using Sweby’s flow limiting method; and the vertical transport (advection and diffusion) is solved with Gross and Crank-Nicolson’s beta methods. Meshes of different resolutions are employed in the solution of the multi-physics problem. The resulting numerical scheme is semi-implicit, computationally efficient, stable and provides second order accuracy in space and in time. The equation systems resulting of the discretization, in finite differences, of the flow and 3D transport are of large scale, linear, sparse and symmetric positive definite (SPD). In the 2D case, the systems are linear, but the equation systems for the transport equation are not symmetric. Therefore, for the solution of SPD equation systems and of the non-symmetric systems we employ, respectively, the methods of Krylov’s sub-space of the conjugate gradient and of the generalized minimum residue. In the case of the solution of 3-diagonal systems, Thomas algorithm and Cholesky algorithm are used. The parallel solution was obtained through two approaches. In data decomposition or partitioning, operation and data are distributed among the processes available and are solved in parallel. In domain decomposition the solution of the global problem is obtained combining the solutions of the local sub-problems. In particular, in this work, Schwarz additive domain decomposition method is used as solution method and as preconditioner. In order to maximize the computation/communication relation, since the computational efficiency of the parallel solution depends directly of the load balancing and of the minimization of the communication between processes, graph-partitioning algorithms were used to obtain the sub-problems or part of the data locally. The resulting parallel computational model is computationally efficient and of high numerical quality.

Neste trabalho desenvolveu-se e implementou-se um modelo computacional paralelo multifísica para a simulação do transporte de substâncias e do escoamento hidrodinâmico, bidimensional (2D) e tridimensional (3D), em corpos de água. Sua motivação está centrada no fato de que as margens e zonas costeiras de rios, lagos, estuários, mares e oceanos são locais de aglomerações de seres humanos, dada a sua importância para as atividades econômica, de transporte e de lazer, causando desequilíbrios a esses ecossistemas. Esse fato impulsiona o desenvolvimento de pesquisas relativas a esta temática. Portanto, o objetivo deste trabalho é o de construir um modelo computacional com alta qualidade numérica, que possibilite simular os comportamentos da hidrodinâmica e do transporte escalar de substâncias em corpos de água com complexa configuração geométrica, visando a contribuir para seu manejo racional. Visto que a ênfase nessa tese são os aspectos numéricos e computacionais dos algoritmos, analisaram-se as características e propriedades numérico-computacionais que as soluções devem contemplar, tais como a estabilidade, a monotonicidade, a positividade e a conservação da massa. As estratégias de soluções enfocam os termos advectivos e difusivos, horizontais e verticais, da equação do transporte. Desse modo, a advecção horizontal é resolvida empregando o método da limitação dos fluxos de Sweby, e o transporte vertical (advecção e difusão) é resolvido com os métodos beta de Gross e de Crank-Nicolson. São empregadas malhas com distintas resoluções para a solução do problema multifísica. O esquema numérico resultante é semi-implícito, computacionalmente eficiente, estável e fornece acurácia espacial e temporal de segunda ordem. Os sistemas de equações resultantes da discretização, em diferenças finitas, das equações do escoamento e do transporte 3D, são de grande porte, lineares, esparsos e simétricos definidos-positivos (SDP). No caso 2D os sistemas são lineares, mas os sistemas de equações para a equação do transporte não são simétricos. Assim, para a solução de sistemas de equações SDP e dos sistemas não simétricos empregam-se, respectivamente, os métodos do subespaço de Krylov do gradiente conjugado e do resíduo mínimo generalizado. No caso da solução dos sistemas 3-diagonal, utiliza-se o algoritmo de Thomas e o algoritmo de Cholesky. A solução paralela foi obtida sob duas abordagens. A decomposição ou particionamento de dados, onde as operações e os dados são distribuídos entre os processos disponíveis e são resolvidos em paralelo. E, a decomposição de domínio, onde obtém-se a solução do problema global combinando as soluções de subproblemas locais. Em particular, emprega-se neste trabalho, o método de decomposição de domínio aditivo de Schwarz, como método de solução, e como pré-condicionador. Para maximizar a relação computação/comunicação, visto que a eficiência computacional da solução paralela depende diretamente do balanceamento de carga e da minimização da comunicação entre os processos, empregou-se algoritmos de particionamento de grafos para obter localmente os subproblemas, ou as partes dos dados. O modelo computacional paralelo resultante mostrou-se computacionalmente eficiente e com alta qualidade numérica. / A multi-physics parallel computational model was developed and implemented for the simulation of substance transport and for the two-dimensional (2D) and threedimensional (3D) hydrodynamic flow in water bodies. The motivation for this work is focused in the fact that the margins and coastal zones of rivers, lakes, estuaries, seas and oceans are places of human agglomeration, because of their importance for economic, transport, and leisure activities causing ecosystem disequilibrium. This fact stimulates the researches related to this topic. Therefore, the goal of this work is to build a computational model of high numerical quality, that allows the simulation of hydrodynamics and of scalar transport of substances behavior in water bodies of complex configuration, aiming at their rational management. Since the focuses of this thesis are the numerical and computational aspects of the algorithms, the main numerical-computational characteristics and properties that the solutions need to fulfill were analyzed. That is: stability, monotonicity, positivity and mass conservation. Solution strategies focus on advective and diffusive terms, horizontal and vertical terms of the transport equation. In this way, horizontal advection is solved using Sweby’s flow limiting method; and the vertical transport (advection and diffusion) is solved with Gross and Crank-Nicolson’s beta methods. Meshes of different resolutions are employed in the solution of the multi-physics problem. The resulting numerical scheme is semi-implicit, computationally efficient, stable and provides second order accuracy in space and in time. The equation systems resulting of the discretization, in finite differences, of the flow and 3D transport are of large scale, linear, sparse and symmetric positive definite (SPD). In the 2D case, the systems are linear, but the equation systems for the transport equation are not symmetric. Therefore, for the solution of SPD equation systems and of the non-symmetric systems we employ, respectively, the methods of Krylov’s sub-space of the conjugate gradient and of the generalized minimum residue. In the case of the solution of 3-diagonal systems, Thomas algorithm and Cholesky algorithm are used. The parallel solution was obtained through two approaches. In data decomposition or partitioning, operation and data are distributed among the processes available and are solved in parallel. In domain decomposition the solution of the global problem is obtained combining the solutions of the local sub-problems. In particular, in this work, Schwarz additive domain decomposition method is used as solution method and as preconditioner. In order to maximize the computation/communication relation, since the computational efficiency of the parallel solution depends directly of the load balancing and of the minimization of the communication between processes, graph-partitioning algorithms were used to obtain the sub-problems or part of the data locally. The resulting parallel computational model is computationally efficient and of high numerical quality.

Orientadores: Akebo Yamakami, Ricardo Coelho Silva / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Elétrica e de Computação / Made available in DSpace on 2018-08-19T06:43:21Z (GMT). No. of bitstreams: 1 Ciappina_JussaraRodrigues_D.pdf: 6752680 bytes, checksum: 3dcb7910cc5098e44936e76045792069 (MD5) Previous issue date: 2011 / Resumo: Este trabalho apresenta dois métodos baseados no método clássico de decomposição de Dantzig-Wolfe e um método heurístico, os quais resolvem problemas com incertezas nos parâmetros utilizando a teoria dos conjuntos fuzzy. O primeiro aborda incertezas somente nos custos, enquanto que, os outros dois abordam incertezas nos custos e nas restrições. Métodos que utilizam técnicas de decomposição são indicados para resolver problemas de grande porte que apresentam uma estrutura especial em uma parte do conjunto das restrições. Um exemplo de problema que apresenta tal estrutura é o problema de fluxo multiproduto. Este problema pode ser modelado através de um grafo, cujos nós representam pontos de oferta, demanda e passagem de produtos que trafegam pelos arcos da rede. O objetivo é determinar o fluxo de cada produto nos arcos, de modo a atender a demanda a um custo mínimo, respeitando as restrições de capacidade dos arcos e as restrições de conservação de fluxo dos nós. Com exceção do terceiro, os demais métodos propostos neste trabalho não se limitam a resolver problemas de fluxo multiproduto fuzzy, também resolvem problemas de programação linear fuzzy que apresentam uma estrutura especial em uma parte do conjunto das restrições / Abstract: In this work we present two methods based in the classical Dantzig-Wolfe decomposition and a heuristic method, which solve problems with uncertainties in the parameters using the theory of fuzzy sets. The first one deals with uncertainties only in costs, while the others two deal with uncertainties in costs and restrictions. Methods using decomposition techniques address problems that have a special structure in the set of restrictions. An example of such a problem that has this structure is the fuzzy multicommodity flow problem. This problem can be modeled by a graph whose nodes represent points of supply, demand and passage of commodities that travel by the arcs of the network. The objective is to determine the flow of each commodity in the arcs, in order to meet demand at a minimal cost while respecting the capacity restrictions of the arcs and the flow conservation restrictions of the nodes. With the exception of the third, the other methods proposed in this work are not limited to solve fuzzy multicommodity flow problems, also solve fuzzy linear programming problems that have a special structure in a part of the set of restrictions / Doutorado / Automação / Doutor em Engenharia Elétrica

Ce travail concerne l’étude d’écoulements incompressibles soumis à la rotation avec un solveur haute précision dans des géométries semi-complexes. La technique numérique mise en œuvre combine des schémas compacts, une méthode de projection multi domaine directe et un traitement efficace de la singularité à l’axe basé sur des conditions de parité dans l’espace de Fourier. Le solveur a été parallélisé avec une approche hybride MPI-OpenMP pour réduire les temps de calcul. Dans un premier temps, les précisions spatiales et temporelles de la méthode numérique et la scalabilité du solveur ont été vérifiées. La capacité du solveur à traiter des écoulements plus complexes a été évaluée en considérant des écoulements de type éclatement tourbillonnaire et un écoulement turbulent en conduite cylindrique. Dans un second temps, plusieurs écoulements typiques des machines tournantes ont été étudiés. Le premier écoulement est un écoulement turbulent incompressible isotherme dans un étage simplifié d’un compresseur haute pression d’une turbine à gaz. Les simulations menées ont mises en évidence l’effet de la rotation sur l’écoulement, notamment sur les instabilités se développant le long des parois et sur les différentes structures cohérentes. Le second cas traité est un écoulement turbulent de jet impactant un disque en rotation avec un fort confinement et transfert thermique. Une attention particulière a été portée sur les champs hydrodynamiques et thermique le long du rotor. Enfin, une étude préliminaire d’un jet turbulent impactant un disque fixe d’épaisseur non nulle dans une configuration moins contrainte avec prise en compte du couplage conduction-convection a été réalisée. / This work deals with the study of rotating incompressible flows with a high accurate solver in semi complex geometries. The numerical method used in this work combines compact schemes, a direct multidomain projection method and an efficient axis treatment based on parity conditions in Fourier space. The use of cylindrical coordinates introduces this mathematical singularity. In order to reduce the calculation time, the solver was parallelized with an hybrid MPI-OpenMP parallelization. First, the spatial and temporal accuracies of the numerical method and the scalability of the solver were checked. Then, the capability of the algorithm to deal with more complex flows was verified. Vortex breakdown flows and turbulent pipe flow were studied. In the second step, typical flows of turbomachineries and rotating systems were considered. The first flow was an incompressible isothermal turbulent flow in a high pressure compressor of gas turbine. The different simulations highlighted the rotation effects on the flows, especially on the instabilities appearing along the walls and the coherent structures. The second considered flow was a turbulent impinging jet on a rotating disk with heat transfer in a small aspect ratio cavity. The hydrodynamic fields and heat transfer near the rotor were analyzed in detail. Finally, a preliminary investigation of an impinging jet on a non-rotating disk in a larger aspect ratio cavity with a coupling between conduction and convection transfer was carried out.

Changes in mortality behavior in 1950s-1970s pose a fluctuation of Czech population's mortality trend and still have an impact on current level of life expectancy. The thesis focuses on mortality process in that period from cross-sectional and longitudinal approach and tries to identify how the mortality changes affected age intervals older than 30 and influenced cohorts which were in these age intervals in observed period. In the thesis Arriaga's decomposition method was used for quantification of mortality changes in particular age intervals and their affect on changes of life expectancy. The obtained results show that mostly age groups older than 60 were more sensitive to mortality fluctuations.

This curriculum project provides instructional activities and materials to be used in teaching decomposition and flowcharting skills to students in grades 5-7. It is indicated and supported by research that the use of a variety of activities geared toward problem solving are of great value to students who will make use of decomposition and flowcharting. This project is specifically designed for teachers in grades 5-7 who wish to prepare their students for programming activities on the microcomputer. The criteria for selection of activities and materials are included. This project includes activities categorized under one of the following headings: sequencing activities, selecting information activities, flowcharting activities, and working with computer activities. Finally, the materials and activities used in this project are listed by source so that teachers might find other references that could be helpful in teaching decomposition and flowcharting to students in grades 5-7.

This master's thesis deals with solving fractional-order ordinary differential equations by the Adomian decomposition method. A part of the work is therefore devoted to the theory of equations containing differential operators of non-integer order, especially the Caputo operator. The next part is devoted to the Adomian decomposition method itself, its properties and implementation in the case of Chen system. The work also deals with bifurcation analysis of this system, both for integer and non-integer case. One of the objectives is to clarify the discrepancy in the literature concerning the fractional-order Chen system, where experiments based on the use of the Adomian decomposition method give different results for certain input parameters compared with numerical methods. The clarification of this discrepancy is based on recent theoretical knowledge in the field of fractional-order differential equations and their systems. The conclusions are supported by numerical experiments, own code implementing the Adomian decomposition method on the Chen system was used.