We present a method for stable buildup of snow on surfaces of arbitrary topology and geometric complexity. This is achieved by tracing quantities of snow, so-called snow packages, through a dynamic wind field. Dual compact level sets are used to represent geometry as well as accumulated snow. The level sets have also proven to be well suited for the internal boundaries for our Navier-Stokes solver, which produces a wind field that changes according to snow buildup. Our method is different from previous work in that all the addition of snow is done by local operations, avoiding computationally expensive global refinement procedures. The main contribution of this work is a dual level set method for particle interaction with level sets.

In this work, based on the Law of conservation of mass and the Newton s second Law, we deducted the equations of Navier-Stokes for the uncompressing fluids. Then we have developed in a detailed way the varying form of these equations. For this formulation, through the Galerkin method, we have proved the existence of weak solutions in Sobolev spaces in a bounded domain whit regular boundary. / Neste trabalho, a partir da lei da conservação da massa e da segunda lei de Newton, deduzimos as equações de Navier-Stokes para um fluido incompress nível. Em seguida fazemos, de forma detalhada, a formulação variacional destas equações. Para esta formulação, através do método de Galerkin, provamos a existência de solução fraca em espaços de Sobolev em um domínio limitado cuja fronteira seja regular.

No description available.

Orientador: Marcelo Martins dos Santos / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Matematica, Estatística e Computação Científica / Made available in DSpace on 2018-08-14T14:51:49Z (GMT). No. of bitstreams: 1 Teixeira_EdsonJose_M.pdf: 638751 bytes, checksum: 1d26a9bbc1ee3ba6c4ee45e29c14c45e (MD5) Previous issue date: 2009 / Resumo: Este trabalho consiste de uma exposição detalhada do resultado provado no artigo "Global weak solutions to 1D compressible isentropic Navier-Stokes equations with density-dependent viscosity" de S. Jiang, Z. P. Xin e P. Zhang (Methods Appl. Anal. - 2005), sobre a existência e unicidade de solução fraca para o sistema de Navier-Stokes unidimensional de um fluido isentrópico compressível com viscosidade dependente da densidade e com fronteira livre em coordenadas lagrangianas, ?t +?2ux = 0 0 < x < 1, t > 0 ut + (P(?))x = (?µ (?)ux)x 0 < x < 1, t > 0 onde ?, u; P(?) e µ(?) são a densidade, velocidade, pressão e viscosidade do fluido, e exigiremos que este fluido satisfaça a condição de fronteira (-P(?) + (?µ(?)ux)= 0. Trataremos do caso particular onde consideramos P(?) = A ?? e µ( ?) = B?a; onde A, B > 0,? > 1 e 0 < a < 1 são constantes. Acrescentaremos uma condicão inicial (?0,u0). / Abstract: The present work makes a well-detailed exposition about the main results given in the paper "Global weak solutions to 1D compressible isentropic Navier-Stokes equations with densitydependent viscosity" by S. Jiang, Z. P. Xin and P. Zhang (Methods Appl. Anal. - 2005). The problem in this paper has a free boundary but in lagrangian coordinates the equations are the following, ?t +?2ux = 0 0 < x < 1, t > 0 ut + (P(?))x = (?µ (?)ux)x 0 < x < 1, t > 0 and the boundary becomes the fixed points x = 0 and x = 1; Here ?, u; P(?) and µ(?) are, respectively, the density, velocity, pressure and the viscosity of the fluid. The boundary condition, at x = 0 and x = 1, is given by (-P(?) + (?µ(?)ux)= 0. Although the pressure and viscosity may have more general forms, to be more specific, the authors consider only the special case P(?) = A ?? e µ( ?) = B?a, with A; B > 0,? > 1 and 0 <a< 1 being constants. An initial condition (?0,u0) is also given at time t = 0. / Mestrado / Analise, Equações Diferenciais Parciais / Mestre em Matemática

Orientador: Marko A. Rojas Medar / Tese (doutorado) - Universidade Estadual de Campinas, Instituto de Matematica, Estatistica e Computação Cientifica / Made available in DSpace on 2018-07-23T21:34:49Z (GMT). No. of bitstreams: 1 OrtegaTorres_ElvaEliana_D.pdf: 5966361 bytes, checksum: f08e9f1f1c04a07735054d7ca6ea4112 (MD5) Previous issue date: 1998 / Resumo: Não informado / Abstract: Not informed / Doutorado / Doutor em Matemática Aplicada

Um problema relevante na modelagem de escoamentos com movimento de corpo rígido consiste na consideração de forças externas bem como forças que o fluido exerce sobre o próprio corpo rígido. Um outro problema importante refere-se à descrição da trajetória dos corpos rígidos. Essas duas questões são o objeto de estudo deste trabalho. No sentido de solucioná-las, foram implementadas duas extensões ao modelador de movimentos do sistema Freeflow-3D. A primeira incorpora um tipo de movimento definido por forças externas e a segunda um tipo de movimento definido por interpolação linear por partes. / A relevant problem in the modeling of the flow of fluids with rigid body movements consists in the consideration of externai forces, as well as forces that the fluids apply in the rigid body. Another important problem in this domain refers to the description of the path of the rigid bodies. Studying these two issues was the goal of this work. Towards solving these problems, two extensions were implemented in the movement modeler of the Freeflow-3D system. The first one, adds to the modeler a new type of movement, defined by externai forces. The second, adds another type of movement defined by piecewise linear interpolation.

Moment methods derived from the kinetic theory of gases can be used for the prediction of continuum and non-equilibrium flows and offer numerical advantages over other methods, such as the Navier-Stokes model. Models developed in this fashion are described by first-order hyperbolic partial differential equations (PDEs) with stiff local relaxation source terms. The application of discontinuous-Galerkin (DG) methods for the solution of such models has many benefits. Of particular interest is the third-order accurate, coupled space-time discontinuous-Galerkin Hancock (DGH) method. This scheme is accurate, as well as highly efficient on large-scale distributed-memory computers. The current study outlines a general implementation of the DGH method used for the parallel solution of moment methods in one, two, and three dimensions on modern distributed clusters. An algorithm for adaptive mesh refinement (AMR) was developed alongside the implementation of the scheme, and is used to achieve even higher accuracy and efficiency. Many different first-order hyperbolic and hyperbolic-relaxation PDEs are solved to demonstrate the robustness of the scheme. First, a linear convection-relaxation equation is solved to verify the order of accuracy of the scheme in three dimensions. Next, some classical compressible Euler problems are solved in one, two, and three dimensions to demonstrate the scheme's ability to capture discontinuities and strong shocks, as well as the efficacy of the implemented AMR. A special case, Ringleb's flow, is also solved in two-dimensions to verify the order of accuracy of the scheme for non-linear PDEs on curved meshes. Following this, the shallow water equations are solved in two dimensions. Afterwards, the ten-moment (Gaussian) closure is applied to two-dimensional Stokes flow past a cylinder, showing the abilities of both the closure and scheme to accurately compute classical viscous solutions. Finally, the one-dimensional fourteen-moment closure is solved.

In this thesis a method for studying the asymptotic behavior of solutions to dissipative partial differential equations is developed, motivated by the study of the compressible Navier-Stokes equations in the past works of Hoff and Zumbrun,1995, Hoff and Zumbrun, 1997. In its most basic form, this method allows one to compute n^th order approximations in terms of Hermite functions of solutions of the heat equation having n^th order moments. The main advantage is that these approximations can be efficiently computed, and are often given explicitly in terms of elementary functions. It is shown how this method can be extended to increasingly complicated systems, leading the way toward the asymptotic analysis of the compressible Navier-Stokes equations. A number of challenges must be overcome to apply this method to the compressible Navier-Stokes system. For technical reasons, the analysis is carried out on the divergence and curl of the velocity field, and hence a means of recovering the velocity field from these quantities is established first. The linear part of the evolution is then studied, and an extended version of the artificial viscosity decomposition previously developed (Kawashima, Hoff and Zumbrun1995) is introduced. This decomposition is in terms of the heat and combined heat-wave operators, and hence general estimates on their evolution in weighted L^p spaces are obtained. A modified compressible Navier-Stokes system is then introduced which captures the dominant behavior of the linear evolution and possesses similar nonlinear terms. Solutions to this modified system are proven to exist in weighted spaces, showing that solutions initially having a certain number of moments possess this same number of moments for all time. An analysis of the asymptotic behavior of the modified compressible Navier-Stokes system is then carried out, and it is shown that the method developed herein extends and unifies the approach of Hoff and Zumbrun with that of Gallay and Wayne, 2002a, Gallay and Wayne, 2002b, where it was originally developed to study the behavior of the incompressible Navier-Stokes equations. The thesis is concluded with a discussion of how the results obtained for the modified compressible Navier-Stokes system pave the way for an analysis of the true compressible Navier-Stokes system, the generalization of this asymptotic analysis to arbitrary order, and with a comparison of this asymptotic analysis to that found in the recent work of Kagei and Okita, 2017.

Using the Navier-Stokes equations along with a continuity equation, a one-dimensional model is developed to describe the spin coating process of thin polymeric films. The resulting model is a system of a parabolic partial differential equation coupled with an integral equation as well as with an ordinary differential equation describing the motion of a moving boundary. Viscosity and diffusivity are allowed to be varied in the model. To be able to perform the finite element approximation of the model equations, the moving boundary is fixed. Then the finite element method is applied along with the so called Method of Lines resulting in a semi-discrete problem, a large system of ordinary differential equations which is then solved with MATLAB. We present an existence and uniqueness result what concerns the semi-discrete solutions. Finally, we illustrate numerically the behavior of the solutions to our model.

We deal with the numerical solution of the Navier-Stokes equations describing a motion of viscous compressible flows. The governing equations are discretized with the aid of discontinuous Galerkin finite element method which is based on a discontinuous piecewise polynomial approximation. The discretizations leads to a large nonlinear algebraic system. In order to solve this system efficiently, we develop the so-called p-multigrid solution strategy which employ as a projec- tion and a restriction operators the L2 -projection in the spaces of polynomial functions on each element separately. The p-multigrid technique is studied, deve- loped and implemented in the code ADGFEM. The computational performance of the method is presented.