Mathematical Aspects of Navier Stoke's Equations

The Navier-Stokes differential equations describe the motion of fluids which are incompressible. The three-dimensional Navier-Stokes equations misbehave very badly although they are relatively simple-looking. The solutions could wind up being extremely unstable even with nice, smooth, reasonably harmless initial conditions. A mathematical understanding of the outrageous behaviour of these equations would dramatically alter the
field of fluid mechanics. This paper describes why the three-dimensional Navier-Stokes equations are not solvable, i.e., the equations cannot be used to model turbulence, which is
a three-dimensional phenomenon.

This work is on the aerodynamics design of the body and frame of a motor tricycle using SolidWorks 2011 modeling system. Its Computational Fluid Dynamics (CFD) feature was used to run simulation tests at a target speed of 150 km/h to evaluate the aerodynamic performance of the tricycle. Simulation results presented shows that lift and drag forces are diminished considerably and that 170 km/h is the maximum speed to be traveled by the tricycle for a smooth and stable ride. Also, the design reduced drastically the effects of lift and drag forces, increased the tricycle’s stability, traction and performance as well as minimized the weight of the tricycle as a result of the use of high performance to mass ratio materials such as carbon fiber for the body and alloy steel for the frame and rims.

This PhD thesis focuses on numerical and analytical methods for simulating the dynamics of volcanic ash plumes.
The study starts from the fundamental balance laws for a multiphase gas– particle mixture, reviewing the existing models and developing a new set of Partial Differential Equations (PDEs), well suited for modeling multiphase dispersed turbulence. In particular, a new model generalizing the equilibrium–Eulerian model to two-way coupled compressible flows is developed.
The PDEs associated to the four-way Eulerian-Eulerian model is studied, in- vestigating the existence of weak solutions fulfilling the energy inequalities of the PDEs. In particular, the convergence of sequences of smooth solutions to such a set of weak solutions is showed.
Having explored the well-posedness of multiphase systems, the three-dimensional compressible equilibrium–Eulerian model is discretized and numerically solved by using the OpenFOAM® numerical infrastructure. The new solver is called ASHEE, and it is verified and validated against a number of well understood benchmarks and experiments. It demonstrates to be capable to capture the key phenomena involved in the dynamics of volcanic ash plumes. Those are: turbulence, mixing, heat transfer, compressibility, preferential concentration of particles, plume entrainment.
The numerical solver is tested by taking advantage of the newest High Perfor- mance Computing infrastructure currently available.
Thus, ASHEE is used to simulate two volcanic plumes in realistic volcanological conditions. The influence of model configuration on the numerical solution is analyzed. In particular, a parametric analysis is performed, based on: 1) the kinematic decoupling model; 2) the subgrid scale model for turbulence; 3) the discretization resolution.
In a one-dimensional and steady-state approximation, the multiphase flow model is used to derive a model for volcanic plumes in a calm, stratified atmosphere. The corresponding Ordinary Differential Equations (ODEs) are written in a compact, dimensionless formulation. The six non-dimensional parameters characterizing a multiphase plume are then written. The ODEs is studied both numerically and analytically. Different regimes are analyzed, extracting the first integral of motion and asymptotic solutions. An asymptotic analytical solution approximating the model in the general regime is derived and compared with numerical results. Such a solution is coupled with an electromagnetic model providing the infrared intensity emitted by a volcanic ash plume. Key vent parameters are then retrieved by means of inversion techniques applied to infrared images measured during a real volcanic eruption.

In parts 1 - 3 the method was described for solving cases of unidirectional flow. It was, obviously, an introduction to the concept of the method. In part 4 it is shown how this method can be used to solve an exemple of 2 dimensional flow. The square box problem with no convection is solved with this method in order to show its viability to be applied in cases of multidirectional fow

The motion of fluids which are incompressible could be described by the Navier-Stokes differential equations. However, the
three-dimensional Navier-Stokes equations for modelling turbulence misbehave very badly although they are relatively simple-looking. The solutions could wind up being extremely unstable even with nice, smooth, reasonably harmless initial conditions. A mathematical understanding of the outrageous behaviour of these equations would greatly affect the field of fluid mechanics. In this paper, which had been published in an international journal in 2010, a reasoned, practical approach towards resolving the issue is adopted and a practical, statistical kind of mathematical solution is proposed.

The paper's focus is the calculation of unsteady incompressible 2D flows past airfoils. In the framework of the primitive variable Navier–Stokes equations, the initial and boundary conditions must be assigned so as to be compatible, to assure the correct prediction of the flow evolution. This requirement, typical of all incompressible flows, viscous or inviscid, is often violated when modelling the flow past immersed bodies impulsively started from rest. Its fulfillment can however be restored by means of a procedure enforcing compatibility, consisting in a pre-processing of the initial velocity field, here described in detail. Numerical solutions for an impulsively started multiple airfoil have been obtained using a finite element incremental projection method. The spatial discretization chosen for the velocity and pressure are of different order to satisfy the inf–sup condition and obtain a smooth pressure field. Results are provided to illustrate the effect of employing or not the compatibility procedure, and are found in good agreement with those obtained with a non-primitive variable solver. In addition, we introduce a post-processing procedure to evaluate an alternative pressure field which is found to be more accurate than the one resulting from the projection method. This is achieved by considering an appropriate ‘unsplit’ version of the momentum equation, where the velocity solution of the projection method is substituted. Copyright © 2004 John Wiley & Sons, Ltd.

This article is focused on the choice of a suitable turbulence model for simulations of an industrial pump's intake, from the perspective of accuracy and, partially, also the CPU time. Twelve steady-state and transient simulations were made on a fine computational mesh, using turbulence models such as: the shear stress transport (SST), the scale-adaptive simulation (SAS), the Reynolds stress model, the explicit algebraic Reynolds-stress model, the detached eddy simulation and the large eddy simulation (LES). The curvature-correction (CC) option was assessed for the SST and SAS turbulence models. The results were compared with the LES and with published experimental results. Although all the models could predict the main floor vortex, there were still some substantial differences. It can be able to conclude that it is better to use either the SST-CC turbulence model, due to its low-computational resources and far better results than the SST model, or the SAS-CC turbulence model, since its predictions are quite similar to the LES results. In the final step, good agreement with experimental results was shown for a longer simulation with the SAS-CC turbulence model.

The motion of a fluid around a circular cylinder presents interesting phenomena including flow separation, wake and turbulence. The physics of these are enshrined in the continuity equation and the Navier-Stokes Equations (NSE). Therefore, their studies are important in mathematics and physics. They also have engineering applications. These studies can either be carried out experimentally, computationally, or theoretically. Theoretical studies of a cylinder flow using classical potential flow theory (CPT) have some gaps when compared to experiments. Attempting to bridge these gaps, this article introduces refined potential flow theory (RPT) and employs it on a stationary circular cylinder incompressible crossflow at Reynolds number 3,900. It leverages experimental observations, physical deductions and some agreements between CPT and experiments in the theoretical development. This results in the incompressible Eulerian Kwasu function which is a quasi-irrotational stream function that satisfies the governing equations and boundary conditions. It captures vorticity, boundary layer, shed wake vortices, three-dimensional effects, and static unsteadiness. The Lagrangian form of the function is exploited for the flow pathlines that are used to incorporate dynamic unsteadiness. A gravity analogy is used to predict the separation, transition, and reattachment points. The analogy introduces the perifocal frame of fluid motion. The forces are obtained in this frame with a change of variable. The drag prediction is within the error bound of measured data. The RPT pressure distribution, separation point and Strouhal number are also within acceptable ranges. Energy spectra analyses of the wake velocity display Kolmogorov’s Five-Thirds law of homogeneous isotropic turbulence.

Applications of understanding are extracted from the Clay Mathematics Institute's paper on the Navier-Stokes equations which seeks to find a force, pressure, and initial conditions in an infinitely differentiable non-divergent vector space of an imcompressible fluid.  A proof for equation (9) of that paper is produced.

The two-dimensional simulation described in part 4 is considerably more complex than the three first parts, that threated one-dimensional flow. However, it is limited to the presence of linear therms only, which is a very limited universe. The most important cases that need to be simulated with numerical methods ar the ones in which the non-linear therms are present, for in those cases no analytical solution is possible. In this part 5, it will be shown how to use this method when the convective therms ARE present.

Mathematical simulation of an external 3-D flow past a short solid cylinder show some astonishing behaviour. The author invite researcher to verify the existence of such fenomena and, if found true, to review all that have been stated for 2-D flows.

Stokes proposed a very ingenious method for describing 2-D flows, which uses one single equation that represents 3 different ones. However, this method is completelly useless to describe 3-D flows. Those have been a challenge to mechanical engineers. So far there is no method described in the literature that allows analytical solution to 3-D flows. Numerical methods use routines that are extremely expensive, in terms of computer resources. This article proposes a method that allows analytical solution to 3-D flows - if the convective terms are not present - and present some advatages for using in numerical simulators. Because it is too long, it has been divided into parts. This is part one, an introduction and the description of its use only in the simplest situation of all: no convective terms and recttangular coordinates.

As a continuation of the description of the proposed numeric method shown in parts 1 and 2, this part 3 shows how to use this method and its limits when both, convective and viscous transport of movement acts together.

PART 4 IS NOW DISPOSABLE!

In this paper, we propose a first order action functional for a large class of systems that generalize the relativistic perfect fluids in the Kaehler parametrization to noncommutative spacetimes. We calculate the equations of motion for the fluid potentials and the energy-momentum tensor in the first order in the noncommutative parameter. The density current does not receive any noncommutative corrections and it is conserved under the action of the commutative generators Pμ but the energy-momentum tensor is not. Therefore, we determine the set of constraints under which the energy-momentum tensor is divergenceless. Another set of constraints on the fluid potentials is obtained from the requirement of the invariance of the action under the generalization of the volume preserving transformations of the noncommutative spacetime. We show that the proposed action describes noncommutative fluid models by casting the energy-momentum tensor in the familiar fluid form and identifying the corresponding energy and momentum densities. In the commutative limit, they are identical to the corresponding quantities of the relativistic perfect fluids. The energy-momentum tensor contains a dissipative term that is due to the noncommutative spacetime and vanishes in the commutative limit. Finally, we particularize the theory to the case when the complex fluid potentials are characterized by a function K(z,z¯) that is a deformation of the complex plane and show that this model has important common features with the commutative fluid such as infinitely many conserved currents and a conserved axial current that in the commutative case is associated to the topologically conserved linking number.

In this paper we address the problem of the quantization of the perfect relativistic fluids formulated in terms of the Kaehler parametrization. This fluid model describes a large set of interesting systems such as the power law energy density fluids, Chaplygin gas, etc. In order to maintain the generality of the model, we apply the BRST method in the reduced phase space in which the fluid degrees of freedom are just the fluid potentials and the fluid current is classically resolved in terms of them. We determine the physical states in this setting, the time evolution and the path integral formulation.