One of the key open problems in the field of theoretical fluid mechanics concerns the possibility of the singularity formation in solutions of the 3D Navier-Stokes system in finite time. This phenomenon is associated with the behaviour of the enstrophy, which is an L2 norm of the vorticity and must become unbounded if such a singularity occurs. Although there is no blow-up in the 2D Navier-Stokes equation, we would like to investigate how much enstrophy can a planar incompressible flow in a bounded domain produce given certain initial enstrophy. We address this issue by formulating an optimization problem in which the time derivative of the enstrophy serves as the objective functional and solve it using tools of the optimization theory and calculus of variations. We propose an efficient computational approach which is based on the iterative steepest-ascent procedure. In addition, we introduce an easy-to-implement method of computing the gradient of the objective functional. Finally, we present computational results addressing the key question of this project and provide numerical evidence that the maximum enstrophy growth exhibits the scaling dE/dt ~ C*E*E for C>0 and very small E. All computations are performed using the Chebyshev spectral method. / Thesis / Master of Science (MSc) / For many decades, scientists have been investigating fundamental aspects of the Navier-Stokes equation, a central mathematical model arising in fluid mechanics. Although the equation is widely used by engineers to describe numerous flow phenomena, it is still an open question whether the Navier-Stokes system always admits physically meaningful solutions. To address this issue, we want to explore its mathematical aspects deeper by analyzing the behaviour of the enstrophy, which is a quantity associated with the vorticity of the flow and a convenient measure of the regularity of the solution. In this study, we consider a planar and incompressible flow bounded by solid walls. Using basic tools of mathematical analysis and optimization theory, we propose a computational method enabling us to find out how much enstrophy can such a flow produce instantaneously. We present numerical evidence that this instantaneous growth of enstrophy has a well-defined asymptotic behavior, which is consistent with physical assumptions.

The Lagrangian trajectory code TRACMASS is extended to the atmosphere to examine the tropi- cal Hadley Cells using fields from the ERA-Interim reanalysis dataset. The analysis is made using both pressure, temperature and specific humidity as vertical coordinates. By letting a trajectory represent a mass transport and tracing millions of trajectories in a domain between the latitudes 15°N and 15°S, the mass stream function based on trajectories is obtained (Lagrangian stream function). By separating the trajectories into classes depending on their starting point and des- tination (“North-to-North”, “North-to-South”, “South-to-North” and “South-to-South”), the mass stream function is decomposed into four paths. This can not be done if the stream function is cal- culated directly from the velocity fields (Eulerian stream function). Using this technique, the mass transports recirculating within the cells are compared to the mass transports between the cells, giving further insight to the structure of the Hadley Circulations. The magnitudes of the mass stream functions are presented by converting the volume flux unit Sverdrup into a mass flux unit. It is found that the recirculating transports of the northern and southern cells are 473 Sv and 508 Sv respectively. The inter-hemispheric mass transports are 126 Sv northward and 125 Sv southward. It is also found that far from all trajectories follow paths sim- ilar to the stream lines, since the stream lines are zonal and temporal means and the particle trajectories chaotic.

Wavelike spiral asymmetries are an intriguing aspect of Tropical Cyclone dynamics. Previous work hypothesized that some of them are Vortex Rossby Waves propagating on the radial gradient of mean–flow relative vorticity. In the Intermittently Forced Vortex Rossby Wave theory, intermittent convection near the eyewall wind maximum excites them so that they propagate wave energy outward and converge angular momentum inward. The waves’ energy is absorbed as the perturbation vorticity becomes filamented near the outer critical radii where their Doppler–shifted frequencies and radial group velocities approaches zero. This process may initiate outer wind maxima by weakening the mean–flow just inward from the critical radius. The waves are confined to a relatively narrow annular waveguide because of their slow tangential phase velocity and the narrow interval between the Rossby wave cut–off frequency, where the radial wavenumber is locally zero, and the zero frequency, where it is locally infinite.

The linear Barotropic Non-Divergent simulation of a vortex on a beta plane is consistent with Willoughby’s earlier shallow-water divergent results in that it produced an unbounded accelerating westward and poleward motion without an asymptotic limit. However, Montgomery’s work which yielded finite linear drift speeds for his completely cyclonic vortex was inconsistent with ours. The nonlinearly-forced streamfunction exhibited a beta-gyre like structure, but with opposite polarity phase to the linear gyres. Utilization of the linear model with time-dependent, but otherwise beta-like, forcing revealed increasing magnitude and phase reversal in the neighborhood of a low cyclonic frequency. Here, the mean bounded vortex has an outer waveguide that supports Vortex Rossby Wave propagation that is faster than the mean flow and confined to a very narrow band of frequencies between zero and the Vortex Rossby Wave cutoff. The low frequency waves constitute the beta-gyre mode described previously by Willoughby.