Related papers: On the importance of numerical integration details…
Computational modelling of metal-electrolyte reactions is central to the understanding and prediction of a wide range of physical phenomena, yet this is often challenging owing to the presence of numerical oscillations that arise due to…
The choice of a time integration scheme is a crucial aspect of any transient fluid simulation, and Smoothed-Particle Hydrodynamics (SPH) is no exception. The influence of the time integration scheme on energy balance is here addressed. To…
The Discrete Element Method is widely employed for simulating granular flows, but conventional integration techniques may produce unphysical results for simulations with static friction when particle size ratios exceed $R \approx 3$. These…
A new symplectic time-reversible algorithm for numerical integration of the equations of motion in magnetic liquids is proposed. It is tested and applied to molecular dynamics simulations of a Heisenberg spin fluid. We show that the…
Accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and…
We present a reversible and symplectic algorithm called ROLL, for integrating the equations of motion in molecular dynamics simulations of simple fluids on a hypersphere $\mathcal{S}^d$ of arbitrary dimension $d$. It is derived in the…
A new approach is developed to integrate numerically the equations of motion for systems of interacting rigid polyatomic molecules. With the aid of a leapfrog framework, we directly involve principal angular velocities into the integration,…
Viscoelastic flows occur widely, and numerical simulations of them are important for a range of industrial applications. Simulations of viscoelastic flows are more challenging than their Newtonian counterparts due to the presence of…
The Stoermer-Verlet-leapfrog group of integrators commonly used in molecular dynamics simulations has long become a textbook subject and seems to have been studied exhaustively. There are, however, a few striking effects in performance of…
Simulation of many-particle system evolution by molecular dynamics takes to decrease integration step to provide numerical scheme stability on the sufficiently large time interval. It leads to a significant increase of the volume of…
Numerical simulation of compressible fluid flows is performed using the Euler equations. They include the scalar advection equation for the density, the vector advection equation for the velocity and a given pressure dependence on the…
A new methodology is developed to integrate numerically the equations of motion for classical many-body systems in molecular dynamics simulations. Its distinguishable feature is the possibility to preserve, independently on the size of the…
Preserving scalar boundedness is important for numerical schemes used in turbulent compressible multi-component flow simulations to prevent unphysical results and unstable simulations. However, ensuring scalar boundedness for high-order,…
Accounting for the Reynolds number is critical in numerical simulations of turbulence, particularly for subsonic flow. For Smoothed Particle Hydrodynamics (SPH) with constant artificial viscosity coefficient alpha, it is shown that the…
In this paper, we consider mathematical modeling and numerical simulation of non-isothermal compressible multi-component diffuse-interface two-phase flows with realistic equations of state. A general model with general reference velocity is…
Many applications of computational fluid dynamics require multiple simulations of a flow under different input conditions. In this paper, a numerical algorithm is developed to efficiently determine a set of such simulations in which the…
Scale-resolving flow simulations often feature several million [thousand] spatial [temporal] discrete degrees of freedom. When storing or re-using these data, e.g., to subsequently train some sort of data-based surrogate or compute…
Time-symmetric integration schemes share with symplectic schemes the property that their energy errors show a much better behavior than is the case for generic integration schemes. Allowing adaptive time steps typically leads to a loss of…
Despite its numerical challenges, finite element method is used to compute viscous fluid flow. A consensus on the cause of numerical problems has been reached; however, general algorithms---allowing a robust and accurate simulation for any…
Atmospheric systems incorporating thermal dynamics must be stable with respect to both energy and entropy. While energy conservation can be enforced via the preservation of the skew-symmetric structure of the Hamiltonian form of the…