A posteriori error estimates of hp spectral element method for parabolic optimal control problems

<abstract><p>In this paper, we investigate the spectral element approximation for the optimal control problem of parabolic equation, and present a hp spectral element approximation scheme for the parabolic optimal control problem. For improve the accuracy of the algorithm and construct an adaptive finite element approximation. Under the Scott-Zhang type quasi-interpolation operator, a $ L^2(H^1)-L^2(L^2) $ posteriori error estimates of the hp spectral element approximated solutions for both the state variables and the control variable are obtained. Adopting two auxiliary equations and stability results, a $ L^2(L^2)-L^2(L^2) $ posteriori error estimates are derived for the hp spectral element approximation of optimal parabolic control problem.</p></abstract>

On Parameter Identification for Reaction-Dominated Pore-Scale Reactive Transport Using Modified Bee Colony Algorithm

Parameter identification is an important research topic with a variety of applications in industrial and environmental problems. Usually, a functional has to be minimized in conjunction with parameter identification; thus, there is a certain similarity between the parameter identification and optimization. A number of rigorous and efficient algorithms for optimization problems were developed in recent decades for the case of a convex functional. In the case of a non-convex functional, the metaheuristic algorithms dominate. This paper discusses an optimization method called modified bee colony algorithm (MBC), which is a modification of the standard bees algorithm (SBA). The SBA is inspired by a particular intelligent behavior of honeybee swarms. The algorithm is adapted for the parameter identification of reaction-dominated pore-scale transport when a non-convex functional has to be minimized. The algorithm is first checked by solving a few benchmark problems, namely finding the minima for Shekel, Rosenbrock, Himmelblau and Rastrigin functions. A statistical analysis was carried out to compare the performance of MBC with the SBA and the artificial bee colony (ABC) algorithm. Next, MBC is applied to identify the three parameters in the Langmuir isotherm, which is used to describe the considered reaction. Here, 2D periodic porous media were considered. The simulation results show that the MBC algorithm can be successfully used for identifying admissible sets for the reaction parameters in reaction-dominated transport characterized by low Pecklet and high Damkholer numbers. Finite element approximation in space and implicit time discretization are exploited to solve the direct problem.

Multi-Incidence Holographic Profilometry for Large Gradient Surfaces with Sub-Micron Focusing Accuracy

Surface reconstruction for micro-samples with large discontinuities using digital holography is a challenge. To overcome this problem, multi-incidence digital holographic profilometry (MIDHP) has been proposed. MIDHP relies on the numerical generation of the longitudinal scanning function (LSF) for reconstructing the topography of the sample with large depth and high axial resolution. Nevertheless, the method is unable to reconstruct surfaces with large gradients due to the need of: (i) high precision focusing that manual adjustment cannot fulfill and (ii) preserving the functionality of the LSF that requires capturing and processing many digital holograms. In this work, we propose a novel MIDHP method to solve these limitations. First, an autofocusing algorithm based on the comparison of shapes obtained by the LSF and the thin tilted element approximation is proposed. It is proven that this autofocusing algorithm is capable to deliver in-focus plane localization with submicron resolution. Second, we propose that wavefield summation for the generation of the LSF is carried out in Fourier space. It is shown that this scheme enables a significant reduction of arithmetic operations and can minimize the number of Fourier transforms needed. Hence, a fast generation of the LSF is possible without compromising its accuracy. The functionality of MIDHP for measuring surfaces with large gradients is supported by numerical and experimental results.

NUMERICAL SOLUTION OF THE PROBLEM FOR POISSON’S EQUATION WITH THE USE OF DAUBECHIES WAVELET DISCRETE-CONTINUAL FINITE ELEMENT METHOD

Numerical solution of the problem for Poisson’s equation with the use of Daubechies wavelet discrete continual finite element method (specific version of wavelet-based discrete-continual finite element method) is under consideration in the distinctive paper. The operational initial continual and discrete-continual formulations of the problem are given, several aspects of finite element approximation are considered. Some information about the numerical implementation and an example of analysis are presented.

Implicit-explicit multistep formulations for finite element discretisations using continuous interior penalty

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.

A mixed elasticity formulation for fluid-poroelastic structure interaction

We develop a mixed finite element method for the coupled problem arising in the interaction between a free fluid governed by the Stokes equations and flow in deformable porous medium modeled by the Biot system of poroelasticity. Mass conservation, balance of stress, and the Beavers-Joseph-Saffman condition are imposed on the interface. We consider a fully mixed Biot formulation based on a weakly symmetric stress-displacement-rotation elasticity system and Darcy velocity-pressure flow formulation. A velocity-pressure formulation is used for the Stokes equations. The interface conditions are incorporated through the introduction of the traces of the structure velocity and the Darcy pressure as Lagrange multipliers. Existence and uniqueness of a solution are established for the continuous weak formulation. Stability and error estimates are derived for the semi-discrete continuous-in-time mixed finite element approximation. Numerical experiments are presented to verify the theoretical results and illustrate the robustness of the method with respect to the physical parameters.

Nonlinearly Preconditioned FETI Solver for Substructured Formulations of Nonlinear Problems

We consider the finite element approximation of the solution to elliptic partial differential equations such as the ones encountered in (quasi)-static mechanics, in transient mechanics with implicit time integration, or in thermal diffusion. We propose a new nonlinear version of preconditioning, dedicated to nonlinear substructured and condensed formulations with dual approach, i.e., nonlinear analogues to the Finite Element Tearing and Interconnecting (FETI) solver. By increasing the importance of local nonlinear operations, this new technique reduces communications between processors throughout the parallel solving process. Moreover, the tangent systems produced at each step still have the exact shape of classically preconditioned linear FETI problems, which makes the tractability of the implementation barely modified. The efficiency of this new preconditioner is illustrated on two academic test cases, namely a water diffusion problem and a nonlinear thermal behavior.

Numerical method for solving the piecewise constant source inverse problem of an elliptic equation from a partial boundary observation data

We present results of numerical investigation of the source term recovery in a boundary value problem for an elliptic equation. An additional information about the solution is considered as its normal derivative taken on a part of the boundary. Such source inverse problem is related with inverse gravimetry problem of determining an inhomogeneity from gravitational potential anomalies on the Earth’s surface. We propose an iterative method for numerical recovery of the source term on the base of minimization of the observation residual by a gradient type method. The numerical implementation is based on finite element approximation using the FEniCS scientific computing platform and the dolfin-adjoint package. The capabilities of the developed computational algorithm are illustrated by results of numerical solutions of two dimensional test problems.