Manuel Solano

Tonatiuh Sánchez-Vizuet, Manuel Solano

In axisymmetric fusion reactors, the equilibrium magnetic configuration can be expressed in terms of the solution to a semi-linear elliptic equation known as the Grad-Shafranov equation, the solution of which determines the poloidal component of the magnetic field. When the geometry of the confinement region is known, the problem becomes an interior Dirichlet boundary value problem. We propose a high order solver based on the Hybridizable Discontinuous Galerkin method. The resulting algorithm (1) provides high order of convergence for the flux function and its gradient, (2) incorporates a novel method for handling piecewise smooth geometries by extension from polygonal meshes, (3) can handle geometries with non-smooth boundaries and x-points, and (4) deals with the semi-linearity through an accelerated two-grid fixed-point iteration. The effectiveness of the algorithm is verified with computations for cases where analytic solutions are known on configurations similar to those of actual devices (ITER with single null and double null divertor, NSTX, ASDEX upgrade, and Field Reversed Configurations).

Weifeng Qiu, Manuel Solano, Patrick Vega

We generalize the technique of [Solving Dirichlet boundary-value problems on curved domains by extensions from subdomains, SIAM J. Sci. Comput. 34, pp. A497--A519 (2012)] to elliptic problems with mixed boundary conditions and elliptic interface problems involving a non-polygonal interface. We study first the treatment of the Neumann boundary data since it is crucial to understand the applicability of the technique to curved interfaces. We provide numerical results showing that, in order to obtain optimal high order convergence, it is desirable to construct the computational domain by interpolating the boundary/interface using piecewise linear segments. In this case the distance of the computational domain to the exact boundary is only $O(h^2)$.

Esteban Henríquez, Manuel Solano

We analyze a high order unfitted hybridizable discontinuous Galerkin (HDG) method for an optimal control problem governed by a convection-diffusion equation posed in a domain with piecewise-wise $\mathcal{C}^2$ boundary $\partial Ω$. The computational domain $Ω_h$ does not necessarily fit $Ω$ and the Transfer Path Method (TPM) is used to transfer the boundary data from $\partial Ω$ to $\partial Ω_h$ through segments of direction $\boldsymbol{m}$. Under closeness conditions between $\partial Ω_h$ and $\partial Ω$ and on the transfer vector $\boldsymbol{m}$, we prove optimal order of convergence in the $L^2$-norm for all variables of the state and adjoint problems. We also show numerical examples to complement the theory.