Mathematics and Statistics

Permanent URI for this collection


Recent Submissions

Now showing 1 - 2 of 2
  • Item
    (2022-03-22) Pham, Dinh Dong; Cesmelioglu, Aycil; Cheng, Eddie; Horvath, Tamas; Schmidt, Darrell; Shillor, Meir
    The main focus of this thesis is on finding highly accurate and robust numerical methods to solve a complex flow and transport problem governed by the fully-coupled time-dependent Stokes-Darcy-transport equations. This problem has many applications one of which is groundwater contamination by pollutants transported via surface/subsurface flow. It consists of two main ingredients; the time-dependent Stokes-Darcy equations describing the flow, and the time-dependent advection-diffusion equation for the transport of chemicals via this flow. Therefore, the first part of this thesis is dedicated to studying the time-dependent Stokes-Darcy problem that describes the free flow and porous media flow on two different parts of a domain and their interaction at the common interface. We introduce a hybridized discontinuous Galerkin (HDG) method which provides exact mass conservation and pressure robustness and handles the interface conditions via facet unknowns. We prove well-posedness and a priori error estimates in the energy norm, and provide numerical experiments that show optimal convergence and robustness of the method with respect to the problem parameters. The second part deals with the time-dependent advection-diffusion equation where we again use an HDG method for the spatial discretization. We show the existence and uniqueness of the semi-discrete transport problem and prove a priori error estimates in the energy norm. A number of numerical experiments are presented for different boundary conditions and we observe optimal rates of convergence in each case. Combining the two parts by a sequential algorithm, we solve the fully coupled time-dependent Stokes-Darcy-transport problem. The coupling of the flow and transport is introduced by the dependence of the fluid viscosity and source/sink terms on the concentration and by the dependence of the dispersion/diffusion tensor in the porous media domain on the advective fluid velocity. Our sequential algorithm employs a linearizing decoupling strategy based on the backward Euler time-stepping where the Stokes-Darcy and the transport equations are solved sequentially by time-lagging the concentration. The well-posedness results and a priori error estimates for the velocity and the concentration in the energy norm are presented and numerical examples demonstrating optimal convergence and mass conservation are provided.
  • Item
    (2021-12-06) Al-Khawaja, Thanaa Ali Kadhim A; Shillor, Meir; Spagnuolo, Anna Maria; Ogunyemi, Theophilus; Andrews, Kevin
    This study presents three different mathematical versions of the HANDY (Human And Nature DYnamics) model for the socioeconomic dynamics of a large stratified society. The basic model was introduced in the ground breaking publications of Motesharrei (dissertation 2014) and Motesharrei et. al. (2016). The original model consists of a nonlinear system of four ordinary differential equations (ODEs) which describe the development, in time, of a ’very simple’ society consisting of two populations: the Elite (rich) and Commoners (workers). Included also are the use of natural (renewable and nonrenewable) resources and the accumulation of human wealth. The model’s solutions depict the dynamics of these variables. Motesharrei’s main impetus and interest was to use the model as a tool for evaluating the conditions that contribute to the flourishing, sustainability, or collapse of complex societies. This dissertation expands the basic HANDY model and studies its mathematical properties and those of its three extensions. It establishes the existence of solutions to the models, as well as their uniqueness, boundedness and positivity. Furthermore, it investigates the stability of the systems’ steady states, which describe the long-time behavior of the societies. It also presents a number of qualitatively different computer simulations, providing insights into potential behaviors of the societies described by these models. The main contributions of this work are the mathematical analysis of the basic HANDY model, its three expansions and their analysis, and computer simulations. The first extension, the HANDY-SM model, includes social mobility. Rich individuals may go bankrupt and become workers, and some workers may become rich. It also allows for two different aspects of inequality, through variations in salaries and the wealth structure. The second extension, the HANDY-MC-I model, includes the Middle Class population, making the model more practical when applied to modern societies. It expands the system into five ODEs, and allows for social mobility among the three populations. Finally, in the third extension, HANDY-MC-II, two variables describe the natural resources: the renewable resources (wood, solar and wind energies), and nonrenewable resources (coal, oil, gas). This particular extension makes the model more realistic, but it also adds considerable complexity since it consists of six nonlinear coupled ODEs. The model simulations depict the consequences of having three different populations with different income status, two natural resources, and unequal contributions to wealth structure. Analysis of the models’ steady states shows that the state is stable when the populations and wealth die out but nature (the resources) is at its equilibrium. The model has also asymptotically stable, nonzero steady states to which the populations, the resources and the wealth converge in the long-time limit. The simulations also show the existence of periodic solutions in which the populations, the natural resources and wealth undergo large oscillations, indicating cycles of ‘boom and bust.’ Finally, the simulations demonstrate that the models may have chaotic solutions, pointing to a high level of unpredictability. This dissertation describes three increasingly more complex HANDY models. It paves the way and raises mathematically interesting topics for their further study. In particular, the uniqueness of the solutions, and the questions of the existence of periodic solutions, limit cycles and chaos, remain unresolved, yet. Furthermore, it suggests the possibility of tailoring such models to existing societies, and then using them as tools for evaluation of the potential outcomes of various policy decision