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)

7
votes
2answers
1k views

Physics of simple collisions

I'm building a physics simulator for a graphics course, and so far I have it implementing gravitational and Coulomb forces. I want to add collisions next, but I'm not exactly sure how to go about ...
7
votes
5answers
812 views

Church-Turing hypothesis as a fundamental law of physics

The Church-Turing hypothesis says one can not build a computing device which has more computing power (in terms of computability) than the abstract model of Turing machine. So, there is something in ...
3
votes
1answer
191 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 ρ, ...
13
votes
1answer
3k views

Solving Schrödinger's equation with Crank-Nicolson method

I am trying to numerically solve Schrödinger's equation with Cayley's expansion ($\hbar=1$) $$\psi(x,t+\Delta t)=e^{-i H\Delta t}\psi(x,t)\approx\frac{1-\frac{1}{2}i H\Delta t}{1+\frac{1}{2}i H\Delta ...
2
votes
1answer
600 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
1answer
67 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. ...
4
votes
2answers
1k views

Radial Schrodinger equation with inverse power law potential

Recently I read a paper about solving radial Schrodinger equation with inverse power law potential. Consider the radial Schrodinger equation(simply set $\mu=\hbar=1$): ...
3
votes
2answers
180 views

Simulate the universe?

Alright, Lets assume that I have a computer with limited calculation speed (1-4GHz) but unlimited parallel processing capability and unlimited memory capacity to go with it. Under this assumption ...
15
votes
5answers
969 views

Home-made lattice calculation?

The topic of Lattice QCD or Lattice gauge theory or even Lattice field theory is quite old now. And the main reason for the interest in the topic is the ability to calculate nonperturbative stuff on a ...
19
votes
5answers
3k views

What do theoretical physicists need from computer scientists?

I recently co-authored a paper (not online yet unfortunately) with some chemists that essentially provided answers to the question, "What do chemists need from computer scientists?" This included the ...
8
votes
3answers
2k views

How are physics and computer science getting united? [closed]

How is theoretical computer science getting united with physics? Phenomena like Quantum Computing uses Quantum Mechanics to be able to compute things, how are computers helping not just to model our ...
4
votes
2answers
241 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 ...
4
votes
2answers
2k views

On numerically solving the Schrödinger equation

I just read a paper 'A pocket calculator determination of energy eigenvalues' by J Killingbeck (1979). Link: http://iopscience.iop.org/0305-4470/10/6/001 I have some questions about section 2. Why ...
7
votes
2answers
5k views

Solving one dimensional Schrodinger equation with finite difference method

Consider the one-dimensional Schrodinger equation $$-\frac{1}{2}D^2 \psi(x)+V(x)\psi(x)=E\psi(x)$$ where $D^2=\dfrac{d^2}{dx^2},V(x)=-\dfrac{1}{|x|}$. I want to calculate the ground state ...
6
votes
4answers
4k views

Why is the canonical ($NVT$) ensemble often used for (classical) molecular dynamics (MD) simulations?

Molecular dynamics (MD) simulation is a common approach to the (classical) many-body problem. It relies on integration of Newton's equations of motion to simulate the trajectories of many (e.g., ...
5
votes
1answer
110 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
246 views

Discretizing the Wave Equation in polar coordinates

I want to discretize the wave equation $$\frac{1}{c^2}\frac{\partial^2\psi\left(\vec{r},t\right)}{\partial t^2}=\triangle\psi\left(\vec{r},t\right)$$ in polar coordinates. I find the following ...
5
votes
2answers
375 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. ...
7
votes
1answer
133 views

Minimal Extension of Wave Equation to Include Dispersion

Let's say you are modeling some process with the wave equation $\frac{1}{c^{2}}\frac{\partial^{2}\psi}{\partial t^{2}} = \nabla^{2}\psi$. You wish to improve your model by including dispersive ...
3
votes
1answer
93 views

Is there an algorithm for N body simulations in General Relativity [duplicate]

I am new to general relativity but have a background in computer science. Why is it so hard to do n-body simulations in GR? For example, there could be a program which takes the properties (mass, ...
1
vote
1answer
619 views

Calculating Kramers-Kronig using Mathematica [closed]

First of all, I know there is a Mathematica group in beta, but I don't think the problem of the following is directly a Mathematica issue. I am trying to calculate the change of the refractive index ...
7
votes
1answer
410 views

Numerical schemes, time integration algorithms and energy conservation

What does it mean when someone says a numerical scheme or a time integration algorithm is "energy conserving". How can a numerical scheme "gain" or "lose" or "conserve" energy apart from the numerical ...
5
votes
1answer
130 views

How close to the critical point is sufficient close for measuring critical exponents?

I am learning Monte Carlo and just manage to simulate a phase transition by computing the heat capacity or the susceptibility. I wish I can also compute critical exponents.To this purpose, I have read ...
4
votes
1answer
325 views

How to solve the heat equation for compound materials with different heat conductivities numerically?

I'm solving the heat equation with time dependent boundary conditions numerically in a 2D system using the ADI scheme. For the purpose of this question, let's assume a constant heat conductivity and ...
3
votes
1answer
89 views

Calculating a two-dimensional orbital path with infinite granularity (non-Euler integration)

For a game I am making, I am trying to calculate the position of an orbiting object around one or more bodies. I have successfully implemented this gravity simulation by calculating the force, then ...
3
votes
0answers
145 views

Difference between a Fixed Point and a Limit Point in implementations of the Renormalization Group (RNG) in Large Eddy Simulation (LES) model

In the introduction of this paper, it is explained that and how the application of a dynamic subrid scale model for turbulence into a large eddy simulation (LES) model corresponds to doing one ...
2
votes
0answers
96 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
2answers
812 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
134 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
384 views

Hysteresis Curve and how to implement it using Preisach model or other models

As a homework I need to draw a hysteresis curve (preferably an interactive one) using Matlab or any other programming language. The problem is I have trouble finding a good algorithm to do so. I need ...
1
vote
1answer
75 views

Getting nonphysical results when solving for the index of refraction of a slab?

I'm trying to computationally find the refractive indices (real and imaginary) for a thin slab suspended in air (so the only indices to deal with are air and my material's). I've experimentally taken ...
1
vote
2answers
168 views

What is a good model for computing water dropping on a surface?

Before introducing my question I would like to outline the fact that I'm a coder, so I can be wrong when using some kind of terminology. What is a good model for computing the flow, the shape and all ...
1
vote
3answers
920 views

Calculate kinematics of body movement from the set of spatial coordinates

Short intro I have a set of 3D (x,y,z) spatial coordinates of arm movement obtained using motion capture system. The example set of such coordinates looks like this (rounded up): ...
-2
votes
1answer
140 views

Can you suggest me some good resource for learning Mathematica for physics? [closed]

I am into HEP. I know nothing in Mathematica and have to start from basics. HEP requires a lot of computation for Feynman diagrams etc could you please suggest some good resource to start learning ...