# Questions tagged [simulations]

The tag has no usage guidance.

1,077 questions
Filter by
Sorted by
Tagged with
47 views

### Metropolis-Hastings and underlying Markov process

I FOUND A MISTAKE IN THE POST, TO BE EDITD OR DELETED SHORLTY I tried to understand the workings of the Metropolis-Hasting algorithm. As far as I can understand, it allows to draw samples from an ...
• 407
10 views

8 views

### Create a plausible vibration signal from a sliding RMS window

So lets say I have observed many vibration signals and build a sliding window of the RMS of said vibration signals as well as other statistical measurements such as the variance and nth order ...
• 141
19 views

### Adding gravity to (Elcott07 "Stable, Circulation-Preserving, Simplicial Fluids")

The paper is quite popular, and I thought to implement it as an exercise for DEC. I added an initial vortex force, and indeed it's preserved throughout the simulation due to incompressibility and ...
• 101
16 views

### Are there any resources for simulating the Kuramoto model/solving the system of coupled ODEs for the oscillators?

I want to simulate the Kuramoto model and I am following a review by Strogatz to read up the theory. Are there resources that explain how to simulate the model and solve the coupled ODEs, or the mean ...
5 views

### What formula should I use to calculate the total energy of a linear polymer chain?

I am using the Bead-spring model to model a polymer chain. Suppose, the polymer has 3 monomers in its chain: $$A_1-A_2-A_3$$ And, I can use the following formula to calculate the total Lennard-Jones ...
• 109
28 views

### Where should I apply the harmonic spring function in the case of polymer simulation?

I need to simulate the off-lattice movement of a polymer chain in a 3D space using Monte Carlo simulation. Note this is not a simulation of polymer growth; rather this is about polymer movement/motion....
• 109
30 views

### $n$-link pendulum solution methods

I want to simulate a $n$-link pendulum without using the Lagrangian method. The reason is that the formula for the Lagrangian becomes large and depends on the number of the links $n$. Therefore, I ...
• 159
18 views

### Polymer simulation

I am absolutely new to Polymer simulation. I am trying to understand the simulation by analyzing source code written by others. Can anyone tell me what are the differences between the following three ...
• 109
48 views

### Which forces act on a ball in the air (including air friction) including their equation? [closed]

I am trying to write a ballistic throw simulation with the forces which act on the ball I take into consideration being the gravitational force, and drag. I know that the Magnus effect and another ...
• 103
1 vote
62 views

24 views

### How to take care of singularities or discontinuities in a numerical Time-Dependent Schrodinger Equation simulation using the Split-Operator method?

I am currently trying to numerically simulate a two-dimensional Hydrogen atom using a split-operator method (elaborated below). I have tested that the code produces a good result for an isotropic ...
• 15
5k views

### Is angular momentum just a convenience?

I'm wondering whether angular momentum is just a convenience that I could hypothetically solve any mechanics problems without ever using the concept of angular momentum. I came up with this question ...
• 627
24 views

### Do most climate models use the absolute temperature or the temperature anomaly under the hood?

Climate models for the climate impact of anthropogenic global warming almost never report the climate impacts of the average global temperature reaching some absolute threshold. Instead, they almost ...
• 41.3k
1 vote
37 views

### How can I model convection in a simple and effective way?

I have a project in which we're trying to model heating of a room. We are modelling the room in 3D. The room has realistic walls, windows and a door, which will loose heat through them. To balance the ...
• 11
14 views

### Is there a way to estimate the force of this impact?

https://www.youtube.com/watch?v=UKuPuRRNa5Q&ab_channel=BigCatPublicSafetyAct Let's assume that the tiger weighs 500 lbs, is there a way to estimate the force absorbed by its shoulders when it hit ...
39 views

### When should smartphone be declared as static using accelerometer reading? [closed]

I wish to know when should I declare an smart phone static by reading its 3 axis accelerometer reading. I have notice that modulo of all 3 axes gives 9.8 when ever smart phone is static. What does it ...
• 237
1 vote
34 views

### Reference request: Numerical techniques, Monte-Carlo (MC), Density Matrix Renormalization Group (DMRG), Dynamical Mean-Field Theory (DMFT)

I am an undergraduate student and my previous learning in physics is more on theory instead of numerics. I would be very grateful if you can point me to good introductory lecture notes/lecture videos ...
857 views

### How can one detect a device (phone) falling using its 3-axis accelerometer?

I am working on an Android app to use its 3-axis inbuilt accelerometer to find if the phone falls. I googled and found that if the following value is less than 2: $$a = \sqrt{x^2 + y^2 + z^2}$$ (where ...
• 237
51 views

### Pointers for dealing with instability in fluid simulation?

I have been working on a 2D fluid simulator, and I have fallen victim to the problem of jumping in before doing enough research. I naively implemented a finite difference form of the continuity and ...
46 views

### Boundary conditions and normalization of wavefunction in time-dependent Schrodinger equation

I came across this confusion when trying to numerically solve the time-dependent Schrodinger equation. The short version of my question is: when solving for time-dependent wavefunctions $\psi(t,x)$ on ...
18 views

### Calculate optimal new vector for a virtual spaceship

I'm working on a 3D game as a demo for some unrelated technology. In it I have a spaceship that I move in the game engine by applying a force along a vector $(x,y,z)$. So I'm happily flying along ...
• 251
73 views

### Simulation algorithm question

I have been looking at various simulation algorithms, like Barnes-Hut and a few others, and the Laplace operator, the wave equation, etc., and I was wondering why it's necessary to compute the mass/...
30 views

### Rate Dispersion is the distribution of rates on a rate spectrum. What are the quantities physically associated in a Rate Spectrum?

I am working on a project to decipher the "origin of rate dispersion" arriving because of Heterogeneity in our system using correlation function analysis in Python. Rate Dispersion, non-...
32 views

### How does (air) pressure balance under constraints?

I'm playing around with building a very basic simulation of (air) pressure in a grid environment. The idea is that I have this 2d (or 3d) grid of cells. Each cell has a certain atmosphere. Between the ...
103 views

### Meaning of 'thermalization' in Markov Chain Monte Carlo simulations

In performing MCMC simulations, it is standard practice to 'equilibriate' or 'thermalize' the system and then discard the initial data before useful sampling is done. My question is about the concept ...
78 views

### Making Lorentz equation dimensionless for simulations

$$\frac{dv}{dt}=\frac{e}{m}(E+\frac{v}{c}\times B)$$ I'm making a numerical simulation and I need to make the equation dimensionless. I'm having trouble in doing so. I was ...
• 23
32 views

### Difference between the softening parameter in a Plummer model of a star cluster vs. the softening parameter in a softened gravitational potential

I am doing a gravitational N-body simulation of a star cluster. I set up my initial conditions using the Plummer model, which has potential $$\Phi_P = -\frac{G M}{\sqrt{r^2 + a^2}},\tag{1}$$ where $M$ ...
25 views

### Mixing position and force based simulation

I have read quite a bit on rigid body physics engines and feel I have a reasonable understanding of how systems evolve with force based simulation. What I am kind of confused on still is how you mix ...
• 784
100 views

### Bogoliubov-de Gennes equation and how to diagonalize

Can someone please share how you would go on to solve the Bogoliubov-de Gennes equation for a Bose-Einstein condensate with a potential. In particular, how would you deal with the diagonalization when ...
• 303
1 vote
140 views

### Poisson-Nearst-Planck equations with normal distribution as initial condition

I am having a hard time trying to think of this model: Imagine you have a normal distribution of +q charges in 2D $$\rho(r,t=0)=\rho_0 e^{-\frac{r^2}{\sigma_0^2}}$$ where $\sigma_0$ is the width of ...
• 11
1 vote
38 views

### Absorbing regions in Physics

How would you define an absorbing region (at the walls of a finite numerical simulation) in physics and how would you implement it? Can it be simulated by a high potential, low negative potential at ...
• 303
18 views

### What is the reason that in DPD (LAMMPS) we can define a protein (or any material with a flexible geometry) without the material loosing its shape?

I am new to LAMMPS/DPD (Dissipative particle dynamics), but as I read through DPD model which is a coarse graining model, there are only three main forces between beads (beads are the spheres ...
20 views

### LAMMPS: How to enter a simple potential for Dissipative Particle Dynamics (DPD)?

I apologize if this is off-topic and if so would be grateful for someone to point me to a relevant source of information. Dissipative Particle Dynamics (DPD) is a quick simulation technique that ...
1 vote
20 views

### Difference Between a Time-Domain and a Simulation Domain Acoustic Pressure Simulation

I would appreciate it if someone would be able to clarify my doubts. I'm currently working on acoustic simulations in COMSOL and is currently confused on the information provided by the simulations ...
20 views

### Transparent Boundary Conditions for Finite Difference ADI PR 2D TDSE solution

I want to put (non-dirichlet) boundary conditions inside the code I wrote to solve the 2dim TDSE using the alternating direction implicit Peaceman - Rachford method.  (1 + iB\Delta t/2 ) \psi^{n+1/2}...
• 101
14 views

### Trinocular optical ray tracing reference

I am new to the optical tracing, was looking for ray tracing refence that shows what kind of lens used in trinocular optical train. Most what is out there is an explnation of microscope theory and ...
1 vote
37 views

### Conservation of energy in a numerical simulation of a particle in magnetic field

I'm currently working on my physics project in Mathematica. My goal is to simulate a 3d trajectory of a charged particle, moving in constant magnetic field generated by a loop (with the electric ...
• 13
25 views

### Conceptualizing time series data of fluctuating sizes of particle aggregates

I am working with simulation data (a time series of positions) of aggregating particles. I want to look at the overall distribution of aggregate size. A colleague calculated the number of aggregates ...
• 448
1 vote
In HMC (good intro here) the Hamiltonian is defined as \begin{align} H(p,\theta) &= -\log p(q,\theta) \\ &= -\log p(q) - \log p(\theta). \end{align} Here $p$ is the momentum coordinate and \$\...