Highly accurate numerical computation of implicitly defined volumes using the Laplace-Beltrami operator

This paper introduces a novel method for the efficient and accurate computation of the volume of a domain whose boundary is given by an orientable hypersurface which is implicitly given as the iso-contour of a sufficiently smooth level-set function. After spatial discretization, local approximation of the hypersurface and application of the Gaussian divergence theorem, the volume integrals are transformed to surface integrals. Application of the surface divergence theorem allows for a further reduction to line integrals which are advantageous for numerical quadrature. We discuss the theoretical foundations and provide details of the numerical algorithm. Finally, we present numerical results for convex and non-convex hypersurfaces embedded in cuboidal domains, showing both high accuracy and thrid- to fourth-order convergence in space.Comment: 25 pages, 17 figures, 3 table

On the Numerical Evaluation of One-Loop Amplitudes: the Gluonic Case

We develop an algorithm of polynomial complexity for evaluating one-loop amplitudes with an arbitrary number of external particles. The algorithm is implemented in the Rocket program. Starting from particle vertices given by Feynman rules, tree amplitudes are constructed using recursive relations. The tree amplitudes are then used to build one-loop amplitudes using an integer dimension on-shell cut method. As a first application we considered only three and four gluon vertices calculating the pure gluonic one-loop amplitudes for arbitrary external helicity or polarization states. We compare our numerical results to analytical results in the literature, analyze the time behavior of the algorithm and the accuracy of the results, and give explicit results for fixed phase space points for up to twenty external gluons.Comment: 22 pages, 9 figures; v2: references added, version accepted for publicatio

A moving control volume approach to computing hydrodynamic forces and torques on immersed bodies

We present a moving control volume (CV) approach to computing hydrodynamic forces and torques on complex geometries. The method requires surface and volumetric integrals over a simple and regular Cartesian box that moves with an arbitrary velocity to enclose the body at all times. The moving box is aligned with Cartesian grid faces, which makes the integral evaluation straightforward in an immersed boundary (IB) framework. Discontinuous and noisy derivatives of velocity and pressure at the fluid-structure interface are avoided and far-field (smooth) velocity and pressure information is used. We re-visit the approach to compute hydrodynamic forces and torques through force/torque balance equation in a Lagrangian frame that some of us took in a prior work (Bhalla et al., J Comp Phys, 2013). We prove the equivalence of the two approaches for IB methods, thanks to the use of Peskin's delta functions. Both approaches are able to suppress spurious force oscillations and are in excellent agreement, as expected theoretically. Test cases ranging from Stokes to high Reynolds number regimes are considered. We discuss regridding issues for the moving CV method in an adaptive mesh refinement (AMR) context. The proposed moving CV method is not limited to a specific IB method and can also be used, for example, with embedded boundary methods