计算物理
Fine-scale simulation of complex systems governed by multiscale partial differential equations (PDEs) is computationally expensive and various multiscale methods have been developed for addressing such problems. In addition, it is…
In the present study, the multiphase volume distribution problem, where there can be an arbitrary number of phases, is addressed using a consistent and conservative volume distribution algorithm. The proposed algorithm satisfies the…
An implementation of a lattice-based approach for computing the topological skyrmion charge is provided for the open source micromagnetics code MuMax3. Its accuracy with respect to an existing method based on finite difference derivatives…
A problem of the erroneous duality gap caused by the presence of symmetries is solved in this paper utilizing point group theory. The optimization problems are first divided into two classes based on their predisposition to suffer from this…
The computational cost of fluid simulations increases rapidly with grid resolution. This has given a hard limit on the ability of simulations to accurately resolve small scale features of complex flows. Here we use a machine learning…
The Monte Carlo method is a powerful technique for computing thermodynamic magnetic states of otherwise unsolvable spin Hamiltonians, but the method becomes computationally prohibitive with increasing number of spins and the simulation of…
The integration of the equations of motion of N interacting particles, represents a classical problem in many branches of physics and chemistry. The direct N-body problem is at the heart of simulations studying Coulomb Crystals. We present…
We study a model of two-dimensional classical dimers on the square lattice with strong geometric constraints (there is exactly one bond with the nearest point for every point in the lattice). This model corresponds to the quantum dimer…
Transcranial ultrasound therapy is increasingly used for the non-invasive treatment of brain disorders. However, conventional numerical wave solvers are currently too computationally expensive to be used online during treatments to predict…
Monte Carlo methods provide detailed and accurate results for radiation transport simulations. Unfortunately, the high computational cost of these methods limits its usage in real-time applications. Moreover, existing computer codes do not…
We develop a numerical solver for three-dimensional wave propagation in coupled poroelastic-elastic media, based on a high-order discontinuous Galerkin (DG) method, with the Biot poroelastic wave equation formulated as a first order…
Strongly interacting quantum systems described by non-stoquastic Hamiltonians exhibit rich low-temperature physics. Yet, their study poses a formidable challenge, even for state-of-the-art numerical techniques. Here, we investigate…
Density functional theory (DFT) is an efficient instrument for describing a wide range of nanoscale phenomena: wetting transition, capillary condensation, adsorption, etc. In this paper, we suggest a method for obtaining the equilibrium…
In this paper, a diffuse-interface lattice Boltzmann method (DI-LBM) is developed for fluid-particle interaction problems. In this method, the sharp interface between the fluid and solid is replaced by a thin but nonzero thickness…
A PCA-based, machine learning version of the SPH method is proposed. In the present scheme, the smoothing tensor is computed to have their eigenvalues proportional to the covariance's principal components, using a modified octree data…
The Kirchhoff integral is a fundamental integral in scattering theory, appearing in both the Kirchhoff approximation, as well as the small slope approximation. In this work, a functional Taylor series approximation to…
We present new approaches for solving constrained multicomponent nonlinear Schr\"odinger equations in arbitrary dimensions. The idea is to introduce an artificial time and solve an extended damped second order dynamic system whose…
Quantum Monte Carlo (QMC) forces have been studied extensively in recent decades because of their importance with spectroscopic observables and geometry optimization. Here we benchmark the accuracy and statistical cost of QMC forces. The…
Magnetization dynamics in magnetic materials is modeled by the Landau-Lifshitz-Gilbert (LLG) equation. In the LLG equation, the length of magnetization is conserved and the system energy is dissipative. Implicit and semi-implicit schemes…
We present a high-energy neutrino event generator, called LeptonInjector, alongside an event weighter, called LeptonWeighter. Both are designed for large-volume Cherenkov neutrino telescopes such as IceCube. The neutrino event generator…