Subdiscipline of theoretical physics which consists in the study of physics problems using numerical algorithms. PLEASE NOTE that questions about computational methods and/or programming are OFF-TOPIC on Physics.SE.

learn more… | top users | synonyms (1)

4
votes
1answer
490 views

Philae lander simulation off by factor of ~3

I'm trying to simulate the Philae landing by writing a program to compute the position of the lander vs time. According to various mission websites, the orbiter will match its orbit to the rotation of ...
2
votes
0answers
82 views

Modelling gravitational potential of a galaxy

I am interested in modelling the gravitational potential of a disc-shaped galaxy with radius $R$, i.e. solving the 2D Poisson equation numerically by Gauss-Seidel relaxation: $$\nabla^2 \phi = 4\pi G ...
0
votes
0answers
31 views

Optimize locations of vertices and edge length using physics without simulation

I have a directed graph with set of vertices and edges. Each vertices don't have a position coordinate to be displayed on 2D computer screen. I want to add unique position coordinate to each vertices ...
-2
votes
1answer
57 views

Software for calculating the trajectory of a body in a coordinate system

I'm looking for a program, which would simulate the path of motion of a body in a coordinate system, given the force acting upon the body. I'd type in the initial conditions such as the velocity, ...
4
votes
3answers
328 views

Good source for numerical simulations of Wigner function?

I'm interested in simulating the time evolution of a Wigner function for a harmonic oscillator (and possibly some other potentials) and I can't seem to find a good resource for that. My background in ...
5
votes
1answer
134 views

Least action principle — numerical simulation strangeness

I'm trying to get some experience with the least action principle, and for this I chose a simple 1-dimensional problem of a particle moving in some field. The least action principle would then look ...
5
votes
1answer
313 views

Ray tracing in General Relativity

I would like to find out what one would see at the Schwarzschild radius of a massive non-rotating black hole, if the black hole is surrounded by a bright ring. For that, I would place the observer at ...
2
votes
0answers
210 views

How to compute matsubara frequency summation over computer?

Matsubara frequency sum usually takes the following form: $S_\eta = \frac{1}{\beta}\sum_{i\omega_n} g(i\omega_n).$ But in my problem, $g(i\omega_n)$ is a lengthy expression, which can not be ...
1
vote
1answer
145 views

Numerical error with simulation of electric charge in homogeneous magnetic field [closed]

So, I am trying to make an 2D animation of electric charge in homogeneous magnetic field which is perpendicular to charge's velocity. I've got the "circular" motion but the problem is that the speed ...
1
vote
2answers
1k views

Gas viscosity at high pressure, high temperature

EDIT 1 PER COMMENTS I am wanting to model nitrogen gas viscosity as a function of pressure and temperature OR learn of an existing equation that models nitrogen viscosity for the pressure and ...
1
vote
1answer
387 views

Turning a finite difference equation into code (2d Schrodinger equation)

I am trying to convert the following finite difference equations into code (taken from the bottom of page 12 of this thesis by Maike Schulte Numerical Solution of the Schrodinger Equation on Unbounded ...
2
votes
2answers
177 views

Kernel normalization in Smoothed Particle Hydrodynamcs

I am trying to write some code for a fluid simulation. I have read about Smoothed Particle Hydrodynamics (SPH) and my question is related to the properties of the smoothing kernel. How do I set the ...
0
votes
1answer
427 views

Physical simulation in python [closed]

Is there a "standard" python package used to aid in physical simulations? What is the most popular? edit: perhaps I should have worded this question differently. Something more to the effect: Could ...
0
votes
1answer
242 views

Position based dynamics constraint scaling factor

Reading through Müller et al., Position Based Dynamics, 2007 I got lost when passing from equation (5) $$\Delta p = \frac{C(p)}{|\nabla_pC(p)|^2}\nabla_pC(p)$$ to equation (6) (and applying the ...
1
vote
0answers
58 views

Does Planck's constant imply limits to computing *results*

... I don't mean quantum effects limiting hardware fabrication sizing. Such small scales have for some time exhibited issues. Rather, along the lines of imagining the smallest possible divisions of ...
1
vote
1answer
72 views

How do I set up the tridiagonal matrix for a heat diffusion with layers of different thermal diffusivity?

I have Scala code that recreates the Crank-Nicolson solutions for the diffusion equations, and matches 'Excel for Scientists and Engineers' (Joe Billo, Wiley). However, I would like to be able to ...
0
votes
0answers
83 views

Examples of application of detour matrices in physics?

Are there any good examples of application of detour matrices in physics?
1
vote
1answer
177 views

A puzzle of thermalization in simulating the 3D XY-model

I am learning the classical Monte Carlo simulation. When I simulate the 3D XY-model $$ \beta H = -\beta J \sum_{<i,j>} cos(\theta_i-\theta_j) $$ where $\beta$ is the inverse of the temperature ...
3
votes
0answers
69 views

Textbooks on algorithms for the perturbative calculation of High energy physics

For the perturbative calculation of High energy physics, I have known some packages such as FeynArts, FeynCalc, MadGraph, CompHEP, GiNaC, and so on. But I am wondering whether there exists a textbook ...
1
vote
0answers
41 views

Computer coding of Perdew Burke and Ernzerho (PBE) method

I have been assigned the task to write a computer code to implement the Perdew Burke and Ernzerhof (PBE) method. Does anyone know a good reference which can make the coding of this method easier?
0
votes
1answer
103 views

Broken Thin Lens Algorithm [closed]

We all know the thin lens equation. For $o$ being a horizontal object distance and $f$ being the focal length, the horizontal image distance $i$ is described by: $$\frac{1}{f} = \frac{1}{o} + ...
0
votes
1answer
225 views

Monte carlo simulation for continuous spin model (e.g. XY or Heisenberg model)

Unlike the Ising model, the XY model and the Heisenberg model have a continuous spectrum. So one need discretize them for a numerical simulation. But how to make sure the discretization procedure ...
1
vote
3answers
147 views

Why does the Metropolis algorithm allow changes even for ∆E > 0?

In the Metropolis Monte Carlo algorithm, why can you accept changes even for ∆E > 0 (provided that a random number is less than a given probability ratio, e.g. exp(-β∆E))?
3
votes
0answers
86 views

Non-linear Wave Equation - Numerical Methods

Motivation: I'm working with a highly non-linear spherical wave-like equation (second order PDE). The equation can be written on the form $$\ddot{u} = f(t, \dot{u},\dot{u}',u',u'')$$ where ...
2
votes
1answer
81 views

Question about Metropolis Monte Carlo in the case of equal energies

If configuration A is equal to configuration B in a Metropolis Monte Carlo method, do you still do the attempted update?
2
votes
0answers
37 views

Coefficients and Parameters for contracted Gaussian basis sets

This is a repost from Chemistry.stackexchange in the hopes that someone here will be able to help me. Any help at all would be greatly appreciated. As far as I understand, an STO-NG contracted ...
2
votes
2answers
210 views

How do you measure numerically the central charge of a system?

Let's say that you are doing some Monte-Carlo simulations of a statistical system on a lattice and you observe scale invariance, meaning that you are at a conformal point. Can you get a numerical ...
2
votes
1answer
343 views

Unexpected eccentricity in moon orbit simulation

I've used Mathematica's NDSolve function to calculate the orbit of the Moon around the Earth. I used the following initial positions (perigee): ...
3
votes
1answer
93 views

Complex semi-definite programming

I'm doing some calculations and I want to simulate them in python or matlab (or whatever). However I use hermitian matrices and I don't really manage to find a library which enables me to calculate ...
2
votes
2answers
136 views

Modeling a potential well

I attempted to simulate the interaction of a moving particle and a potential well in Mathematica. The particle should experience a force of -$1/r^2$, if the equation for the potential well is -$1/r$. ...
1
vote
1answer
171 views

Is there a normalized form of the Euler equation discretized with finite volumes?

I want to calculate a flux on my fpga using the Euler equations with the finite volume method. Unfortunately the values of the state variables differ a lot. For example the pressure has a value of ...
1
vote
0answers
125 views

How to numerically solve a complex equation? [closed]

I want to know that if you are given a very complex equation g(x)=A(T). How could you solve for x, which is a function of variable T. To be more specific, I encounter a polylogarithmic function I need ...
1
vote
0answers
156 views

XCrysDen Structure file [closed]

Does anyone know how to directly convert a .cif file to a Xcryden structure file(.xsf) ? I know how to extract the lattice vectors and the atom positions from a .cif file, but don't whether the .cif ...
3
votes
2answers
250 views

Boundary Element Method or Boundary Integral Method Computational Aspects

I have to solve a Helmholtz equation inside a simply connected domain. I know that in general the boundary integral can be written as, $$\phi(x)=\int_V G(x,x') \rho(x')\ d^3x'+\int_S ...
4
votes
0answers
106 views

Nature as a calculator [closed]

I wondered whether physical experiments can find the outcome of otherwise intractable mathematical problems. For example, solitons or other non-perturbative effects in QFT cannot be seen at any ...
0
votes
0answers
51 views

Poisson equation solver with specific boundary condition

I want to solve 2d Poisson equation with this Boundary conditions below $$u(-5,y) = 0 , \\ \frac{\partial u(x,y)}{ \partial x} = 0 \,\,{\rm at}\,\, x = -5\\ u(x,-5) = u(x,5)$$ Now my question is ...
1
vote
0answers
137 views

Lattice Gas Cellular Automata - HPP model square lattice

In the HPP model of LGCA, a square lattice is used and there is only one collision configuration as mentioned in figure (taken from the book Lattice Gas Cellular Automata and Lattice Boltzmann models ...
2
votes
1answer
255 views

Simulation of fluid flow using Euler equation

I have been looking on Euler's equations for a while and can't grasp one thing. Suppose we have initial system state with volumes of fluid "hanging" in air (time is frozen and equal to zero), each of ...
2
votes
0answers
106 views

Doing numerics in physics [closed]

Soon, I am going to write my master thesis in theoretical physics. I assume there, and later on in my career, I will have to do more serious numerics than I did up to this point. That's why I want to ...
1
vote
1answer
222 views

Ground state Phase Diagram of Bose-Hubbard Model

The Hamiltonian of Bose-Hubbard model reads as $$H=-J\sum\limits_{<i,j>}b_i^{\dagger}b_j+h.c.+\frac{U}{2}\sum\limits_{i}n_i(n_i-1)-\mu\sum\limits_in_i~~~~~~~~~(1)$$ For this we plot phase ...
0
votes
0answers
37 views

Verlet timestep

I am using velocity Verlet in molecular dynamics. Is just a simple question, if a simulation using time-step femtosecond, in velocity Verlet is just necessary $dt = 10^{-15}$ to use femtosecond?
1
vote
1answer
71 views

Calculate Euler equations of fluid dynamics without division?

I'm working on the calculation of the Euler equations with the finite volume method. Unfortunately I'm not allowed to do a division. So I'm wondering if there's a form which does not need a division. ...
2
votes
1answer
313 views

Numerical Ising Model: Swendsen–Wang algorithm, Percolation theory?

When you look at the original paper of Swendsen and Wang in 1987: "Nonuniversal critical dynamics in Monte Carlo simulations" it is somewhat mentioned that the proposed algorithm uses percolation ...
1
vote
2answers
82 views

Difference in calculated and simulated ellipsies

My task here is to determine orbit parameters, using current values: $\mu=GM$ - standard gravitational parameter $r$ - distance to the object with Mass $M$ $v$ - speed of the object in the point $r$ ...
3
votes
1answer
258 views

Euler Equations, Sod shock tube & conservation

Conservation of momentum? I am considering the Euler equations in conservative form and solving the Sod shock tube problem I have written a Godunov finite volume type solver. It solves for density ρ, ...
2
votes
0answers
63 views

Euler equation with single state variables [duplicate]

I want to simulate a flow using the Euler equations. For this reason I'm wondering if there's a modified version of the Euler equations. At the moment my formula looks like this: $$ ...
0
votes
0answers
117 views

Need help simulating a fan

I'm trying to write some code to approximately simulate a number fans in an arbitrarily shaped room. I'm not real strong in physics, and could use some help. If it helps I'm not interested in fine ...
0
votes
2answers
123 views

Energy of damped harmonic oscillator begins to increase with very large Q in numerical integration

I have numerically integrated the (reduced) homogeneous equation of a damped harmonic oscillator in order to see how the error propagates. $$\frac{d^2 X}{d\phi^2} + ...
9
votes
2answers
2k views

How do I calculate the Reynolds number in multiphase flows?

I am modeling a gas flowing through a liquid. How do I calculate the Reynolds number in multiphase flows? And, at what Reynolds number should I consider the flow to be turbulent? The problem is of a ...
3
votes
1answer
117 views

Update velocity or position first in computation?

I am trying to make a simulation of a vibrating string. The string is divided into $n$ points, and each point along the string is acted upon by a force due to the positioning of its neighbors. I ...