Slashdot Mirror


Phase Change in Fluids Simulated

brendotroy writes "After decades of work by the physics and computer science communities, scientists at the University of Rochester have finally created a mathematical model that will allow scientists to simulate and understand phase changes. This discovery 'could have an impact on everything from decaffeinating coffee to improving fuel cell efficiency in automobiles of the future.'"

12 of 181 comments (clear)

  1. Re:Too dense by techno-vampire · · Score: 5, Informative
    Can it show why lakes don't freeze from the bottom up as water approaches 0 Celsius?

    No, because that has to do with an entirely different, well-understood phenominum. Unlike most substances, water gets less dense when it gets near its freezing point instead of continuing to get denser. When it freezes, it gets even less dense. (This is caused by something called "hydrogen bonds," but I'm not going to go into that.) Thus, ice is slightly less dense than the water surrounding it, making it float.

    --
    Good, inexpensive web hosting
  2. Phase Change and Complexity by Quirk · · Score: 3, Informative
    The wiki page gives a general introduction to phase change. My limited exposure to phase change has come from trying to fathom the various ideas put forth under the banner of complexity. The Santa Fe Institute is the home base for Complexity Theory. A search on the Institutes site turns up a plethora of articles on phase change. One of the godfathers of complexity theory, Stuart Kauffman makes many references to the idea of phase change as it applies to his ideas of the origins of life and open, non-equilibrium systems. The International Society for Complexity, Information, and Design (ISCID) also puts up some interesting material.

    The above gives an introduction to phase change as it is considered in terms of Complexity Theory. Approaching phase change through complexity theory, even for an outsider like myself, gives insight into how far reaching are the results of insight into phase change.

    --
    "Academicians are more likely to share each other's toothbrush than each other's nomenclature."
    Cohen
    1. Re:Phase Change and Complexity by wass · · Score: 5, Informative
      The article says :
      "This is the first time a computer program could simulate a phase transition because the computers would always bog down at what's known as the 'critical slowdown.' We figured out a way to perform a kind of end-run around that critical point slowdown and the results allow us to calculate certain critical point properties for the first time."
      There have been previous methods to look at systems arbitrarily close to the critical point in phase transitions, and the article is misleading when it says this area has been off-limits to computers. I haven't read the actual Physical Review Letters article, but it appears the authors have come up with a novel algorithm, perhaps more ideally suited for fluids, to overcome critical slowdown. But this is not the first such algorithm, and there has been loads of prior computer simulation of phase transitions and critical phenomena to boot.

      For example, in the Ising Model or the Potts Model, one can examine system parameters arbitrarily close to the critical point, in finite time, using a Cluster Algorithm. This page gives some information on how the cluster algorithm. The page has a java applet graphically depicting the system for a variety of algorithms.

      Just for completeness, here's an Ising model applet that I wrote, which doesn't just have a system animation, but allows you to calculate and plot data (specific heat, magnetization, etc) as the system passes through the critical point. This applet uses the Metropolis algorithm for time advancement, hence it is subject to critical slowdown. In that respect, the applet is flawed because close to the critical point I don't generate enough Metropolis iterations to ensure the subsequent frame is sufficiently thermally indepdent from the previous state. However, the cluster algorithm would remove these limitations. This applet has actually been used in graduate physics classes at Johns Hopkins to demonstrate magnetic phase transitions.

      And also for completeness, here's a Potts model applet, but it doesn't acquire data for plotting like the Ising model. The Potts applet actually uses the Microcanonical ensemble, whereby the energy of the system is conserved, but the Ising applet uses the Canonical ensemble, where the system is in contact with a heat bath at some settable temperature.

      And in case anyone's curious, these applets (except for the first one) are part of the Java Virtual Physics Lab , which contains a few different physics java simulations I wrote to help with conceptual understanding.

      --

      make world, not war

    2. Re:Phase Change and Complexity by Anonymous Coward · · Score: 1, Informative

      General simulation methods are a very active area of research, actually. We are doing atom-by-atom simulation, but your suggestion that "in the worst case scenario you can do a very computationally intensive simulation" indicates that you're not familiar with just how computationally intensive these simulations are.

      * Because the potential energy (Monte Carlo) / forces (molecular dynamics) use pairwise summation of potentials, performance scales as the square of the number of atoms; performance becomes prohibitive with large system sizes. (A system that I ran recently took 10 hours to do a sufficient number of iterations on a single unit cell in a crystal; if I were to scale that up to eight unit cells, it would take 640 hours.)

      * In a typical molecular dynamics simulation, timesteps are on the order of picoseconds or femtoseconds, with overall simulation time topping out in the nanoseconds -- these timescales are necessary to capture the atomic interactions correctly. If correlation times are macroscopic, then the amount of simulation time required to capture them would be prohibitive. Second-order phase transition is certainly not the only problem that has these issues: one of the research projects in our lab involves measuring transport properties of relatively large molecules.

      For these reasons, there is a lot of ongoing research into algorithms for making larger simulations and more difficult systems more accessible.

  3. Slashdotted by Doomedsnowball · · Score: 5, Informative

    MEDIA CONTACT: Jonathan Sherwood (585) 273-4726, jonathan.sherwood@rochester.edu
    January 6, 2006
    Phase Change in Fluids Finally Simulated After Decades of Effort
    Eldred Chimowitz and Yonathan Shapir

    Everyone knows what happens to water when it boils--everyone, that is, except computers. Modeling the transformation process of matter moving from one phase to another, such as from liquid to gas, has been all but impossible near the critical point. This is due to the increasingly complex way molecules behave as they approach the change from one phase to another. Researchers at the University of Rochester, however, have now created a mathematical model that will allow scientists to simulate and understand phase changes, which could have an impact on everything from decaffeinating coffee to improving fuel cell efficiency in automobiles of the future. The findings have been published in Physical Review Letters.

    "This problem has baffled scientists for decades," says Yonathan Shapir, professor of physics and chemical engineering at the University of Rochester, and co-author of the paper. "This is the first time a computer program could simulate a phase transition because the computers would always bog down at what's known as the 'critical slowdown.' We figured out a way to perform a kind of end-run around that critical point slowdown and the results allow us to calculate certain critical point properties for the first time."

    "Critical slowdown" is a phenomenon that happens as matter moves from one phase to another near the critical point. As molecules in a gas, for instance, are cooled, they lose some of their motion, but are still moving around and bumping into each other. As the temperature drops to where the gas will change into a liquid, the molecules' motion becomes correlated, or connected, across larger and larger distances. That correlation is a bit like deciding where to go to dinner--quick and easy with two people, but takes forever for a group of 20 to take action. The broadening correlation dramatically increases the time it takes for the gas to reach an overall equilibrium, and that directly leads to an increase in computing time required, approaching infinity and bogging down as the gas crosses the point of phase change.

    To illustrate the effect, imagine a perfectly pure and still lake. If you drop a pebble into this lake, its ripples would spread outward, dissipating until the lake had returned to a calm equilibrium again. But, if you were to take this impossibly perfect lake just barely above the critical point and drop your pebble, the ripples would remain as ripples much longer--likely bouncing off the distant shores. This imaginary lake would take seemingly forever to return to its calm equilibrium again.

    The research team of Shapir, Eldred Chimowitz, professor in the Department of Chemical Engineering, and physics graduate student Subhranil De created a novel approach to tackle the phase-change process. They devised a computational model consisting of two separate reservoirs of fluid at equilibrium and near the critical point threshold. One reservoir was slightly more pressurized than its neighbor. The reservoirs were opened to each other and the pressure difference caused the fluids to mix. The team let the simulation run until the entire system reached thermodynamic equilibrium. By watching the rate that equilibrium returned, the team was able to calculate the behavior at the critical point. Their simulation findings match predictions and experimental results, including very precise measurements performed in microgravity on the Space Shuttle.

    "In principle, it's a difficult calculation," says Chimowitz. "Fluid systems require a different class of models than the common lattice models used by researchers who have studied dynamic critical behavior. These different classes give rise to different dynamic critical exponents and we found them, for the first time, in real fluid systems."

    The best known examples of phase changes are perhaps water to ice and

    --
    7h3$3 4r3n'7 7h3 Ðr01Ð$ ¥0 4r3 £00|{1n9 f0r. M0v3 4£0n9. --OB1
  4. Re:The Japanese have already done this by Kagura · · Score: 2, Informative

    NOT work safe, and disgusting. It's hentai, with terrible body mutilation in graphic, well-drawn detail.

  5. Not just plain phase changes by chthonicdaemon · · Score: 3, Informative

    I would like to point out that the article is not about plain phase changes, but rather about phase changes near the critical point , where liquid and gas phases become indistinguishable. Predicting ideal phase change behaviour has been done, but the critical point poses some unique challenges.

    --
    Languages aren't inherently fast -- implementations are efficient
  6. Re:Great but... by wass · · Score: 4, Informative
    In short... this does nothing for our "understanding" of phase changes.

    Wrong, I highly suggest you take a Phase Transitions and Critical Phenomena class if you want to see the utility of methods such as this. As I noted in another post, though, this isn't the first method to allow computer simulations of points arbitrarily close to criticality, there have been other algorithms (eg Cluster algorithm) to allow this too. But every new algorithm to get past critical slowdown is very useful.

    What we've learned in the past several decades in critical phenomena is how parameters change close to the critical point. For example, look up Critical Exponents and Scaling. What is very interesting is that critical exponents are unique to a universality class. So if you are able to take a new system and show that it boils down (no pun intended) to a previously-studied universality class, you can know instantly how various parameters will scale and change as a function of temperature, magnetic field, etc, close to the critical point.

    And to give you an example of this, look at Superconductivity. It was originally discovered by Onnes in 1911, but it took 46 years until the BCS Theory was adequately able to explain how Cooper Pairs form and how resistanceless supercurrent can flow quantum-mechanically. Such a theory is referred to as 'microscopic', meaning it deals with the fundamental physics involved, specifically the electron-phonon-electron interaction and how the Fermi sea is unstable to Pair condensation.

    However - alot of work was done prior to BCS dealing with 'macroscopic' theory, whereby certain laws were able to be formed (eg London equations for classical electrodynamics of a superconductor), we just didn't understand how or why they were valid.

    One such important example is the Ginzburg-Landau theory (Landau won the Nobel Prize decades ago, Ginzburg just got it a couple of years ago), which extends Landau's Theory of 2nd-Order Phase Transitions to use a complex order parameter, which can vary in space. This yields the Ginzburg-Landau equations, which describe VERY WELL the behavior of a superconductor close to the transition point. It was using these equations that Josephson was able to come up with the concept of the Josephson Effect (earning him a Nobel Prize). And Abrikosov was able to come up wit the idea of Type II superconductors and vortices (he also won a Nobel Prize for this work). And after the BCS theory was understood, Gor'kov was able to show that the Ginzburg-Landau equations are a limiting case of the BCS theory close to the critical point.

    However, the point of all this is that it was shown, before the microscopic BCS interactions were understood, scientists were able to do ALOT of things with the Ginzburg-Landau equations. What makes these so great is that they are able to approximate quantum mechanics decently, which the London equations were unable to do. And the best part is that scientists today (myself included) still use Ginzburg-Landau equations to model superconductors. It's just that much easier to use these equations for many interactions than the lower-level BCS theory. But amazingly, these equations were known BEFORE the BCS theory!

    So back to your comment, such study of critical phenomena teaches us a great deal about systems in criticality, even if the methods involved are decoupled from the microscopic physics. Especially if one can determine the universality class of an unknown system. And for very complicated systems, critical exponents will be difficult to determine analytically and must be solved numerically. Hence the importance of simulations and algorithms such as this.

    --

    make world, not war

  7. Re:Too dense by sadtrev · · Score: 2, Informative

    Water freezes at the surface of the lake because that's where it's coldest.

    The surface is in contact with cold air that takes away the the water's latent heat. In macroscopic terms this is by evaporation and by conduction across the thermal boundary layer of the air. Changes in the temperature of the ground are more gradual and will slow down freezing by supplying heat to the water at the bottom.

  8. More on water and ice density by jdfox · · Score: 2, Informative
    The physics of phase changes in water is just amazing. The water and ice we use everyday behaves that way, but amorphous ice is denser than liquid-phase water. Wikipedia has a good, concise primer on the physics of water.

    A friend of mine learned the hard way about how water expands as it freezes and its density drops. She put water inside glass Christmas ornaments, then put them in the freezer with the idea of floating pretty ornament-icecubes in her Christmas party punchbowl. She didn't leave any room for the expansion inside the ornaments though. So just as she was readying the hors d'oeuvres for the oven, she heard small explosions in her freezer, and cautiously opened the door to find ice and thin shards of glass all over everything.

    She didn't see the funny side of it at the time, but she does now. :)

  9. Re:just like non-alcoholic beer by cblguy · · Score: 2, Informative
    Spoken from a Budweiser-Miller-Coors experience, obviously.

    Go try a Fuller's London Porter. It has flavors of coffee and chocolate. And yes, it's beer.

    Coffee flavored beer is called a stout (and the lighter flavored version is a porter). Forget Guinness and Murphy's, they're far too watery. I'm talking a real stout. If you see one, grab a Great Divide Oak Aged Yeti Imperial Stout - it's a drink to behold... But warning, once you travel down the dark path of real stouts, you'll never drink a BMC beer again...

  10. The paper by Anonymous Coward · · Score: 1, Informative
    A. Chen, E.H. Chimowitz, S. De, and Y. Shapir, Phys. Rev. Lett. 95, 255701 (2005)

    The liquid-gas system is expected to exhibit distinct dynamic behavior in the fluid's critical region (model H). We present molecular dynamics simulations of a Lennard-Jones fluid model starting from specially designed, near-equilibrium, initial conditions. By following the fluid's relaxation towards equilibrium, we calculate the requisite transport coefficients in the critical region. The results yield the scaling behavior of the thermal diffusion coefficient DT~xi^-1.023+/-0.018 (xi is the correlation length) and a nonconventional divergent heat conductivity, all of which are in accord with mode-coupling and renormalization group predictions, as well as some experimental data.