Applied Mathematics Commons^™

Stabilized Weak-Gradient Discontinuous Finite Elements With Optimal Error Estimates For Second-Order Elliptic Pdes, Aymen Laadhari

Mathematical Modelling and Numerical Simulation with Applications

This work introduces an accurate finite element approach employing a new stabilized discrete weak gradient, designed for second-order elliptic problems on arbitrary conforming meshes. We formulate the approach within a discontinuous Galerkin framework and derive a consistent and coercive bilinear form. Appropriate error analysis on a model problem confirms optimal convergence. Building on the core analysis, we extend the method to more challenging settings, including time-dependent heterogeneous scenarios and a biophysically realistic optimal-control model of photobleaching in the budding yeast cell. We further illustrate the versatility of the weak-gradient construction by applying it to an unsteady level-set equation relevant to …

Investigations Of Conjugate Heat Transfer And Fluid Flow In Partitioned Porous Cavity Using Darcy-Forchheimer Model: Finite Element-Based Computations, Nasir Yasin, Shafee Ahmad, Muhammad Umair, Zahir Shah, Narcisa Vrinceanu, Ghadah Alhawael

Mathematics & Statistics Faculty Publications

The conjugate heat transfer and fluid flow has vast applications in thermal engineering, particularly for cooling in thermal devices, and automobile engines. This study investigates conjugate heat transfer in 2D enclosures, featuring thin solid fins attached to a porous bottom wall. The porous medium is considered isotropic and homogeneous by the Darcy-Forchheimer model, with fluid phases in local thermal equilibrium. The boundary conditions at the porous fluid interface ensure continuity of the velocities, stresses, temperature, and heat flux. The phenomenon is mathematically modelled by obtaining a set of partial differential equations. The finite element method (FEM) is used to perform …

Mathematical Modeling: Finite Element Analysis And Computations Arising In Fluid Dynamics And Biological Applications, Jorge Reyes

UNLV Theses, Dissertations, Professional Papers, and Capstones

It is often the case when attempting to capture real word phenomena that the resulting mathematical model is too difficult and even not feasible to be solved analytically. As a result, a computational approach is required and there exists many different methods to numerically solve models described by systems of partial differential equations. The Finite Element Method is one of them and it was pursued herein.This dissertation focuses on the finite element analysis and corresponding numerical computations of several different models. The first part consists of a study on two different fluid flow models: the main governing model of fluid …

Implementation Of A Least Squares Method To A Navier-Stokes Solver, Jada P. Lytch, Taylor Boatwright, Ja'nya Breeden

Rose-Hulman Undergraduate Mathematics Journal

The Navier-Stokes equations are used to model fluid flow. Examples include fluid structure interactions in the heart, climate and weather modeling, and flow simulations in computer gaming and entertainment. The equations date back to the 1800s, but research and development of numerical approximation algorithms continues to be an active area. To numerically solve the Navier-Stokes equations we implement a least squares finite element algorithm based on work by Roland Glowinski and colleagues. We use the deal.II academic library , the C++ language, and the Linux operating system to implement the solver. We investigate convergence rates and apply the least squares …

Simulating Dislocation Densities With Finite Element Analysis, Ja'nya Breeden, Dow Drake, Saurabh Puri

REU Final Reports

A one-dimensional set of nonlinear time-dependent partial differential equations developed by Acharya (2010) is studied to observe how differing levels of applied strain affect dislocation walls. The framework of this model consists of a convective and diffusive term which is used to develop a linear system of equations to test two methods of the finite element method. The linear system of partial differential equations is used to determine whether the standard or Discontinuous Galerkin method will be used. The Discontinuous Galerkin method is implemented to discretize the continuum model and the results of simulations involving zero and non-zero applied strain …

Smooth Global Approximation For Continuous Data Assimilation, Kenneth R. Brown

Theses and Dissertations

This thesis develops the finite element method, constructs local approximation operators, and bounds their error. Global approximation operators are then constructed with a partition of unity. Finally, an application of these operators to data assimilation of the two-dimensional Navier-Stokes equations is presented, showing convergence of an algorithm in all Sobolev topologies.

Decoupled Finite Element Methods For General Steady Two-Dimensional Boussinesq Equations, Lioba Boveleth

Masters Theses

"This work presents two kinds of decoupled finite element methods for the steady natural convection problem in two dimensions. Firstly, the standard Galerkin finite element method is derived in detail stating algorithms needed for the realization in MATLAB. A numerical example verifies the error convergence. Secondly, using iteration, the Boussinesq equations are decoupled into the Navier-Stokes equations and a parabolic problem. The resulting problems are solved either in parallel or sequentially. Finally, the same numerical example as before is used to confirm the convergence and analyze the methods in terms of iteration performance. In addition to a higher flexibility and …

Trefftz Finite Elements On Curvilinear Polygons, Akash Anand, Jeffrey S. Ovall, Samuel E. Reynolds, Steffen Weisser

Mathematics and Statistics Faculty Publications and Presentations

We present a Trefftz-type finite element method on meshes consisting of curvilinear polygons. Local basis functions are computed using integral equation techniques that allow for the efficient and accurate evaluation of quantities needed in the formation of local stiffness matrices. To define our local finite element spaces in the presence of curved edges, we must also properly define what it means for a function defined on a curved edge to be "polynomial" of a given degree on that edge. We consider two natural choices, before settling on the one that yields the inclusion of complete polynomial spaces in our local …

Simulating Backward Wave Propagation In Metamaterial With Radial Basis Functions, Jichun Li, Bing Nan

Mathematical Sciences Faculty Research

In this paper we revisit the radial basis function (RBF) meshless method and implement it to solve the time-dependent Maxwell's equations in metamaterials. Numerical simulations of the backward wave propagation phenomena are successfully achieved using this RBF method with properly chosen free parameters in the RBFs.

Numerical Analysis And Fluid Flow Modeling Of Incompressible Navier-Stokes Equations, Tahj Hill

UNLV Theses, Dissertations, Professional Papers, and Capstones

The Navier-Stokes equations (NSE) are an essential set of partial differential equations for governing the motion of fluids. In this paper, we will study the NSE for an incompressible flow, one which density ρ = ρ0 is constant.

First, we will present the derivation of the NSE and discuss solutions and boundary conditions for the equations. We will then discuss the Reynolds number, a dimensionless number that is important in the observations of fluid flow patterns. We will study the NSE at various Reynolds numbers, and use the Reynolds number to write the NSE in a nondimensional form.

We will …

Numerical Solution Of Mhd Bioconvection In A Porous Square Cavity Due To Oxytactic Microorganisms, Chandra S. Balla, Kishan Naikoti

Applications and Applied Mathematics: An International Journal (AAM)

The present paper investigates the magnetohydrodynamic (MHD) bioconvection flow in a porous square cavity filled with oxytactic microorganism. The bioconvection flow and heat transfer in porous media is formulated using Darcy model of Boussinesq approximation. Finite element method based on Galerkin weighted residual scheme is used to solve the governing partial differential equations. The computational numerical results are illustrated in the form of streamlines, isotherms, isoconcentrations of oxygen and microorganisms, average Nusselt number and average Sherwood number. In the present study the effects of key parameters such as bioconvection Rayleigh number (Rb), Rayleigh number (Ra), Peclet number (Pe) magnetic field …

Spectral Discretization Errors In Filtered Subspace Iteration, Jay Gopalakrishnan, Luka Grubišić, Jeffrey S. Ovall

Mathematics and Statistics Faculty Publications and Presentations

We consider filtered subspace iteration for approximating a cluster of eigenvalues (and its associated eigenspace) of a (possibly unbounded) selfadjoint operator in a Hilbert space. The algorithm is motivated by a quadrature approximation of an operator-valued contour integral of the resolvent. Resolvents on infinite dimensional spaces are discretized in computable finite-dimensional spaces before the algorithm is applied. This study focuses on how such discretizations result in errors in the eigenspace approximations computed by the algorithm. The computed eigenspace is then used to obtain approximations of the eigenvalue cluster. Bounds for the Hausdorff distance between the computed and exact eigenvalue clusters …

Full Field Computing For Elastic Pulse Dispersion In Inhomogeneous Bars, A. Berezovski, R. Kolman, M. Berezovski, D. Gabriel, V. Adamek

Publications

In the paper, the finite element method and the finite volume method are used in parallel for the simulation of a pulse propagation in periodically layered composites beyond the validity of homogenization methods. The direct numerical integration of a pulse propagation demonstrates dispersion effects and dynamic stress redistribution in physical space on example of a one-dimensional layered bar. Results of numerical simulations are compared with analytical solution constructed specifically for the considered problem. Analytical solution as well as numerical computations show the strong influence of the composition of constituents on the dispersion of a pulse in a heterogeneous bar and …

Finite Element Solution Of The Two-Dimensional Incompressible Navier-Stokes Equations Using Matlab, Endalew G. Tsega, V. K. Katiyar

Applications and Applied Mathematics: An International Journal (AAM)

The Navier–Stokes equations are fundamental in fluid mechanics. The finite element method has become a popular method for the solution of the Navier-Stokes equations. In this paper, the Galerkin finite element method was used to solve the Navier-Stokes equations for two-dimensional steady flow of Newtonian and incompressible fluid with no body forces using MATLAB. The method was applied to the lid-driven cavity problem. The eight-noded rectangular element was used for the formulation of element equations. The velocity components were located at all of 8 nodes and the pressure variable is located at 4 corner of the element. From location of …

The Auxiliary Space Preconditioner For The De Rham Complex, Jay Gopalakrishnan, Martin Neumüller, Panayot S. Vassilevski

Portland Institute for Computational Science Publications

We generalize the construction and analysis of auxiliary space preconditioners to the n-dimensional finite element subcomplex of the de Rham complex. These preconditioners are based on a generalization of a decomposition of Sobolev space functions into a regular part and a potential. A discrete version is easily established using the tools of finite element exterior calculus. We then discuss the four-dimensional de Rham complex in detail. By identifying forms in four dimensions (4D) with simple proxies, form operations are written out in terms of familiar algebraic operations on matrices, vectors, and scalars. This provides the basis for our implementation of …

Derivation Of The Hellinger-Reissner Variational Form Of The Linear Elasticity Equations, And A Finite Element Discretization, Bram Fouts

REU Final Reports

In this paper we are going to derive the linear elasticity equations in the Strong Form to the Hellinger Reissner Form. We find a suitable solution to solve our stress tensor. Then we will use finite element discretization from. We will run tests on a unit cube and multiple other shapes, which are described at the end. We view the different magnitudes of the displacement vector of each shape.

A New Formulation Of Time Boundary Integral Equation For Acoustic Wave Scattering In The Presence Of A Uniform Mean Flow, Fang Q. Hu, Michelle E. Pizzo, Douglas M. Nark

Mathematics & Statistics Faculty Publications

It has been well-known that under the assumption of a constant uniform mean flow, the acoustic wave propagation equation can be formulated as a boundary integral equation, in both the time domain and the frequency domain. Compared with solving partial differential equations, numerical methods based on the boundary integral equation have the advantage of a reduced spatial dimension and, hence, requiring only a surface mesh. However, the constant uniform mean flow assumption, while convenient for formulating the integral equation, does not satisfy the solid wall boundary condition wherever the body surface is not aligned with the uniform mean flow. In …

Effect Of Solid Volume Fraction On Forced Convective Flow Of Nanofluid Through Direct Absorption Solar Collector, Salma Parvin, Rehena Nasrin, M. A. Alim

Applications and Applied Mathematics: An International Journal (AAM)

The present work numerically investigates the heat transfer performance and entropy generation of forced convection through a direct absorption solar collector. The working fluid is Cu-water nanofluid. The simulations focus specifically on the effect of solid volume fraction of nanoparticle on the mean Nusselt number, total entropy generation, Bejan number and collector efficiency. Also Isotherms, heat function and entropy generation are presented for various solid volume fraction. The governing partial differential equations are solved using penalty finite element method with Galerkins weighted residual technique. The results show that the mean Nusselt number and mean entropy generation increases as the volume …

Prandtl Number Effect On Assisted Convective Heat Transfer Through A Solar Collector, Rehen Nasrin, Salma Parvin, M. A. Alim

Applications and Applied Mathematics: An International Journal (AAM)

Numerical study of the influence of Prandtl number on forced convective heat transfer through a riser pipe of a flat plate solar collector is done. The working fluid is Al2O3/water nanofluid. By Finite Element Method the governing partial differential equations are solved. The effect of the Prandtl number on the temperature and velocity field has been depicted. Comprehensive average Nusselt number, average bulk temperature, mean velocity, mid-height temperature inside the pipe, mean output temperature and collector efficiency are presented for the governing parameter mentioned above. Nu increases by 16% with the variation of Pr from 4.6 to 6.6 using nanofluid. …

On The Coupling Of Dpg And Bem, Thomas Führer, Norbert Heuer, Michael Karkulik

Mathematics and Statistics Faculty Publications and Presentations

We develop and analyze strategies to couple the discontinuous Petrov-Galerkin method with optimal test functions to (i) least-squares boundary elements and (ii) various variants of standard Galerkin boundary elements. An essential feature of our method is that, despite the use of boundary integral equations, optimal test functions have to be computed only locally. We apply our findings to a standard transmission problem in full space and present numerical experiments to validate our theory.

Multiphysics Modeling To Enhance Understanding Of Microwave Heating Of Food In Domestic Ovens, Jiajia Chen

Department of Agricultural and Biological Systems Engineering: Dissertations, Theses, and Student Research

Nonuniform heating is the biggest issue in the microwave heating of prepared meals. Multiphysics based models are promising tools to improve microwave heating uniformity by properly designing the food product. However, limited availability of accurate temperature-dependent material properties, inadequate model prediction accuracy, and high computational power and complexity in model development are three gaps that greatly limited the application of these models in the food industry.

To fill in the gaps, firstly, we developed a multitemperature calibration protocol to measure temperature-dependent dielectric properties (dielectric constant and loss factor). The temperature-dependent dielectric and thermal (thermal conductivity and specific heat capacity) properties …

Numerical Simulations Of Traffic Flow Models, Puneet Lakhanpal

UNLV Theses, Dissertations, Professional Papers, and Capstones

Traffic flow has been considered to be a continuum flow of a compressible liquid having a certain density profile and an associated velocity, depending upon density, position and time. Several one-equation and two-equation macroscopic continuum flow models have been developed which utilize the fluid dynamics continuity equation and help us find analytical solutions with simplified initial and boundary conditions. In this thesis, the one-equation Lighthill Witham and Richards (LWR) model combined with the Greenshield's model, is used for finding analytical and numerical solutions for four problems: Linear Advection, Red Traffic Light turning into Green, Stationary Shock and Shock Moving towards …

Comparison Of Mesh And Meshless Methods For Partial Differential Equations Of Galerkin Form, Wallace F. Atterberry

UNLV Theses, Dissertations, Professional Papers, and Capstones

There are two purposes of this research project. The first purpose is to compare two types of Galerkin methods: The finite element mesh method and moving least sqaures meshless Galerkin (EFG) method. The second purpose of this project is to determine if a hybrid between the mesh and meshless method is beneficial.

This manuscript will be divided into three main parts. The first part is chapter one which develops the finite element method. The second part (Chapter two) will be developing the meshless method. The last part will provide a method for combining the mesh and meshless methods for a …

Symmetric Nonconforming Mixed Finite Elements For Linear Elasticity, Jay Gopalakrishnan, Johnny Guzmán

Mathematics and Statistics Faculty Publications and Presentations

We present a family of mixed methods for linear elasticity that yield exactly symmetric, but only weakly conforming, stress approximations. The method is presented in both two and three dimensions (on triangular and tetrahedral meshes). The method is efficiently implementable by hybridization. The degrees of freedom of the Lagrange multipliers, which approximate the displacements at the faces, solve a symmetric positive-definite system. The design and analysis of this method is motivated by a new set of unisolvent degrees of freedom for symmetric polynomial matrices. These new degrees of freedom are also used to give a new simple calculation of the …

Commuting Smoothed Projectors In Weighted Norms With An Application To Axisymmetric Maxwell Equations, Jay Gopalakrishnan, Minah Oh

Mathematics and Statistics Faculty Publications and Presentations

We construct finite element projectors that can be applied to functions with low regularity. These projectors are continuous in a weighted norm arising naturally when modeling devices with axial symmetry. They have important commuting diagram properties needed for finite element analysis. As an application, we use the projectors to prove quasioptimal convergence for the edge finite element approximation of the axisymmetric time-harmonic Maxwell equations on nonsmooth domains. Supplementary numerical investigations on convergence deterioration at high wavenumbers and near Maxwell eigenvalues and are also reported.

A Second Elasticity Element Using The Matrix Bubble, Jay Gopalakrishnan, Johnny Guzmán

Mathematics and Statistics Faculty Publications and Presentations

We presented a family of finite elements that use a polynomial space augmented by certain matrix bubbles in Cockburn et al. (2010) A new elasticity element made for enforcing weak stress symmetry. Math. Comput., 79, 1331–1349 . In this sequel we exhibit a second family of elements that use the same matrix bubble. This second element uses a stress space smaller than the first while maintaining the same space for rotations (which are the Lagrange multipliers corresponding to a weak symmetry constraint). The space of displacements is of one degree less than the first method. The analysis, while similar to …

Finite Element Analysis In Porous Media For Incompressible Flow Of Contamination From Nuclear Waste, Abbas Al-Bayati, Saad A. Manaa, Ekhlass S. Ahmed

Applications and Applied Mathematics: An International Journal (AAM)

A non-linear parabolic system is used to describe incompressible nuclear waste disposal contamination in porous media, in which both molecular diffusion and dispersion are considered. The Galerkin method is applied for the pressure equation. For the brine, radionuclide and heat, a kind of partial upwind finite element scheme is constructed. Examples are included to demonstrate certain aspects of the theory and illustrate the capabilities of the kind of partial upwind finite element approach.

A Class Of Discontinuous Petrov–Galerkin Methods. Part Iv: The Optimal Test Norm And Time-Harmonic Wave Propagation In 1d., Jeffrey Zitelli, Leszek Demkowicz, Jay Gopalakrishnan, D. Pardo, V. M. Calo

Mathematics and Statistics Faculty Publications and Presentations

The phase error, or the pollution effect in the finite element solution of wave propagation problems, is a well known phenomenon that must be confronted when solving problems in the high-frequency range. This paper presents a new method with no phase errors for one-dimensional (1D) time-harmonic wave propagation problems using new ideas that hold promise for the multidimensional case. The method is constructed within the framework of the discontinuous Petrov–Galerkin (DPG) method with optimal test functions. We have previously shown that such methods select solutions that are the best possible approximations in an energy norm dual to any selected test …

A Class Of Discontinuous Petrov–Galerkin Methods. Part I: The Transport Equation, Leszek Demkowicz, Jay Gopalakrishnan

Mathematics and Statistics Faculty Publications and Presentations

Considering a simple model transport problem, we present a new finite element method. While the new method fits in the class of discontinuous Galerkin (DG) methods, it differs from standard DG and streamline diffusion methods, in that it uses a space of discontinuous trial functions tailored for stability. The new method, unlike the older approaches, yields optimal estimates for the primal variable in both the element size h and polynomial degree p, and outperforms the standard upwind DG method.

Unified Hybridization Of Discontinuous Galerkin, Mixed, And Continuous Galerkin Methods For Second Order Elliptic Problems, Bernardo Cockburn, Jay Gopalakrishnan, Raytcho Lazarov

Mathematics and Statistics Faculty Publications and Presentations

We introduce a unifying framework for hybridization of finite element methods for second order elliptic problems. The methods fitting in the framework are a general class of mixed-dual finite element methods including hybridized mixed, continu- ous Galerkin, non-conforming and a new, wide class of hybridizable discontinuous Galerkin methods. The distinctive feature of the methods in this framework is that the only globally coupled degrees of freedom are those of an approximation of the solution defined only on the boundaries of the elements. Since the associated matrix is sparse, symmetric and positive definite, these methods can be efficiently implemented. Moreover, the …