Taxila LBM: a parallel,modular lattice Boltzmann framework for simulating pore-scale flow in porous media

The lattice Boltzmann method is a popular tool for pore-scale simulation of flow. This is likely due to the ease of including complex geometries such as porous media and representing multiphase and multifluid flows. Many advancements, including multiple relaxation times, increased isotropy, and others have improved the accuracy and physical fidelity of the method. Additionally, the lattice Bolzmann method is computationally very efficient, thanks to the explicit nature of the algorithm and relatively large amount of local work. The combination of many algorithmic options and efficiency means that a software framework enabling the usage and comparison of these advancements on computers from laptops to large clusters has much to offer. In this paper, we introduce Taxila LBM, an open-source software framework for lattice Boltzmann simulations. We discuss the design of the framework and lay out the features available, including both methods in the literature and a few new enhancements which generalize methods to complex geometries. We discuss the trade-off of accuracy and performance in various methods, noting how the Taxila LBM makes it easy to perform these comparisons for real problems. And finally, we demonstrate a few common applications in pore-scale simulation, including the characterization of permeability of a Berea sandstone and analysis of multifluid flow in heterogenous micromodels.     

Dedication: Dr. Scott W. Nixon (1943–2012)

The availability of suspended sediments will be a dominant factor influencing the stability of tidal wetlands as sea levels rise. Watershed-derived sediments are a critical source of material supporting accretion in many tidal wetlands, and recent declines in wetland extent in several large river delta systems have been attributed in part to declines in sediment delivery. Little attention has been given, however, to changes in sediment supply outside of large river deltas. In this study, significant declines in suspended sediment concentrations (SSCs) over time were observed for 25 of 61 rivers examined that drain to the East and Gulf Coasts of the USA. Declines in fluvial SSC were significantly correlated with increasing water retention behind dams, indicating that human activities play a role in declining sediment delivery. There was a regional pattern to changes in fluvial sediment, and declines in SSCs were also significantly related to rates of relative sea level rise (RSLR) along the coast, such that wetlands experiencing greater RSLR also tend to be receiving less fluvial sediment. Tidal wetlands in the Mid-Atlantic, Mississippi River Delta, and Texas Gulf especially may become increasingly vulnerable due to rapid RSLR and reductions in sediment. These results also indicate that past rates of marsh accretion may not be indicative of potential future accretion due to changes in sediment availability. Declining watershed sediment delivery to the coastal zone will limit the ability of tidal marshes to keep pace with rising sea levels in some coastal systems.     

A 2000-yr reconstruction of air temperature in the Great Basin of the United States with specific reference to the Medieval Climatic Anomaly

A sediment core representing the past two millennia was recovered from Stella Lake in the Snake Range of the central Great Basin in Nevada. The core was analyzed for sub-fossil chironomids and sediment organic content. A quantitative reconstruction of mean July air temperature (MJAT) was developed using a regional training set and a chironomid-based WA-PLS inference model (r²_jack = 0.55, RMSEP = 0.9°C). The chironomid-based MJAT reconstruction suggests that the interval between AD 900 and AD 1300, corresponding to the Medieval Climate Anomaly (MCA), was characterized by MJAT elevated 1.0°C above the subsequent Little Ice Age (LIA), but likely not as warm as recent conditions. Comparison of the Stella Lake temperature reconstruction to previously published paleoclimate records from this region indicates that the temperature fluctuations inferred to have occurred at Stella Lake between AD 900 and AD 1300 correspond to regional records documenting hydroclimate variability during the MCA interval. The Stella Lake record provides evidence that elevated summer temperature contributed to the increased aridity that characterized the western United States during the MCA.     

A temperature-induced reversible transformation between paratacamite and herbertsmithite

The crystal chemistry of paratacamite has been re-evaluated by studying a crystal from the holotype specimen BM86958 of composition Cu_3.71Zn_0.29(OH)₆Cl₂ using single-crystal X-ray diffraction at 100, 200, 300, 353, 393 and 423 K. At 300 K paratacamite has space group $R\bar{3}$ with unit-cell parameters a 13.644 and c 14.035 Å and exhibits a pronounced subcell, a′ = ½a and c′ = c, analogous to that of the closely related mineral herbertsmithite, Cu₃Zn(OH)₆Cl₂. Between 353 and 393 K, paratacamite undergoes a reversible phase transformation to the herbertsmithite-like substructure, space group $R\bar{3}m$ , unit-cell parameters a 6.839 and c 14.072 Å (393 K). The transformation is characterised by a gradual reduction in intensity of superlattice reflections, which are absent at 393 and 443 K. On cooling from 443 to 300 K at ~10 K min^?1, the superlattice reflections reappear and the refined structures ( $R\bar{3}$ ) of the initial and recovered 300 K states are almost identical. The complete reversibility of the transformation establishes that paratacamite of composition Cu_3.71Zn_0.29(OH)₆Cl₂ is thermodynamically stable at ambient temperatures. The nature of the rhombic distortion of the M(2)O₆ octahedron is discussed by considering two possibilities that are dependent upon the nature of cation substitution in the interlayer sites.     

Strategically Locating Wildlife Crossing Structures for Florida Panthers Using Maximal Covering Approaches

Crossing structures are an effective method for mitigating habitat fragmentation and reducing wildlife‐vehicle collisions, although high construction costs limit the number that can be implemented in practice. Therefore, optimizing the placement of crossing structures in road networks is suggested as a strategic conservation planning method. This research explores two approaches for using the maximal covering location problem (MCLP) to determine optimal sites to install new wildlife crossing structures. The first approach is based on records of traffic mortality, while the second uses animal tracking data for the species of interest. The objective of the first is to cover the maximum number of collision sites, given a specified number of proposed structures to build, while the second covers as many animal tracking locations as possible under a similar scenario. These two approaches were used to locate potential wildlife crossing structures for endangered Florida panthers (Puma concolor coryi) in Collier, Lee, and Hendry Counties, Florida, a population whose survival is threatened by excessive traffic mortality. Historical traffic mortality records and an extensive radio‐tracking dataset were used in the analyses. Although the two approaches largely select different sites for crossing structures, both models highlight key locations in the landscape where these structures can remedy traffic mortality and habitat fragmentation. These applications demonstrate how the MCLP can serve as a useful conservation planning tool when traffic mortality or animal tracking data are available to researchers.     

Evolution of the mass function of dark matter haloes

We use a high-resolution ΛCDM numerical simulation to calculate the mass function of dark matter haloes down to the scale of dwarf galaxies, back to a redshift of 15, in a  50 h ⁻¹ Mpc  volume containing 80 million particles. Our low-redshift results allow us to probe low-σ density fluctuations significantly beyond the range of previous cosmological simulations. The Sheth & Tormen mass function provides an excellent match to all of our data except for redshifts of 10 and higher, where it overpredicts halo numbers increasingly with redshift, reaching roughly 50 per cent for the  10¹⁰–10¹¹ M_⊙  haloes sampled at redshift 15. Our results confirm previous findings that the simulated halo mass function can be described solely by the variance of the mass distribution, and thus has no explicit redshift dependence. We provide an empirical fit to our data that corrects for the overprediction of extremely rare objects by the Sheth & Tormen mass function. This overprediction has implications for studies that use the number densities of similarly rare objects as cosmological probes. For example, the number density of high-redshift  ( z ≃ 6) QSOs  , which are thought to be hosted by haloes at 5σ peaks in the fluctuation field, are likely to be overpredicted by at least a factor of 50 per cent. We test the sensitivity of our results to force accuracy, starting redshift and halo-finding algorithm.     

Seasonal change in the deep atmosphere of Uranus

We present high-resolution radio maps of Uranus, made from data collected in 1994 at wavelengths of 2 and 6 cm, which show large-scale changes occurring deep and rapidly in the troposphere. Brightness features in these maps are significantly different from those observed throughout the 1980's. These differences are not due to the changing viewing geometry, but result from atmospheric changes in the 5 to 50 bar region. All the observations show strong latitudinal variations in absorber abundance and/or temperature, causing the South Pole to appear brighter than lower latitudes. The transition between bright pole and darker latitudes is always near −45°, but between 1989 and 1994 the contrast between the regions increased significantly. This suggests that the large-scale circulation in the upper 50 bars of the uranian Southern Hemisphere changed. Older, disk-averaged microwave observations have suggested that seasonal variability occurs, but these new maps are the first to provide detailed timing and location information which can be used to test dynamical models.     

Understanding geomorphic responses to environmental change: a 19 000‐year case study from semi‐arid central Texas,USA

Understanding the impact of past climatic changes on landscape stability is crucial in order to predict and mitigate the effects of future changes. However, in arid and semi‐arid environments, reconstructions are often hampered by a poor understanding of the relationship between sediment deposition and climate. We present here data from central Texas, a region that is sensitive to environmental change but has received relatively little attention. The study integrates a chronology of 29 optically stimulated luminescence (OSL) ages from six sedimentary sites in a range of depositional contexts with a 19 000 a climate record derived from pollen extracted from the adjacent Boriack Bog. By comparing the two records, we aimed to assess the relationship between climate change and geomorphic activity. Data show that extensive aeolian and colluvial deposition occurred during the mid to late Holocene, with sedimentation generally increasing during more arid phases. However, a number of depositional events on slopes were associated with moister episodes, and sediment was also deposited in a summit setting immediately after phases of increased precipitation. Linkages between climate change and geomorphic response were therefore complex, being controlled by both sediment supply and transport energy. The climatic fluctuations identified in the Boriack Bog record highlight the sensitivity of central Texas to environmental change, while disparities in conditions recorded between it and other palaeoenvironmental sites in the southern USA emphasise the need for further work in order to enhance understanding of landscape sedimentary response to climate. Copyright © 2010 John Wiley & Sons, Ltd.