1
vote
1answer
52 views

Rigid body collision, 3 circles in contact

I'm working on a 2d physics simulation. It's a continuous time simulation, that is, it uses swept shapes over the time-frame and geometrical/vector 'analysis' to determine most immediate time of ...
0
votes
1answer
55 views

Failure of a simple stat mech simulation

so I did a simple simulation that I thought would yield a Boltzmann distribution, but it failed to, and I was wondering if anyone has insight into why it failed. Ok, so I had a simple discrete system ...
1
vote
0answers
41 views

How does the Lennard Jones Potential changes for interaction between particles of different sizes?

I am interested in incorporating a Lennard-Jones potential in a simulation. When the interaction only involves the same type of particle, with same characteristics, we can use reduced units, scaling ...
1
vote
0answers
57 views

Gerstner Wave Formula to Vertex information $(x,y,z)$

I have created a program, that draws a plane mesh, and allows to me edit or animate all the separate vertices along it. The idea behind it, is for me to be able to create a gridMesh based wave ...
1
vote
0answers
40 views

Find out ground sates for large 2D classical spin model

Reaching the ground state of a large 2D classical spin model (e.g. classical Heisenberg model) might be a relatively difficult task while using conventional "flip/reject" Monte Carlo method. The ...
4
votes
1answer
75 views

Can computers accurately model all of the details (to the subatomic level) of macro objects in collisions?

Frequently when trying to solve cosmology questions physicists turn to computer simulations of the universe (albeit massively simplified) in order to verify or disprove their hypotheses. This got me ...
0
votes
2answers
95 views

Grand canonical ensemble with interaction, simulation doubts

First of all, this is a follow-up of my first question. The idea is the same, every state would consist of some particles in a line (with an energy associated) and the particles can not be nearer than ...
0
votes
1answer
231 views

Acceptance probability 2D Ising Model

Disclaimer: I just found a possible solution - eventhough i don't really understand, whats wrong with my prior approach. Edit: I just tried to calculate it from scratch and found the following: $E ...
1
vote
1answer
108 views

Simulation of an one dimensional driven diffusive system [closed]

I'm currently writing a simulation in python with scipy and matplotlib to reproduce an one dimensional driven diffusive system described in this paper from M.R. Evans et al. The system consists of ...
7
votes
1answer
566 views

Riemann invariants…Any physical interpretation?

I am really new to the CFD simulation, and started some simple algorithms recently. I then got introduced to the Riemann Invariants. Can any one provide some physical interpretation? Also, why is ...
2
votes
0answers
132 views

simple MHD simulation code for (self) education and play with [closed]

I would like some super simple computational code for solving magnetohydrodynamics problems. High accuracy nor performance is not my concern. I wan't it just to visually explore qualitative behavior ...
0
votes
3answers
169 views

Meaning of “Gradient with respect to coordinates of particle” in SPH

I'm currently trying to implement a simple SPH simulation based on a variety of papers. However as I'm not a trained physicist nor mathematician I have a small issue with the following notation and ...
6
votes
1answer
505 views

What are the limitations of Smoothed-Particle Hydrodynamics?

I've been excited by some of the possibilities of Smoothed-Particle Hydrodynamics (SPH). I have seen some very exciting demonstrations of their use in 3D graphics, but I am wondering how well the ...
1
vote
1answer
108 views

Deriving the change in the Helmholtz free energy in the context of the free energy perturbation method

I am reading Free Energy Calculations: Theory and Applications in Chemistry and Biology by Chipot and Pohorille. At the beginning of the text (page 19, for example), the authors define the Helmholtz ...
2
votes
1answer
109 views

Mutual Interaction of $N$-Particles in a Cartesian Plane

I am making a simulation of $N$-Particles in a cartesian plane and need help with understanding the basics. At anytime, in my particle system, I will have $N$ number of particles. I am treating the ...
1
vote
1answer
256 views

Help: 3D visualization of magnetic field around moving point charge

The diagrams I'm able to find online only show the concentric field lines in the particle's plane, perpendicular to the motion, which of course generalizes to the cylindrical shape around a conductor. ...
1
vote
2answers
831 views

How do I integrate the Poisson equation to determine the electric potential along a particular direction (e.g., $z$)?

This question is a sequel of sorts to my earlier (resolved) question about a recent paper. In the paper, the authors performed molecular dynamics (MD) simulations of parallel-plate supercapacitors, ...
2
votes
1answer
339 views

How is acceleration computed in the leap-frog algorithm in computer simulations of particles?

I am reading the book Computer Simulation of Liquids by Allen and Tildesley (here is another link). On page 80, the authors describe the leap-frog algorithm, which is used extensively in molecular ...
1
vote
2answers
226 views

How to decide convergence of the ensemble average in a Monte Carlo process

I am simulating electromagnetic scattering off a rough surface. The usual process is to do a Monte Carlo simulation, which is briefly described as follows. Generate a randomly rough surface, and ...
5
votes
2answers
298 views

In condensed matter simulations, how is particle number density computed in practice?

I have been reading a recent paper. In it, the authors performed molecular dynamics (MD) simulations of parallel-plate supercapacitors, in which liquid resides between the parallel-plate electrodes. ...
2
votes
0answers
314 views

Comparison of different ab-initio codes

One may find on the web a lot of different computational packages to perform "ab-initio" calculations of electron structure of the solids. Usually, the documentation is not quite transparent about the ...
1
vote
1answer
222 views

Numerical algorithms to generate a random wavefunction from a thermal ensemble

I am seeking an algorithm to generate a random wavefunction = $\sum {c_i |\varphi _i\rangle }$ from a thermal ensemble, whose density matrix $\rho \sim e^{-\beta H}$, without the need to diagonalize ...
5
votes
3answers
336 views

How to include random force in the simulation (Classical Molecular Dynamics)

I need to implement a random force in my code according to the fluctuation dissipation theorem. I have a Gaussian distribution function ready width average 0 and distribution 1 and I know I need to ...
5
votes
3answers
233 views

Simulate a physical impact of objects made of finite, small elements

I want to simulate an impact between two bodies according to gravity, and eventually considering other forces to stick matter together. I'd like to use python to do this, but I am open to ...