Moment methods from the kinetic theory of gases exist as an alternative to the Navier-Stokes model. Models in this family are described by first-order hyperbolic PDEs with local relaxation. They provide a natural treatment for non-equilibrium effects and expand the regime for which the model is physically applicable past the Navier-Stokes level (when the continuum assumption breaks down). Discontinuous-Galerkin (DG) methods are very well suited for distributed parallel solution of first-order PDEs. This is because the optimal locality of the method minimizes needed communication between computational processes. One highly efficient, coupled space-time DG method that achieves third-order accuracy in both space and time while using only linear elements is the discontinuous-Galerkin Hancock (DGH) scheme, which was specifically designed for the efficient solution of PDEs resulting from moment closures. Third-order accuracy is obtained through the use of a technique originally proposed by Hancock. The combination of moment methods with the DGH discretization leads to a very efficient numerical treatment for viscous compressible gas flows that is accurate both in and out of local thermodynamic equilibrium. This thesis describe the first-ever implementation of this scheme for the solution of moment equations on large-scale distributed-memory computers. This implementation uses solution-directed automatic mesh refinement to increase accuracy while reducing cost. A linear hyperbolic-relaxation equation is used to verify the order of accuracy of the scheme. Next a supersonic compressible Euler case is used to demonstrate the mesh refinement as well as the scheme’s ability to capture sharp discontinuities. Third, a moment-closure is then used to compute a viscous mixing layer. This serves to demonstrate the ability of the first-order PDEs and the DG scheme to efficiently compute viscous solutions. A moment-closure is used to compute the solution for Stokes flow past a circular cylinder. This case reinforces the hyperbolic PDEs’ ability to accurately predict viscous phenomena. As this case is very low speed, it also demonstrates the numerical technique’s ability to accurately solve problems that are ill-conditioned due to the extremely low Mach number. Finally, the parallel efficiency of the scheme is evaluated on Canada’s largest supercomputer. It may be surprising to some that viscous flow behaviour can be accurately predicted by first-order PDEs. However, the applicability of hyperbolic moment methods to both continuum and non-equilibrium gas flows is now well established. Such a first-order treatment brings many physical and computational advantages to gas flow prediction.

In this paper a continuous Galerkin multiscale method (CGMM) and a discontinuous Galerkin multiscale method (DGMM) are proposed, both based on the variational multiscale method for solving partial differential equations numerically. The solution is decoupled into a coarse and a fine scale contribution, where the fine-scale contribution is computed on patches with localized right hand side. Numerical experiments are presented where exponential decay of the error is observed when increasing the size of the patches for both CGMM and DGMM. DGMM gives much better accuracy when the same size of the patches are used.

No description available.

Near shore hydrodynamics has been an important research area dealing with coastal processes. The nearshore coastal region is the region between the shoreline and a fictive offshore limit which usually is defined as the limit where the depth becomes so large that it no longer influences the waves. This spatially limited but highly energetic zone is where water waves shoal, break and transmit energy to the shoreline and are governed by highly dispersive and non-linear effects. An accurate understanding of this phenomena is extremely useful, especially in emergency situations during hurricanes and storms. While the shallow water assumption is valid in regions where the characteristic wavelength exceeds a typical depth by orders of magnitude, Boussinesq-type equations have been used to model near-shore wave motion. Unfortunately these equations are complex system of coupled non-linear and dispersive differential equations that have made the developement of numerical approximations extremely challenging. In this dissertation, a local discontinuous Galerkin method for Boussinesq-Green Naghdi Equations is presented and validated against experimental results. Currently Green-Naghdi equations have many variants. We develop a numerical method in one horizontal dimension for the Green-Naghdi equations based on rotational characteristics in the velocity field. Stability criterion is also established for the linearized Green-Naghdi equations and a careful proof of linear stability of the numerical method is carried out. Verification is done against a linearized standing wave problem in flat bathymetry and h,p (denoted by K in this thesis) error rates are plotted. The numerical method is validated with experimental data from dispersive and non-linear test cases. / text

<p>Domain decomposition method provides a solution for the very large electromagnetic</p><p>system which are impossible for single domain methods. Discontinuous Galerkin</p><p>(DG) method can be viewed as an extreme version of the domain decomposition,</p><p>i.e., each element is regarded as one subdomain. The whole system is solved element</p><p>by element, thus the inversion of the large global system matrix is no longer necessary,</p><p>and much larger system can be solved with the DG method compared to the</p><p>continuous Galerkin (CG) method.</p><p>In this work, the DG method is implemented on a subdomain level, that is, each subdomain contains multiple elements. The numerical flux only applies on the</p><p>interfaces between adjacent subdomains. The subodmain level DG method divides</p><p>the original large global system into a few smaller ones, which are easier to solve,</p><p>and it also provides the possibility of parallelization. Compared to the conventional</p><p>element level DG method, the subdomain level DG has the advantage of less total</p><p>DoFs and fexibility in interface choice. In addition, the implicit time stepping is </p><p>relatively much easier for the subdomain level DG, and the total CPU time can be</p><p>much less for the electrically small or multiscale problems.</p><p>The hybrid of elements are employed to reduce the total DoF of the system.</p><p>Low-order tetrahedrons are used to catch the geometry ne parts and high-order</p><p>hexahedrons are used to discretize the homogeneous and/or geometry coarse parts.</p><p>In addition, the non-conformal mesh not only allow dierent kinds of elements but</p><p>also sharp change of the element size, therefore the DoF can be further decreased.</p><p>The DGTD method in this research is based on the EB scheme to replace the</p><p>previous EH scheme. Dierent from the requirement of mixed order basis functions</p><p>for the led variables E and H in the EH scheme, the EB scheme can suppress the</p><p>spurious modes with same order of basis functions for E and B. One order lower in</p><p>the basis functions in B brings great benets because the DoFs can be signicantly</p><p>reduced, especially for the tetrahedrons parts.</p><p>With the basis functions for both E and B, the EB scheme upwind </p><p>ux and</p><p>EB scheme Maxwellian PML, the eigen-analysis and numerical results shows the</p><p>eectiveness of the proposed DGTD method, and multiscale problems are solved</p><p>eciently combined with the implicit-explicit hybrid time stepping scheme and multiple</p><p>kinds of elements.</p><p>The EB scheme DGTD method is further developed to allow arbitrary anisotropic</p><p>media via new anisotropic EB scheme upwind </p><p>ux and anisotropic EB scheme</p><p>Maxwellian PML. The anisotropic M-PML is long time stable and absorb the outgoing</p><p>wave eectively. A new TF/SF boundary condition is brought forward to</p><p>simulate the half space case. The negative refraction in YVO4 bicrystal is simulated</p><p>with the anisotropic DGTD and half space TF/SF condition for the rst time with</p><p>numerical methods.</p> / Dissertation

Coastal regions are the most densely populated regions of the world. The populations of these regions continue to grow which has created a high demand for water that stresses existing water resources. Coastal aquifers provide a source of water for coastal populations and are generally part of a larger system where freshwater aquifers are hydraulically connected with a saline surface-water body. They are characterized by salinity variations in space and time, sharp freshwater/saltwater interfaces which can lead to dramatic density differences, and complex groundwater chemistry. Mismanagement of coastal aquifers can lead to saltwater intrusion, the displacement of fresh water by saline water in the freshwater regions of the aquifers, making them unusable as a freshwater source. Saltwater intrusion is of significant interest to water resource managers and efficient simulators are needed to assist them. Numerical simulation of saltwater intrusion requires solving a system of flow and transport equations coupled through a density equation of state. The scale of the problem domain, irregular geometry and heterogeneity can require significant computational resources. Also, modeling sharp transition zones and accurate flow velocities pose numerical challenges. Discontinuous Galerkin (DG) finite element methods (FEM) have been shown to be well suited for modeling flow and transport in porous media but a fully coupled DG formulation has not been applied to the variable density flow and transport model. DG methods have many desirable characteristics in the areas of numerical stability, mesh and polynomial approximation adaptivity and the use of non-conforming meshes. These properties are especially desirable when working with complex geometries over large scales and when coupling multi-physics models (e.g. surface water and groundwater flow models). In this dissertation, we investigate a new combined local discontinuous Galerkin (LDG) and non-symmetric, interior penalty Galerkin (NIPG) formulation for the non-linear coupled flow and solute transport equations that model saltwater intrusion. Our main goal is the formulation and numerical implementation of a robust, efficient, tightly-coupled combined LDG/NIPG formulation within the Department of Defense (DoD) Proteus Computational Mechanics Toolkit modeling framework. We conduct an extensive and systematic code and model verification (using established benchmark problems and proven convergence rates) and model validation (using experimental data) to verify accomplishment of this goal. Lastly, we analyze the accuracy and conservation properties of the numerical model. More specifically, we derive an a priori error estimate for the coupled system and conduct a flow/transport model compatibility analysis to prove conservation properties. / text

A look into an implementation of turbulence model into the ADGFEM code for viscous flow. Discretization, theory background and development of the method will be carried during this thesis. Also some numerical examples of the application of the code will be provided. 1

abstract: This dissertation describes a process for interface capturing via an arbitrary-order, nearly quadrature free, discontinuous Galerkin (DG) scheme for the conservative level set method (Olsson et al., 2005, 2008). The DG numerical method is utilized to solve both advection and reinitialization, and executed on a refined level set grid (Herrmann, 2008) for effective use of processing power. Computation is executed in parallel utilizing both CPU and GPU architectures to make the method feasible at high order. Finally, a sparse data structure is implemented to take full advantage of parallelism on the GPU, where performance relies on well-managed memory operations. With solution variables projected into a kth order polynomial basis, a k+1 order convergence rate is found for both advection and reinitialization tests using the method of manufactured solutions. Other standard test cases, such as Zalesak's disk and deformation of columns and spheres in periodic vortices are also performed, showing several orders of magnitude improvement over traditional WENO level set methods. These tests also show the impact of reinitialization, which often increases shape and volume errors as a result of level set scalar trapping by normal vectors calculated from the local level set field. Accelerating advection via GPU hardware is found to provide a 30x speedup factor comparing a 2.0GHz Intel Xeon E5-2620 CPU in serial vs. a Nvidia Tesla K20 GPU, with speedup factors increasing with polynomial degree until shared memory is filled. A similar algorithm is implemented for reinitialization, which relies on heavier use of shared and global memory and as a result fills them more quickly and produces smaller speedups of 18x. / Dissertation/Thesis / Doctoral Dissertation Aerospace Engineering 2015

This thesis studies the numerical solution of non-linear convection-diffusion problems using the space- time discontinuous Galerkin method, which perfectly suits the space as well as time local adaptation. We aim to develop a posteriori error estimates reflecting the spatial, temporal, and algebraic errors. These estimates are based on the measurement of the residuals in dual norms. We derive these estimates and numerically verify their properties. Finally, we derive an adaptive algorithm and apply it to the numerical simulation of non-stationary viscous compressible flows. Powered by TCPDF (www.tcpdf.org)

As traditional communication technology we use in our day-to-day life reaches its limitations, the international community searches for new methods to communicate information. One such novel approach is the so-called molecular communication system. During the last few decades, molecular communication systems become more and more popular. The main difference between traditional communication and molecular communication systems is that in the latter, information transfer occurs through chemical means, most often between microorganisms. This process already happens all around us naturally, for example, in the human body. Even though the molecular communication topic is attractive to researchers, and a lot of theoretical results are available - one cannot claim the same about the practical use of molecular communication. As for experimental results, a few studies have been done on the macroscale, but investigations at the micro- and nanoscale ranges are still lacking because they are a challenging task. In this work, a self-contained introduction of the underlying theory of molecular communication is provided, which includes knowledge from different areas such as biology, chemistry, communication theory, and applied mathematics. Two numerical methods are implemented for three well-studied partial differential equations of the MC field where advection, diffusion, and the reaction are taken into account. Numerical results for test cases in one and three dimensions are presented and discussed in detail. Conclusions and essential analytical and numerical future directions are then drawn.