A two-dimensional discontinuous Galerkin finite element method algorithm in the time domain was developed for calculation of the radar cross-section of an arbitrary object. The algorithm was formed using local nodal basis functions in each element and coupling them via numerical upwind flux. Both transverse electric and transverse magnetic polarization, as well as three different dispersive material models, were handled. The computational domain was effectively truncated with low reflections using the uniaxial perfectly matched layer method. Two different time stepping methods were used, low-storage explicit Runge-Kutta and Leap-Frog, to allow for flexibility in the time step and application of a stabilization method. The algorithm was verified with geometries, which have analytical expressions, and an existing validated code. The algorithm was also compared to an existing algorithm, which utilized the continuous finite element method with implicit time stepping, and showed outstanding performance regarding computation time and memory allocation. Since the developed algorithm had explicit time stepping could no general conclusions favoring any of the methods beyond these specific algorithms be made. The results still encouraged continued development of the DGFEM algorithm, where the expansion into three dimensions and optimizations could be explored further.

Indiana University-Purdue University Indianapolis (IUPUI) / The interplay between biomechanics and blood perfusion in the optic nerve head (ONH) has a critical role in ocular pathologies, especially glaucomatous optic neuropathy. Elucidating the complex interactions of ONH perfusion and tissue structure in health and disease using current imaging methodologies is difficult, and mathematical modeling provides an approach to address these limitations. The biophysical phenomena governing the ONH physiology occur at different scales in time and space and porous media theory provides an ideal framework to model them. We critically review fundamentals of porous media theory, paying particular attention to the assumptions leading to a continuum biphasic model for the phenomenological description of fluid flow through biological tissues exhibiting viscoelastic behavior. The resulting system of equations is solved via a numerical method based on a novel hybridizable discontinuous Galerkin finite element discretization that allows accurate approximations of stresses and discharge velocities, in addition to solid displacement and fluid pressure. The model is used to theoretically investigate the influence of tissue viscoelasticity on the blood perfusion of the lamina cribrosa in the ONH. Our results suggest that changes in viscoelastic properties of the lamina may compromise tissue perfusion in response to sudden variations of intraocular pressure, possibly leading to optic disc hemorrhages.

No description available.

No description available.

No description available.

Truncation errors and computational cost are obstacles that still hinder large-scale applications of the Computational Fluid Dynamics method. The discontinuous Galerkin method is one of the high-order schemes utilized extensively in recent years, which is locally conservative, stable, and high-order accurate. Besides that, it can handle complex geometries and irregular meshes with hanging nodes. In this document, the nondimensional compressible Euler equations and Reynolds- Averaged Navier-Stokes equations are discretized by discontinuous Galerkin methods with a two-equations turbulence model on both structured and unstructured meshes. The traditional equation of state for an ideal gas model is substituted by a multispecies thermodynamics model in order to complete the governing equations. An approximate Riemann solver is used for computing the convective flux, and the diffusive flux is approximated with some internal penalty based schemes. The temporal discretization of the partial differential equations is either performed explicitly with the aid of Rung-Kutta methods or with semi-implicit methods. Inspired by the artificial viscosity diffusion based limiter for shock-capturing method, which has been extensively studied, a novel and robust technique based on the introduction of mass diffusion to the species governing equations to guarantee that the species mass fractions remain positive has been thoroughly investigated. This contact-surface-capturing method is conservative and a high order of accuracy can be maintained for the discontinuous Galerkin method. For each time step of the algorithm, any trouble cell is first caught by the contact-surface discontinuity detector. Then some amount of mass diffusions are added to the governing equations to change the gas mixtures and arrive at an equilibrium point satisfying some conditions. The species properties are reasonable without any oscillations. Computations are performed for many steady and unsteady flow problems. For general non-mixing fluid flows, the classical air-helium shock bubble interaction problem is the central test case for the high-order discontinuous Galerkin method with a mass diffusion based limiter chosen. The computed results are compared with experimental, exact, and empirical data to validate the fluid dynamic solver.

No description available.

No description available.

No description available.

No description available.