Spelling suggestions: "subject:"discretization schemes"" "subject:"iscretization schemes""
1 |
Symplectic integration of constrained Hamiltonian systemsLeimkuhler, Benedict, Reich, Sebastian January 1994 (has links)
A Hamiltonian system in potential form (formula in the original abstract) subject to smooth constraints on q can be viewed as a Hamiltonian system on a manifold, but numerical computations must be performed in Rn. In this paper methods which reduce "Hamiltonian differential algebraic equations" to ODEs in Euclidean space are examined. The authors study the construction of canonical parameterizations or local charts as well as methods based on the construction of ODE systems in the space in which the constraint manifold is embedded which preserve the constraint manifold as an invariant manifold. In each case, a Hamiltonian system of ordinary differential equations is produced. The stability of the constraint invariants and the behavior of the original Hamiltonian along solutions are investigated both numerically and analytically.
|
2 |
Efeitos numéricos na simulação de escoamentos gás-sólido em leito fuidizado borbulhante utilizando a teoria cinética dos escoamentos granulares /Souza, Meire Pereira de. January 2009 (has links)
Orientador: Hélio Aparecido Navarro / Banca: Vicente Luiz Scalon / Banca: Luben Cabezas Gomez / Resumo: No presente trabalho desenvolve-se um estudo de modelagem matemática e simulação numérica do escoamento bifásico gás-sólido em leito fluidizado borbulhante. Utiliza-se o modelo Euleriano de duas fases separadas formulando o tensor das tensões da fase sólida através da teoria cinética dos escoamentos granulares. As simulações numéricas são realizadas através do código fonte MFIX (Multiphase Flow with Interphase eXchanges) desenvolvido no NETL (National Energy Technology Laboratory). Os resultados de simulação numérica são avaliados por meio da análise da influência dos seguintes parâmetros: malha computacional e esquemas de discretização dos termos convectivos das equações de conservação. Com base nos estudos teóricos e resultados obtidos durante o trabalho conclui-se que esquemas de primeira, tais como FOUP são altamente difusivos, já os resultados obtidos utilizando o esquema de alta ordem, Superbee, produziu resultados de melhor qualidade para as malhas testadas neste trabalho. Além disso, os resultados mostraram-se bastante dependentes do tamanho da malha computacional. / Abstract: In the present work is described a mathematical model and numerical simulation of gas-solid two-phase flow in a bubbling fluidized bed. It is used the Eulerian gas-solid two-fluid model and the solid phase stress tensor is modeled considering the kinetic theory of granular flows. The numerical simulations were developed using the MFIX (Multiphase Flow with Interphase eXchanges) code developed in NETL (National Energy Technology Laboratory). The numerical diffusion is analyzed considering a single bubbling detachment and its movement process in a two-dimensional bubbling fluidized bed using the bubble shape as a metric for results description. The influence of computacional grid it is also analyzed. It is concluded that SuperBee scheme produces the better results and analysis about estimating uncertainty in grid refinement should be studied. / Mestre
|
3 |
Étude et modélisation des équations différentielles stochastiques / High weak order discretization schemes for stochastic differential equationRey, Clément 04 December 2015 (has links)
Durant les dernières décennies, l'essor des moyens technologiques et particulièrement informatiques a permis l'émergence de la mise en œuvre de méthodes numériques pour l'approximation d'Equations Différentielles Stochastiques (EDS) ainsi que pour l'estimation de leurs paramètres. Cette thèse aborde ces deux aspects et s'intéresse plus spécifiquement à l'efficacité de ces méthodes. La première partie sera consacrée à l'approximation d'EDS par schéma numérique tandis que la deuxième partie traite l'estimation de paramètres. Dans un premier temps, nous étudions des schémas d'approximation pour les EDSs. On suppose que ces schémas sont définis sur une grille de temps de taille $n$. On dira que le schéma $X^n$ converge faiblement vers la diffusion $X$ avec ordre $h in mathbb{N}$ si pour tout $T>0$, $vert mathbb{E}[f(X_T)-f(X_T^n)] vertleqslant C_f /n^h$. Jusqu'à maintenant, sauf dans certains cas particulier (schémas d'Euler et de Ninomiya Victoir), les recherches sur le sujet imposent que $C_f$ dépende de la norme infini de $f$ mais aussi de ses dérivées. En d'autres termes $C_f =C sum_{vert alpha vert leqslant q} Vert partial_{alpha} f Vert_{ infty}$. Notre objectif est de montrer que si le schéma converge faiblement avec ordre $h$ pour un tel $C_f$, alors, sous des hypothèses de non dégénérescence et de régularité des coefficients, on peut obtenir le même résultat avec $C_f=C Vert f Vert_{infty}$. Ainsi, on prouve qu'il est possible d'estimer $mathbb{E}[f(X_T)]$ pour $f$ mesurable et bornée. On dit alors que le schéma converge en variation totale vers la diffusion avec ordre $h$. On prouve aussi qu'il est possible d'approximer la densité de $X_T$ et ses dérivées par celle $X_T^n$. Afin d'obtenir ce résultat, nous emploierons une méthode de calcul de Malliavin adaptatif basée sur les variables aléatoires utilisées dans le schéma. L'intérêt de notre approche repose sur le fait que l'on ne traite pas le cas d'un schéma particulier. Ainsi notre résultat s'applique aussi bien aux schémas d'Euler ($h=1$) que de Ninomiya Victoir ($h=2$) mais aussi à un ensemble générique de schémas. De plus les variables aléatoires utilisées dans le schéma n'ont pas de lois de probabilité imposées mais appartiennent à un ensemble de lois ce qui conduit à considérer notre résultat comme un principe d'invariance. On illustrera également ce résultat dans le cas d'un schéma d'ordre 3 pour les EDSs unidimensionnelles. La deuxième partie de cette thèse traite le sujet de l'estimation des paramètres d'une EDS. Ici, on va se placer dans le cas particulier de l'Estimateur du Maximum de Vraisemblance (EMV) des paramètres qui apparaissent dans le modèle matriciel de Wishart. Ce processus est la version multi-dimensionnelle du processus de Cox Ingersoll Ross (CIR) et a pour particularité la présence de la fonction racine carrée dans le coefficient de diffusion. Ainsi ce modèle permet de généraliser le modèle d'Heston au cas d'une covariance locale. Dans cette thèse nous construisons l'EMV des paramètres du Wishart. On donne également la vitesse de convergence et la loi limite pour le cas ergodique ainsi que pour certains cas non ergodiques. Afin de prouver ces convergences, nous emploierons diverses méthodes, en l'occurrence : les théorèmes ergodiques, des méthodes de changement de temps, ou l'étude de la transformée de Laplace jointe du Wishart et de sa moyenne. De plus, dans dernière cette étude, on étend le domaine de définition de cette transformée jointe / The development of technology and computer science in the last decades, has led the emergence of numerical methods for the approximation of Stochastic Differential Equations (SDE) and for the estimation of their parameters. This thesis treats both of these two aspects. In particular, we study the effectiveness of those methods. The first part will be devoted to SDE's approximation by numerical schemes while the second part will deal with the estimation of the parameters of the Wishart process. First, we focus on approximation schemes for SDE's. We will treat schemes which are defined on a time grid with size $n$. We say that the scheme $ X^n $ converges weakly to the diffusion $ X $, with order $ h in mathbb{N} $, if for every $ T> 0 $, $ vert mathbb{E} [f (X_T) -f (X_T^n)]vert leqslant C_f / h^n $. Until now, except in some particular cases (Euler and Victoir Ninomiya schemes), researches on this topic require that $ C_f$ depends on the supremum norm of $ f $ as well as its derivatives. In other words $C_f =C sum_{vert alpha vert leqslant q} Vert partial_{alpha} f Vert_{ infty}$. Our goal is to show that, if the scheme converges weakly with order $ h $ for such $C_f$, then, under non degeneracy and regularity assumptions, we can obtain the same result with $ C_f=C Vert f Vert_{infty}$. We are thus able to estimate $mathbb{E} [f (X_T)]$ for a bounded and measurable function $f$. We will say that the scheme converges for the total variation distance, with rate $h$. We will also prove that the density of $X^n_T$ and its derivatives converge toward the ones of $X_T$. The proof of those results relies on a variant of the Malliavin calculus based on the noise of the random variable involved in the scheme. The great benefit of our approach is that it does not treat the case of a particular scheme and it can be used for many schemes. For instance, our result applies to both Euler $(h = 1)$ and Ninomiya Victoir $(h = 2)$ schemes. Furthermore, the random variables used in this set of schemes do not have a particular distribution law but belong to a set of laws. This leads to consider our result as an invariance principle as well. Finally, we will also illustrate this result for a third weak order scheme for one dimensional SDE's. The second part of this thesis deals with the topic of SDE's parameter estimation. More particularly, we will study the Maximum Likelihood Estimator (MLE) of the parameters that appear in the matrix model of Wishart. This process is the multi-dimensional version of the Cox Ingersoll Ross (CIR) process. Its specificity relies on the square root term which appears in the diffusion coefficient. Using those processes, it is possible to generalize the Heston model for the case of a local covariance. This thesis provides the calculation of the EMV of the parameters of the Wishart process. It also gives the speed of convergence and the limit laws for the ergodic cases and for some non-ergodic case. In order to obtain those results, we will use various methods, namely: the ergodic theorems, time change methods or the study of the joint Laplace transform of the Wishart process together with its average process. Moreover, in this latter study, we extend the domain of definition of this joint Laplace transform
|
4 |
Efeitos numéricos na simulação de escoamentos gás-sólido em leito fuidizado borbulhante utilizando a teoria cinética dos escoamentos granularesSouza, Meire Pereira de [UNESP] 12 January 2009 (has links) (PDF)
Made available in DSpace on 2014-06-11T19:26:18Z (GMT). No. of bitstreams: 0
Previous issue date: 2009-01-12Bitstream added on 2014-06-13T20:54:32Z : No. of bitstreams: 1
souza_mp_me_bauru.pdf: 1121017 bytes, checksum: 8ed8fe65e5cd46111ec16064bc42b009 (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / No presente trabalho desenvolve-se um estudo de modelagem matemática e simulação numérica do escoamento bifásico gás-sólido em leito fluidizado borbulhante. Utiliza-se o modelo Euleriano de duas fases separadas formulando o tensor das tensões da fase sólida através da teoria cinética dos escoamentos granulares. As simulações numéricas são realizadas através do código fonte MFIX (Multiphase Flow with Interphase eXchanges) desenvolvido no NETL (National Energy Technology Laboratory). Os resultados de simulação numérica são avaliados por meio da análise da influência dos seguintes parâmetros: malha computacional e esquemas de discretização dos termos convectivos das equações de conservação. Com base nos estudos teóricos e resultados obtidos durante o trabalho conclui-se que esquemas de primeira, tais como FOUP são altamente difusivos, já os resultados obtidos utilizando o esquema de alta ordem, Superbee, produziu resultados de melhor qualidade para as malhas testadas neste trabalho. Além disso, os resultados mostraram-se bastante dependentes do tamanho da malha computacional. / In the present work is described a mathematical model and numerical simulation of gas-solid two-phase flow in a bubbling fluidized bed. It is used the Eulerian gas-solid two-fluid model and the solid phase stress tensor is modeled considering the kinetic theory of granular flows. The numerical simulations were developed using the MFIX (Multiphase Flow with Interphase eXchanges) code developed in NETL (National Energy Technology Laboratory). The numerical diffusion is analyzed considering a single bubbling detachment and its movement process in a two-dimensional bubbling fluidized bed using the bubble shape as a metric for results description. The influence of computacional grid it is also analyzed. It is concluded that SuperBee scheme produces the better results and analysis about estimating uncertainty in grid refinement should be studied.
|
5 |
Modelagem euleriana do escoamento gás-sólido em leito fluidizado circulante: análise da influência de parâmetros físicos e numéricos nos resultados de simulação / Eulerian modeling of the gas-solid flow in a circulating fluidized bed: analysis of the physical and numerical parameters influence in the simulation resultsSilva, Renato César da 03 February 2006 (has links)
No presente trabalho desenvolve-se um estudo de modelagem matemática e simulação numérica do escoamento bifásico gás-sólido na coluna ascendente de um leito fluidizado circulante. Utiliza-se o modelo euleriano de duas fases separadas considerando dois procedimentos diferentes para a modelagem do tensor das tensões da fase sólida: modelo tradicional e a teoria cinética dos escoamentos granulares (TCEG). As simulações numéricas são conduzidas com a utilização do código MFIX que é um software livre e disponível na rede (Internet). Os resultados da simulação numérica são avaliados por meio da análise da influência dos seguintes parâmetros: malha computacional, correlações para o computo do tensor das tensões da fase sólida e esquemas de discretização dos termos advectivos. Também se desenvolve estudo de caracterização de estruturas coerentes - \"clusters\". De forma complementar foram realizadas duas análises teóricas compreendendo: uma análise da influência das diversas correlações utilizadas na TCEG para o computo da viscosidade dinâmica do sólido; e uma análise enfocando o emprego de diversos esquemas de discretização para os termos advectivos presentes nas equações de conservação (Foup, Muscl, Van Leer, Minmod e Superbee). De todos os estudos e resultados apresentados no trabalho conclui-se que os escoamentos gás-sólido em leitos fluidizados circulantes são muito complexos, sendo necessário a realização de futuras pesquisas para uma melhor compreensão dos fenômenos físicos inerentes a esses escoamentos. / In the present work is described a mathematical model and numerical study simulation of the gas-solid flow in the riser of a circulating fluidized bed. It is used the two fluids eulerian model considering two different procedures for the solid phase stress tensor modeling: the traditional model and the kinetic theory of granular flows (KTGF). The numerical simulation results are evaluated through the influence analysis of the following parameters: computational mesh, correlations for computing the solid phase stress tensor and the discretization of the advective terms. It is also presented a study concerning the characterization coherent structures - \"clusters\". Complementing the above studies were accomplished two theoretical analyses comprehending: an influence analysis of several correlations used in the KTGF for computing the dynamic viscosity of the solid phase; and an analysis concerning several discretization schemes for the advective terms present in the conservative equations. Considering the developed studies and the obtained results it is concluded that the gas-solid flows in circulating fluidized beds are very complex, being necessary future research works for a better comprehension of the inherent physical phenomena to these flows.
|
6 |
Modelagem euleriana do escoamento gás-sólido em leito fluidizado circulante: análise da influência de parâmetros físicos e numéricos nos resultados de simulação / Eulerian modeling of the gas-solid flow in a circulating fluidized bed: analysis of the physical and numerical parameters influence in the simulation resultsRenato César da Silva 03 February 2006 (has links)
No presente trabalho desenvolve-se um estudo de modelagem matemática e simulação numérica do escoamento bifásico gás-sólido na coluna ascendente de um leito fluidizado circulante. Utiliza-se o modelo euleriano de duas fases separadas considerando dois procedimentos diferentes para a modelagem do tensor das tensões da fase sólida: modelo tradicional e a teoria cinética dos escoamentos granulares (TCEG). As simulações numéricas são conduzidas com a utilização do código MFIX que é um software livre e disponível na rede (Internet). Os resultados da simulação numérica são avaliados por meio da análise da influência dos seguintes parâmetros: malha computacional, correlações para o computo do tensor das tensões da fase sólida e esquemas de discretização dos termos advectivos. Também se desenvolve estudo de caracterização de estruturas coerentes - \"clusters\". De forma complementar foram realizadas duas análises teóricas compreendendo: uma análise da influência das diversas correlações utilizadas na TCEG para o computo da viscosidade dinâmica do sólido; e uma análise enfocando o emprego de diversos esquemas de discretização para os termos advectivos presentes nas equações de conservação (Foup, Muscl, Van Leer, Minmod e Superbee). De todos os estudos e resultados apresentados no trabalho conclui-se que os escoamentos gás-sólido em leitos fluidizados circulantes são muito complexos, sendo necessário a realização de futuras pesquisas para uma melhor compreensão dos fenômenos físicos inerentes a esses escoamentos. / In the present work is described a mathematical model and numerical study simulation of the gas-solid flow in the riser of a circulating fluidized bed. It is used the two fluids eulerian model considering two different procedures for the solid phase stress tensor modeling: the traditional model and the kinetic theory of granular flows (KTGF). The numerical simulation results are evaluated through the influence analysis of the following parameters: computational mesh, correlations for computing the solid phase stress tensor and the discretization of the advective terms. It is also presented a study concerning the characterization coherent structures - \"clusters\". Complementing the above studies were accomplished two theoretical analyses comprehending: an influence analysis of several correlations used in the KTGF for computing the dynamic viscosity of the solid phase; and an analysis concerning several discretization schemes for the advective terms present in the conservative equations. Considering the developed studies and the obtained results it is concluded that the gas-solid flows in circulating fluidized beds are very complex, being necessary future research works for a better comprehension of the inherent physical phenomena to these flows.
|
Page generated in 0.1236 seconds