Search

In this paper, we develop a framework for the discretization of a mixed formulation of quasi-reversibility solutions to ill-posed problems with respect to Poisson's equations. By carefully choosing test and trial spaces a formulation that is stable in a certain residual norm is obtained. Numerical stability and optimal convergence are established based on the conditional stability property of the problem. Tikhonov regularisation is necessary for high order polynomial approximation, , but its weak consistency may be tuned to allow for optimal convergence. For low order elements a simple numerical scheme with optimal convergence is obtained without stabilization. We also provide a guideline for feasible pairs of finite element spaces that satisfy suitable stability and consistency assumptions. Numerical experiments are provided to illustrate the theoretical results.

In this paper we continue the work on implicit-explicit (IMEX) time discretizations for the incompressible Oseen equations that we started in \cite{BGG23} (E. Burman, D. Garg, J. Guzm\`an, {\emph{Implicit-explicit time discretization for Oseen's equation at high Reynolds number with application to fractional step methods}}, SIAM J. Numer. Anal., 61, 2859--2886, 2023). The pressure velocity coupling and the viscous terms are treated implicitly, while the convection term is treated explicitly using extrapolation. Herein we focus on the implicit-explicit Crank-Nicolson method for time discretization. For the discretization in space we consider finite element methods with stabilization on the gradient jumps. The stabilizing terms ensures inf-sup stability for equal order interpolation and robustness at high Reynolds number. Under suitable Courant conditions we prove stability of the implicit-explicit Crank-Nicolson scheme in this regime. The stabilization allows us to prove error estimates of order $O(h^{k+\frac12} + \tau^2)$. Here $h$ is the mesh parameter, $k$ the polynomial order and $\tau$ the time step. Finally we discuss some fractional step methods that are implied by the IMEX scheme. Numerical examples are reported comparing the different methods when applied to the Navier-Stokes' equations., Comment: 33 pages

We consider a stable unique continuation problem for the wave equation where the initial data is lacking and the solution is reconstructed using measurements in some subset of the bulk domain. Typically fairly sophisticated space-time methods have been used in previous work to obtain stable and accurate solutions to this reconstruction problem. Here we propose to solve the problem using a standard discontinuous Galerkin method for the temporal discretization and continuous finite elements for the space discretization. Error estimates are established under a geometric control condition. We also investigate two preconditioning strategies which can be used to solve the arising globally coupled space-time system by means of simple time-stepping procedures. Our numerical experiments test the performance of these strategies and highlight the importance of the geometric control condition for reconstructing the solution beyond the data domain.

We consider a parabolic-parabolic interface problem and construct a loosely coupled prediction-correction scheme based on the Robin-Robin splitting method analyzed in [J. Numer. Math., 31(1):59--77, 2023]. We show that the errors of the correction step converge at $\mathcal O((\Delta t)^2)$, under suitable convergence rate assumptions on the discrete time derivative of the prediction step, where $\Delta t$ stands for the time-step length. Numerical results are shown to support our analysis and the assumptions.

We consider a loosely coupled, non-iterative Robin-Robin coupling method proposed and analyzed in [J. Numer. Math., 31(1):59--77, 2023] for a parabolic-parabolic interface problem and prove estimates for the discrete time derivatives of the scalar field in different norms. When the interface is flat and perpendicular to two of the edges of the domain we prove error estimates in the $H^2$-norm. Such estimates are key ingredients to analyze a defect correction method for the parabolic-parabolic interface problem. Numerical results are shown to support our findings.

In this paper, we consider the unique continuation problem for the Schr\"odinger equations. We prove a H\"older type conditional stability estimate and build up a parameterized stabilized finite element scheme adaptive to the \textit{a priori} knowledge of the solution, achieving error estimates in interior domains with convergence up to continuous stability. The approximability of the scheme to solutions with only $H^1$-regularity is studied and the convergence rate for solutions with regularity higher than $H^1$ is also shown. Comparisons in terms of different parameterization for different regularities will be illustrated with respect to the convergence and condition numbers of the linear systems. Finally, numerical experiments will be given to illustrate the theory.

This paper discusses the effect of rotation on the boundary layer in high Reynolds number flow over a ridge using a numerical method based on stabilised finite elements that captures steady solutions up to Reynolds number of order $10^6$. The results are validated against boundary layer computations in shallow flows and for deep flows against experimental observations reported in Machicoane et al. (Phys. Rev. Fluids, 2018). In all cases considered the boundary layer remains attached, even at large Reynolds numbers, provided the Rossby number of the flow is sufficiently small. At any fixed Rossby number the flow detaches at sufficiently high Reynolds number to form a steady recirculating region in the lee of the ridge. At even higher Reynolds numbers no steady flow is found. This disappearance of steady solutions closely reproduces the transition to unsteadiness seen in the laboratory.

We consider finite element approximations of unique continuation problems subject to elliptic equations in the case where the normal derivative of the exact solution is known to reside in some finite dimensional space. To give quantitative error estimates we prove Lipschitz stability of the unique continuation problem in the global H1-norm. This stability is then leveraged to derive optimal a posteriori and a priori error estimates for a primal-dual stabilised finite method.

Introducing a coupling framework reminiscent of FETI methods, but here on abstract form, we establish conditions for stability and minimal requirements for well-posedness on the continuous level, as well as conditions on local solvers for the approximation of subproblems. We then discuss stability of the resulting Lagrange multiplier methods and show stability under a mesh conditions between the local discretizations and the mortar space. If this condition is not satisfied we show how a stabilization, acting only on the multiplier can be used to achieve stability. The design of preconditioners of the Schur complement system is discussed in the unstabilized case. Finally we discuss some applications that enter the framework.

We consider finite element approximations of ill-posed elliptic problems with conditional stability. The notion of {\emph{optimal error estimates}} is defined including both convergence with respect to mesh parameter and perturbations in data. The rate of convergence is determined by the conditional stability of the underlying continuous problem and the polynomial order of the finite element approximation space. A proof is given that no finite element approximation can converge at a better rate than that given by the definition, justifying the concept. A recently introduced class of finite element methods with weakly consistent regularisation is recalled and the associated error estimates are shown to be quasi optimal in the sense of our definition.

We study unique continuation over an interface using a stabilized unfitted finite element method tailored to the conditional stability of the problem. The interface is approximated using an isoparametric transformation of the background mesh and the corresponding geometrical error is included in our error analysis. To counter possible destabilizing effects caused by non-conformity of the discretization and cope with the interface conditions, we introduce adapted regularization terms. This allows to derive error estimates based on conditional stability. The necessity and effectiveness of the regularization is illustrated in numerical experiments. We also explore numerically the effect of the heterogeneity in the coefficients on the ability to reconstruct the solution outside the data domain. For Helmholtz equations we find that a jump in the flux impacts the stability of the problem significantly less than the size of the wavenumber., Comment: 39 pages, 9 figures

In this paper, we provide a theoretical analysis of the recently introduced weakly adversarial networks (WAN) method, used to approximate partial differential equations in high dimensions. We address the existence and stability of the solution, as well as approximation bounds. More precisely, we prove the existence of discrete solutions, intended in a suitable weak sense, for which we prove a quasi-best approximation estimate similar to Cea's lemma, a result commonly found in finite element methods. We also propose two new stabilized WAN-based formulas that avoid the need for direct normalization. Furthermore, we analyze the method's effectiveness for the Dirichlet boundary problem that employs the implicit representation of the geometry. The key requirement for achieving the best approximation outcome is to ensure that the space for the test network satisfies a specific condition, known as the inf-sup condition, essentially requiring that the test network set is sufficiently large when compared to the trial space. The method's accuracy, however, is only determined by the space of the trial network. We also devise a pseudo-time XNODE neural network class for static PDE problems, yielding significantly faster convergence results than the classical DNN network., Comment: 29 pages, 8 figures

We consider a unique continuation problem where the Dirichlet trace of the solution is known to have finite dimension. We prove Lipschitz stability of the unique continuation problem and design a finite element method that exploits the finite dimensionality to enhance stability. Optimal a priori and a posteriori error estimates are shown for the method. The extension to problems where the trace is not in a finite dimensional space, but can be approximated to high accuracy using finite dimensional functions is discussed. Finally, the theory is illustrated in some numerical examples.

The Poisson--Boltzmann equation is widely used to model electrostatics in molecular systems. Available software packages solve it using finite difference, finite element, and boundary element methods, where the latter is attractive due to the accurate representation of the molecular surface and partial charges, and exact enforcement of the boundary conditions at infinity. However, the boundary element method is limited to linear equations and piecewise constant variations of the material properties. In this work, we present a scheme that couples finite and boundary elements for the Poisson--Boltzmann equation, where the finite element method is applied in a confined {\it solute} region, and the boundary element method in the external {\it solvent} region. As a proof-of-concept exercise, we use the simplest methods available: Johnson--N\'ed\'elec coupling with mass matrix and diagonal preconditioning, implemented using the Bempp-cl and FEniCSx libraries via their Python interfaces. We showcase our implementation by computing the polar component of the solvation free energy of a set of molecules using a constant and a Gaussian-varying permittivity. We validate our implementation against the finite difference code APBS (to 0.5\%), and show scaling from protein G B1 (955 atoms) up to immunoglobulin G (20\,148 atoms). For small problems, the coupled method was efficient, outperforming a purely boundary integral approach. For Gaussian-varying permittivities, which are beyond the applicability of boundary elements alone, we were able to run medium to large sized problems on a single workstation. Development of better preconditioning techniques and the use of distributed memory parallelism for larger systems remains an area for future work. We hope this work will serve as inspiration for future developments for molecular electrostatics with implicit solvent models.

We revisit the celebrated Kohn-Vogelius penalty method and discuss how to use it for the unique continuation problem where data is given in the bulk of the domain. We then show that the primal-dual mixed finite element methods for the elliptic Cauchy problem introduced in \cite{BLO18} (\emph{E. Burman, M. Larson, L. Oksanen, Primal-dual mixed finite element methods for the elliptic Cauchy problem, SIAM J. Num. Anal., 56 (6), 2018}) can be interpreted as a Kohn-Vogelius penalty method and modify it to allow for unique continuation using data in the bulk. We prove that the resulting linear system is invertible for all data. Then we show that by introducing a singularly perturbed Robin condition on the discrete level sufficient regularization is obtained so that error estimates can be shown using conditional stability. Finally we show how the method can be used for the identification of the diffusivity coefficient in a second order elliptic operator with partial data. Some numerical examples are presented showing the performance of the method for unique continuation and for impedance computed tomography with partial data., Comment: Submitted to the conference proceedings of the Peruvian conference on Scientific Computing, October 2023

In this article, we design and analyze an arbitrary-order stabilized finite element method to approximate the unique continuation problem for laminar steady flow described by the linearized incompressible Navier--Stokes equation. We derive quantitative local error estimates for the velocity, which account for noise level and polynomial degree, using the stability of the continuous problem in the form of a conditional stability estimate. Numerical examples illustrate the performances of the method with respect to the polynomial order and perturbations in the data. We observe that the higher order polynomials may be efficient for ill-posed problems, but are also more sensitive for problems with poor stability due to the ill-conditioning of the system., Comment: 26 PAGES

We introduce an arbitrary order, stabilized finite element method for solving a unique continuation problem subject to the time-harmonic elastic wave equation with variable coefficients. Based on conditional stability estimates we prove convergence rates for the proposed method which take into account the noise level and the polynomial degree. A series of numerical experiments corroborates our theoretical results and explores additional aspects, e.g. how the quality of the reconstruction depends on the geometry of the involved domains. We find that certain convexity properties are crucial to obtain a good recovery of the wave displacement outside the data domain and that higher polynomial orders can be more efficient but also more sensitive to the ill-conditioned nature of the problem., Comment: 36 pages, 10 figures

In this note we design a cut finite element method for a low order divergence free element applied to a boundary value problem subject to Stokes' equations. For the imposition of Dirichlet boundary conditions we consider either Nitsche's method or a stabilized Lagrange multiplier method. In both cases the normal component of the velocity is constrained using a multiplier, different from the standard pressure approximation. The divergence of the approximate velocities is pointwise zero over the whole mesh domain, and we derive optimal error estimates for the velocity and pressures, where the error constant is independent of how the physical domain intersects the computational mesh, and of the regularity of the pressure multiplier imposing the divergence free condition.

In this paper we will review recent advances in the application of the augmented Lagrange multiplier method as a general approach for generating multiplier--free stabilised methods. We first show how the method generates Galerkin/Least Squares type schemes for equality constraints and then how it can be extended to develop new stabilised methods for inequality constraints. Application to several different problems in computational mechanics is given.

We develop a discrete extension operator for trimmed spline spaces consisting of piecewise polynomial functions of degree $p$ with $k$ continuous derivatives. The construction is based on polynomial extension from neighboring elements together with projection back into the spline space. We prove stability and approximation results for the extension operator. Finally, we illustrate how we can use the extension operator to construct a stable cut isogeometric method for an elliptic model problem.

In this paper we propose, analyze, and test numerically a pressure-robust stabilized finite element for a linearized problem in incompressible fluid mechanics, namely, the steady Oseen equation with low viscosity. Stabilization terms are defined by jumps of different combinations of derivatives for the convective term over the element faces of the triangulation of the domain. With the help of these stabilizing terms, and the fact the finite element space is assumed to provide a point-wise divergence-free velocity, an $\mathcal O\big(h^{k+\frac12}\big)$ error estimate in the $L^2$-norm is proved for the method (in the convection-dominated regime), and optimal order estimates in the remaining norms of the error. Numerical results supporting the theoretical findings are provided., Comment: 3 figures, 23 pages

In this note, we develop a new stabilization mechanism for cut finite element methods that generalizes previous approaches of ghost penalty type in two ways: (1) The quantity that is stabilized and (2) The choice of elements that are connected in the stabilization. In particular, we can stabilize functionals of the discrete function such as finite element degrees of freedom. We subsequently show that the kernel of our ghost penalty operator defines a finite element space based on discrete extensions in the spirit of those introduced in Burman, E.; Hansbo, P. and Larson, M. G., CutFEM Based on Extended Finite Element Spaces, arXiv2101.10052, 2021., Comment: 31 pages, 12 figures

We present a loosely coupled, non-iterative time-splitting scheme based on Robin-Robin coupling conditions. We apply a novel unified analysis for this scheme applied to both a Parabolic/Parabolic coupled system and a Parabolic/Hyperbolic coupled system. We show for both systems that the scheme is stable, and the error converges as $\mathcal{O}\big(\Delta t \sqrt{T +\log{\frac{1}{\Delta t}}}\big)$, where $\Delta t$ is the time step

We consider the null controllability problem for the wave equation, and analyse a stabilized finite element method formulated on a global, unstructured spacetime mesh. We prove error estimates for the approximate control given by the computational method. The proofs are based on the regularity properties of the control given by the Hilbert Uniqueness Method, together with the stability properties of the numerical scheme. Numerical experiments illustrate the results.

In this paper we discuss a hybridised method for FEM-BEM coupling. The coupling from both sides use a Nitsche type approach to couple to the trace variable. This leads to a formulation that is robust and flexible with respect to approximation spaces and can easily be combined as a building block with other hybridised methods. Energy error norm estimates and the convergence of Jacobi iterations are proved and the performance of the method is illustrated on some computational examples.

In this paper we consider the semi-discretization in space of a first order scalar transport equation. For the space discretization we use standard continuous finite elements. To obtain stability we add a penalty on the jump of the gradient over element faces. We recall some global error estimates for smooth and rough solutions and then prove a new local error estimate for the transient linear transport equation. In particular we show that in the stabilized method the effect of non-smooth features in the solution decay exponentially from the space time zone where the solution is rough so that smooth features will be transported unperturbed. Locally the $L^2$-norm error converges with the expected order $O(h^{k+\frac12})$. We then illustrate the results numerically. In particular we show the good local accuracy in the smooth zone of the stabilized method and that the standard Galerkin fails to approximate a solution that is smooth at the final time if discontinuities have been present in the solution at some time during the evolution.

We propose two different discrete formulations for the weak imposition of the Neumann boundary conditions of the Darcy flow. The Raviart-Thomas mixed finite element on both triangular and quadrilateral meshes is considered for both methods. One is a consistent discretization depending on a weighting parameter scaling as $\mathcal O(h^{-1})$, while the other is a penalty-type formulation obtained as the discretization of a perturbation of the original problem and relies on a parameter scaling as $\mathcal O(h^{-k-1})$, $k$ being the order of the Raviart-Thomas space. We rigorously prove that both methods are stable and result in optimal convergent numerical schemes with respect to appropriate mesh-dependent norms, although the chosen norms do not scale as the usual $L^2$-norm. However, we are still able to recover the optimal a priori $L^2$-error estimates for the velocity field, respectively, for high-order and the lowest-order Raviart-Thomas discretizations, for the first and second numerical schemes. Finally, some numerical examples validating the theory are exhibited.

We present a new approach for the mechanically consistent modelling and simulation of fluid-structure interactions with contact. The fundamental idea consists of combining a relaxed contact formulation with the modelling of seepage through a porous layer of co-dimension 1 during contact. For the latter, a Darcy model is considered in a thin porous layer attached to a solid boundary in the limit of infinitesimal thickness. In combination with a relaxation of the contact conditions the computational model is both mechanically consistent and simple to implement. We analyse the approach in detailed numerical studies with both thick- and thin-walled solids, within a fully Eulerian and an immersed approach for the fluid-structure interaction and using fitted and unfitted finite element discretisations.

In the present work we show some results on the effect of the Smagorinsky model on the stability of the associated perturbation equation. We show that in the presence of a spectral gap, such that the flow can be decomposed in a large scale with moderate gradient and a small amplitude fine scale with arbitratry gradient, the Smagorinsky model admits stability estimates for perturbations, with exponential growth depending only on the large scale gradient. We then show in the context of stabilized finite element methods that the same result carries over to the approximation and that in this context, for suitably chosen finite element spaces the Smagorinsky model acts as a stabilizer yielding close to optimal error estimates in the $L^2$-norm for smooth flows in the pre-asymptotic high Reynolds number regime.

We develop a general framework for construction and analysis of discrete extension operators with application to unfitted finite element approximation of partial differential equations. In unfitted methods so called cut elements intersected by the boundary occur and these elements must in general by stabilized in some way. Discrete extension operators provides such a stabilization by modification of the finite element space close to the boundary. More precisely, the finite element space is extended from the stable interior elements over the boundary in a stable way which also guarantees optimal approximation properties. Our framework is applicable to all standard nodal based finite elements of various order and regularity. We develop an abstract theory for elliptic problems and associated parabolic time dependent partial differential equations and derive a priori error estimates. We finally apply this to some examples of partial differential equations of different order including the interface problems, the biharmonic operator and the sixth order triharmonic operator., Comment: 32 pages, 9 figures

We consider a finite element method with symmetric stabilisation for the discretisation of the transient convection--diffusion equation. For the time-discretisation we consider either the second order backwards differentiation formula or the Crank-Nicolson method. Both the convection term and the associated stabilisation are treated explicitly using an extrapolated approximate solution. We prove stability of the method and the $\tau^2 + h^{p+{\frac12}}$ error estimates for the $L^2$-norm under either the standard hyperbolic CFL condition, when piecewise affine ($p=1$) approximation is used, or in the case of finite element approximation of order $p \ge 1$, a stronger, so-called $4/3$-CFL, i.e. $\tau \leq C h^{4/3}$. The theory is illustrated with some numerical examples.

In this note we develop a fully explicit cut finite element method for the wave equation. The method is based on using a standard leap frog scheme combined with an extension operator that defines the nodal values outside of the domain in terms of the nodal values inside the domain. We show that the mass matrix associated with the extended finite element space can be lumped leading to a fully explicit scheme. We derive stability estimates for the method and provide optimal order a priori error estimates. Finally, we present some illustrating numerical examples., Comment: 29 pages, 15figures

In this paper we discuss a level set approach for the identification of an unknown boundary in a computational domain. The problem takes the form of a Bernoulli problem where only the Dirichlet datum is known on the boundary that is to be identified, but additional information on the Neumann condition is available on the known part of the boundary. The approach uses a classical constrained optimization problem, where a cost functional is minimized with respect to the unknown boundary, the position of which is defined implicitly by a level set function. To solve the optimization problem a steepest descent algorithm using shape derivatives is applied. In each iteration the cut finite element method is used to obtain high accuracy approximations of the pde-model constraint for a given level set configuration without re-meshing. We consider three different shape derivatives. First the classical one, derived using the continuous optimization problem (optimize then discretize). Then the functional is first discretized using the CutFEM method and the shape derivative is evaluated on the finite element functional (discretize then optimize). Finally we consider a third approach, also using a discretized functional. In this case we do not perturb the domain, but consider a so-called boundary value correction method, where a small correction to the boundary position may be included in the weak boundary condition. Using this correction the shape derivative may be obtained by perturbing a distance parameter in the discrete variational formulation. The theoretical discussion is illustrated with a series of numerical examples showing that all three approaches produce similar result on the proposed Bernoulli problem., Comment: 28 pages; 7 figures

We derive a residual based a-posteriori error estimate for the outer normal flux of approximations to {the diffusion problem with variable coefficient}. By analyzing the solution of the adjoint problem, we show that error indicators in the bulk may be defined to be of higher order than those close to the boundary, which lead to more economic meshes. The theory is illustrated with some numerical examples., Comment: 27 pages, 13 figures, 3 tables

Discretization of Navier-Stokes' equations using pressure-robust finite element methods is considered for the high Reynolds number regime. To counter oscillations due to dominating convection we add a stabilization based on a bulk term in the form of a residual-based least squares stabilization of the vorticity equation supplemented by a penalty term on (certain components of) the gradient jump over the elements faces. Since the stabilization is based on the vorticity equation, it is independent of the pressure gradients, which makes it pressure-robust. Thus, we prove pressure-independent error estimates in the linearized case, known as Oseen's problem. In fact, we prove an $O(h^{k+\frac12})$ error estimate in the $L^2$-norm that is known to be the best that can be expected for this type of problem. Numerical examples are provided that, in addition to confirming the theoretical results, show that the present method compares favorably to the classical residual-based SUPG stabilization.

We consider a fully discrete loosely coupled scheme for incompressible fluid-structure interaction based on the time semi-discrete splitting method introduced in {\emph{[Burman, Durst \& Guzm\'an, arXiv:1911.06760]}}. The splittling method uses a Robin-Robin type coupling that allows for a segregated solution of the solid and the fluid systems, without inner iterations. For the discretisation in space we consider piecewise affine continuous finite elements for all the fields and ensure the inf-sup condition by using a Brezzi-Pitk\"aranta type pressure stabilization. The interfacial fluid-stresses are evaluated in a variationally consistent fashion, that is shown to admit an equivalent Lagrange multiplier formulation. We prove that the method is unconditionally stable and robust with respect to the amount of added-mass in the system. Furthermore, we provide an error estimate that shows the error in the natural energy norm for the system is $\mathcal O\big(\sqrt{T}(\sqrt{\Delta t} + h)\big)$ where $T$ is the final time, $\Delta t$ the time-step length and $h$ the space discretization parameter.

We show error estimates for a cut finite element approximation of a second order elliptic problem with mixed boundary conditions. The error estimates are of low regularity type where we consider the case when the exact solution $u \in H^s$ with $s\in (1,3/2]$. For Nitsche type methods this case requires special handling of the terms involving the normal flux of the exact solution at the the boundary. For Dirichlet boundary conditions the estimates are optimal, whereas in the case of mixed Dirichlet-Neumann boundary conditions they are suboptimal by a logarithmic factor., Comment: 21 pages, 4 figures

We consider the numerical approximation of the ill-posed data assimilation problem for stationary convection-diffusion equations and extend our previous analysis in [Numer. Math. 144, 451--477, 2020] to the convection-dominated regime. Slightly adjusting the stabilized finite element method proposed for dominant diffusion, we draw upon a local error analysis to obtain quasi-optimal convergence along the characteristics of the convective field through the data set. The weight function multiplying the discrete solution is taken to be Lipschitz and a corresponding super approximation result (discrete commutator property) is proven. The effect of data perturbations is included in the analysis and we conclude the paper with some numerical experiments., Comment: 25 pages, 16 figures; final published version

We develop a robust cut finite element method for a model of diffusion in fractured media consisting of a bulk domain with embedded cracks. The crack has its own pressure field and can cut through the bulk mesh in a very general fashion. Starting from a common background bulk mesh, that covers the domain, finite element spaces are constructed for the interface and bulk subdomains leading to efficient computations of the coupling terms. The crack pressure field also uses the bulk mesh for its representation. The interface conditions are a generalized form of conditions of Robin type previously considered in the literature which allows the modeling of a range of flow regimes across the fracture. The method is robust in the following way: 1. Stability of the formulation in the full range of parameter choices; and 2. Not sensitive to the location of the interface in the background mesh. We derive an optimal order a priori error estimate and present illustrating numerical examples.

We consider boundary element methods where the Calder\'on projector is used for the system matrix and boundary conditions are weakly imposed using a particular variational boundary operator designed using techniques from augmented Lagrangian methods. Regardless of the boundary conditions, both the primal trace variable and the flux are approximated. We focus on the imposition of Dirichlet conditions on the Helmholtz equation, and extend the analysis of the Laplace problem from \emph{Boundary element methods with weakly imposed boundary conditions} to this case. The theory is illustrated by a series of numerical examples., Comment: 24 pages

In this note we will explore some applications of the recently constructed piecewise affine, $H^1$-conforming element that fits in a discrete de Rham complex [Christiansen and Hu, Generalized finite element systems for smooth differential forms and Stokes' problem. Numer. Math. 140 (2018)]. In particular we show how the element leads to locking free methods for incompressible elasticity and viscosity robust methods for the Brinkman model.

In this note we introduce a mixed dimensional Stokes-Darcy coupling where a $d$ dimensional Stokes' flow is coupled to a Darcy model on the $d-1$ dimensional boundary of the domain. The porous layer introduces tangential creeping flow along the boundary and allows for the modelling of boundary flow due to surface roughness. This leads to a new model of flow in fracture networks with reservoirs in an impenetrable bulk matrix. Exploiting this modelling capability, we then formulate a fluid-structure interaction method with contact, where the porous layer allows for mechanically consistent contact and release. Physical seepage in the contact zone due to rough surfaces is modelled by the porous layer. Some numerical examples are reported, both on the Stokes'-Darcy coupling alone and on the fluid-structure interaction with contact in the porous boundary layer.

We consider $C^1$-continuous approximations of the Kirchhoff plate problem in combination with a mesh dependent augmented Lagrangian method on a simply supported Signorini boundary.

We consider a stabilized finite element method based on a spacetime formulation, where the equations are solved on a global (unstructured) spacetime mesh. A unique continuation problem for the wave equation is considered, where data is known in an interior subset of spacetime. For this problem, we consider a primal-dual discrete formulation of the continuum problem with the addition of stabilization terms that are designed with the goal of minimizing the numerical errors. We prove error estimates using the stability properties of the numerical scheme and a continuum observability estimate, based on the sharp geometric control condition by Bardos, Lebeau and Rauch. The order of convergence for our numerical scheme is optimal with respect to stability properties of the continuum problem and the interpolation errors of approximating with polynomial spaces. Numerical examples are provided that illustrate the methodology.