Laplacian Eigenfunctions

You can get very interesting result when plotting different linear combination of eigenfunctions of the Laplacian on a given (convex) shape. Here I use a square plate [0,1] by [0,1] and show the zero iso-value of the some of such linear combination. As time evolves, the linear combination changes.

Since the Eigenfunctions are symmetric, I only plot them on the [0, 0.5] interval. In that case, the eigenfunction are of the form

With (n1,n2) such that

The associated eigenvalue l= 340*pi^2 is of multiplicity 4 (12,14) (18,4) (14,12) (4,18). Below, same thing with l = 5525 which is of multiplicity 12.

Those isocontours are called nodal lines. By the unique continuation property, those nodal lives consist of curves in the interior of the domain (her a square, but the result is valid for any convex shape) which are infinetly differentiable. You may also note, but that is a well known theorem, that when nodal lines interesect a boundary (except may be at corner points), then they form egal angle. A single nodal line intersects the boundary with right angles, two intersect it at 60° angles, and so forth. This also explain why traction/compression cracks always cut boundaries of their domain at right angles, but that's another story.

Well, my point here was not simply to show funny stuff, but also to point out that using gnuplot and image magick is really simple to create animations.

Comments

Shear waves, medecine and brain

Yesterday evening, too bored by what TV was proposing to me, I decided to watch a conference of Mathias Fink , a french researcher working on multidisciplinary application of waves. Specially shear waves. Here is a brief summary of his talk. In solids, waves have two principal components: compression waves (P-waves for primary) moving in the direction of propagation, and shear waves (S-waves, for secondary) that make ripples in the plane orthogonal to that direction. Since compression waves propagate in the direction of propagation, they move faster than shear waves. Usually ultrasound equipment in medicine only use compressional waves. But since human tissues have a high bulk modulus, the P-wave speed is relatively constant (around 1580 m/s). Human tissues are very stiff if you apply isotropic constraints on them (like pressure of water). However M. Fink and his colleagues proposed a new way to investigate human tissues by first sending a strong compressional wave in ...

We're not playing dices !

Software programming and design is an amazingly complex task. Specially when it concerns numerical applications, that generally require optimization to get results in a reasonable time. For that purpose the basic pattern is usually to write the code as simple as possible, debug it and when it works to begin the optimization process. This asks for some nerves, and patience, two things I usually don't have in real life , but for programming yes. This design pattern makes sense because we are doing things really sequentially. You first wrote some c++ class and then add some feature progressively, and when you come up with essentially different concept (I mean a concept that should be well separated from the first one), you write an other class and so on... Usually the class-writing process follows the solving process you have in mind. For instance you have to make some initialization on your model first, then you declare some variables for the computation, the computation then f...

My Desk is an Optimal Container

When you do informatics codes, sooner or later you are concerned with data containers and their associated complexities for inserting, getting, sorting or whatsoever. Some of them are lists, vectors, map, set etc... Sometimes you use a container for some stuff and you get unexpected low behavior. Then you go into your manager's office and ask for complexity of what you are doing. By the way, he convinces you that you are the last of the idiots of computer scientists (thing I take as an honor), and you just change your container. But the people who are working on containers (have you ever looked to the STL (standard template library) to see how containers are implemented ? Not sure those people are human...) should try to investigate how my desk is working. It is very messy, there is a lot of stuff on it (even food sometimes) but I always find my items in a constant time . Well, almost always. Notably when someone put it in order (presumably the housekeeper or some terribly i...

New Horizons

Search This Blog