Wednesday, April 28, 2010

Oil leak, volcano eruption, and MODIS imagery

The Moderate Resolution Imaging Spectroradiometer (MODIS) is a sensor developed by Raytheon that rides on NASA’s earth observing satellites Terra and Aqua. Terra is now 10 years old (launched 1999) far exceeding its expected design life of 5 years. Aqua was launched in 2002.

The Earth Observing System (EOS) acquires images of all natural (or not) phenomena on planet Earth (e.g. fires, volcano eruptions, sand storms, smoke, pollution, etc.) NASA makes images from the satellites available to the public within hours after collection. As they say ‘it is almost real time imagery.’ For example, here is one of the images collected today (4/28/10) around 9:30 UTC. Seems rainy in Crete!

Europe_3_03.2010118.terra.1km

 

Some of the images are handpicked and posted with explanation in a gallery maintained by NASA. The image of the day in that gallery is an image of the oil leak from the site where an oil platform sank in the Gulf of Mexico. The latest estimate is that 5,000 barrels of oil are leaking from the site per day (up from 1,000 barrels per day.)

image04282010_250m

The incident happened close to the south-most part of the oil slick which travels towards Mississippi Delta.

Another cool set of images were collected a few days ago and they show the Eyjafjallajokull Volcano in Iceland. In the one shown below Iceland is covered by ashes and steam. The volcano is located in the southern part of the island, close to the coast.

Iceland.A2010110.1155.1km

And of course, here is an image of the ash plume over Europe taken on 4/16/10 at 10:45 UTC. The eruption started on 4/14/10 so this image shows the airspace over Northern Europe while almost all commercial flights were cancelled.

Germany.A2010106.1045.2km

England, most of Germany, most of Holland, and Poland are clearly covered by the ashes from the volcano.

Tuesday, April 20, 2010

Google posts censorship demands by country

A couple of days ago Google posted on the web the censorship demands made by countries to Google and YouTube between July 1, 2009 and December 31, 2009. The website URL is http://www.google.com/governmentrequests/.

It is a pretty interesting tool. If you click the little box with the number of data or removal requests per country then Google will show how it has complied (or not) with the request and what the break down is.

No data on China. Google a little red clickable box indicating that ‘Chinese officials consider censorship demands as state secrets, so we cannot disclose that information at this time.

Interestingly, Brazil is number one on both counts (data and removal requests.) The reason for that seems to be a a social networking site called Orkut that is owned and operated by Google. In Brazil Google complied in 291/3663 request and 218/291 were related to Orkut (presumably taking down offensive material.)

Interesting tool!

Monday, April 19, 2010

World Cup 2010: The Aerodynamics Of Jabulani

Jabulani is the name of the official ball of the World Cup in South Africa made by Adidas (like all balls since 1970.) Jabulani means ‘to celebrate’ in Zulu. The ball was unveiled last December in Cape Town.

Adidas posted a few of videos about the ball. In the videos they mention that the design of the ball took 6 years (!!) and that is made of eight thermally bonded 3D panels and is perfectly round. In one of the videos they show the whole manufacturing process which is pretty amazing for a ball.

In a couple of other videos they talk about the aerodynamic design of the ball.

In one of them they initially show flow past a smooth sphere in a wind tunnel.

Flow past sphere_

In the video frame shown here the streamlines are clearly visible as the air flows from left to right. Because of viscous effects (friction, thermal conduction, etc) the flow separates over the back of the ball creating a recirculating flow (i.e. turbulent wake) downstream of body. This separated flow is responsible for the pressure (or form) drag which contributes greatly to the total drag of the ball thus affecting its flight performance characteristics.

In the video, the separation points seem to be either ahead of the maximum thickness of the sphere or very close to it. That indicates a Reynolds number (Re) that is less than 5e5 (for a smooth sphere.) For such an Re the separation points form at about 82 degrees measured from the stagnation point (front center of the sphere.)

The dimensionless Reynolds number (named after its discoverer, Osborne Reynolds) is a measure of the ratio of inertia forces to viscous forces in a flow. It is given by Re = rVL/m, where V is the free stream velocity, r is the free stream density, L is the characteristic length of the body immersed in the fluid (in this case the diameter), and m is the free stream viscosity coefficient.

Assuming that a standard size 5 smooth ball (about 0.23 m in diameter) would fly in standard air (constant r = 1.2 kgr/m3) at 90 km/hr (25 m/sec) after a decent kick the Re would equal:

Re = 1.2 * 25 *0.23 /1002e-6 = 6.9e+3

This is well within the regime of Re < 5e+5. Above that the separation points would move further back and thus pressure drag would actually reduce. That is because the boundary layer (a thin layer of fluid very close to the surface of the body that the viscous effects matter) would transition from laminar to turbulent which separates at about 125 degrees (measured from the forward stagnation point.) This results in a thinner turbulent wake with higher pressures and a lower drag coefficient for the transition point (i.e Re ~ 5e+5 for the sphere.)

However, it would be impossible to achieve something like in football (i.e. the velocity of the ball in flight would have to be a few thousands km/hr!) Therefore another technique is needed to push the stagnation points further back and reduce the turbulent wake and thus the pressure drag.

Enter the ridges!

Jabulani_2

Later in the video the ball replaces the sphere in the wind tunnel but unfortunately they do not show the flow pattern around the ball, presumably to protect proprietary information. The engineering manager for Adidas, Dr. Tim Lucas, explains that the designers spent a lot of time on the design of the ridges on the ball and ‘they play a very important role so that the ball is stable in flight.

Presumably, the ridges break up the laminar flow within the boundary layer and transition it to turbulent at much lower Re thereby reducing the form drag and improving flight performance. The ‘rough’ surface of Jabulani should lower the critical Reynolds number and possibly make sudden drop in drag coefficient attainable, which would be a pain for goalkeepers!

Finally, these effects would affect the bending of the ball. If the critical Re is low enough to be in the region relevant to football then you could have Re larger than the critical one. In that case a kick that adds backspin to the ball would cause separation on the bottom side to occur earlier than on the top side. Therefore the bottom side will incur higher pressure and thus generate a higher upward force (lift.) Higher lift would result in longer balls and possibly flatter trajectories. More bad news for the goalkeepers!!

Thursday, April 15, 2010

MATLAB: Generate Carousel I

In earlier posts I described how to generate a rotating cylinder using MATLAB and how to define to coloring scheme of the rotating cylinder. The reason for these exercises was that I wanted to write code that would create a carousel for entertainment and because I can. I selected MATLAB because among other great attributes it is perfect for prototyping.

I am not done with the design yet but so far it consists of 3 elements: Cylinder, cone, and circle. The cylinder was described in the aforementioned posts. For the other two elements:

CONE:

The roof of the carousel is a cone described using the following parametric equations:

x = (h-u)/hrcostheta

y = (h-u)/hrsintheta

z = u

where, h is the height of the cone, r is the radius of the base,  u in [0,h] and theta in [0,2pi).

In the MATLAB codespace the following entry was used to describe a cone:

rCc = 15; %radius
thetaCc = linspace(0, 2*pi, N);
hCc = linspace(4, 7, M); %height

[UCc, hCc_] = meshgrid(thetaCc,hCc);

[XCc,YCc,ZCc] = gkcone(rCc,hCc_,UCc);

hhCc = surf(XCc,YCc,ZCc);
set(hhCc,'CData',colorCylinder);

=============================================================

function [x,y,z] = gkcone(radius,height,theta)

step = (max(max(height)) - height) ./max(max(height));

x = step .*radius .*cos( theta );
y = step .*radius .*sin( theta );
z = height;

CIRCLE:

Parametric equations were used to describe the circle too. The equations are the same as the cylinder but with a zero z-axis component.

I created 4 cylinders, 4 circles, and a cone with the following result:

The color scheme of the circles is solid white or solid gray. These where generated using:

color_ = ones(M,N+1,3); % WHITE

color_ = ones(M,N+1,3) * 0.9; % GRAY

Sunday, April 11, 2010

FIFA World Ranking Top 20: April 2010

Rank (previous) Team Points Comment
1 (1) Spain 1602  
2 (2) Brazil 1589  
3 (3) Netherlands 1261  
4 (6) Portugal 1214  
5 (4) Italy 1183 Cameroon managed a scoreless tie in Italy.
6 (5) Germany 1157 Lost to Argentina at home.
7 (8) England 1120 Won 3-1 over Egypt at home.
8 (7) France 1171 Lost to Spain at home. No blame for this one!
9 (9) Argentina 1075  
10 (11) Croatia 1051  
11 (10) Greece 1032 Slipped one place after the loss to Senegal.
12 (12) Russia 1009
13 (14) Chile 974  
14 (17) Egypt 967 Hm! Something is wrong with this ranking.
15 (13) Serbia 950  
16 (18) USA 948 Lost to Holland and gained!
17 (15) Mexico 921  
18 (19) Uruguay 902  
19 (23) Australia 898 Won a Continental Qualifier at home against the powerhouse of Indonesia!
20 (20) Cameroon 887  
20 (15) Switzerland 887 Lost to Uruguay at home.

No major changes. The next game for Greece is a World Cup friendly against Korea DPR at home on May 25.

Full table at the FIFA World Ranking website.

Wednesday, April 7, 2010

112th Element: Copernicium

For some reason this caught my attention. It is not exactly news since I remember reading about the discovery of element 112 by researchers in Germany years ago, while I was studying at the University of Washington in 1996.

The interesting aspect here is that the scientific community officially accepted the name Copernicium with the symbol Cn after 14 years of trying to resolve a dispute over who produced the element first. The Germans won! Needless to say that it was named after Nicolaus Copernicus that most regard as the father of modern astronomy.

Also see articles:

And only a few days ago it was reported that a team of Russian scientists synthesized element 117 (aka Ununseptium). Element 118 was discovered before 117.

Also see article:

So, given that there are currently 118 elements in the periodic table it might take a bit of time to name them all and give them a symbol!

Monday, April 5, 2010

Saturn's Inner Moon: Mimas

This is an image of Mimas, one of Saturn's moons (it has over 50 of them.) It was constructed by putting together 6 pictures taken by the Cassini spacecraft as it flew as close as 9500 km on February 13, 2010. Mimas has a diameter of 196 km and the image clearly shows a crater (called Herschel Crater after the British astronomer that discovered it) that has a diameter of 140 km. It is amazing that whatever left such a crater after impact did not destroy Mimas completely.



The image information from the NASA JPL webpage can be found in the catalog page for PIA12570.

BTW, the moon was named after the Titan Mimas from the Greek Mythology!