This work is a numerical study of Burgers' equation with Neumann boundary conditions. The goal is to determine the long term behavior of solutions. We develop and test two separate finite element and Galerkin schemes and then use those schemes to compute the response to various initial conditions and Reynolds numbers. It is known that for sufficiently small initial data, all steady state solutions of Burgers' equation with Neumann boundary conditions are constant. The goal here is to investigate the case where initial data is large. Our numerical results indicate that for certain initial data the solution of Burgers' equation can approach non-constant functions as time goes to infinity. In addition, the numerical results raise some interesting questions about steady state solutions of nonlinear systems. / Master of Science

The Burgers equation, and related partial differential equations (PDEs), can be numerically challenging for small values of the viscosity parameter. For example, these equations can develop discontinuous solutions (or solutions with large gradients) from smooth initial data. Aside from numerical stability issues, standard numerical methods can also give rise to spurious oscillations near these discontinuities. In this study, we consider an equivalent form of the Burgers equation given by Constantin and Iyer, whose solution can be written as the expected value of a stochastic differential equation. This equivalence is used to develop a numerical method for approximating solutions to Burgers equation. Our preliminary analysis of the algorithm reveals that it is a natural generalization of the method of characteristics and that it produces approximate solutions that actually improve as the viscosity parameter vanishes. We present three examples that compare our algorithm to a recently published reference method as well as the vanishing viscosity/entropy solution for decreasing values of the viscosity. / Master of Science / Burgers equation is a Partial Differential Equation (PDE) used to model how fluids evolve in time based on some initial condition and viscosity parameter. This viscosity parameter helps describe how the energy in a fluid dissipates. When studying partial differential equations, it is often hard to find a closed form solution to the problem, so we often approximate the solution with numerical methods. As our viscosity parameter approaches 0, many numerical methods develop problems and may no longer accurately compute the solution. Using random variables, we develop an approximation algorithm and test our numerical method on various types of initial conditions with small viscosity coefficients.

The problem of establishing local existence and uniqueness of solutions to systems of differential equations is well understood and has a long history. However, the problem of proving global existence and uniqueness is more difficult and fails even for some very simple ordinary differential equations. It is still not known if the 3D Navier-Stokes equation have global unique solutions and this open problem is one of the Millennium Prize Problems. However, many of these mathematical models are extremely useful in the understanding of complex physical systems. For years people have considered methods for modifying these equations in order to obtain models that still capture the observed fundamental physics, but for which one can rigorously establish global results. In this thesis we focus on a taming method to achieve this goal and apply taming to modeling and numerical problems. The method is also applied to a class of nonlinear differential equations with conservative nonlinearities and to Burgers’ Equation with Neumann boundary conditions. Numerical results are presented to illustrate the ideas. / Master of Science

We study the Korteweg-de Vries-Burgers equation. With a deep investigation into the spectral and smoothing properties of the linearized system, it is shown by applying Banach Contraction Principle and Gronwall's Inequality to the integral equation based on the variation of parameters formula and explicit representation of the operator semigroup associated with the linearized equation that, under appropriate assumption appropriate assumption on initial states w(x, 0), the nonlinear system is well-posed and its solutions decay exponentially to the mean value of the initial state in H1(O, 1) as t -> +". / Ph. D.

The main objective of this research work is to apply the discontinuous Galerkin method to a classical partial differential equation to investigate the properties of the numerical solution and compare the numerical solution to the analytical solution by using discontinuous Galerkin method. This scheme is applied to 1-D non-linear conservation equation (Burgers equation) in which the governing differential equation is simplified model of the inviscid Navier-stokes equations. In this work three cases are studied. They are sinusoidal wave profile, initial shock discontinuity and initial linear distribution. A grid and time step refinement is performed. Riemann fluxes at each element interfaces are calculated. This scheme is applied to forward differentiation method (Euler's method) and to second order Runge-kutta method of this work.

Optimal control of PDEs has a crucial place in many parts of sciences and industry. Over the last decade, there have been a great deal in, especially, control problems of elliptic problems. Optimal control problems of Burgers equation that is as a simplifed model for turbulence and in shock waves were recently investigated both theoretically and numerically. In this thesis, we analyze the space-time simultaneous discretization of control problem for Burgers equation. In literature, there have been two approaches for discretization of optimization problems: optimize-then-discretize and discretize-then-optimize. In the first part, we follow optimize-then-discretize appoproach. It is shown that both distributed and boundary time dependent control problem can be transformed into an elliptic pde. Numerical results obtained with adaptive and non-adaptive elliptic solvers of COMSOL Multiphysics are presented for both the unconstrained and the control constrained cases. As for second part, we consider discretize-then-optimize approach. Discrete adjoint concept is covered. Optimality conditions, KKT-system, lead to a saadle point problem. We investigate the numerical treatment for the obtained saddle point system. Both direct solvers and iterative methods are considered. For iterative mehods, preconditioners are needed. The structures of preconditioners for both distributed and boundary control problems are covered. Additionally, an a priori error analysis for the distributed control problem is given. We present the numerical results at the end of each chapter.

CSIR (INDIA), IFCPAR / The physics of turbulence is the study of the chaotic and irregular behaviour in driven fluids. It is ubiquitous in cosmic, terrestrial and laboratory environments. To describe the properties of a simple incompressible fluid it is sufficient to know its velocity at all points in space and as a function of time. The equation of motion for the velocity of such a fluid is the incompressible Navier–Stokes equation. In more complicated cases, for example if the temperature of the fluid also fluctuates in space and time, the Navier–Stokes equation must be supplemented by additional equations. Incompressible fluid turbulence is the study of solutions of the Navier–Stokes equation at very high Reynolds numbers, Re, the dimensionless control parameter for this problem. The chaotic nature of these solutions leads us to characterise them by their statistical properties. For example, statistical properties of fluid turbulence are characterised often by structure functions of velocity. For intermediate range of length scales, that is the inertial range, these structure functions show multiscaling. Most studies concentrate on equal-time structure functions which describe the equal-time statistical properties of the turbulent fluid. Dynamic properties can be measured by more general time-dependent structure functions. A major challenge in the field of fluid turbulence is to understand the multiscaling properties of both the equal-time and time-dependent structure functions of velocity starting from the Navier–Stokes equation. In this thesis we use numerical and analytical techniques to study scaling and multiscaling of equal-time and time-dependent structure functions in turbulence not only in fluids but also in advection of passive-scalars and passive vectors, and in randomly forced Burgers equation.

abstract: Solution methods for certain linear and nonlinear evolution equations are presented in this dissertation. Emphasis is placed mainly on the analytical treatment of nonautonomous differential equations, which are challenging to solve despite the existent numerical and symbolic computational software programs available. Ideas from the transformation theory are adopted allowing one to solve the problems under consideration from a non-traditional perspective. First, the Cauchy initial value problem is considered for a class of nonautonomous and inhomogeneous linear diffusion-type equation on the entire real line. Explicit transformations are used to reduce the equations under study to their corresponding standard forms emphasizing on natural relations with certain Riccati(and/or Ermakov)-type systems. These relations give solvability results for the Cauchy problem of the parabolic equation considered. The superposition principle allows to solve formally this problem from an unconventional point of view. An eigenfunction expansion approach is also considered for this general evolution equation. Examples considered to corroborate the efficacy of the proposed solution methods include the Fokker-Planck equation, the Black-Scholes model and the one-factor Gaussian Hull-White model. The results obtained in the first part are used to solve the Cauchy initial value problem for certain inhomogeneous Burgers-type equation. The connection between linear (the Diffusion-type) and nonlinear (Burgers-type) parabolic equations is stress in order to establish a strong commutative relation. Traveling wave solutions of a nonautonomous Burgers equation are also investigated. Finally, it is constructed explicitly the minimum-uncertainty squeezed states for quantum harmonic oscillators. They are derived by the action of corresponding maximal kinematical invariance group on the standard ground state solution. It is shown that the product of the variances attains the required minimum value only at the instances that one variance is a minimum and the other is a maximum, when the squeezing of one of the variances occurs. Such explicit construction is possible due to the relation between the diffusion-type equation studied in the first part and the time-dependent Schrodinger equation. A modication of the radiation field operators for squeezed photons in a perfect cavity is also suggested with the help of a nonstandard solution of Heisenberg's equation of motion. / Dissertation/Thesis / Ph.D. Applied Mathematics for the Life and Social Sciences 2013

en.txt After an introduction into current sheet physics, with emphasis to solar physics we showed some formulations of finite element method. We implemented and evaluated new discontinuous finite element with Taylor basis and introduced deal.II library with an example of burgers equation. While the program is dimension independent, we compared our solution with a one-dimensional analytical solution. Finally, using previously derived LSFEM formulation, we solved simple current sheet problem using deal.II. Stránka 1

The aim of this project is to study the viscous Burgers' equation for the case where the viscosity is constant, but also when it contains a jump in viscosity. In the first case where the viscosity is constant, Burgers' is simply solved on a singular domain. For the case with jump in viscosity, Burgers' is solved on multiple domains with different viscosity. The different domains are then connected by applying inner boundary conditions at an interface in order to produce a singular solution. The inner boundary conditions are imposed using three different methods; simultaneous approximation term (SAT), projection and hybrid method, where the hybrid method is a combination of both the SAT and projection method. These methods are used in combination with a stable and high-order accurate summation by parts (SBP) finite difference approximation in MATLAB. The three methods are then compared to each other with respect to the least square error and the corresponding convergence rate to determine which method is the most preferable to use. The methods resulting in the highest convergence rates are the projection and the hybrid methods. These methods manage to live up to the expected convergence rates for all operators with different orders of accuracy and are therefore both good methods to use. However, the best method to use is the projection method since it is much simpler to implement than the two other methods but still achieves just as good convergence rates as the hybrid method.