Adaptive optics (AO) systems have become an essential part of ground-based telescopes and enable diffraction-limited imaging at near-IR and mid-IR wavelengths. For several key science applications the required wavefront quality is higher than what current systems can deliver. For instance obtaining high quality diffraction-limited images at visible wavelengths requires residual wavefront errors to be well below 100 nm RMS. High contrast imaging of exoplanets and disks around nearby stars requires high accuracy control of low-order modes that dominate atmospheric turbulence and scatter light at small angles where exoplanets are likely to be found. Imaging planets using a high contrast corona graphic camera, as is the case for the Spectro-Polarimetric High-contrast Exoplanet Research (SPHERE) on the Very Large Telescope (VLT), and the Gemini Planet Imager (GPI), requires even greater wavefront control accuracy. My dissertation develops a highly sensitive non-linear Curvature Wavefront Sensor (nlCWFS) that can deliver diffraction-limited (λ/D) images, in the visible, by approaching the theoretical sensitivity limit imposed by fundamental physics. The nlCWFS is derived from the successful curvature wavefront sensing concept but uses a non-linear reconstructor in order to maintain sensitivity to low spatial frequencies. The nlCWFS sensitivity makes it optimal for extreme AO and visible AO systems because it utilizes the full spatial coherence of the pupil plane as opposed to conventional sensors such as the Shack-Hartmann Wavefront Sensor (SHWFS) which operate at the atmospheric seeing limit (λ/r₀). The difference is equivalent to a gain of (D/r₀)² in sensitivity, for the lowest order mode, which translates to the nlCWFS requiring that many fewer photons. When background limited the nlCWFS sensitivity scales as D⁴, a combination of D² gain due to the diffraction limit and D² gain due to telescope's collecting power. Whereas conventional wavefront sensors only benefit from the D² gain due to the telescope's collecting power. For a 6.5 m telescope, at 0.5 μm, and seeing of 0.5", the nlCWFS can deliver for low order modes the same wavefront measurement accuracy as the SHWFS with 1000 times fewer photons. This is especially significant for upcoming extremely large telescopes such as the Giant Magellan Telescope (GMT) which has a 25.4 m aperture, the Thirty Meter Telescope (TMT) and the European Extremely Large Telescope (E-ELT) which has a 39 m aperture.

In the present work we study the numerical solution of shallow water equations. We introduce a vectorial notation of equations laws of conservation from which we derive the shallow water equations (SWE). There is the simplify its derivation, notation and the most important features. The original contribution is to derive equations for shallow water without the using of Leibniz's formula. There we report the finite volume method with the numerical flow of Vijayasundaram type for SWE. We present a description of the linear reconstruction, quadratic reconstruction and ENO reconstruction and their using for increasing of order accuracy. We demonstrate using of linear reconstruction in finite volume method of second order accuracy. This method is programmed in Octave language and used for solving of two problems. We apply the method of the ADER type for the shallow water equations. This method was originally designed for the Euler's equation.

A finite-volume procedure, comprising a gradient-reconstruction technique and a multidimensional limiter, has been proposed for upwind algorithms on unstructured grids. The high-resolution strategy, with its inherent dependence on a wide computational stencil, does not suffer from a catastrophic loss of accuracy on a grid with poor connectivity as reported recently is the case with many unstructured-grid limiting procedures. The continuously-differentiable limiter is shown to be effective for strong discontinuities, even on a grid which is composed of highly-distorted triangles, without adversely affecting convergence to steady state. Numerical experiments involving transient computations of two-dimensional scalar convection to steady-state solutions of Euler and Navier-Stokes equations demonstrate the capabilities of the new procedure.

A high-order central essentially non-oscillatory (CENO) finite-volume scheme in combination with a block-based adaptive mesh refinement (AMR) algorithm is proposed for solution of hyperbolic and elliptic systems of conservation laws on body- fitted multi-block mesh. The spatial discretization of the hyperbolic (inviscid) terms is based on a hybrid solution reconstruction procedure that combines an unlimited high-order k-exact least-squares reconstruction technique following from a fixed central stencil with a monotonicity preserving limited piecewise linear reconstruction algorithm. The limited reconstruction is applied to computational cells with under-resolved solution content and the unlimited k-exact reconstruction procedure is used for cells in which the solution is fully resolved. Switching in the hybrid procedure is determined by a solution smoothness indicator. The hybrid approach avoids the complexity associated with other ENO schemes that require reconstruction on multiple stencils and therefore, would seem very well suited for extension to unstructured meshes. The high-order elliptic (viscous) fluxes are computed based on a k-order accurate average gradient derived from a (k+1)-order accurate reconstruction. A novel h-refinement criterion based on the solution smoothness indicator is used to direct the steady and unsteady refinement of the AMR mesh. The predictive capabilities of the proposed high-order AMR scheme are demonstrated for the Euler and Navier-Stokes equations governing two-dimensional compressible gaseous flows as well as for advection-diffusion problems characterized by the full range of Peclet numbers, Pe. The ability of the scheme to accurately represent solutions with smooth extrema and yet robustly handle under-resolved and/or non-smooth solution content (i.e., shocks and other discontinuities) is shown for a range of problems. Moreover, the ability to perform mesh refinement in regions of smooth but under-resolved and/or non-smooth solution content to achieve the desired resolution is also demonstrated.

A high-order central essentially non-oscillatory (CENO) finite-volume scheme in combination with a block-based adaptive mesh refinement (AMR) algorithm is proposed for solution of hyperbolic and elliptic systems of conservation laws on body- fitted multi-block mesh. The spatial discretization of the hyperbolic (inviscid) terms is based on a hybrid solution reconstruction procedure that combines an unlimited high-order k-exact least-squares reconstruction technique following from a fixed central stencil with a monotonicity preserving limited piecewise linear reconstruction algorithm. The limited reconstruction is applied to computational cells with under-resolved solution content and the unlimited k-exact reconstruction procedure is used for cells in which the solution is fully resolved. Switching in the hybrid procedure is determined by a solution smoothness indicator. The hybrid approach avoids the complexity associated with other ENO schemes that require reconstruction on multiple stencils and therefore, would seem very well suited for extension to unstructured meshes. The high-order elliptic (viscous) fluxes are computed based on a k-order accurate average gradient derived from a (k+1)-order accurate reconstruction. A novel h-refinement criterion based on the solution smoothness indicator is used to direct the steady and unsteady refinement of the AMR mesh. The predictive capabilities of the proposed high-order AMR scheme are demonstrated for the Euler and Navier-Stokes equations governing two-dimensional compressible gaseous flows as well as for advection-diffusion problems characterized by the full range of Peclet numbers, Pe. The ability of the scheme to accurately represent solutions with smooth extrema and yet robustly handle under-resolved and/or non-smooth solution content (i.e., shocks and other discontinuities) is shown for a range of problems. Moreover, the ability to perform mesh refinement in regions of smooth but under-resolved and/or non-smooth solution content to achieve the desired resolution is also demonstrated.