Related papers: An improved current potential method for fast comp…

In many imaging applications where segmented features (e.g. blood vessels) are further used for other numerical simulations (e.g. finite element analysis), the obtained surfaces do not have fine resolutions suitable for the task. Increasing…

Finding an $\epsilon$-stationary point of a nonconvex function with a Lipschitz continuous Hessian is a central problem in optimization. Regularized Newton methods are a classical tool and have been studied extensively, yet they still face…

This work uses high-order discontinuous Galerkin discretization techniques as a generic, parameter-free, and reliable tool to accurately predict transitional and turbulent flows through medical devices. Flows through medical devices are…

With the advent of neoclassically optimised stellarators, optimising stellarators for turbulent transport is an important next step. The reduction of ion-temperature-gradient-driven turbulence has been achieved via shaping of the magnetic…

As an alternative to solving of Poisson equation in Particle-in-Cell methods, a new construction of current density exactly satisfying continuity equation in finite differences is developed. This procedure called density decomposition is…

Stark deceleration enables the production of cold and dense molecular beams with applications in trapping, collisional studies, and precision measurement. Improving the efficiency of Stark deceleration, and hence the achievable molecular…

We present a novel method for minimizing the effects of radiative depolarization in electron storage rings by use of vertical orbit bumps in the arcs. Electron polarization is directly characterized by the RMS of the so-called spin orbit…

Resolving numerically Vlasov-Poisson equations for initially cold systems can be reduced to following the evolution of a three-dimensional sheet evolving in six-dimensional phase-space. We describe a public parallel numerical algorithm…

Inversions of rotational splittings have shown that there exists at the base of the solar convection zone a region called the tachocline in which high radial gradients of the rotation rate occur. The usual linear regularization methods tend…

We present recent improvements in the simulation of regolith sampling processes in microgravity using the numerical particle method smooth particle hydrodynamics (SPH). We use an elastic-plastic soil constitutive model for large deformation…

The capacitance of capacitive energy storage devices can not be directly measured, but can be estimated from the input and output signals expressed in the time or frequency domains. Here the time-domain voltage-charge relationship in…

The Ryabinkin-Kohut-Staroverov (RKS) and Kanungo-Zimmerman-Gavini (KZG) methods offer two approaches to find exchange-correlation (XC) potentials from ground state densities. The RKS method utilizes the one- and two-particle reduced density…

Single Helical Axis (SHAx) configurations are emerging as the natural state for high current reversed field pinch (RFP) plasmas. These states feature the presence of transport barriers in the core plasma. Here we present a method for…

One aim of dimensionality reduction is to discover the main factors that explain the data, and as such is paramount to many applications. When working with high dimensional data, autoencoders offer a simple yet effective approach to learn…

Stencil computation is one of the most important kernels in various scientific and engineering applications. A variety of work has focused on vectorization techniques, aiming at exploiting the in-core data parallelism. Briefly, they either…

In this work we consider the problem of optimizing a stellarator subject to hard constraints on the design variables and physics properties of the equilibrium. We survey current numerical methods for handling these constraints, and…

'Coil design' is an inverse problem in which arrangements of wire are designed to generate a prescribed magnetic field when energized with electric current. The design of gradient and shim coils for magnetic resonance imaging (MRI) are…

We consider a class of nonconvex nonsmooth multicomposite optimization problems where the objective function consists of a Tikhonov regularizer and a composition of multiple nonconvex nonsmooth component functions. Such optimization…

Weighted compact nonlinear schemes (WCNS) [Deng and Zhang, JCP 165(2000): 22-44] were developed to improve the performance of the compact high-order nonlinear schemes (CNS) by utilizing the weighting technique originally designed for WENO…

We have developed a novel Monte Carlo method for simulating the dynamical evolution of stellar systems in arbitrary geometry. The orbits of stars are followed in a smooth potential represented by a basis-set expansion and perturbed after…