Unruled Notebook

Entries from June 2008

Whale Vel

June 29, 2008 · 6 Comments

Just then we finished cutting a picture of a whale shark from the recent Young World magazine and pasting it in a special notebook. And as I was helping the wife in her kitchen work by talking to her soothingly about life, along comes the kid and declares

Kid: I am going underwater to take pictures of whale shark. And I am not taking you.

Me: (peeved) OK, let me go to Pazhani and see Lord Muruga. He also has a Vel (Tamil for javelin or spear and pronounced as whale) with him and let me see that one.

Kid: My whale is big. Bigger than your murugar vel.

Me: Lord Muruga can get bigger whenever he chooses to. So will his Vel, along the way.

Kid: My whale will be underwater.

Me: If Lord Muruga throws it, this vel will also go underwater.

Kid: My whale can travel. It can go even to Australia.

Me: (where does this kid learn all this) OK, I can ask muruga to drop it real fast so that my Vel also can travel to Australia underwater.

Kid: (by now almost shrieking, not realizing I am also on the verge) My whale is a mammal. It can have babies…

Me: !@#^&

(the last I overheard, the kid was busy instructing all of mom’s paint brushes that she is taking along with her underwater to stay inside her back-pack as they could get wet and dilute their “color”…)

Categories: Muse
Tagged: , , , ,

Bidisperse Porous Media

June 28, 2008 · 2 Comments

ResearchBlogging.orgA porous medium whose solid matrix itself is another porous medium can be considered as a bi-disperse  or bi-dispersed porous medium (BDPM). Z. Q. Chen et al. in their 2000 paper defined such a BDPM to be composed of clusters of large particles that are agglomerates of small particles. So the macro-porous medium is made of a solid matrix region which in turn is made of a micro-porous medium. In between such micro-porous clusters there are definite macro-pores. Likewise, the micro-porous medium could contain micro or even smaller sized pores. Usually, these micro and macro pores are connected and the same fluid pervades through all of them.

Some applications of such BDPM modeling are in PEM Fuel cells, porous wicks in heat pipes, modeling granular mixtures and chemical reactors.

Accompanying two pictures, one a schematic (modified from Fig. 1 in this 2007 ijhmt paper), the other a possible BDPM reality, should help visualizing the concept.

Another useful explanation for the BDPM from the 2006, 3rd edition of Convection in Porous Medium by Nield and Bejan, is

A BDPM … considered as a porous medium in which fractures or tunnels have been introduced.

The porous medium in the above definition is our micro-porous medium and the fractures or tunnels are the macro-pores that runs through them, resulting in a BDPM.

For studying the heat transfer in such BDPM, it could be assumed the micro-porous medium can be represented by a single representative temperature, implicitly assuming local thermal equilibrium between the fluid and the solid in that micro-porous medium.

Mass diffusion in such BDPM have been studied in late 1980s in the chemical engineering literature, where chemical reactors can be modeled as BDPM.

Standard (good) research questions are how to determine the effective permeability and effective thermal conductivity of such a BDPM. Use of the terminology BDPM dates back to 1990s in works concerning porous medium of dual porosity. Because of the fine structure embedded in a coarser structure of similar characteristics and function, fractals are used to model such BDPM, although BDPM offers only a two step reduction when compared to standard fractal objects. Fractal models for BDPM were formulated by Yu and Cheng in 2002. Two experiments were also attempted (published in 2000) to determine the thermal conductivity of the BDPM. More recently (2003 onwards), Nield and Kuznetsov have worked on modeling convection in such BDPM.

A reference list of research literature mostly on convection in BDPM is provided below, which should be updated periodically. If you happen to notice any relevant paper, kindly put a link (if available) or reference in the comments.

NIELD, D., & KUZNETSOV, A. (2008). Natural convection about a vertical plate embedded in a bidisperse porous medium International Journal of Heat and Mass Transfer, 51 (7-8), 1658-1664 DOI: 10.1016/j.ijheatmasstransfer.2007.07.011

References

  • Daniel Arnost, Petr Schneider, Effective diffusivities from dynamic diffusion cell: the general moments analysis, Chemical Engineering Science Volume 49, Issue 3, , 1994, Pages 393-399. [ doi:10.1016/0009-2509(94)87010-1 ]
  • Burghardt, A.; Rogut, J. & Gotkowska, J., Diffusion coefficients in bidisperse porous structures, Chemical Engineering Science, 1988, 43, 2463-2476
  • Chen, Z. Q.; Cheng, P. & Hsu, C. T., A theoretical and experimental study on stagnant thermal conductivity of bi-dispersed porous media, International Communications in Heat and Mass Transfer, 2000, 27, 601-610 [ doi:10.1016/S0735-1933(00)00142-1 ]
  • Z. Q. Chen, P. Cheng, T. S. Zhao, An experimental study of two phase flow and boiling heat transfer in bi-dispersed porous channels, International Communications in Heat and Mass Transfer Volume 27, Issue 3, , April 2000, Pages 293-302. [ doi:10.1016/S0735-1933(00)00110-X ]
  • Markicevic, B. & Djilali, N., Two-scale modeling in porous media: Relative permeability predictions, Physics of Fluids, 2006, 18, 033101 [http://link.aip.org/link/?PHFLE6/18/033101/1 ] [ DOI:10.1063/1.2174877 ]
  • Moutsopoulos, K. & Koch, D., Hydrodynamic and boundary-layer dispersion in bidisperse porous media, Journal of Fluid Mechanics, 1999, 385, 359-79
  • Nield, D. A. & Bejan, A., Convection in Porous Media, Springer-Verlag, Newyork, 2006 (pages 25-26 and 99-100)
  • Nield, D.A. and Kuznetsov, A.V., Natural convection about a vertical plate embedded in a bidisperse porous medium, International Journal of Heat and Mass Transfer, 51(7), p.1658-1664, Apr 2008 [doi:10.1016/j.ijheatmasstransfer.2007.07.011 ]
  • Nield, D. & Kuznetsov, A., The effect of combined vertical and horizontal heterogeneity on the onset of convection in a bidisperse porous medium, International Journal of Heat and Mass Transfer, 2007, 50, 3329-3339 [doi:10.1016/j.ijheatmasstransfer.2007.01.027]
  • Nield, D. & Kuznetsov, A., The onset of convection in a bidisperse porous medium, International Journal of Heat and Mass Transfer, 2006, 49, 3068-3074 [doi:10.1016/j.ijheatmasstransfer.2006.02.008]
  • Nield, D. & Kuznetsov, A., Ingham, D. & Pop, I. (ed.), Heat transfer in bidisperse porous media, Pergamon, 2005, 34-59
  • Nield, D. & Kuznetsov, A., A two-velocity two-temperature model for a bi-dispersed porous medium: Forced convection in a channel, Transport in Porous Media, 2005, 59, 325-339 [Abstract] [ DOI: 10.1007/s11242-004-1685-y ]
  • Nield, D. & Kuznetsov, A., Thermally developing forced convection in a channel occupied by a porous medium saturated by a non-Newtonian fluid, International Journal of Heat and Mass Transfer, 2005, 48, 1214-1218
  • Nield, D. & Kuznetsov, A., Forced convection in a bi-disperse porous medium channel: a conjugate problem, International Journal of Heat and Mass Transfer, 2004, 47, 5375-5380
  • Kuznetsov, A. and Nield, D., Thermally Developing Forced Convection in a Bidisperse Porous Medium, Journal of Porous Media, 2006, 5, 393-402. [DOI: 10.1615/JPorMedia.v9.i5]
  • Wang, G.; Johannessen, E.; Kleijn, C. R.; de Leeuw, S. W. & Coppens, M., Optimizing transport in nanostructured catalysts: A computational study
    Chemical Engineering Science, 2007, 62, 5110-5116
  • Yu, B. & Cheng, P., A fractal permeability model for bi-dispersed porous media, International Journal of Heat and Mass Transfer, 2002, 45, 2983-2993 [ doi:10.1016/S0017-9310(02)00014-5 ]
  • Yu, B. & Cheng, P., Fractal models for the effective thermal conductivity of bidispersed porous media, Journal of Thermophysics and Heat Transfer, 2002, 16, 22-29

Some more related work

  • Mitra Dadvar, Muhammad Sahimi, The effective diffusivities in porous media with and without nonlinear reactions, Chemical Engineering Science Volume 62, Issue 5, , March 2007, Pages 1466-1476. [ doi:10.1016/j.ces.2006.12.002]
  • C. R. Ethier, “Flow through mixed fibrous porous materials,” AIChE J. 37, 1227 (1991). [Inspec] [ISI]
  • Ashi Ofir, Snir Dor, Larisa Grinis, Arie Zaban, Thomas Dittrich, and Juan Bisquert, Porosity dependence of electron percolation in nanoporous TiO[sub 2] layers,  J. Chem. Phys. 128, 064703 (2008), DOI:10.1063/1.2837807 [ http://link.aip.org/link/?JCPSA6/128/064703/1 ]
  • G. A. Heeter, A. I. Liapis, Estimation of pore diameter for intraparticle fluid flow in bidisperse porous chromatographic particles, Journal of Chromatography A Volume 761, Issues 1-2, , 14 February 1997, Pages 35-40. [ doi:10.1016/S0021-9673(96)00791-1 ]
  • A. -R. A. Khaled, K. Vafai, The role of porous media in modeling flow and heat transfer in biological tissues, International Journal of Heat and Mass Transfer Volume 46, Issue 26, , December 2003, Pages 4989-5003. [link] [ doi:10.1016/S0017-9310(03)00301-6 ]
  • R. Pitchumani and B. Ramakrishnan , A fractal geometry model for evaluating permeabilities of porous preforms used in liquid composite molding. Int. J. Heat Mass Transfer 42 (1999), pp. 2219–2232. Abstract [ doi:10.1016/S0017-9310(98)00261-0 ]
  • B.M. Yu and J.H. Li, Some fractal characters of porous media, Fractals 9 (2001), pp. 365–372.
  • Li-Wu Fan, Ya-Cai Hu, Tian Tian, Zi-Tao Yu, The prediction of effective thermal conductivities perpendicular to the fibres of wood using a fractal model and an improved transient measurement technique, International Journal of Heat and Mass Transfer Volume 49, Issues 21-22, , October 2006, Pages 4116-4123. [ doi:10.1016/j.ijheatmasstransfer.2006.03.027]
  • Boming Yu, Mingqing Zou, Yongjin Feng, Permeability of fractal porous media by Monte Carlo simulations, International Journal of Heat and Mass Transfer Volume 48, Issue 13, , June 2005, Pages 2787-2794. [ doi:10.1016/j.ijheatmasstransfer.2005.02.008 ]
  • Mingqing Zou, Boming Yu, Yongjin Feng, Peng Xu, A Monte Carlo method for simulating fractal surfaces, Physica A: Statistical Mechanics and its Applications Volume 386, Issue 1, , 1 December 2007, Pages 176-186. [ doi:10.1016/j.physa.2007.07.058 ]
  • Jinsui Wu, Boming Yu, A fractal resistance model for flow through porous media, International Journal of Heat and Mass Transfer Volume 50, Issues 19-20, , September 2007, Pages 3925-3932. [ doi:10.1016/j.ijheatmasstransfer.2007.02.009 ]
  • Peng Xu, Boming Yu, Meijuan Yun, Mingqing Zou, Heat conduction in fractal tree-like branched networks, International Journal of Heat and Mass Transfer Volume 49, Issues 19-20, , September 2006, Pages 3746-3751.[ doi:10.1016/j.ijheatmasstransfer.2006.01.033 ]

Categories: Porous Medium · Read List · Research Notes
Tagged: , , , , , , ,

Optimization and Genetic Algorithms

June 26, 2008 · 3 Comments

When we say we want to optimize some situation or a solution for a problem, we usually mean we are looking for a way to achieve a combination of better of many (or at least two) worlds. We want only better of these worlds (or problem parameters) because achieving the best in each one of them is seen to restrict the other worlds (or problem parameters) from becoming its best.

To give a quick example, to maximize his potential energy, Humpty Dumpty needs to sit (and remain) on the wall. To maximize his kinetic energy, he needs to have a great fall and be at a place just before actually hitting the ground. Obviously the best position for maximizing one parameter restricts the other from becoming maximized. To have forever a better of both worlds of potential and kinetic energy, Humpty Dumpty should have a fall and should not touch the ground forever.

Without contemplating further on how to fall upward and turning back to our optimization, these restrictions on the problem parameters are what are known as constraints in an optimization problem. An optimization without the related constraint(s) doesn’t necessarily optimize. There is an excellent introductory article by Prof. A. K. Mallik in the June 2008 issue of Resonance (see reference [1] below) that explains optimization problems in elementary geometry, without much calculus – a necessary mathematical tool for several optimization methods.

So if we want the best of one particular world in our optimization problem solution, then it may not be possible to achieve simultaneously the best in other worlds. Such multi variate optimization problem are abundant in all the engineering disciplines. In fact, to engineer means to optimize something for a given situation with the given (restricted) resources. Close to our context, given the finiteness of time (per day or month or in one’s lifetime), how to do science, research, learn, teach, write, crib, sleep, blog, family, dream and vegetate, all for as much time as possible, is a multi parameter optimization problem that has as many solutions as we desire it to be. And that is where effective solution search techniques enter.

In engineering we encounter an optimization problem for which we anticipate multiple solutions based on restrictive ranges of real life parameters that affect the problem. More than one such solution nevertheless might satisfy the optimization constraints and yield a reasonably correct (or required) global solution within an accepted margin of discrepancy. To find which of these solutions is better (uniqueness based on restrictions) would not be possible through perhaps a linear solution methodology.

Genetic Algorithm is a search technique, proposed by John Henry Holland in the 1970s, that finds for us the better (optimized) solution by iteratively combining an initial pool of solutions, each successive iteration bettering in some way the solution from the previous iteration., in a sequence of “genetic evolution” process. The final satisfactory solution set or “genetically evolved population”, fits the imposed optimization constraints for the problem.

Genetic Algorithms involve an “explorative logic” which ensures that a large number of solutions, marginally better or worse are considered, while avoiding convergence on local maxima. For instance, while seeking a globally optimized solution, if only a single constraint of “checking whether a variable become greater than a value” is used to verify the arrival of the correct solution, it is possible that such a solution satisfies only a local maximum for that variable when it reaches a maximum value. Obviously, in a problem involving constraints for more than one variable this local maximum solution need not be the correct one.

In genetic algorithm each solution is considered as a genome so the initial solution set (randomly generated and need not contain the optimized solution) is a population of genomes differing unique enough in some minor changes in the parameter set considered to generate the solutions. These solutions are then iteratively changed to better solutions by doing iteratively two things. First the parameters governing the solution are changed slightly by an exhaustive sequence (random and systematic) covering the entire range for all the parameters involved. To facilitate this process usually the solution population is represented as a genome sequence or a computer byte with multiple bits. These binary representation of the “individual” apportions specific number of bits for parameters involved, depending on the problem.

Each new generation of “genome population” thus generated are checked with an “objective or fitness function” for their validity or “fertility”. Depending on the fitness values, pairs of individuals from this solution set called “parents” are chosen for “breeding” using operators borrowed form natural genetics. The level of fitness of an individual dictates its chances of reproducing and surviving through generations. Thus GA starts with a randomly chosen population and refines them over generations.

Breeding is done in two ways. Mimicking the evolution process as an algorithm, mutation introduces some randomness to the solution “generations”, say, by randomly changing some parameter over a range. Crossover usually combines two parent solutions and produces two more children or offspring solutions that carry “genetic information” or parameter combinations from the parents. In the binary representation mentioned above, crossover can be understood to involve swapping of some of the bits of the two selected parents from a specified bit position. Mutation randomly alters a bit in the representation of the offspring depending on mutation probability thus aiding the “explorative logic”. A new population is created thus, which retains the best individuals of the previous generation and replaces the rest of the individuals by the offspring.

To perceive the analogue, a sample of chromosomal mutations out of the five possible ones and the concept of crossover is shown in the accompanying pictures. Source of these pictures are [here] and [here], carrying further explanations. Since the performance of GA is dependent on the choice of various GA parameters like Crossover Probability, Mutation Probability etc., prediction of the choice of parameters for ensuring best results is seldom possible.

The solution subset that deviates away from the desired solution (deduced from their deviation from the objective function) are “faded out” of subsequent iterations. By fading out we mean the solution subsets are not immediately discarded in the next iteration itself. They are carried over for enough future iterations to ensure their evolutionary refuse towards the required solution – and allowed a “natural extinction” in the simulated “survival of the fittest” competition.

The above sequence of operations is repeated for many iterations, smoking a reasonably high end computer depending on the problem parameters or the size of the binary representation of the sample solution. After a few or million such generations GA converges to the best solution. At this simulation instant almost all the individuals in the solution population represent the same solution, revealing identical values for their fitness function. Thus the population is unlikely to evolve further.

To summarize the above explanation: There are three important components in writing a successful genetic algorithm (1) objective function should be properly defined to catch the optimum you are looking for (2) genetic representation or generic solution represented as a “genome” should be properly defined and used (3) definition of the genetic operators that govern the mutation or crossover of one solution into another that is closer to the desired solution.

Shall stop the basics here. Use the reference links for further reading. In the next installment, I shall explain an application of GA to designing “optimum” phase change material heat sinks for electronics cooling.

Further Reading

1. Optimization Problems in Elementary Geometry, A. K. Mallik, Resonance – Journal of science education, June 2008, Vol. 13, Num. 6, pp. 561 – 582. [link to pdf and June 2008 issue]

2. For web resources, check the wikipedia write-up for Genetic Algorithm and the external links given there.

3. John Henry Holland, Adaptation in Natural and Artificial Systems: 2nd edition. MIT Press. 1992

4. Global Optimization Algorithms – Theory and Application (link is a pdf document) by Thomas Weise. A free book on genetic algorithm.

5.  When Least is Best by Paul J. Nahin. A good book on optimization solutions without using calculus.

Categories: Maths · Science Notes · Thermal Sciences
Tagged: , , , , , , , , ,

Recent Carnatic Music Recommendations

June 23, 2008 · 7 Comments

Here are some observations on three Carnatic musicians I listened recently. You will perhaps be tempted to skip the post but I assure I at least don’t pass off these opinions as erudite promulgations from a holier than thou. On the other hand, if you haven’t yet attempted to enjoy this South Indian Classical Music form, you are missing something good in life. Same holds true for any classical music form and music in general.

Sanjay Subramanian is one of my contemporary favorites. Reasons abound on why it is so the important ones being sincerity, confidence, creativity and the willingness to experiment and stretch himself and the music. Of course my qualifying remarks are superfluous if you already have listened to him at least once.

In a recent concert I attended at the TTD Devasthanam in T. Nagar, Chennai, he gave expositions on manIrangU and sAmA, two relatively uncommon ragas in the present day concert repertoire. Manirangu if not handled deftly, would sound like a madhyamAvathi or srI. sAmA is another raga that could get daunting. In Shankar’s film Anniyan, there was a song (listen here) tuned by Harris Jayaraj that had shades of sAmA (another one was in nAtai). Even a well tread Mohanam reveals its folkish melodies when Sanjay sings it.There was also a cute Annamacharya krithi that starts Palukute, set in AbhEri. In a recent album of Sanjay this song was preceded with a raga exposition in AbhEri that reminded of rare melodious phrases one gets to hear only when this rAga is delivered in a nAgaswaram. In fact that recent album by Sanjay titled Chaturshram contains four apt rare gems in bEgadA, AbhEri, nAtaikurinji, sindhubhairavi ragas, delivered in good steam.

Every time I listen to, accompanist Nagai Murali amazes me with his sweetest violin tone, even while he follows/shadows/repeats the briga rich creative solos of Sanjay. And whenever a rare raga is exposited by the main artist – as done in the above concert the other day – I always get the feeling Murali does the exposition even better, when he gets his turn. Arun Prakash on the mridangam is a musician’s favorite. Teams like this usually deliver irrespective of the kirtanas (songs) and concert location.

Another Sanjay concert I listened to is the unedited CD version of one of his Dec 2007 concert. After a standard set of srI, varAli, mOhanam etc., there was an exposition of a rare mElakarthA raga shUlini. I wouldn’t rate the effort because I am hearing the raga for the first time. There is no reference point in the (recorded) hoary past of Carnatic Music to which I can compare this effort (a good thing). Sanjay is one contemporary musician who has a penchant for trying out rare mElakarthAs in concerts as a sub main or main piece. Over the years I have listened to him elaborate some rarities like nAtakapriya, gAyakapriya, natabhairavi, jyOthiswarUpini, chalanAta, kOsalam, sucharithra, rIshabapriya, kOmalangi (a janya raga) and now, shUlini.

There is a dogmatic “purist” group of seasoned listeners for whom a Carnatic concert must be from a selection out of about thirty ragas – a drop in the ocean of Carnatic music with a million possibilities. Familiarity to these ragas and some of their songs over the years and how well it is regurgitated in a concert breeds their approval and credit for an artist. Nothing fundamentally wrong with this expectation. Such familiarity also breeds their form of punditry and eventually, art stagnation. Artists, like any other creator, are very insecure. They need constant encouragement for practicing the art form their way. Only then the boundaries of art could be stretched and art itself flourish. In this context, as a listener, one should be ready to go to the next concert, even when such experiment flop. Fortunately, there is enough support and recognition for artists like Sanjay and their experiments, amongst a slice of concert goers.

Now for two more non-run-of-the-mill singers.

Prassala Ponnammal is an aged gem whose recent concert I listened to had a remarkably soothing swarajathi (rAvE HimagirikumAri in tOdi) in the right ambulating and lilting tempo. I couldn’t come out of the reverie for a long time. Such songs identified as swarajathis are set to slow and very slow tempos with long sententious and intricate melodies (swara korvais) accompanied with suitable lyrics. Syama Sastri, one of the Carnatic musical trinities, has excelled in composing some remarkable swarajathis. It would take about 10 to 15 minutes to sing just the song. It takes a certain confidence and composure to sing a swarajathi as a sub main piece of a concert, as Prassala Ponnammal did that day. In the past Semmangudi Srinivasa Iyer used to perform swarajathis (ambA kAmAkshi in particular) with elan, a tradition T. M. Krishna attempts to imbibe nowadays in his concerts. In Prassala Ponnammal’s concert there was also an AlApanai (exposition) of Arabhi that was good. There were two nEravals at this concert, a rarity nowadays and both were intelligently done.

Vasundhra Rajagopal, notwithstanding her Srirangam nativity, is another lady I respect for the music. I discovered her late in her career and the first time I listened to her two years back at the Music Academy morning concert she gave good AlApanais (expositions) of relatively uncommon ragas like mandAri, husEni and dEvagAndhAri. It is since ages I hear an AlApanai for husEni – the only recorded versions I have are by Ariyakudi Ramanuja Iyengar delivered in the 1960s. There was an RTP (ragam thanam pallavi) for dEvagAndhAri too. Since then I listened to her thematic concert on Arunachala Kavi’s compositions set to tune by herself in reasonably eclectic ragas. I missed attending (and regret till this day) her thematic concert on Dec 31, 2007, on Paasurappadi Ramayanam of Peria Vachan Pillai. The concert also had accompanying commentaries by Velukkudi Krishnan (any recordings of this concert available?).

Rummaging through Landmark store’s dirty bottom boxes, last week I was able to relieve them of her 2001 dated CD on dEvagAndhAri for a paltry price. This CD also contained initial discourses on the raga by Carnatic music theorists in which there was a specific mention that dEvagAndhAri is not taken for full expositions as in a RTP.

After such punditries came an excellent forty five minute long RTP in dEvagAndhAri by Vasundhra.

Categories: Carnatic Music · Muse
Tagged: , , , ,

Quotes

June 22, 2008 · Leave a Comment

1) There is only one basic human right, the right to do as you damn well please. And with it comes the only basic human duty, the duty to take the consequences.

2) I like to think of anything stupid I’ve done as a “learning experience.” It makes me feel less stupid.

3) You know your children are growing up when they stop asking you where they came from and refuse to tell you where they’re going.

- P. J. O’ Rourke [ more from him]

Categories: Asides · Quotes