Related papers: Simulating isothermal aging of snow
The melting point of ice Ih, as well as the temperature of maximum density (TMD) in the liquid phase, has been computed using the path integral Monte Carlo method. Two new models are introduced; TIP4PQ_D2O and TIP4PQ_T2O which are…
Context. Cometary dust particles are subjected to various forces after being lifted off the nucleus. These forces define the dynamics of dust, trajectories, alignment, and fragmentation, which, in turn, have a significant effect on the…
An open source software package for simulating thermal neutron propagation in geometry is presented. In this system, neutron propagation can be treated by either the particle transport method or the ray-tracing method. Supported by an…
A diathermal wall between two heat baths at different temperatures can be mimicked by a layer of independent spin pairs with some internal energy and where each spin $\sigma_a$ is flipped by thermostat $a$ ($a=1,2$). The transition rates…
We investigate the effect of heating by luminosity sources in a simulation of clustered star formation. Our heating method involves a simplified continuum radiative transfer method that calculates the dust temperature. The gas temperature…
We study an induced dynamics in the space of energy of single-spin-flip Monte Carlo algorithm. The method gives an efficient reweighting technique. This dynamics is shown to have relaxation times proportional to the specific heat. Thus, it…
The basic problem in equilibrium statistical mechanics is to compute phase space average, in which Monte Carlo method plays a very important role. We begin with a review of nonlocal algorithms for Markov chain Monte Carlo simulation in…
We present a non-isothermal mesoscopic model for investigation of the phase transition dynamics of thermoresponsive polymers. Since this model conserves energy in the simulations, it is able to correctly capture not only the transient…
We examine the molecular dynamics of crystal growth in the presence of surface melting and surface impurities, and from this propose a detailed microscopic model for the growth of ice from the vapor phase. Our model naturally accounts for…
The effect of edge on wetting and layering transitions of a three-dimensional spin-1/2 Ising model is investigated, in the presence of longitudinal and surface magnetic fields, using mean field (MF) theory and Monte Carlo (MC) simulations.…
Polar precipitation archived in ice caps contains information on past temperature conditions. Such information can be retrieved by measuring the water isotopic signals of $\delta{}^{18}\mathrm{O}$ and $\delta\mathrm{D}$ in ice cores. These…
The influence of nearest-neighbor diffusion on the decay of a metastable low-coverage phase (monolayer adsorption) in a square lattice-gas model of electrochemical metal deposition is investigated by kinetic Monte Carlo simulations. The…
Freezing is a fundamental physical phenomenon that has been studied over many decades; yet the role played by surfaces in determining nucleation has remained elusive. Here we report direct computational evidence of surface induced…
We introduce and test a new and highly efficient method for treating the thermal and radiative effects influencing the energy equation in SPH simulations of star formation. The method uses the density, temperature and gravitational…
In this text, we implement a monte carlo algorithm to study thermodynamic properties of ice. Our program, written in Python, is open-sourced and available at https://github.com/AKnightWing/ColdAsIce. We develop a novel scheme to compute the…
The objective of this work is to investigate the coupling of fluid dynamics, heat transfer and mass transfer during the impact and evaporation of droplets on a heated solid substrate. A laser-based thermoreflectance method is used to…
The Kinetic Monte Carlo (KMC) method has become an important tool for examination of phenomena like surface diffusion and thin film growth because of its ability to carry out simulations for time scales that are relevant to experiments. But…
Iceberg meltwater is a critical freshwater flux from the cryosphere to the oceans. Global climate simulations therefore require simple and accurate parameterisations of iceberg melting. Iceberg shape is an important but often neglected…
It has been shown by means of impurity diffusion simulation that ion-implanted boron redistribution at the annealing temperatures 800^{\circ}C and lower is governed by the long-range migration of nonequilibrium impurity interstitials…
An in-house numerical model has been developed to study the temperature destratification induced by bubble plumes in reservoirs. The mean flow of the mixed fluid phase is solved using the Reynolds-Averaged Navier-Stokes Equations, while the…