The tag has no wiki summary.

learn more… | top users | synonyms

2
votes
0answers
47 views

Precession of Mercury (Python simulation)

I was trying to simulate the precession of Mercury based on the perturbed solution, and my questions about its implementation in python can be seen here: ...
4
votes
5answers
269 views

Resolving General relativity and Newtonian mechanics to a computer [on hold]

I know this is considered an old subject long ridiculed by many as the folly of layman. But I work in the field of computer simulation, specifically in producing fully functional 3D interactive ...
0
votes
0answers
17 views

Graphene plasmonics, Integro-differential eigenvalue equation [closed]

In my research around plasmonics in 2-Dimensional nanostructures like graphene, I encounter an eigenvalue integro-differential equation of the form: ...
1
vote
0answers
35 views

Is there a numerical calculation for black hole - neutron star merges?

I think it might be even a more frequent event as the black hole merges. It would be similar to the black hole merges in the gravitational wave spectrum, but it could have a very clear neutrino and ...
1
vote
1answer
30 views

Meaning of imaginary eigenvalues in advection equation

I'm reading the book Fluid Simulation for Computer Graphics (Amazon link), and stumped by the following sentence in Chapter 3: what's happening is that the eigenvalues of the Jacobian generated ...
0
votes
2answers
61 views

What is the fastest algorithm to solve a Nonlinear second order differential equation numerically

I am trying to solve a second order non linear differential equation in one variable. Using RK4 I am getting good accuracy and is working fine. But the problem is my range is very high so it will ...
3
votes
2answers
56 views

Numerical modelling of a step function in time in a hydrodynamic system. (Runge Kutta fourth order)

So I'm trying to model a hydrodynamic system that introduces a sudden "jump" in the value of a function at a specific time. The system is solved with a Runge-Kutta fourth order method. I have a ...
1
vote
0answers
15 views

Modeling simple laser induced population transfer via adiabatic passage in python [migrated]

I'm trying to model adiabatic passage between two levels in a three-level atom interacting with two laser fields using Scipy and Numpy.. I'm not sure if my model is wrong due to my incorrectly ...
3
votes
1answer
46 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, ...
0
votes
0answers
33 views

Derivative of time with respect to proper time

When wants to solve the Schwarzschild-Two-Body-Problem with the Runge-Kutta-Method, the second derivative of the time $t$ with respect to the proper time of the moving particle $\tau$ is needed. How ...
0
votes
1answer
73 views

How to solve highly oscillating differential equation [closed]

The equation looks like: $$x''(t)+bx'(t)+c x(t)+dx^3(t)=0.$$ This is the motion of a particle in a potential $cx^2/2+dx^4/4$ with friction force $bx'$. In my case, the friction term is very small and ...
0
votes
0answers
93 views

How to solve the Dirac equation numerically?

The effective Hamiltonian for my system is: $$ H=\nu_{F} \sigma\cdot\left(q-By\hat x\right) $$ where $\sigma$ and $q$ are the Pauli matrices and the momentum operator respectively and $\nu_{F}$ and ...
1
vote
1answer
52 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 ...
0
votes
0answers
27 views

Electromagnetic boundary conditions for modelling symmetrical geometry

I stumbled upon this article: http://www.comsol.com/blogs/exploiting-symmetry-simplify-magnetic-field-modeling/ Since the article does not contain any mathematical formulations, I was wondering how ...
0
votes
1answer
53 views

Using the Fourier transform to find the natural frequencies of coupled oscillators

How can I find the natural frequencies of a system consisting of a pair of coupled oscillators using Fourier transforms? The System consists of two masses and three springs. One of the springs ...
3
votes
2answers
497 views

Numerical solution to Schrödinger equation - eigenvalues

This is my first question on here. I'm trying to numerically solve the Schrödinger equation for the Woods-Saxon Potential and find the energy eigenvalues and eigenfunctions but I am confused about how ...
1
vote
1answer
33 views

Numerical construction of phase space for a dynamical system

Suppose I have a standard, deterministic dynamical system. For concreteness I'll assume it's a two variable system of the form, $$ \dot x_1 = f(x_1,x_2; \theta_1)\\ \dot x_2 = g(x_1,x_2; \theta_2) $$ ...
0
votes
1answer
25 views

Is it realistic for soundwaves under water to “sink” or “float”?

I'm studying soundwaves under water and I had a numerical problem that I was asking about. http://stackoverflow.com/questions/28904017/are-my-matlab-iterations-working Now I wonder if you can tell ...
1
vote
0answers
55 views

Numerical Solution of 1D Boltzmann Transport Equation

I need to solve the one-dimensional Boltzmann transport equation in a semiconductor numerically, and I want to take a deterministic approach toward the problem (i.e. not use Monte-Carlo or similar ...
2
votes
1answer
58 views

Estimating divergence of set of vectors

I have a set of points where directions and intensities of a flow are given (in 3D). Is it possible to estimate the divergence of the flow defined by those vectors? I only need a rough estimate and I ...
0
votes
0answers
32 views

What does invariant exactly mean and how does it get the invariant?

I have read many journal about simulation of regularized long wave. In numerical test section,many researcher use invariant of mass,momentum and energy to check accuracy of their method but i found ...
2
votes
2answers
124 views

Hamiltonian mechanics really useful for numerical integration? Lagrangian can become 1st-order

(I'm talking about the classical mechanics.) Many texts say that Euler-Lagrange equations are difficult to treat numerically because they are second-order ODEs, ${f_i(\boldsymbol{q, \dot{q}, ...
1
vote
1answer
89 views

Hartree Fock equations

I don't understand how the Hartree Fock equations define an iterative method! For this discussion, I am referring to the HF equations as described here: click me! Basically if you guess a bunch of ...
0
votes
0answers
18 views

What is the transfer function in fft beam propagation for unpolarized light?

What is the transfer function in fft beam propagation for unpolarized light ? How to construct the fft beam propagation ? This is for homework. For coherent light the beam propagation is E(x,z) ...
0
votes
0answers
39 views

Why are functional representations of systems important in physics or computational physics?

This was an addendum to a previous question I asked, but I figured I should make it it's own discussion. Assuming I am able derive a functional representation for any dynamical system (dissipative, ...
0
votes
1answer
84 views

Leapfrog method in Particle-in-cell

Recently I wanted to write a 3D electromagnetic Particle-in-cell code with c++.I know that I should use leapfrog method.For example,when I calculates the position and velocity of particles,i should ...
0
votes
0answers
13 views

What is the phase-amplitude numerical method?

What is the phase-amplitude numerical method? I heard its used to calculate long range interactions numerically, but I cannot find any papers discussing its method of implementation.
0
votes
0answers
69 views

Stability of Mathieu's equation and parameteric resonance

I am given the following equation (Mathieu's equation) in my subject of Numerical Analysis : $$ \frac{d^2 x}{dt^2}=-\omega^2(1+\epsilon\cos(t))x $$ I am supposed to find those frequencies $\omega$ ...
1
vote
0answers
177 views

Problems while numerically computing band structure using k.p theory

I want to use k.p theory to numerically compute the band structure of a bulk semiconductor. The band I like to include are the lowest conduction band (cb), the heavy-hole (hh), the light-hole (lh) and ...
0
votes
0answers
35 views

Differential Equations for Two-Dimensional N-Body Simulation

So I recently asked a question about improving the stability of a two-dimensional orbital simulation (as before I was using Euler integration). I was told to use Runge-Kutta 4 for a more stable ...
3
votes
1answer
71 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 ...
4
votes
1answer
207 views

Integral of Sedov's self-similar solution to the spherical blast wave problem

I'm studying the Taylor-Sedov self-similar solution to the problem of a strong explosion in a homogenoeus atmosphere. The problem is discussed in Landau & Lifschitz VI (in the 2nd edition it's ...
0
votes
0answers
16 views

How to conserve power at the reflection boundary of perfect conductor?

How to conserve energy during the reflection ? I am taking the abs value of $\sum E_{field}$ along spacial direction and plotting at each time step. Then it is normalized to the input E filed. ...
1
vote
1answer
136 views

Numerical solution for Friedmann equations

My problem today is to solve the Friedmann equations, for those who aren't familiar with them, here they are (in my specific case): $$ \left ( \frac{\dot{a}}{a^2} \right )^2 = \frac{\rho_1}{a^4} - ...
1
vote
1answer
204 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: ...
0
votes
0answers
71 views

Boundary conditions of stream function

I have to do an problem about solving numerically the flow that goes under an airfoil. The airfoil has a flap deployed downwards and I need to solve the mesh that it's under the airfoil. I have drawn ...
3
votes
3answers
527 views

What's the difference between “numerical methods” & “mathematical analysis” as said by Feynman in his lectures?

While reading his lectures, I came to these lines: On the basis of Newton's second law of motion,which gives the relation between the acceleration of any body & the force acting on it,any ...
1
vote
1answer
73 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: $$ ...
7
votes
2answers
237 views

Arrhenius Fit: Linear or exponential form?

I have a seemingly easy question about performing an Arrhenius fit to the equation $$y = A \times \exp \left( -\frac{E_A}{RT} \right)$$ I can either fit this in the exponential form using a ...
0
votes
0answers
20 views

Ising Monte-Carlo and Three point functions

I'm looking for literature on the calculation of three points function in the 2d Ising Model using numerical methods, especially around the critical point. By $Z_2$ symmetry, three spin insertions is ...
5
votes
2answers
82 views

Is it possible to propagate a relativistic system of particles in time using Verlet?

The Verlet algorithm and its derivations are very popular methods to integrate Newton's equations of motion in time and obtain a trajectory for a system with $N$ particles. I work with classical ...
4
votes
1answer
444 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
57 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 ...
4
votes
1answer
171 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 ...
0
votes
0answers
25 views

Adams-Moulton and BDF methods

1.What are the differences between Adams-Moulton and BDF methods. Which one is better and which one computes the solution faster? I think Adams-Moulton is a better method as it can get to the ...
1
vote
1answer
106 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
626 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 ...
3
votes
0answers
66 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
2answers
174 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 ...
1
vote
1answer
100 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 ...