Related papers: $\varphi$enics: Vainshtein screening with the fini…
A new field of numerical astrophysics is introduced which addresses the solution of large, multidimensional structural or slowly-evolving problems (rotating stars, interacting binaries, thick advective accretion disks, four dimensional…
We develop a cut finite element method for the Bernoulli free boundary problem. The free boundary, represented by an approximate signed distance function on a fixed background mesh, is allowed to intersect elements in an arbitrary fashion.…
Many problems of theoretical and practical interest involve finding a convex or concave function. For instance, optimization problems such as finding the projection on the convex functions in $H^k(\Omega)$, or some problems in economics. In…
A generalized geometric method is developed for constructing exact solutions of gravitational field equations in Einstein theory and generalizations. First, we apply the formalism of nonholonomic frame deformations (formally considered for…
Discrete ordinate ($S_N$) and filtered spherical harmonics ($FP_N$) based schemes have been proven to be robust and accurate in solving the Boltzmann transport equation but they have their own strengths and weaknesses in different physical…
The chameleon is a theorised scalar field that couples to matter and possess a screening mechanism, which weakens observational constraints from experiments performed in regions of higher matter density. One consequence of this screening…
Chameleon scalar fields can screen their associated fifth forces from detection by changing their mass with the local density. These models are an archetypal example of a screening mechanism, and have become an important target for both…
We introduce \texttt{featom}, an open source code that implements a high-order finite element solver for the radial Schr\"odinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or…
Effective field theories featuring light scalar fields play a pivotal role in addressing fundamental questions in (astro)particle physics and cosmology. However, such theories often confront hierarchy problems in the absence of a symmetry.…
We consider linear systems arising from the use of the finite element method for solving scalar linear elliptic problems. Our main result is that these linear systems, which are symmetric and positive semidefinite, are well approximated by…
We present a methodology combining neural networks with physical principle constraints in the form of partial differential equations (PDEs). The approach allows to train neural networks while respecting the PDEs as a strong constraint in…
In this paper, we give a detailed account of the algorithm outlined in [1] for Feynman integral reduction and $\varepsilon$-factorised differential equations. The algorithm consists of two steps. In the first step, we use a new geometric…
Here I develop the simplest method in order to evaluate whether or not the Vainshtein mechanism can operate for a given set of parameters in a given solution. The method is based on the formulation of the mechanism in terms of the…
A high order cut finite element method is formulated for solving the elastic wave equation. Both a single domain problem and an interface problem are treated. The boundary or interface are allowed to cut through the background mesh. To…
In this work, we present an efficient numerical implementation of the finite element method for modal analysis that leverages various symmetry operations, including spatial symmetry in point groups and space-time symmetry in…
We introduce the concept of data-driven finite element methods. These are finite-element discretizations of partial differential equations (PDEs) that resolve quantities of interest with striking accuracy, regardless of the underlying mesh…
Finite element simulations have been used to solve various partial differential equations (PDEs) that model physical, chemical, and biological phenomena. The resulting discretized solutions to PDEs often do not satisfy requisite physical…
We present an upgraded version of \textsc{MG-MAMPOSSt}, an extension of the \textsc{MAMPOSSt} algorithm that performs Bayesian fits of models of mass and velocity anisotropy profiles to the distribution of tracers in projected phase space,…
A finite element methodology for large classes of variational boundary value problems is defined which involves discretizing two linear operators: (1) the differential operator defining the spatial boundary value problem; and (2) a Riesz…
Viable modifications of gravity on cosmological scales predominantly rely on screening mechanisms to recover Einstein's Theory of General Relativity in the Solar System, where it has been well tested. A parametrisation of the effects of…