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)

1
vote
2answers
212 views

Incompressible Navier-Stokes pressure solve in simulations

I am a complete newcomer when it comes to fluid simulations. I'm currently working through some tutorials to understand the idea of of the discretized Navier-Stokes equations for numerical ...
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$ ...
1
vote
1answer
79 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
1answer
207 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 ...
1
vote
2answers
137 views

Discrete approximation of charge density

Given the electric potential $\Phi(r)$ and the Poisson's equation: $$ \nabla^2 \Phi(r) = - 4\pi \rho(r)$$ Consider the 2-dimensional case and let's say that I want to discretize this using a square ...
1
vote
1answer
153 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 ...
1
vote
1answer
135 views

reversible cellular automata

Let's suppose a cellular automaton has a value $b(r,t)$ belongs to $Q$ at site $r$ and time $t$, where $Q$ is the set of possible states at each site. Let $N(r, t)$ be the values of the states of all ...
1
vote
2answers
403 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 ...
1
vote
1answer
448 views

Does the universe follow a general physical and computational principle of evolution?

Evolution is a principle in biology, whereby organisms evolve their ability to replicate and multiply in numbers over successive generations. From a computational point of view, the organisms employ a ...
1
vote
1answer
646 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 ...
1
vote
1answer
57 views

Conversion of satellite coordinates from ITRF to J2000

I have coordinates of various satellites in two coordinate systems: Cartesian coordinates in the international terrestrial reference frame (ITRF) RA / DEC in J2000 epoch, as derived from plate ...
1
vote
1answer
40 views

Single fluid simulation algorithm

Every fluid simulation I've seen implicitly involves a fluid that isn't mentioned, usually air. This is clear because the simulation usually has an ambient temperature, and always an ambient velocity. ...
1
vote
1answer
27 views

Why will a Fourier transform of a pickup placed at the middle of a string only result in odd harmonics?

My physics teacher explained a lot of information about strings but he never went in depth with my question. What happens to the even harmonics and why are the odd ones the only ones picked up?
1
vote
1answer
78 views

How do I calculate the effect of a prograde, retrograde, radial and anti-radial burn on the orbital elements of a two-dimensional orbit?

My relevant knowledge: This question works in two dimensions, with a two-body problem (see bottom for context). To my understanding, a two-dimensional orbit has the following orbital elements (3D ...
1
vote
1answer
29 views

Expression for spontaneous emission in solids

I have yet to find a reliable source for the expression of the rate of spontaneous emission in a solid. Can anyone confirm if the following is correct? The basic ingredients of the calculation are ...
1
vote
1answer
72 views

What is the velocity of Sun due to Earth alone?

If Earth and Sun were in a isolated system, will the Earth's motion around Sun will be similar? What will be Sun's and Earth's velocity when Earth is at its aphelion? Please note that it's not a ...
1
vote
1answer
66 views

Formulating a symplectic integrator for a non-local Hamiltonian

I recently asked two questions, Q. [1] and Q. [2], regarding reformulating non-local Lagrangians as Hamiltonians. In these questions, the Hamiltonian is formulated as an integral because of it's ...
1
vote
1answer
265 views

Can I use Runge-Kutta to solve these equations?

Edit: I'm going to give some more background and derivation to show how I got to these equations. I am basically following the derivation that is found in the appendix of the following paper: ...
1
vote
1answer
120 views

Heat equation with heat radiation and heat transfer

If I want to calculate steady temperature distribution on a one-dimensional stick, and I need to consider both the heat radiation and heat transfer, then my equation will be in the form: $$ ...
1
vote
1answer
139 views

Thermal conduction through 1D rod with changing thermal diffusivity?

I've been working with solving partial differential equations in 1D numerically using the finite difference scheme. My case is for example this: I've got a rod that's insulated around the sides such ...
1
vote
1answer
672 views

Why is the ground state energy of the Heisenberg XXZ Model unbounded for some values of $J$?

At the moment, I'm looking at numerically studying the Heisenberg XXZ model. The Hamiltonian is given below: $$ H=\sum_{j=1}^{N-1}\left(J S_j^z ...
1
vote
1answer
178 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 ...
1
vote
1answer
223 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 ...
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. ...
1
vote
1answer
154 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
1answer
117 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 ...
1
vote
1answer
95 views

what is two species Bose Hubbard Model? what chage it takes on Fock basis representation?

The Hamiltonian will look like \begin{equation} \ H= -t^a\sum_{<i,j>}^{}(a_i^\dagger a_i + h.c.)-t^b\sum_{<i,j>}^{}(b_i^\dagger b_i + ...
1
vote
1answer
759 views

Magnetic field of permanent magnet [closed]

Hi I am looking for some code to (computationally) calculate the magnetic field of a rectangular permanent magnet in 3D space. Preferably I could edit it to add in my own dimensions. Does anyone know ...
1
vote
1answer
136 views

Parameter determining argon phase

Currently I am working a molecular simulation to determine phases of an argon NPT ensemble using Lennard Jones potential. Mainly I use the radial distribution function to determine solid, liquid, or ...
1
vote
1answer
353 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
322 views

Where does Computer Science background students fit in Theoretical Physics [closed]

I am basically an Electronics student - background in computer science (that's where I want to work). I applied for an internship in USA in a research institute where the group is focused in ...
1
vote
1answer
314 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 ...
1
vote
1answer
185 views

What software can I use to simulate supersonic gas particle flow? [closed]

Please please help me out with this. I am trying to find a software/program that I could use the simulate the scattering effect of light when it strikes gas particles moving at supersonic speeds? ...
1
vote
0answers
20 views

Vortex in a Pipe / Solutions or approximations?

I want to study a vortex flow in a pipe; in other words helical flow. A couple introductory texts on fluid mechanics (e.g. Shames, Mechanics of Fluids) describe the solution to vortex flow in 2D: ...
1
vote
0answers
16 views

Efficient ways to solve complicated kinematics problems

I'm curious about the range of approaches one can take to perform motion/kinematic simulations on many 3D objects. Currently I am familiar with the iterative process of doing such a simulation, that ...
1
vote
0answers
20 views

Monte Carlo simulation of Breit-Wigner

I am using a Monte Carlo approach to integrate a scattering amplitude which has a propagator-squared term like $f(s)=\frac{1}{(s - 4)^2 + 1}$, between $s=0$ and $s=20$. Using a uniform random ...
1
vote
0answers
44 views

About multi dimensional numerical integration [closed]

I want to calculate 4 dimensional phase space integration.If anybody help me with a general 4 dimensional integration that will be OK. I can apply trapezoidal or Simpson's rule in lower dimensional ...
1
vote
0answers
36 views

Model of the solar system orbital periods [closed]

I've written a program to model the solar system. I am now testing its accuracy. I've checked conservation of angular momentum, and it is perfectly conserved. I've also plotted the path of the ...
1
vote
0answers
20 views

Package for flavour physics computation? [closed]

Which softwares one usually uses for flavour physics related calculations? I saw one called SUSY_FLAVOR. I haven't used it yet though. What in case of other BSM models other than SUSY? Is there any ...
1
vote
0answers
47 views

Software engineer vs physics engineer vs computational physics [closed]

Don't know if this question is off topic, but i think it is the best place to get an answer. I am finishing the 3rd year of a degree in physics engineering. I enjoyed the first year, but in the ...
1
vote
1answer
22 views

How to choose basis functions that contribute most efficiently per term? [closed]

I would like to approximate some positive, scalar function $f(x,y) > 0$ on a 2D field of finite size i.e. $x\in [a,b],y\in[c,d]$ I am familiar with the set of basis functions used in the Fourier ...
1
vote
0answers
67 views

Experiments to resolve dillema between continuity and discrete

Which experiments/experimental methods are suggested to resolve an alternative about the structure of our universe space and time - is it continuous or is it discrete in a very small scale, especially ...
1
vote
1answer
58 views

Modeling transient heat transfer between two 1-D materials

So I have been working on creating a 1-D model of a material with a layer of heat resistant material stuck to it. For this I have been trying to use the implicit finite-differential method. The ...
1
vote
0answers
17 views

Coupled Diff Equation from Bose Einstein distribution [closed]

I am a student doing physics hons and have had very little experience in programming. This semester we are supposed to do a computational project in thermodynamics. I have to solve these two coupled ...
1
vote
0answers
71 views

Obtaining propagating solutions for Schrodinger equation from known bound states (in 2 and 3 dimensions)?

If I found all the bound states for a certaing potential in 2 or 3 dimensions (numerically), can I immediately obtain some information about the propagating solutions for the same potential (such as ...
1
vote
1answer
26 views

Adding fluctuations to a hydrodynamic simulation to trigger instabilities

I'm doing a 2D combustion hydrodynamic simulation and there's a hydrodynamic instability that should be triggered because of the particular physical properties of my system. The key to the instability ...
1
vote
0answers
25 views

Minimal possible energy to make calculation [duplicate]

Are there any restrictions to energy that must be consumed to make calculations? Can calculations theoretically be done without consuming it?
1
vote
0answers
34 views

Computational method for finding edge states?

I am actually interested to learn how to calculate edge states in 1D topological systems using computational methods, Q. can anyone tell me which method is best suited and easy to calculate edge ...
1
vote
0answers
99 views

What are the differences between the Jetphox, Pythia and Herwig event generators?

I know Jetphox is a parton-level event NLO generator program. But I want to know more about other generator programs such as Pythia and Herwig. What are the differences? I am undergraduate student so ...
1
vote
0answers
94 views

Selecting physical solutions in numerical eigenvalue problems

I try to solve a certain time-independent Schrodinger equation numerically, using the method of finite differences. My boundary conditions are such that the finite difference method gives me an ...