Permafrost creates a polygonal landscape, irregularity that makes simulating thawing’s impact on climate change a challenge requiring advanced algorithms and high-performance computers. (Photo: Konstanze Piel, Alfred Wegener Institute.)

Permafrost – soil that stays frozen year-round – comprises about one-quarter of the Northern Hemisphere’s land. For now.

Climate changes are turning permafrost into mush over an increasing area. Besides creating soggy soil, this melting generates warming that drives further climate change.

No one is sure, however, how much impact this feedback loop has. Finding the answer demands high-performance computing, advanced algorithms and, as it turns out, some unexpected interactions.

That includes the intersecting careers of two Department of Energy Computational Science Graduate Fellowship (DOE CSGF) recipients. The first was Glenn Hammond, now a computational geohydrologist at Sandia National Laboratories in New Mexico. As a graduate student, Hammond served a 1999 practicum with Peter Lichtner, a multiphase flow and reactive-transport modeling expert at Los Alamos National Laboratory. Hammond’s project involved building PFLOTRAN, a parallel-processing version of the existing FLOTRAN code that models how groundwater reacts with the soil and rock it flows through. Parallel processing parcels pieces of a computational task out to multiple computer processors that work simultaneously to quickly complete the job.

“When you take existing code and try to parallelize it, it can be a bit of a challenge,” Hammond says. “As a new student learning parallelization techniques and new code, I started developing my own code from scratch, modeled very similar to FLOTRAN.”

Later, Richard Mills, who recently joined Intel Corp. after working as a computational scientist at Oak Ridge National Laboratory, spent his 2003 DOE CSGF practicum working under Lichtner to continue Hammond’s PFLOTRAN work. Mills implemented the code’s initial solver for groundwater flow.

Many other scientists also contributed to PFLOTRAN over the years and, Mills says, “it turned into a fairly comprehensive code for modeling surface and subsurface hydrology.” In other words, it simulates water and the associated biogeochemical reactions that occur as it moves over or through soil.

Coding the Equations

Although the code is designed to run in parallel on high-performance computers, Mills says the team wanted to make PFLOTRAN flexible enough to run on both the biggest supercomputers and on laptops.

Hammond and Mills grew up in a C and C++ coding world, but most scientists who study hydrology use an older language, FORTRAN. The two former fellows followed that tradition for PFLOTRAN but enlisted modern FORTRAN features to add a degree of object-oriented structure more characteristic of C or C++ code.

(Visited 1,916 times, 1 visits today)

Page: 1 2

Published by

Recent Posts

Optimized for discovery

Stefan Wild came to Argonne National Laboratory in the summer of 2006 for what was… Read More

October, 2021

Happy returns

Sandia’s Humberto Silva III applies uncertainty quantification to model smooth atmospheric reentry from space. Read More

September, 2021

Crafting community

The Department of Energy Computational Science Graduate Fellowship celebrates 30 years of cultivating leaders and… Read More

August, 2021

New podcast features labs and fellows

Series celebrates DOE Computational Science Fellowship's 30th anniversary. Read More

July, 2021

Extreme 3-D

Argonne’s Advanced Photon Source upgrades to large-sample 3-D imaging beyond the depth of field –… Read More

June, 2021

Targeting tumors

Cancer biology gets the supercomputing treatment on Oak Ridge's Summit and Lawrence Livermore's Sierra. Read More

April, 2021