Researcher profile

Neela Nataraj

Neela Nataraj contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

7 published item(s)

preprint2022arXiv

Lowest-order equivalent nonstandard finite element methods for biharmonic plates

The popular (piecewise) quadratic schemes for the biharmonic equation based on triangles are the nonconforming Morley finite element, the discontinuous Galerkin, the $C^0$ interior penalty, and the WOPSIP schemes. Those methods are modified in their right-hand side $F\in H^{-2}(Ω)$ replaced by $F\circ (JI_{\rm M}) $ and then are quasi-optimal in their respective discrete norms. The smoother $JI_{\rm M}$ is defined for a piecewise smooth input function by a (generalized) Morley interpolation $I_{\rm M}$ followed by a companion operator $J$. An abstract framework for the error analysis in the energy, weaker and piecewise Sobolev norms for the schemes is outlined and applied to the biharmonic equation. Three errors are also equivalent in some particular discrete norm from [Carstensen, Gallistl, Nataraj: Comparison results of nonstandard $P_2$ finite element methods for the biharmonic problem, ESAIM Math. Model. Numer. Anal. (2015)] without data oscillations. This paper extends the work [Veeser, Zanotti: Quasi-optimal nonconforming methods for symmetric elliptic problems, SIAM J. Numer. Anal. 56 (2018)] to the discontinuous Galerkin scheme and adds error estimates in weaker and piecewise Sobolev norms.

preprint2022arXiv

Stability of mixed FEMs for non-selfadjoint indefinite second-order linear elliptic PDEs

For a well-posed non-selfadjoint indefinite second-order linear elliptic PDE with general coefficients $\mathbf A, \mathbf b,γ$ in $L^\infty$ and symmetric and uniformly positive definite coefficient matrix $\mathbf A$, this paper proves that mixed finite element problems are uniquely solvable and the discrete solutions are uniformly bounded, whenever the underlying shape-regular triangulation is sufficiently fine. This applies to the Raviart-Thomas (RT) and Brezzi-Douglas-Marini (BDM) finite element families of any order and in any space dimension and leads to the best-approximation estimate in $H(div)\times L^2$ as well as in in $L^2\times L^2$ up to oscillations. This generalises earlier contributions for piecewise Lipschitz continuous coefficients to $L^\infty$ coefficients. The compactness argument of Schatz and Wang for the displacement-oriented problem does not apply immediately to the mixed formulation in $H(div)\times L^2$. But it allows the uniform approximation of some $L^2$ contributions and can be combined with a recent $L^2$ best-approximation result from the medius analysis. This technique circumvents any regularity assumption and the application of a Fortin interpolation operator.

preprint2021arXiv

A priori and a posteriori error analysis of the Crouzeix-Raviart and Morley FEM with original and modified righthand sides

This article on nonconforming schemes for $m$ harmonic problems simultaneously treats the Crouzeix-Raviart ($m=1$) and the Morley finite elements ($m=2$) for the original and for modified right-hand side $F$ in the dual space $V^*:=H^{-m}(Ω)$ to the energy space $V:=H^{m}_0(Ω)$. The smoother $J:V_{\rm nc} \to V$ in this paper is a companion operator, that is a linear and bounded right-inverse to the nonconforming interpolation operator $I_{\rm nc} :V\to V_{\rm nc}$, and modifies the discrete right-hand side $F_h:=F\circ J \in V_{\rm nc}^*$. The best-approximation property of the modified scheme from Veeser et al. (2018) is recovered and complemented with an analysis of the convergence rates in weaker Sobolev norms. Examples with oscillating data show that the original method may fail to enjoy the best-approximation property but can also be better than the modified scheme. The a~posteriori analysis of this paper concerns data oscillations of various types in a class of right-hand sides $F\in V^*$. The reliable error estimates involve explicit constants and can be recommended for explicit error control of the piecewise energy norm. The efficiency follows solely up to data oscillations and examples illustrate this can be problematic.

preprint2020arXiv

Convergence analysis of a numerical scheme for a tumour growth model

We consider a one--spatial dimensional tumour growth model [2, 3, 4] that consists of three dependent variables of space and time: volume fraction of tumour cells, velocity of tumour cells, and nutrient concentration. The model variables satisfy a coupled system of semilinear advection equation (hyperbolic), simplified linear Stokes equation (elliptic), and semilinear diffusion equation (parabolic) with appropriate conditions on the time-dependent boundary, which is governed by an ordinary differential equation. We employ a reformulation of the model defined in a larger, fixed time-space domain to overcome some theoretical difficulties related to the time-dependent boundary. This reformulation reduces the complexity of the model by removing the need to explicitly track the time-dependent boundary, but nonlinearities in the equations, noncoercive operators in the simplified Stokes equation, and interdependence between the unknown variables still challenge the proof of suitable a priori estimates. A numerical scheme that employs a finite volume method for the hyperbolic equation, a finite element method for the elliptic equation, and a backward Euler in time--mass lumped finite element in space method for the parabolic equation is developed. We establish the existence of a time interval $(0,T_{\ast})$ over which, using compactness techniques, we can extract a convergent subsequence of the numerical approximations. The limit of any such convergent subsequence is proved to be a weak solution of the continuous model in an appropriate sense, which we call a threshold solution. Numerical tests and justifications that confirm the theoretical findings conclude the paper.

preprint2020arXiv

Discontinuous Galerkin Finite Element Methods for the Landau-de Gennes Minimization Problem of Liquid Crystals

We consider a system of second order non-linear elliptic partial differential equations that models the equilibrium configurations of a two dimensional planar bistable nematic liquid crystal device. Discontinuous Galerkin finite element methods are used to approximate the solutions of this nonlinear problem with non-homogeneous Dirichlet boundary conditions. A discrete inf-sup condition demonstrates the stability of the discontinuous Galerkin discretization of a well-posed linear problem. We then establish the existence and local uniqueness of the discrete solution of the non-linear problem. An a priori error estimate in the energy norm is derived and a best approximation property is demonstrated. Further, we prove the quadratic convergence of Newton's iterates along with complementary numerical experiments.

preprint2020arXiv

Hessian discretisation method for fourth order semi-linear elliptic equations: applications to the von Kármán and Navier--Stokes models

This paper deals with the Hessian discretisation method (HDM) for fourth order semi-linear elliptic equations with a trilinear nonlinearity. The HDM provides a generic framework for the convergence analysis of several numerical methods, such as, the conforming and non-conforming finite element methods (ncFEMs) and methods based on gradient recovery (GR) operators. The Adini ncFEM and GR method, a specific scheme that is based on cheap, local reconstructions of higher-order derivatives from piecewise linear functions, are analysed for the first time for fourth order semi-linear elliptic equations with trilinear nonlinearity. Four properties namely, the coercivity, consistency, limit-conformity and compactness enable the convergence analysis in HDM framework that does not require any regularity of the exact solution. Two important problems in applications namely, the Navier--Stokes equations in stream function vorticity formulation and the von Kármán equations of plate bending are discussed. Results of numerical experiments are presented for the Morley ncFEM and GR method.

preprint2020arXiv

Morley Finite Element Method for the von Kármán Obstacle Problem

This paper focusses on the von Kármán equations for the moderately large deformation of a very thin plate with the convex obstacle constraint leading to a coupled system of semilinear fourth-order obstacle problem and motivates its nonconforming Morley finite element approximation. The first part establishes the well-posedness of the von Kármán obstacle problem and also discusses the uniqueness of the solution under an a priori and an a posteriori smallness condition on the data. The second part of the article discusses the regularity result of Frehse from 1971 and combines it with the regularity of the solution on a polygonal domain. The third part of the article shows an a priori error estimate for optimal convergence rates for the Morley finite element approximation to the von Kármán obstacle problem for small data. The article concludes with numerical results that illustrates the requirement of smallness assumption on the data for optimal convergence rate.