Questions tagged [numerics]

DO NOT USE THIS TAG just because the question contains numerics!

1
vote
0answers
621 views

How to calculate spacing in k-space after DFT of spatial signal

Note: In the following i'm going to use t as a time variable [s] and x as a spatial variable [m]. I'm currently working with the Discrete Fourier Transform (DFT), in order to get frequency ...
2
votes
1answer
196 views

rounding when doing calculation in physics

I came across a calculation for projectile motion in my physics book as follow: $$r = \sqrt{x^{2} + y^{2}} = \sqrt{(4.5\,m)^{2} + (-1.2\,m)^{2}} = 4.7\,m$$ If I round the results in (1), (2) and (3),...
1
vote
1answer
227 views

Uncertainties of coefficients in linear fit [closed]

I'm trying to calculate the uncertainty of the coefficients in a linear fit, but I'm not sure how to do it. on the one hand, I can extract the error with matlab statistical tools but on the other hand ...
1
vote
0answers
109 views

Operator splitting for Maxwell's equations

In this book I've read about operator splitting method for solving differential equations. I for example I have an equation in the form $$\dot \omega (t) = A \omega(t)$$, where $A$ is some operator ...
1
vote
0answers
54 views

Numerical error in implementation of universal variables [closed]

I am trying to implement (in Python for now) low thrust orbit propagation for spacecraft using universal variables. For a given central body with the gravitational parameter $\mu$ and an orbit with ...
1
vote
0answers
47 views

Capturing superfluid condensation with exact diagonalization

Doing exact diagonalization on bosonic systems is tricky, because the possibility of multiple occupancy means that even the single-site Hilbert space is infinite-dimensional. It's my understanding ...
4
votes
1answer
103 views

How to conserve momentum when numerically-integrating the path of a charged particle through a magnetic field

I have the python script attached below, which is intended to track the trajectory of a charged particle in a static, uniform magnetic field. It is very simple to calculate the instantaneous force ...
0
votes
1answer
991 views

Which procedure is correct? [closed]

A problem is given in my textbook pg.no-191 as Example 5.10 A solenoid has a core of a material with relative permeability $\mu_r=400$. The windings of the solenoid are insulated from the core and ...
4
votes
1answer
502 views

Geodesics in Kerr

I'm interested in plotting the trajectories of null geodesics near an uncharged rotating black hole (described by the Kerr solution) which involves a system of first order differential equations. Kerr ...
1
vote
1answer
112 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
124 views

Evaluate singularity in numerical iterative physical optics integral

My goal is to implement a code to numerically compute the radar cross section (RCS) of an open-ended PEC cavity by the approach of "iterative physical optics" [1]. The idea is to iteratively compute ...
7
votes
1answer
2k 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, ...
6
votes
1answer
1k 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 ...
0
votes
0answers
518 views

Numerical solution of the BCS gap equation, with Coulomb potential

I'm interested in the excitons condensation, which is described by an equation which is very similar to the standard BCS gap equation. I will be referring to this article: http://arxiv.org/abs/cond-...
49
votes
2answers
2k views

Can we infer the existence of periodic solutions to the three-body problem from numerical evidence?

I recently found out about the discovery of 13 beautiful periodic solutions to the three-body problem, described in the paper Three Classes of Newtonian Three-Body Planar Periodic Orbits. Milovan ...
2
votes
1answer
218 views

Simulating a 1-dimensional wave on (a segment of) an infinite line

I'm trying to numerically simulate a 1-dimensional with a chain of linked harmonic oscillators as described here (the result can be seen here). The simulation behaves like a wave on a finite line ...
7
votes
1answer
1k 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 ...
0
votes
0answers
134 views

numerical diagonalisation comparison of complex Hermition matrix and real matrix

we are doing K space calculation(transnational invariant basis) of Hamiltonian, and trying to compare with exact method for same nos of sites and particles, we are getting big difference in ground ...
3
votes
0answers
199 views

Monte Carlo for Random Bond Ising ferromagnet

The set-up: Consider the Ising model on an $L \times L$ square lattice, where the coupling of each bond is chosen to be $+J$ (ferromagnetic) with probability $(1-p)$ and $-J$ (antiferromagnetic) with ...
2
votes
1answer
255 views

Electric current streamlines in induction cooking vessel

I am looking for a plot of the typical streamlines of the electric induced currents ("eddy currents") in a induction cooking vessel. How can one theoretically predict the streamlines? How is it ...
5
votes
2answers
534 views

Solve the angular part of Schrodinger equation numerically

I would like to solve the angular part (the one for what is usually called the $\theta$ angle) of a time-independent 3D Schrodinger equation $$ \frac{\mathrm{d}}{\mathrm{d}x}\left[ (1-x^2) \frac{\...
9
votes
2answers
9k views

Superconducting gap, temperature dependence: how to calculate this integral?

Tinkham (page 63) states that the temperature dependence of the gap energy of a superconductor $\Delta(T)$ can be calculated using the following integral: How can this actually be carried out? I am ...
4
votes
1answer
152 views

neutron transport approximations for nuclear rocket modelling

I'm pretty ignorant regarding neutron and nuclear transport modelling, but i'm interested in trying to pursue it for a particular pet project. It regards modelling of nuclear reactions like those ...
2
votes
2answers
324 views

numerical formulation of Dirac equation plus electromagnetic field

I have the following equations describing the electron field in a (classic) electromagnetic field: $$ c\left(\alpha _i\right.{\cdot (P - q(A + A_b) + \beta mc) \psi = E \psi } $$ where $A_b$ is ...
3
votes
3answers
546 views

Can a wavefunction be solved to any arbitrary precision, given enough computer time?

I learned that the wavefunction for the hydrogen atom can be solved analytically (we did the derivation in class), but that for more complicated atoms it is "impossible" to solve and that only ...
-3
votes
3answers
147 views

How to recognize broken candies from whole ones [closed]

Let's say I have a bag full of sugar candy. Some will be whole, some will be dent, some will be broken (in part, or half, etc). Let's say I have a device with an input box where I empty the bag, and ...
7
votes
1answer
912 views

Interpretation of Stiffness Matrix and Mass Matrix in Finite Element Method

I would like to have a general interpretation of the coefficients of the stiffness matrix that appears in FEM. For instance if we are solving a linear elasticity problem and we modelize the relation ...
1
vote
1answer
984 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
419 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 ...
12
votes
2answers
132 views

Numerical Analysis of Elliptic PDEs

I am looking for an elementary reference regarding issues of stability in numerical analysis of non-linear elliptic PDEs, particularly using the finite difference method (but something more ...
2
votes
1answer
959 views

Transparent boundary condition [closed]

I am interested in the finite-difference beam propagation method and its applications. I try to solve the Helmholtz equation. At first, i would like to solve numerically it for the easiest case, ...
5
votes
1answer
914 views

Numeric method to calculate the charge distribution on a conducting surface?

If I have an arbitrary (closed?) conducting surface and a nearby charge density, is there a simple numeric way of computing the induced charge distribution on the surface?
4
votes
3answers
1k views

Numerical software to manipulate a light beam in its plane wave representation?

Any light field can be expressed as a sum of plane waves. Such an ensemble of plane waves is called the plane wave spectrum of the light field. The plane wave spectrum is the Fourier transform of the ...
15
votes
1answer
976 views

Why isn't the Gear predictor-corrector algorithm for integration of the equations of motion symplectic?

Okumura et al., J. Chem. Phys. 2007 states that the Gear predictor-corrector integration scheme, used in particular in some molecular dynamics packages for the dynamics of rigid bodies using ...
1
vote
1answer
1k views

Physics needed to build a top down billiards game [duplicate]

Possible Duplicate: How are these balls reflected after they hit each other? I was wondering what sort of physics equations would I need in order to build a top down billiards game? I tried using ...
5
votes
1answer
336 views

Smoothed particle hydrodynamics in cosmological N-body simulations

What is the role of smoothed particle hydrodynamics (SPH) in cosmological N-body simulations like the Millenium Run (performed with Gadget-2)?
17
votes
1answer
7k 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
0answers
78 views

Degenerated Anderson Model Simulation

I'm trying to simulate the degenerative Anderson model. So depending on an energy difference first orbital and afterwards spin magnetism occurs. First i try to solve an easier ansatz with a limitation ...
2
votes
2answers
511 views

History of the use of the concept of phase space in engineering

Engineering textbooks constantly use the concept of 'phase space' (see e.g. http://www.cs.cmu.edu/~baraff/sigcourse/notesc.pdf). That is, they think of the state of a mechanical system as a high-...
2
votes
1answer
2k views

Using Fourier Transforms to Solve Systems with springs of high frequency

I'm trying to numerically solve the differential equations of motion in a system with multiple springs of very high frequency. Because the solution is often a combination of rapidly-oscillating sine ...
12
votes
5answers
2k views

Basic mechanics problems, unsolvable by brute-force numerical integration

I'm looking for simple problems in theoretical mechanics that are impossible or unreasonably difficult to solve by means of "brute-force" numerical integration of Newton or Euler-lagrange equations. ...
14
votes
7answers
1k views

Binary Black Hole Solution of General Relativity?

This is rather a technical question for experts in General Relativity. An accessible link would be an accepable answer, although any additional discussion is welcome. GR has well known solutions ...
5
votes
1answer
705 views

Mathematica to help for an Hamiltonian problem

I have an Hamiltonian problem whose 2D phase space exhibit islands of stability (elliptic fixed points). I can calculate the area of these islands in some cases, but for other cases I would like to ...