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

3
votes
0answers
86 views

Nature as a calculator [on hold]

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
26 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 ...
2
votes
0answers
55 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
46 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
13 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
0answers
42 views

Which are the good universities for doing PHD in computational physics? [closed]

Currently i am doing MSC(physics) from Pune university, India. And I want to do PHD in computational physics as it interests me the most. So I want to know which universities will be a good option for ...
1
vote
1answer
44 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. ...
1
vote
2answers
61 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$ ...
0
votes
0answers
23 views

Computational package to find the ground state of a particle in 3D domain

I am developing a numerical algorithm to find the ground state of a Hermitian matrix. Obvious applications are quantum many-body systems and particles in various potentials. I am a little stuck with ...
3
votes
1answer
54 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
49 views

Euler equation with single state variables

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
56 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 ...
8
votes
2answers
105 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
54 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 ...
4
votes
1answer
133 views

Numerical analytic continuation for Green's function

Recently, I happened to hear about the possibility of doing analytic continuation numerically. That sounds attractive for the ubiquitous $\mathrm{i}\omega_n\rightarrow\omega+\mathrm{i}0^+$ procedure, ...
0
votes
0answers
21 views

How to easily calculate lengths (or relative lengths) of paths between symmetry points in BZ

I am trying to easily calculate the length between special kpoints within the BZ of the 32 point groups in a crystal system. I am calculating the lengths in order to scale k point sampling along these ...
7
votes
3answers
171 views

Runge Kutta Method for a Lindblad Equation

I am solving a Lindblad equation for a dissipative Harmonic Oscillator. My Hamiltonian is time dependent, My Lindblad Equation can be written as \begin{equation} ...
4
votes
1answer
66 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 ...
13
votes
3answers
405 views

What limitations are there in measuring physical properties accurately?

In a StackOverflow answer, I attempted to explain why a 32-bit float was perfectly adequate for representing the questioner's weight measurement: Physical properties are inaccurately measured ...
1
vote
1answer
68 views

I'm getting weird autocorrelations when simulating an Ising model below the critical temperature

So I'm simulating an Ising model using Monte Carlo and the Metropolis algorithm. After letting it reach equilibrium, I try to calculate the autocorrelation of the magnetization. As long as the system ...
6
votes
0answers
68 views

Feynman's Infinite Amount of Logic for One Tiny Bit of Space

Watching one of Feynman's lectures, I came across something that puzzled me. What was Feynman referring to when he said the following? What goes on in no matter how tiny a region of space and no ...
7
votes
2answers
210 views

Monte-Carlo and $O(n)$ models for non-integer n

$O(n)$ lattice statistical models can be generalized to non integer values of n, starting from their (expanded and resumed in graphs) partition function: $$Z = \sum_{\text{loop configurations}} n^{\# ...
3
votes
2answers
108 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 ...
0
votes
0answers
65 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 ...
3
votes
1answer
107 views

Epsilon Tensor in FeynCalc

A few days ago I started to use the Mathematica package FeynCalc and one thing confuses me: Assume we have a four-vector $p_\mu$ and we contract it with the epsilon tensor. FeynCalc produces ...
3
votes
0answers
76 views

What is the probability of quantum tunneling occurring in this CPU?

You may have noticed over the last few years that Moore's law is no longer applying to the real world. This observation states that over the history of computing hardware, the number of transistors on ...
2
votes
1answer
86 views

How much computer power would we need to model every physical aspect of the universe [closed]

So I know there are equations that give you the total energy in the Universe(Friedman). So about if we we to simulate every single interaction from neurons to physical forces, strong and weak, ...
3
votes
1answer
41 views

Why is the objective function used in Nudged Elastic Band method reasonable?

In Nudged Elastic Band (NEB) method, which is used to find reaction pathways when both initial and final states are know, an objective function is first constructed and then minimized to find reaction ...
0
votes
1answer
51 views

How to Extend Relaxation Methods for 2D Laplace Equation given in Jackson E&M to 3D?

In Jackson (3 ed) chapter 1.13 an outline is given for using relaxation to solve laplace equation in 2D. The general procedure in 3D involves minimizing the quantity $$I[\Psi]=\frac{1}{2}\int_V ...
3
votes
0answers
109 views

Comparison between Cadabra and other Symbolic Computer Algebra software [closed]

Does anyone has some experience about working with Cadabra and it's (dis)advantage in comparison to other Symbolic Computer Algebra software such as Maple and Mathematica (physics package) in the ...
1
vote
1answer
70 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 ...
1
vote
0answers
35 views

How to get discretization coefficients of matrix $A$ in Finite Volume Method (FVM)? [closed]

First we have Discretization of the Transport Equation $$ \frac{\partial \rho \phi}{\partial t} + \nabla(\rho U \phi) - \nabla (\rho \Gamma_\phi \nabla \phi) = S_\phi (\phi) $$ In Finite Volume ...
0
votes
1answer
66 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
1answer
58 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 ...
3
votes
0answers
374 views

Numerical problem in solving the Bogoliubov de Gennes equations- methods to solve?

I am trying to solve an assignment on solving the Bogoliubov de Gennes equations self-consistently in Matlab. BdG equations in 1-Dimension are as follows:- $$\left(\begin{array}{cc} ...
1
vote
0answers
52 views

Application of wavelets in computational electromagnetics

I am a first-year graduate student (electrical eng.); wavelets is one of the courses that I am taking presently. As part of the course, I plan to explore application of wavelets in computational ...
1
vote
1answer
37 views

2nd - 3rd Year Physics Labs [closed]

I'm looking to propose a new physics laboratory experiment, does anyone have any particularly fond memories of an experiment they carried out whilst at University? I'd like to include a section that ...
3
votes
1answer
151 views

confusion in discrete transform to solve kronig penney matrix equation in fourier space

I have a periodic potential $$V(x) =\sum_{K}e^{iKx}V_{K} =\sum_{n}e^{\iota2\pi nx/a}V_{n} $$ where $K =\frac{2\pi n}a$ is the reciprocal lattice vector and $a$ is the lattice constant and $n =\pm ...
1
vote
0answers
219 views

Numerical solution of two coupled second order differential equations of motion [closed]

Is there a numerical algorithm for solving a pair of coupled second order differential equations? This question arises from a homework problem that I have that involves two dimensional projectile ...
0
votes
0answers
26 views

What jobs can you get in the field of physics? [duplicate]

What good, fun and reasonably well paid jobs can you get in the field of physics?
1
vote
0answers
59 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 ...
6
votes
2answers
167 views

Eigenvalue problem for differential equations in QM

I have a very simple question with regard to numerical methods in physics. I want to solve the eigenvalue problem for a particle moving in an arbitrary potential. Let's take 1D to be concrete. I.e. I ...
1
vote
2answers
110 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
0answers
71 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
1answer
73 views

Motion of rigid body system in absense of work

In the absence of work on the system, is there a closed form equation for the motion of a set of constrained rigid bodies (let's say, using Revolute (ie: simple pivot) constraints)? If the bodies are ...
2
votes
2answers
56 views

Sampling a distribution (from a galaxy model)

I am reading the following article: http://www.kof.zcu.cz/st/dis/schwarzmeier/galaxy_models.html and am currently at section 5.6 (positions of bodies in a galaxy). I am trying to redo the simulations ...
1
vote
0answers
302 views

Is Bremermann's limit redundant with Landauer's limit for all practical purposes?

Trying to understand the physical limits to computation, I notice that among these we have two types of limits that constrain the minimum allowable energy for a computation. Limits that constrain ...
1
vote
1answer
108 views

Smoothed Particle Hydrodynamics, Chemical Reactions and Stiffness [closed]

Dear people from stackexchange physics. I have been using SPH for quite a while to simulate free boundary flow, and just recently we tried to include some kind of (simple) chemical reaction in our ...
4
votes
1answer
132 views

How to find positions of $n$ masses in Newton mechanics?

I ran into a problem while doing research. The problem can be described as: consider the original $n$-body problem, and if we fix the position of them(unknowns), no interaction among them, they don't ...
1
vote
0answers
47 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 ...