Source author record

Artur Palha

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

9 published item(s)

preprint2021arXiv

A mass-, kinetic energy- and helicity-conserving mimetic dual-field discretization for three-dimensional incompressible Navier-Stokes equations, part I: Periodic domains

We introduce a mimetic dual-field discretization which conserves mass, kinetic energy and helicity for three-dimensional incompressible Navier-Stokes equations. The discretization makes use of a conservative dual-field mixed weak formulation where two evolution equations of velocity are employed and dual representations of the solution are sought for each variable. A temporal discretization, which staggers the evolution equations and handles the nonlinearity such that the resulting discrete algebraic systems are linear and decoupled, is constructed. The spatial discretization is mimetic in the sense that the finite dimensional function spaces form a discrete de Rham complex. Conservation of mass, kinetic energy and helicity in the absence of dissipative terms is proven at the discrete level. Proper dissipation rates of kinetic energy and helicity in the viscous case is also proven. Numerical tests supporting the method are provided.

preprint2016arXiv

A mass, energy, enstrophy and vorticity conserving (MEEVC) mimetic spectral element discretization for the 2D incompressible Navier-Stokes equations

In this work we present a mimetic spectral element discretization for the 2D incompressible Navier-Stokes equations that in the limit of vanishing dissipation exactly preserves mass, kinetic energy, enstrophy and total vorticity on unstructured grids. The essential ingredients to achieve this are: (i) a velocity-vorticity formulation in rotational form, (ii) a sequence of function spaces capable of exactly satisfying the divergence free nature of the velocity field, and (iii) a conserving time integrator. Proofs for the exact discrete conservation properties are presented together with numerical test cases on highly irregular grids.

preprint2016arXiv

A mimetic spectral element solver for the Grad-Shafranov equation

In this work we present a robust and accurate arbitrary order solver for the fixed-boundary plasma equilibria in toroidally axisymmetric geometries. To achieve this we apply the mimetic spectral element formulation presented in [56] to the solution of the Grad-Shafranov equation. This approach combines a finite volume discretization with the mixed finite element method. In this way the discrete differential operators ($\nabla$, $\nabla\times$, $\nabla\cdot$) can be represented exactly and metric and all approximation errors are present in the constitutive relations. The result of this formulation is an arbitrary order method even on highly curved meshes. Additionally, the integral of the \reviewerone{toroidal current $J_ϕ$} is exactly equal to the boundary integral of the poloidal field over the plasma boundary. This property can play an important role in the coupling between equilibrium and transport solvers. The proposed solver is tested on a varied set of plasma cross-sections (smooth and with an X-point) and also for a wide range of pressure and \reviewertwo{toroidal magnetic flux profiles}. Equilibria accurate up to machine precision are obtained. Optimal algebraic convergence rates of order $(p+1)$ and geometric convergence rates are shown for Soloviev solutions (including high Shafranov shifts), field-reversed configuration (FRC) solutions and spheromak analytical solutions. The robustness of the method is demonstrated for non-linear test cases, in particular on an equilibrium solution with a pressure pedestal.

preprint2015arXiv

A hybrid Eulerian-Lagrangian flow solver

Currently, Eulerian flow solvers are very efficient in accurately resolving flow structures near solid boundaries. On the other hand, they tend to be diffusive and to dampen high-intensity vortical structures after a short distance away from solid boundaries. The use of high order methods and fine grids, although alleviating this problem, gives rise to large systems of equations that are expensive to solve. Lagrangian solvers, as the regularized vortex particle method, have shown to eliminate (in practice) the diffusion in the wake. As a drawback, the modelling of solid boundaries is less accurate, more complex and costly than with Eulerian solvers (due to the isotropy of its computational elements). Given the drawbacks and advantages of both Eulerian and Lagrangian solvers the combination of both methods, giving rise to a hybrid solver, is advantageous. The main idea behind the hybrid solver presented is the following. In a region close to solid boundaries the flow is solved with an Eulerian solver, where the full Navier-Stokes equations are solved (possibly with an arbitrary turbulence model or DNS, the limitations being the computational power and the physical properties of the flow), outside of that region the flow is solved with a vortex particle method. In this work we present this hybrid scheme and verify it numerically on known 2D benchmark cases: dipole flow, flow around a cylinder and flow around a stalled airfoil. The success in modelling these flow conditions presents this hybrid approach as a promising alternative, bridging the gap between highly resolved and computationally intensive Eulerian CFD simulations and fast but less resolved Lagrangian simulations.

preprint2015arXiv

Mimetic spectral element method for Hamiltonian systems

There is a growing interest in the conservation of invariants when numerically solving a system of ordinary differential equations. Methods that exactly preserve these quantities in time are known as geometric integrators. In this paper we apply the recently developed mimetic framework to the solution of a system of first order ordinary differential equations. Depending on the discrete Hodge-* employed, two classes of arbitrary order time integrators are derived. It is shown that the one based on a canonical Hodge-* results in a symplectic integrator, whereas the one based on a Galerkin Hodge-* results in an energy preserving integrator. A set of numerical tests confirms these theoretical results.

preprint2013arXiv

Mimetic Spectral Element advection

We present a discretization of the linear advection of differential forms on bounded domains. The framework previously established is extended to incorporate the Lie derivative, $\mathcal L$, by means of Cartan's homotopy formula. The method is based on a physics-compatible discretization with spectral accuracy . It will be shown that the derived scheme has spectral convergence with local mass conservation. Artificial dispersion depends on the order of time integration.

preprint2013arXiv

Mixed Mimetic Spectral Element method applied to Darcy's problem

We present a discretization for Darcy's problem using the recently developed Mimetic Spectral Element Method. The gist lies in the exact discrete representation of integral relations. In this paper, an anisotropic flow through a porous medium is considered and a discretization of a full permeability tensor is presented. The performance of the method is evaluated on standard test problems, converging at the same rate as the best possible approximation.

preprint2013arXiv

Physics-compatible discretization techniques on single and dual grids, with application to the Poisson equation of volume forms

This paper introduces the basic concepts for physics-compatible discretization techniques. The paper gives a clear distinction between vectors and forms. Based on the difference between forms and pseudo-forms and the $\star$-operator which switches between the two, a dual grid description and a single grid description are presented. The dual grid method resembles a staggered finite volume method, whereas the single grid approach shows a strong resemblance with a finite element method. Both approaches are compared for the Poisson equation for volume forms.

preprint2011arXiv

Mimetic framework on curvilinear quadrilaterals of arbitrary order

In this paper higher order mimetic discretizations are introduced which are firmly rooted in the geometry in which the variables are defined. The paper shows how basic constructs in differential geometry have a discrete counterpart in algebraic topology. Generic maps which switch between the continuous differential forms and discrete cochains will be discussed and finally a realization of these ideas in terms of mimetic spectral elements is presented, based on projections for which operations at the finite dimensional level commute with operations at the continuous level. The two types of orientation (inner- and outer-orientation) will be introduced at the continuous level, the discrete level and the preservation of orientation will be demonstrated for the new mimetic operators. The one-to-one correspondence between the continuous formulation and the discrete algebraic topological setting, provides a characterization of the oriented discrete boundary of the domain. The Hodge decomposition at the continuous, discrete and finite dimensional level will be presented. It appears to be a main ingredient of the structure in this framework.