1,720,994 research outputs found
A Two-Level Method For The Steady-State Quasigeostrophic Equation
The quasi-geostrophic equations (QGE) are a model of large-scale ocean flows. We consider a pure stream function formulation and cite results for optimal error estimates for finding approximate solutions with the finite element method. We examine both the time dependent and steady-state versions of the equations. Numerical experiments verify the error estimates.
We examine the Argyris finite element and derive the transformation matrix necessary to perform calculations on the reference triangle. We use the Argyris element because it is a high-order, conforming finite element for fourth order problems.
In order to increase computational efficiency, we consider a two-level method to linearize the system of equations. This allows us to solve a small, nonlinear system and then use the result to linearize a larger system.Master of Scienc
Discontinuous Galerkin error estimation for linear symmetric hyperbolic systems
In this manuscript we present an error analysis for the discontinuous Galerkin discretization error of multi-dimensional first-order linear symmetric hyperbolic systems of partial differential equations. We perform a local error analysis by writing the local error as a series and showing that its leading term can be expressed as a linear combination of Legendre polynomials of degree p and p+1p+1. We apply these asymptotic results and solve relatively small local problems to compute efficient and asymptotically exact estimates of the finite element error. We present computational results for several linear hyperbolic systems in one, two and three dimensions
Galerkin methods for parabolic problems
The Galerkin finite element method of lines is one of the most popular and powerful numerical techniques for solving transient partial differential equations of parabolic type. The Galerkin finite element method of lines can be viewed as a separation-of-variables technique combined with a weak finite element formulation to discretize the problem in space. This leads to a stiff system of ordinary differential equations that can be integrated by available off-the-shelf implicit ordinary differential equations (ode) solvers based on implicit time-stepping schemes such as backward differentiation formulas (BDF) or implicit Runge-Kutta (IRK) methods
Going Beyond Counting First Authors in Author Co-citation Analysis
The present study examines one of the fundamental aspects of author co-citation analysis (ACA) - the way co-citation
counts are defined. Co-citation counting provides the data on which all subsequent statistical analyses and mappings
are based, and we compare ACA results based on two different types of co-citation counting - the traditional type that
only counts the first one among a cited work's authors on the one hand and a non-traditional type that takes into
account the first 5 authors of a cited work on the other hand. Results indicate that the picture produced through this non-traditional author co-citation counting contains more coherent author groups and is therefore considerably clearer. However, this picture represents fewer specialties in the research field being studied than that produced through the traditional first-author co-citation counting when the same number of top-ranked authors is selected and analyzed. Reasons for these effects are discussed
Formulations, Issues and Comparison of Car-Following Models
Microscopic simulation software use car-following models to capture the interaction of a vehicle and the preceding vehicle traveling in the same lane. In the literature, much research has been carried out in the field of car-following and traffic stream modeling. Microscopic car-following models have been characterized by using the relationship between a vehicle's desired speed and the distance headway (h) between the lead and follower vehicles. On the other hand, macroscopic traffic stream models describe the motion of a traffic stream by approximating for the flow of a continuous compressible fluid. This research work develops and compares three different formulations of car-following models — speed formulation, molecular acceleration, and fluid acceleration formulation. First, four state-of-the-art car-following models namely, Van Aerde, Greenshields, Greenberg and Pipes models, are selected for developing the three aforementioned formulations. Then a comprehensive car-following behavior encompassing steady-state conditions and two constraints — acceleration and collision avoidance — is presented. Specifically, the variable power vehicle dynamics model proposed by Rakha and Lucic (2002) is utilized for the acceleration constraint.
Subsequently, the thesis describes the issues associated with car-following formulations. Recognizing that many different traffic flow conditions exist, three distinct scenarios are selected for comparison purposes. The results demonstrate that the speed formulation ensures that vehicles typically revert to steady-state conditions when vehicles experience a perturbation from steady-state conditions. On the other hand, both acceleration formulations are unable to converge to steady-state conditions when the system experiences a perturbation from a steady-state.
The thesis also attempts to address the question of capacity drop associated with vehicles accelerating from congested conditions. Specifically, the capacity drop proposition is analyzed for the case of a backward recovery (typical of a signalized intersection) and stationary shockwave (typical of a capacity drop on a freeway). In the case of the backward recovery shockwave, the acceleration constraint results in a temporally and spatially confined capacity drop as vehicles accelerate to their desired steady-state speed. This temporally and spatially confined capacity drop results in what is typically termed the start loss of a signalized phase. Subsequently, vehicles attain steady-state conditions, in the case of the speed and molecular acceleration formulations, at the traffic signal stop bar after the initial five vehicle departures. The analysis also demonstrates that after attaining steady-state conditions the capacity may drop for the initial vehicle departures as a result of traffic stream dispersion. This traffic dispersion capacity drop increases as vehicles travel further downstream. Alternatively, in the case of a stationary bottleneck the aggressiveness of vehicle accelerations plays a major role in defining the capacity drop downstream of a bottleneck. The study demonstrates that any temporal headways that may be lost while vehicles accelerate to steady-state conditions may not be recuperated and thus result in capacity drops downstream of a bottleneck. A typical example of this scenario is the traffic stream flow rate downstream of a stop sign, which is significantly less than the roadway capacity. The reduction in capacity is caused by losses in temporal headways between successive vehicles which are not recuperated. The study also demonstrates that the ability to model such a capacity drop does not require the use of a dual-regime traffic stream model as is proposed in the Highway Capacity Manual (HCM). Instead, the use of a single-regime model captures the observed capacity with the introduction of an acceleration constraint to the car-following system of equations.Master of Scienc
Design, Analysis, and Application of Immersed Finite Element Methods
This dissertation consists of three studies of immersed finite element (IFE) methods for inter- face problems related to partial differential equations (PDEs) with discontinuous coefficients. These three topics together form a continuation of the research in IFE method including the extension to elasticity systems, new breakthroughs to higher degree IFE methods, and its application to inverse problems.
First, we extend the current construction and analysis approach of IFE methods in the literature for scalar elliptic equations to elasticity systems in the vector format. In particular, we construct a group of low-degree IFE functions formed by linear, bilinear, and rotated Q1 polynomials to weakly satisfy the jump conditions of elasticity interface problems. Then we analyze the trace inequalities of these IFE functions and the approximation capabilities of the resulted IFE spaces. Based on these preparations, we develop a partially penalized IFE (PPIFE) scheme and prove its optimal convergence rates.
Secondly, we discuss the limitations of the current approaches of IFE methods when we try to extend them to higher degree IFE methods. Then we develop a new framework to construct and analyze arbitrary p-th degree IFE methods. In this framework, each IFE function is the extension of a p-th degree polynomial from one subelement to the whole interface element by solving a local Cauchy problem on interface elements in which the jump conditions across the interface are employed as the boundary conditions. All the components in the analysis, including existence of IFE functions, the optimal approximation capabilities and the trace inequalities, are all reduced to key properties of the related discrete extension operator. We employ these results to show the optimal convergence of a discontinuous Galerkin IFE (DGIFE) method.
In the last part, we apply the linear IFE methods in the literature together with the shape optimization technique to solve a group of interface inverse problems. In this algorithm, both the governing PDEs and the objective functional for interface inverse problems are discretized optimally by the IFE method regardless of the location of the interface in a chosen mesh. We derive the formulas for the gradients of the objective function in the optimization problem which can be implemented efficiently in the IFE framework through a discrete adjoint method. We demonstrate the properties of the proposed algorithm by applying it to three representative applications.Doctor of PhilosophyInterface problems arise from many science and engineering applications modeling the transmission of some physical quantities between multiple materials. Mathematically, these multiple materials in general are modeled by partial differential equations (PDEs) with discontinuous parameters, which poses challenges to developing efficient and reliable numerical methods and the related theoretical error analysis. The main contributions of this dissertation is on the development of a special finite element method, the so called immersed finite element (IFE) method, to solve the interface problems on a mesh independent of the interface geometry which can be advantageous especially when the interface is moving. Specifically, this dissertation consists of three projects of IFE methods: elasticity interface problems, higher-order IFE methods and interface inverse problems, including their design, analysis, and application
Galerkin Projections Between Finite Element Spaces
Adaptive mesh refinement schemes are used to find accurate low-dimensional approximating spaces when solving elliptic PDEs with Galerkin finite element methods. For nonlinear PDEs, solving the nonlinear problem with Newton's method requires an initial guess of the solution on a refined space, which can be found by interpolating the solution from a previous refinement. Improving the accuracy of the representation of the converged solution computed on a coarse mesh for use as an initial guess on the refined mesh may reduce the number of Newton iterations required for convergence. In this thesis, we present an algorithm to compute an orthogonal L^2 projection between two dimensional finite element spaces constructed from a triangulation of the domain. Furthermore, we present numerical studies that investigate the efficiency of using this algorithm to solve various nonlinear elliptic boundary value problems.Master of Scienc
Immersed Finite Elements for a Second Order Elliptic Operator and Their Applications
This dissertation studies immersed finite elements (IFE) for a second order elliptic operator and their applications to interface problems of related partial differential equations.
We start with the immersed finite element methods for the second order elliptic operator with a discontinuous coefficient associated with the elliptic interface problems. We introduce an energy norm stronger than the one used in [111]. Then we derive an estimate for the IFE interpolation error with this energy norm using patches of interface elements. We prove both the continuity and coercivity of the bilinear form in a partially penalized IFE (PPIFE) method. These properties allow us to derive an error bound for the PPIFE solution in the energy norm under the standard piecewise regularity assumption instead of the more stringent regularity used in [111]. As an important consequence, this new estimation further enables us to show the optimal convergence in the norm which could not be done by the analysis presented in [111].
Then we consider applications of IFEs developed for the second order elliptic operator to wave propagation and diffusion interface problems. The first application is for the time-harmonic wave interface problem that involves the Helmholtz equation with a discontinuous coefficient. We design PPIFE and DGIFE schemes including the higher degree IFEs for Helmholtz interface problems. We present an error analysis for the symmetric linear/bilinear PPIFE methods. Under the standard piecewise regularity assumption for the exact solution, following Schatz's arguments, we derive optimal error bounds for the PPIFE solutions in both an energy norm and the usual norm provided that the mesh size is sufficiently small.
{In the second group of applications, we focus on the error analysis for IFE methods developed for solving typical time-dependent interface problems associated with the second order elliptic operator with a discontinuous coefficient.} For hyperbolic interface problems, which are typical wave propagation interface problems, we reanalyze the fully-discrete PPIFE method in [143]. We derive the optimal error bounds for this PPIFE method for both an energy norm and the norm under the standard piecewise regularity assumption in the space variable of the exact solution. Simulations for standing and travelling waves are presented to corroborate the results of the error analysis. For parabolic interface problems, which are typical diffusion interface problems, we reanalyze the PPIFE methods in [113]. We prove that these PPIFE methods have the optimal convergence not only in an energy norm but also in the usual norm under the standard piecewise regularity.Doctor of PhilosophyThis dissertation studies immersed finite elements (IFE) for a second order elliptic operator and their applications to a few types of interface problems.
We start with the immersed finite element methods for the second order elliptic operator with a discontinuous coefficient associated with the elliptic interface problem. We can show that the IFE methods for the elliptic interface problems converge optimally when the exact solution has lower regularity than that in the previous publications.
Then we consider applications of IFEs developed for the second order elliptic operator to wave propagation and diffusion interface problems. For interface problems of the Helmholtz equation which models time-Harmonic wave propagations, we design IFE schemes, including higher degree schemes, and derive error estimates for a lower degree scheme. For interface problems of the second order hyperbolic equation which models time dependent wave propagations, we derive better error estimates for the IFE methods and provides numerical simulations for both the standing and traveling waves. For interface problems of the parabolic equation which models the time dependent diffusion, we also derive better error estimates for the IFE methods
Variations on the Author
“Variations on the Author” discusses two of Eduardo Coutinho’s recent films (Um Dia na Vida, from 2010, and Últimas Conversas, posthumously released in 2015) and their contribution to the general question of documentary authorship. The director’s filmography is characterized by a consistent yet self-effacing form of authorial self-inscription: Coutinho often features as an interviewer that rather than express opinions propels discourses; an interviewer that is good at listening. This mode of self-inscription characterizes him as an author who is not expressive but who is nonetheless markedly present on the screen. In Um Dia na Vida, however, Coutinho is completely absent form the image, while Últimas Conversas, on the contrary, includes a confessional prologue that moves the director from the margins to the center of his films. This article examines the ways in which these works stand out in the filmography of a director who offers new insights into the notion of cinematic authorship
- …
