Tag Info

Hot answers tagged

15

Take a look at the notes on lectures 1 and 2 of Geometric Numerical Integration found here. Quoting from Lecture 2 A numerical one-step method $y_{n+1} = \Phi_h(y_n)$ is called symplectic if, when applied to a Hamiltonian system, the discrete flow $y \mapsto \Phi_h(y)$ is a symplectic transformation for all sufficiently small step sizes. From your ...


8

Crank-Nicholson method is effectively the average of forward (explicit) Euler $\psi(x,t+dt)=\psi(x,t) - i*H \psi(x,t)*dt$ and backward (implicit) Euler method $\psi(x,t+dt)=\psi(x,t) - i*H \psi(x,t+dt)*dt$ The backward component makes Crank-Nicholson method stable. The forward component makes it more accurate, but prone to oscillations. If you want to ...


6

There are no exact solutions, only approximations and numerical solutions. Don't forget that orbiting black holes will radiate gravitational waves so any solution would have to include those and the corresponding decay of the orbit until the black holes coalesce.


5

We have observed that the underlying level of nature is quantum mechanical. Quantum means "a definite quantity" of something so definite quantities can be counted and so integral numbers play a role : a) in the number of particles , in the number of energy levels characterized by quantum numbers ( i.e. integer numbers).b) There are the fields which are ...


5

"They" are probably talking about symplectic integrators. Most numerical integrators for (partial) differential equations do not specifically consider the energy of the system; they are generic integrators capable of solving any set of DEs, and not all DE's have a concept like "energy". When these are applied to a classical dynamics problem concerning ...


5

For a very recent authoritative review of the numerical approach, see Centrella et. al. http://arxiv.org/abs/1010.5260 For the alternate parameterized post Newtonian approach, see Living Reviews of Relativity http://relativity.livingreviews.org/Articles/subject.html and look for articles number 2007-2, 2006-4 and 2003-6.


4

There is rather nice function in Mathematica 7, which allows one to integrate over an arbitrary complicated region. It is Boole:[True,False]$\to${1,0}. Below is just an example taken from Mathematica Documentation Center. If you have a 2D area defined by the inequality $4 x^4-4 x^2+y^2\leq 0$, you can integrate any function $f(x,y)$ over this domain as ...


3

I will not give you the numerical solution, but I will below explain some analytical simplifications that I believe are required to solve the numerical problem. The strategy is simple: try to express all the parameters of the integral in term of dimensionless variables. To achieve a discussion in term of $\delta = \Delta(T) / \Delta(T=0)$ and $\tau = T / ...


3

It seems like they were able to rigorously prove the existence of N-body choreographies by using interval Krawczyk method to show that a minimum exist to the variational problem solved in the subspace of the full phase space satisfying some symmetry conditions. Following the links given I found this paper where they explain the method. It's not exactly a ...


3

Moshe, do you have access to coursework at the University of Idaho? They have a course listing that has your question in the title of the course. Math WS547 Numerical Analysis of Elliptic PDE's (3 cr) WSU Math 546 However, they don't seem to describe the course itself beyond that. I did find this paper entitled: LECTURES on COMPUTATIONAL NUMERICAL ...


3

There is no simple way. The "standard" way is to solve Poisson equation with proper boundary conditions (constant $\varphi$ at the surface). Out of potential distribution it is easy to extract charge distribution. For simple shapes (infinite plane, sphere, etc) it is possible to solve the problem analytically. For arbitrary shape there is no simple ...


3

Franz Pretorius has worked on this and developed animations. http://prl.aps.org/abstract/PRL/v95/i12/e121101 The field is numerical relativity. Matthew Choptuik also, I believe, has done work on this.


3

According to general relativity, a pair of massive bodies that orbit each other emits gravitational waves - for analogous reasons to the reasons why accelerating charges in electrodynamics emit electromagnetic waves. So there can't be any static solutions resembling binary stars or binary black holes. The solutions have to be non-static and a complicated ...


2

Maybe they do; but Fourier transforms have an inherent flaw which makes them less than useful for such cases. That flaw is that they are transforms of steady state conditions: The initial and final conditions of the system are assumed to be the same; and transients are not considered. The transform for which you seek is the Laplace transform. Laplace is a ...


2

A way of physically thinking about this is that a two body problem in general relativity does not generally have closed orbits. If one of the bodies is very large and the other a small satellite the problem is integrable. The periapsis (perihelion) advance of the small satellite is repeated with each orbit, which makes the problem integrable. If the two ...


2

The reason that these sorts of libraries don't exist is because the particular algorithm that you use to do the calculation will depend upon the exact details of the light field and the input and output planes you are trying to compute. For example, let me outline the the simplest case for this sort of calcualtion: The input light field $g_0$ has a slowly ...


2

The most important thing is conservation of momentum to describe the collisions. This part is actually quite straightforward, but before you get to collisions you should model the motion of single balls. Obviously, you will describe the balls classically and probably not at relativistic speeds (though that would be interesting...) so pretty much all you ...


2

So what should the boundary conditions be for a segment of an infinite line. Let us explore the options ($u$ is displacement, $u'$ is slope and $u''$ is curvature): Fixed Ends: $u(0)=0$, $u(L)=0$ Mirror Ends: $u'=0$, $u'(L)=0$ Free Ends: $u''=0$, $u''(L)=0$ Periodic Displacements: $u(0)=u(L)$, $u'(0)=u'(L)$ Periodic Tensions: $u'(0)=u'(L)$, $u''(0)=u''(L)$ ...


2

Two books to get you started in general computational radiation transport: Computational Methods of Neutron Transport, written by E.E. Lewis, edited by W.F., Jr. Miller, ISBN 0-89448-452-4 Monte Carlo Particle Transport Methods: Neutron and Photon Calculations, written by Ivan Lux and Laszlo Koblinger, ISBN 0-8493-6074-9 The supporting materials for the ...


2

The main problem in your proposed equation is that the electromagnetic equation with the D'Alambertian over the vector potential is not in Hamiltonian form, this means that the separation of solutions in Sturm–Liouville eigenstates of the energy operator is not manifest in the equation. Without that, you cannot find eigenstates of the coupled system. You ...


2

Let me take a slightly different perspective to Ron Maimon and say that the answer depends on whether you're after an exact solution of some mathematical model, or whether you want to calculate the exact physical behaviour. Any method for calculating the physical properties of a system rely upon an approximation. If you choose some model you can certainly ...


1

The shape of the induced eddy currents in induction cooking will depend on the shape of the fluctuating magnetic field and the shape of the cooking vessel. Certainly there are commonalities though and I suspect the diagram you have is similar. Also, in the case of induction cooking, the field is varying rather than the vessel moving and this will cause the ...


1

The mechanism you would use to identify the candies would be one of computer vision (which is off topic for this site try http://stats.stackexchange.com/). However, the exact mathematical method you would/could use to determine the broken candies from the okay ones depends on the shape of these candies. Lets say for arguments sake that they were spherical; ...


1

It's not totally clear to me what you want, but a good prescription for generating an ensemble of thermal states is Steve White's "Minimally Entangled Typical Thermal States" algorithm, where one starts with a high-temperature random state then evolves it in imaginary time, much like the second suggestion you listed. This random ensemble of thermal states ...


1

Since you seem to be using finite differences, you should look at the paper of Hadley, titled 'Transparent Boundary Condition for the Beam Propagation Method' - without any treatment of the boundary values you're automatically assuming a Dirichlet boundary condition. You can incorporate the boundary conditions into your square-root differential operator ...


1

You could try MEEP. It is an opensource tool for light manipulation. It automatically does all the FFTs and has a built-in Gaussian source and it also produces nice plots of 1D, 2D, 3D numerical simulations. It also has a python wrapper.


1

I use pylab numpy scipy matplotlib (and matplotlib.mlab) examples f.i. here and doc and nice ref here or use the online integrated python environment sagenb.org (open account or download package, upload files and start working anywhere) nice plots ... to do digital filters explore the online book The Scientist and Engineer's Guide to Digital Signal ...



Only top voted, non community-wiki answers of a minimum length are eligible