Browsing by Subject "Navier-Stokes"
Now showing 1 - 5 of 5
Results Per Page
Sort Options
Item Discontinuous Galerkin Finite Element Method for the Nonlinear Hyperbolic Problems with Entropy-Based Artificial Viscosity Stabilization(2012-07-16) Zingan, Valentin NikolaevichThis work develops a discontinuous Galerkin finite element discretization of non- linear hyperbolic conservation equations with efficient and robust high order stabilization built on an entropy-based artificial viscosity approximation. The solutions of equations are represented by elementwise polynomials of an arbitrary degree p > 0 which are continuous within each element but discontinuous on the boundaries. The discretization of equations in time is done by means of high order explicit Runge-Kutta methods identified with respective Butcher tableaux. To stabilize a numerical solution in the vicinity of shock waves and simultaneously preserve the smooth parts from smearing, we add some reasonable amount of artificial viscosity in accordance with the physical principle of entropy production in the interior of shock waves. The viscosity coefficient is proportional to the local size of the residual of an entropy equation and is bounded from above by the first-order artificial viscosity defined by a local wave speed. Since the residual of an entropy equation is supposed to be vanishingly small in smooth regions (of the order of the Local Truncation Error) and arbitrarily large in shocks, the entropy viscosity is almost zero everywhere except the shocks, where it reaches the first-order upper bound. One- and two-dimensional benchmark test cases are presented for nonlinear hyperbolic scalar conservation laws and the system of compressible Euler equations. These tests demonstrate the satisfactory stability properties of the method and optimal convergence rates as well. All numerical solutions to the test problems agree well with the reference solutions found in the literature. We conclude that the new method developed in the present work is a valuable alternative to currently existing techniques of viscous stabilization.Item Estimates on higher derivatives for the Navier-Stokes equations and Hölder continuity for integro-differential equations(2012-08) Choi, Kyudong; Vasseur, Alexis F.; Caffarelli, Luis; Figalli, Alessio; Gamba, Irene; Morrison, Philip; Pavlovic, NatasaThis thesis is divided into two independent parts. The first part concerns the 3D Navier-Stokes equations. The second part deals with regularity issues for a family of integro-differential equations. In the first part of this thesis, we consider weak solutions of the 3D Navier-Stokes equations with L² initial data. We prove that ([Nabla superscript alpha])u is locally integrable in space-time for any real [alpha] such that 1 < [alpha] < 3. Up to now, only the second derivative ([Nabla]²)u was known to be locally integrable by standard parabolic regularization. We also present sharp estimates of those quantities in local weak-L[superscript (4/([alpha]+1))]. These estimates depend only on the L² norm of the initial data and on the domain of integration. Moreover, they are valid even for [alpha] ≥ 3 as long as u is smooth. The proof uses a standard approximation of Navier-Stokes from Leray and blow-up techniques. The local study is based on De Giorgi techniques with a new pressure decomposition. To handle the non-locality of fractional Laplacians, Hardy space and Maximal functions are introduced. In the second part of this thesis, we consider non-local integro-differential equations under certain natural assumptions on the kernel, and obtain persistence of Hölder continuity for their solutions. In other words, we prove that a solution stays in C[superscript beta] for all time if its initial data lies in C[superscript beta]. Also, we prove a C[superscript beta]-regularization effect from [mathematical equation] initial data. It provides an alternative proof to the result of Caffarelli, Chan and Vasseur [10], which was based on De Giorgi techniques. This result has an application for a fully non-linear problem, which is used in the field of image processing. In addition, we show Hölder regularity for solutions of drift diffusion equations with supercritical fractional diffusion under the assumption [mathematical equation]on the divergent-free drift velocity. The proof is in the spirit of Kiselev and Nazarov [48] where they established Hölder continuity of the critical surface quasi-geostrophic (SQG) equation by observing the evolution of a dual class of test functions.Item On the linear stability problem for Jeffery-Hamel flows(2015-05) Carlson, William Zechariah; Vishik, Mikhail; Chen, Thomas; Pavlovic, Natasa; Vasseur, Alexis; Bogard, DavidWe study the linear stability of a family of Jeffery-Hamel solutions which satisfy a zero flux condition. With a suitable regularization of these velocity profiles we show that the linearized perturbation equation is well-posed on a weighted L² space with a certain class of radial weights, in the example of a half plane or in the whole plane. We prove that the perturbed Stokes operator of this system is the generator of a strongly continuous analytic semigroup. We also describe some formal asymptotics under which the linear stability problem could be reduced to a one dimensional problem for which we state a formal perturbation theory.Item Results towards a Scalable Multiphase Navier-Stokes Solver for High Reynolds Number Flows(2013-07-19) Thompson, Travis BrandonThe incompressible Navier-Stokes equations have proven formidable for nearly a century. The present difficulties are mathematical and computational in nature; the computational requirements, in particular, are exponentially exacerbated in the presence of high Reynolds number. The issues are further compounded with the introduction of markers or an immiscible fluid intended to be tracked in an ambient high Reynolds number flow; despite the overwhelming pragmatism of problems in this regime, and increasing computational efficacy, even modest problems remain outside the realm of direct approaches. Herein three approaches are presented which embody direct application to problems of this nature. An LES model based on an entropy-viscosity serves to abet the computational resolution requirements imposed by high Reynolds numbers and a one-stage compressive flux, also utilizing an entropy-viscosity, aids in accurate, efficient, conservative transport, free of low order dispersive error, of an immiscible fluid or tracer. Finally, an integral commutator and the theory of anti-dispersive spaces is introduced as a novel theoretical tool for consistency error analysis; in addition the material engenders the construction of error-correction techniques for mass lumping schemes.Item Space-time discontinuous Petrov-Galerkin finite elements for transient fluid mechanics(2016-05) Ellis, Truman Everett; Demkowicz, Leszek; Moser, Robert D; Hughes, Thomas J.R; Dawson, Clint N; Bui, TanInitial mesh design for computational fluid dynamics can be a time-consuming and expensive process. The stability properties and nonlinear convergence of most numerical methods rely on a minimum level of mesh resolution. This means that unless the initial computational mesh is fine enough, convergence can not be guaranteed. Any meshes below this minimum resolution level are termed to be in the ``pre-asymptotic regime.'' This condition implies that meshes need to in some way anticipate the solution before it is known. On top of the minimum requirement that the surface meshes must adequately represent the geometry of the problem under consideration, resolution requirements on the volume mesh make the CFD practitioner's job significantly more time consuming. In contrast to most other numerical methods, the discontinuous Petrov-Galerkin finite element method retains exceptional stability on extremely coarse meshes. DPG is also inherently very adaptive. It is possible to compute the residual error without knowledge of the exact solution, which can be used to robustly drive adaptivity. This results in a very automated technology, as the user can initialize a computation on the coarsest mesh which adequately represents the geometry then step back and let the program solve and adapt iteratively until it resolves the solution features. A common complaint of minimum residual methods by computational fluid dynamics practitioners is that they are not locally conservative. In this thesis, this concern is addressed by developing a locally conservative DPG formulation by augmenting the system with Lagrange multipliers. The resulting DPG formulation is then proved to be robust and shown to produce superior numerical results over standard DPG on a selection of test problems. Adaptive convergence to steady incompressible and compressible Navier-Stokes solutions was explored in Jesse Chan's and Nathan Roberts' dissertations. Space-time offers a natural extension to transient problems as it preserves the stability and adaptivity properties of DPG in the time dimension. Space-time also offers more extensive parallelization capability than problems treated with traditional time stepping as it allows multigrid concurrently in both space and time. A proof of concept space-time DPG formulation is developed for transient convection-diffusion. The robust test norms derived for steady convection-diffusion are extended to the space-time case and proofs of robustness are provided. Numerical results verify the robust behavior and near $L^2$ optimality of the resulting solutions. The space-time formulation for convection-diffusion is then extended to transient incompressible and compressible Navier-Stokes by analogy. Several numerical experiments are performed, but a mathematical analysis is not attempted for these nonlinear problems. Several side topics are explored such as a study of the compressible Navier-Stokes equations under various variable transformations and the development of consistent test norms through the concept of physical entropy.