BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

K. B. Nakshatrala

K. B. Nakshatrala appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
Numerical Analysismath.NAComputational Engineering, Finance, and Sciencemath-phmath.MPcond-mat.mtrl-scicond-mat.softphysics.class-phphysics.flu-dyn

Catalog footprint

What is connected

23works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2022arXiv

How to pose material design problems for flow through porous media applications?: Sensitivity of dissipation rate to medium's permeability holds the key

Recent studies have advocated using the total dissipation rate under topology optimization to realize material designs involving the flow of fluids through porous media. However, these studies decided how to pose the design problem, such as maximizing the total dissipation rate for some situations while minimizing for others, by solving one-dimensional problems and justifying their choices using numerical experiments. The rigor is lacking -- a bottleneck for further scientific advancements to computational material design. This paper provides the missing theoretical justification. We identify four classes of boundary value problems using the adjoint state method and analytically calculate the sensitivity of the total dissipation rate to the permeability field. For two of those classes in which the flow of fluids is pressure-driven, the sensitivity is positive -- the total dissipation rate increases if the medium's permeability increases. While for the other two classes, in which the flow is velocity-driven, the trend is the opposite. These sensitivities provide rigorous answers to the central question: how to pose a material design problem for flow through porous media applications. The impact of our work is multi-fold. First, this study further elevates the role of the dissipation rate in posing well-posed material design problems using topology optimization. Second, besides the theoretical significance, the results benefit computational scientists and practitioners to realize optimal designs. Third, given their simplicity yet far-reaching impact, both the approach and results possess immense pedagogical value.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Material degradation due to moisture and temperature. Part 1: Mathematical model, analysis, and analytical solutions

The mechanical response, serviceability, and load bearing capacity of materials and structural components can be adversely affected due to external stimuli, which include exposure to a corrosive chemical species, high temperatures, temperature fluctuations (i.e., freezing-thawing), cyclic mechanical loading, just to name a few. It is, therefore, of paramount importance in several branches of engineering -- ranging from aerospace engineering, civil engineering to biomedical engineering -- to have a fundamental understanding of degradation of materials, as the materials in these applications are often subjected to adverse environments. As a result of recent advancements in material science, new materials like fiber-reinforced polymers and multi-functional materials that exhibit high ductility have been developed and widely used; for example, as infrastructural materials or in medical devices (e.g., stents). The traditional small-strain approaches of modeling these materials will not be adequate. In this paper, we study degradation of materials due to an exposure to chemical species and temperature under large-strain and large-deformations. In the first part of our research work, we present a consistent mathematical model with firm thermodynamic underpinning. We then obtain semi-analytical solutions of several canonical problems to illustrate the nature of the quasi-static and unsteady behaviors of degrading hyperelastic solids.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mechanics-based solution verification for porous media models

This paper presents a new approach to verify accuracy of computational simulations. We develop mathematical theorems which can serve as robust a posteriori error estimation techniques to identify numerical pollution, check the performance of adaptive meshes, and verify numerical solutions. We demonstrate performance of this methodology on problems from flow thorough porous media. However, one can extend it to other models. We construct mathematical properties such that the solutions to Darcy and Darcy-Brinkman equations satisfy them. The mathematical properties include the total minimum mechanical power, minimum dissipation theorem, reciprocal relation, and maximum principle for the vorticity. All the developed theorems have firm mechanical bases and are independent of numerical methods. So, these can be utilized for solution verification of finite element, finite volume, finite difference, lattice Boltzmann methods and so forth. In particular, we show that, for a given set of boundary conditions, Darcy velocity has the minimum total mechanical power of all the kinematically admissible vector fields. We also show that a similar result holds for Darcy-Brinkman velocity. We then show for a conservative body force, the Darcy and Darcy-Brinkman velocities have the minimum total dissipation among their respective kinematically admissible vector fields. Using numerical examples, we show that the minimum dissipation and total mechanical power theorems can be utilized to identify pollution errors in numerical solutions. The solutions to Darcy and Darcy-Brinkman equations are shown to satisfy a reciprocal relation, which has the potential to identify errors in the numerical implementation of boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Do current lattice Boltzmann methods for diffusion and diffusion-type equations respect maximum principles and the non-negative constraint?

The lattice Boltzmann method (LBM) has established itself as a valid numerical method in computational fluid dynamics. Recently, multiple-relaxation-time LBM has been proposed to simulate anisotropic advection-diffusion processes. The governing differential equations of advective-diffusive systems are known to satisfy maximum principles, comparison principles, the non-negative constraint, and the decay property. In this paper, it will be shown that current single- and multiple-relaxation-time lattice Boltzmann methods fail to preserve these mathematical properties for transient diffusion-type equations. It will also be shown that the discretization of Dirichlet boundary conditions will affect the performance of lattice Boltzmann methods in meeting these mathematical principles. A new way of discretizing the Dirichlet boundary conditions is also proposed. Several benchmark problems have been solved to illustrate the performance of lattice Boltzmann methods and the effect of discretization of boundary conditions with respect to the aforementioned mathematical properties for transient diffusion and advection-diffusion equations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On enforcing maximum principles and achieving element-wise species balance for advection-diffusion-reaction equations under the finite element method

We present a robust computational framework for advective-diffusive-reactive systems that satisfies maximum principles, the non-negative constraint, and element-wise species balance property. The proposed methodology is valid on general computational grids, can handle heterogeneous anisotropic media, and provides accurate numerical solutions even for very high Péclet numbers. The significant contribution of this paper is to incorporate advection (which makes the spatial part of the differential operator non-self-adjoint) into the non-negative computational framework, and overcome numerical challenges associated with advection. We employ low-order mixed finite element formulations based on least-squares formalism, and enforce explicit constraints on the discrete problem to meet the desired properties. The resulting constrained discrete problem belongs to convex quadratic programming for which a unique solution exists. Maximum principles and the non-negative constraint give rise to bound constraints while element-wise species balance gives rise to equality constraints. The resulting convex quadratic programming problems are solved using an interior-point algorithm. Several numerical results pertaining to advection-dominated problems are presented to illustrate the robustness, convergence, and the overall performance of the proposed computational framework.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On mesh restrictions to satisfy comparison principles, maximum principles, and the non-negative constraint: Recent developments and new results

This paper concerns with mesh restrictions that are needed to satisfy several important mathematical properties -- maximum principles, comparison principles, and the non-negative constraint -- for a general linear second-order elliptic partial differential equation. We critically review some recent developments in the field of discrete maximum principles, derive new results, and discuss some possible future research directions in this area. In particular, we derive restrictions for a three-node triangular (T3) element and a four-node quadrilateral (Q4) element to satisfy comparison principles, maximum principles, and the non-negative constraint under the standard single-field Galerkin formulation. Analysis is restricted to uniformly elliptic linear differential operators in divergence form with Dirichlet boundary conditions specified on the entire boundary of the domain. Various versions of maximum principles and comparison principles are discussed in both continuous and discrete settings. In the literature, it is well-known that an acute-angled triangle is sufficient to satisfy the discrete weak maximum principle for pure isotropic diffusion. An iterative algorithm is developed to construct simplicial meshes that preserves discrete maximum principles using existing open source mesh generators. Various numerical examples based on different types of triangulations are presented to show the pros and cons of placing restrictions on a computational mesh. We also quantify local and global mass conservation errors using representative numerical examples, and illustrate the performance of metric-based meshes with respect to mass conservation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A monolithic multi-time-step computational framework for first-order transient systems with disparate scales

Developing robust simulation tools for problems involving multiple mathematical scales has been a subject of great interest in computational mathematics and engineering. A desirable feature to have in a numerical formulation for multiscale transient problems is to be able to employ different time-steps (multi-time-step coupling), and different time integrators and different numerical formulations (mixed methods) in different regions of the computational domain. We present two new monolithic multi-time-step mixed coupling methods for first-order transient systems. We shall employ unsteady advection-diffusion-reaction equation with linear decay as the model problem, which offers several unique challenges in terms of non-self-adjoint spatial operator and rich features in the solutions. We shall employ the dual Schur domain decomposition technique to handle the decomposition of domain into subdomains. Two different methods of enforcing compatibility along the subdomain interface will be used in the time discrete setting. A systematic theoretical analysis (which includes numerical stability, influence of perturbations, bounds on drift along the subdomain interface) will be performed. The first coupling method ensures that there is no drift along the subdomain interface but does not facilitate explicit/implicit coupling. The second coupling method allows explicit/implicit coupling with controlled (but non-zero) drift in the solution along the subdomain interface. Several canonical problems will be solved to numerically verify the theoretical predictions, and to illustrate the overall performance of the proposed coupling methods. Finally, we shall illustrate the robustness of the proposed coupling methods using a multi-time-step transient simulation of a fast bimolecular advective-diffusive-reactive system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On multi-time-step monolithic coupling algorithms for elastodynamics

We present a way of constructing multi-time-step monolithic coupling methods for elastodynamics. The governing equations for constrained multiple subdomains are written in dual Schur form and enforce the continuity of velocities at system time levels. The resulting equations will be in the form of differential-algebraic equations. To crystallize the ideas we shall employ Newmark family of time-stepping schemes. The proposed method can handle multiple subdomains, and allows different time-steps as well as different time stepping schemes from the Newmark family in different subdomains. We shall use the energy method to assess the numerical stability, and quantify the influence of perturbations under the proposed coupling method. We also discuss the conditions under which the proposed method will be energy preserving, and the conditions under which the method will be energy conserving. Several numerical examples are presented to illustrate the accuracy and stability properties of the proposed method. We shall also compare the proposed multi-time-step coupling method with some other similar methods available in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A framework for coupling flow and deformation of the porous solid

In this paper, we consider the flow of an incompressible fluid in a deformable porous solid. We present a mathematical model using the framework offered by the theory of interacting continua. In its most general form, this framework provides a mechanism for capturing multiphase flow, deformation, chemical reactions and thermal processes, as well as interactions between the various physics in a conveniently implemented fashion. To simplify the presentation of the framework, results are presented for a particular model than can be seen as an extension of Darcy's equation (which assumes that the porous solid is rigid) that takes into account elastic deformation of the porous solid. The model also considers the effect of deformation on porosity. We show that using this model one can recover identical results as in the framework proposed by Biot and Terzaghi. Some salient features of the framework are as follows: (a) It is a consistent mixture theory model, and adheres to the laws and principles of continuum thermodynamics, (b) the model is capable of simulating various important phenomena like consolidation and surface subsidence, and (c) the model is amenable to several extensions. We also present numerical coupling algorithms to obtain coupled flow-deformation response. Several representative numerical examples are presented to illustrate the capability of the mathematical model and the performance of the computational framework.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A mixed formulation for a modification to Darcy equation based on Picard linearization and numerical solutions to large-scale realistic problems

In this paper we consider a modification to Darcy equation by taking into account the dependence of viscosity on the pressure. We present a stabilized mixed formulation for the resulting governing equations. Equal-order interpolation for the velocity and pressure is considered, and shown to be stable (which is not the case under the classical mixed formulation). The proposed mixed formulation is tested using a wide variety of numerical examples. The proposed formulation is also implemented in a parallel setting, and the performance of the formulation for large-scale problems is illustrated using a representative problem. Two practical and technologically important problems, one each on enhanced oil recovery and geological carbon-dioxide sequestration, are solved using the proposed formulation. The numerical examples show that the predictions based on Darcy model are qualitatively and quantitatively different from that of the predictions based on the modified Darcy model, which takes into account the dependence of the viscosity on the pressure. In particular, the numerical example on the geological carbon-dioxide sequestration shows that Darcy model over-predicts the leakage into an abandoned well when compared to that of the modified Darcy model. On the other hand, the modified Darcy model predicts higher pressures and higher pressure gradients near the injection well. These predictions have dire consequences in predicting damage and fracture zones, and designing the seal, whose integrity is crucial to the safety of a geological carbon-dioxide sequestration geosystem.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A numerical framework for diffusion-controlled bimolecular-reactive systems to enforce maximum principles and non-negative constraint

We present a novel computational framework for diffusive-reactive systems that satisfies the non-negative constraint and maximum principles on general computational grids. The governing equations for the concentration of reactants and product are written in terms of tensorial diffusion-reaction equations. % We restrict our studies to fast irreversible bimolecular reactions. If one assumes that the reaction is diffusion-limited and all chemical species have the same diffusion coefficient, one can employ a linear transformation to rewrite the governing equations in terms of invariants, which are unaffected by the reaction. This results in two uncoupled tensorial diffusion equations in terms of these invariants, which are solved using a novel non-negative solver for tensorial diffusion-type equations. The concentrations of the reactants and the product are then calculated from invariants using algebraic manipulations. The novel aspect of the proposed computational framework is that it will always produce physically meaningful non-negative values for the concentrations of all chemical species. Several representative numerical examples are presented to illustrate the robustness, convergence, and the numerical performance of the proposed computational framework. We will also compare the proposed framework with other popular formulations. In particular, we will show that the Galerkin formulation (which is the standard single-field formulation) does not produce reliable solutions, and the reason can be attributed to the fact that the single-field formulation does not guarantee non-negative solutions. We will also show that the clipping procedure (which produces non-negative solutions but is considered as a variational crime) does not give accurate results when compared with the proposed computational framework.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A numerical methodology for enforcing maximum principles and the non-negative constraint for transient diffusion equations

Transient diffusion equations arise in many branches of engineering and applied sciences (e.g., heat transfer and mass transfer), and are parabolic partial differential equations. It is well-known that, under certain assumptions on the input data, these equations satisfy important mathematical properties like maximum principles and the non-negative constraint, which have implications in mathematical modeling. However, existing numerical formulations for these types of equations do not, in general, satisfy maximum principles and the non-negative constraint. In this paper, we present a methodology for enforcing maximum principles and the non-negative constraint for transient anisotropic diffusion equation. The method of horizontal lines (also known as the Rothe method) is applied in which the time is discretized first. This results in solving steady anisotropic diffusion equation with decay equation at every discrete time level. The proposed methodology for transient anisotropic diffusion equation will satisfy maximum principles and the non-negative constraint on general computational grids, and with no additional restrictions on the time step. We illustrate the performance and accuracy of the proposed formulation using representative numerical examples. We also perform numerical convergence of the proposed methodology. For comparison, we also present the results from the standard single-field semi-discrete formulation and the results from a popular software package, which all will violate maximum principles and the non-negative constraint.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modification to Darcy model for high pressure and high velocity applications and associated mixed finite element formulations

The Darcy model is based on a plethora of assumptions. One of the most important assumptions is that the Darcy model assumes the drag coefficient to be constant. However, there is irrefutable experimental evidence that viscosities of organic liquids and carbon-dioxide depend on the pressure. Experiments have also shown that the drag varies nonlinearly with respect to the velocity at high flow rates. In important technological applications like enhanced oil recovery and geological carbon-dioxide sequestration, one encounters both high pressures and high flow rates. It should be emphasized that flow characteristics and pressure variation under varying drag are both quantitatively and qualitatively different from that of constant drag. Motivated by experimental evidence, we consider the drag coefficient to depend on both the pressure and velocity. We consider two major modifications to the Darcy model based on the Barus formula and Forchheimer approximation. The proposed modifications to the Darcy model result in nonlinear partial differential equations, which are not amenable to analytical solutions. To this end, we present mixed finite element formulations based on least-squares formalism and variational multiscale formalism for the resulting governing equations. The proposed modifications to the Darcy model and its associated finite element formulations are used to solve realistic problems with relevance to enhanced oil recovery. We also study the competition between the nonlinear dependence of drag on the velocity and the dependence of viscosity on the pressure. To the best of the authors' knowledge such a systematic study has not been performed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the performance of high-order finite elements with respect to maximum principles and the non-negative constraint for diffusion-type equations

The main aim of this paper is to document the performance of $p$-refinement with respect to maximum principles and the non-negative constraint. The model problem is (steady-state) anisotropic diffusion with decay (which is a second-order elliptic partial differential equation). We considered the standard single-field formulation (which is based on the Galerkin formalism) and two least-squares-based mixed formulations. We have employed non-uniform Lagrange polynomials for altering the polynomial order in each element, and we have used $p = 1, ..., 10$. It will be shown that the violation of the non-negative constraint will not vanish with $p$-refinement for anisotropic diffusion. We shall illustrate the performance of $p$-refinement using several representative problems. The intended outcome of the paper is twofold. Firstly, this study will caution the users of high-order approximations about its performance with respect to maximum principles and the non-negative constraint. Secondly, this study will help researchers to develop new methodologies for enforcing maximum principles and the non-negative constraint under high-order approximations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A numerical study of fluids with pressure dependent viscosity flowing through a rigid porous medium

In this paper we consider modifications to Darcy's equation wherein the drag coefficient is a function of pressure, which is a realistic model for technological applications like enhanced oil recovery and geological carbon sequestration. We first outline the approximations behind Darcy's equation and the modifications that we propose to Darcy's equation, and derive the governing equations through a systematic approach using mixture theory. We then propose a stabilized mixed finite element formulation for the modified Darcy's equation. To solve the resulting nonlinear equations we present a solution procedure based on the consistent Newton-Raphson method. We solve representative test problems to illustrate the performance of the proposed stabilized formulation. One of the objectives of this paper is also to show that the dependence of viscosity on the pressure can have a significant effect both on the qualitative and quantitative nature of the solution.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A stabilized mixed formulation for unsteady Brinkman equation based on the method of horizontal lines

In this paper, we present a stabilized mixed formulation for unsteady Brinkman equation. The formulation is systematically derived based on the variational multiscale formalism and the method of horizontal lines. The derivation does not need the assumption that the fine-scale variables do not depend on the time, which is the case with the conventional derivation of multiscale stabilized formulations for transient mixed problems. An expression for the stabilization parameter is obtained in terms of a bubble function, and appropriate bubble functions for various finite elements are also presented. Under the proposed formulation, equal-order interpolation for the velocity and pressure (which is computationally the most convenient) is stable. Representative numerical results are presented to illustrate the performance of the proposed formulation. Spatial and temporal convergence studies are also performed, and the proposed formulation performed well.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Enforcing the non-negativity constraint and maximum principles for diffusion with decay on general computational grids

In this paper, we consider anisotropic diffusion with decay, and the diffusivity coefficient to be a second-order symmetric and positive definite tensor. It is well-known that this particular equation is a second-order elliptic equation, and satisfies a maximum principle under certain regularity assumptions. However, the finite element implementation of the classical Galerkin formulation for both anisotropic and isotropic diffusion with decay does not respect the maximum principle. We first show that the numerical accuracy of the classical Galerkin formulation deteriorates dramatically with increase in the decay coefficient for isotropic medium and violates the discrete maximum principle. However, in the case of isotropic medium, the extent of violation decreases with mesh refinement. We then show that, in the case of anisotropic medium, the classical Galerkin formulation for anisotropic diffusion with decay violates the discrete maximum principle even at lower values of decay coefficient and does not vanish with mesh refinement. We then present a methodology for enforcing maximum principles under the classical Galerkin formulation for anisotropic diffusion with decay on general computational grids using optimization techniques. Representative numerical results (which take into account anisotropy and heterogeneity) are presented to illustrate the performance of the proposed formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the performance of the variational multiscale formulation for subsurface flow and transport in heterogeneous porous media

The following work compares two popular mixed finite elements used to model subsurface flow and transport in heterogeneous porous media; the lowest order Raviart-Thomas element and the variational multiscale stabilized element. Comparison is made based on performance for several problems of engineering relevance that involve highly heterogenous material properties (permeability ratios of up to $1\times10^5$), open flow boundary conditions (pressure driven flows), and large scale domains in two dimensions. Numerical experiments are performed to show the degree to which mass conservation is violated when a flow field computed using either element is used as the advection velocity in a transport model. The results reveal that the variational multiscale element shows considerable mass production or loss for problems that involve flow tangential to layers of differing permeability, but marginal violation of local mass balance for problems of less orthogonality in the permeability. The results are useful in establishing rudimentary estimates of the error produced by using the variational mutliscale element for several different types of problems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the vibrations of lumped parameter systems governed by differential-algebraic equations

In this paper, we consider the vibratory motions of lumped parameter systems wherein the components of the system cannot be described by constitutive expressions for the force in terms of appropriate kinematical quantities. Such physical systems reduce to a system of differential-algebraic equations, which invariably need to be solved numerically. To illustrate the issues with clarity, we consider a simple system in which the dashpot is assumed to contain a "Bingham" fluid for which one cannot describe the force in the dashpot as a function of the velocity. On the other hand, one can express the velocity as a function of the force.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Variational structure of the optimal artificial diffusion method for the advection-diffusion equation

In this research note we provide a variational basis for the optimal artificial diffusion method, which has been a cornerstone in developing many stabilized methods. The optimal artificial diffusion method produces exact nodal solutions when applied to one-dimensional problems with constant coefficients and forcing function. We first present a variational principle for a multi-dimensional advective-diffusive system, and then derive a new stable weak formulation. When applied to one-dimensional problems with constant coefficients and forcing function, this resulting weak formulation will be equivalent to the optimal artificial diffusion method. We present representative numerical results to corroborate our theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Non-negative mixed finite element formulations for a tensorial diffusion equation

We consider the tensorial diffusion equation, and address the discrete maximum-minimum principle of mixed finite element formulations. In particular, we address non-negative solutions (which is a special case of the maximum-minimum principle) of mixed finite element formulations. The discrete maximum-minimum principle is the discrete version of the maximum-minimum principle. In this paper we present two non-negative mixed finite element formulations for tensorial diffusion equations based on constrained optimization techniques (in particular, quadratic programming). These proposed mixed formulations produce non-negative numerical solutions on arbitrary meshes for low-order (i.e., linear, bilinear and trilinear) finite elements. The first formulation is based on the Raviart-Thomas spaces, and is obtained by adding a non-negative constraint to the variational statement of the Raviart-Thomas formulation. The second non-negative formulation based on the variational multiscale formulation. For the former formulation we comment on the affect of adding the non-negative constraint on the local mass balance property of the Raviart-Thomas formulation. We also study the performance of the active set strategy for solving the resulting constrained optimization problems. The overall performance of the proposed formulation is illustrated on three canonical test problems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On dual Schur domain decomposition method for linear first-order transient problems

This paper addresses some numerical and theoretical aspects of dual Schur domain decomposition methods for linear first-order transient partial differential equations. In this work, we consider the trapezoidal family of schemes for integrating the ordinary differential equations (ODEs) for each subdomain and present four different coupling methods, corresponding to different algebraic constraints, for enforcing kinematic continuity on the interface between the subdomains. Method 1 (d-continuity) is based on the conventional approach using continuity of the primary variable and we show that this method is unstable for a lot of commonly used time integrators including the mid-point rule. To alleviate this difficulty, we propose a new Method 2 (Modified d-continuity) and prove its stability for coupling all time integrators in the trapezoidal family (except the forward Euler). Method 3 (v-continuity) is based on enforcing the continuity of the time derivative of the primary variable. However, this constraint introduces a drift in the primary variable on the interface. We present Method 4 (Baumgarte stabilized) which uses Baumgarte stabilization to limit this drift and we derive bounds for the stabilization parameter to ensure stability. Our stability analysis is based on the ``energy'' method, and one of the main contributions of this paper is the extension of the energy method (which was previously introduced in the context of numerical methods for ODEs) to assess the stability of numerical formulations for index-2 differential-algebraic equations (DAEs).

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

On the stability of bubble functions and a stabilized mixed finite element formulation for the Stokes problem

In this paper we investigate the relationship between stabilized and enriched finite element formulations for the Stokes problem. We also present a new stabilized mixed formulation for which the stability parameter is derived purely by the method of weighted residuals. This new formulation allows equal order interpolation for the velocity and pressure fields. Finally, we show by counterexample that a direct equivalence between subgrid-based stabilized finite element methods and Galerkin methods enriched by bubble functions cannot be constructed for quadrilateral and hexahedral elements using standard bubble functions.

0 citations0 reviewsNo code linkedOpen access