The tag has no wiki summary.

learn more… | top users | synonyms

0
votes
1answer
45 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
1answer
256 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
25 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
36 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
50 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
28 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
107 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
81 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
15 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
37 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
53 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
54 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
159 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
27 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
57 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
156 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
14 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
96 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
175 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
57 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
437 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
65 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
223 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
18 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
67 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
435 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
55 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
153 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
23 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
95 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
479 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
61 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
164 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
96 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
82 views

How to numerically solve a complex equation? [closed]

I want to know that if you are given a very complex equation g(x)=A(T). How could you solve for x, which is a function of variable T. To be more specific, I encounter a polylogarithmic function I need ...
3
votes
2answers
191 views

Boundary Element Method or Boundary Integral Method Computational Aspects

I have to solve a Helmholtz equation inside a simply connected domain. I know that in general the boundary integral can be written as, $$\phi(x)=\int_V G(x,x') \rho(x')\ d^3x'+\int_S ...
0
votes
0answers
44 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 ...
1
vote
0answers
86 views

Lattice Gas Cellular Automata - HPP model square lattice

In the HPP model of LGCA, a square lattice is used and there is only one collision configuration as mentioned in figure (taken from the book Lattice Gas Cellular Automata and Lattice Boltzmann models ...
2
votes
1answer
150 views

Simulation of fluid flow using Euler equation

I have been looking on Euler's equations for a while and can't grasp one thing. Suppose we have initial system state with volumes of fluid "hanging" in air (time is frozen and equal to zero), each of ...
2
votes
0answers
84 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
111 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
2answers
72 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
2answers
102 views

Energy of damped harmonic oscillator begins to increase with very large Q in numerical integration

I have numerically integrated the (reduced) homogeneous equation of a damped harmonic oscillator in order to see how the error propagates. $$\frac{d^2 X}{d\phi^2} + ...
0
votes
1answer
54 views

Numerical Computation of Linbald Equation [closed]

Can anybody suggest me a good algorithm for the time evolution of the reduced density matrix using Linbald equation. My Hamiltonian is time dependent. I am aware about Qotoolbox and Qutip. I have ...
0
votes
2answers
254 views

Numerical solving of 2D and 3D Schrodinger equations

I am studying 2D quantum scattering models for my Bachelor's thesis. Somewhat like these: ,with Dirichlet ($\psi \mid_\Gamma = 0$) boundary conditions on the "walls" of the waveguide and the ...
0
votes
0answers
91 views

Difference between normal mode methods and wavenumber integration?

Normal mode and wavenumber integration methods allow the evaluation of an integral transform solution. The text book that I am reading states that the normal mode method evaluates the field as a sum ...
5
votes
1answer
166 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
28 views

Order of Monte Carlo integration and frequency summation

I am currently trying to calculate an integration formula of a linear response function by Monte Carlo method. It is a multiple integration over three 3D vectors, i.e., nine dimensions in all. And ...