Optimising the Termofluids CFD code for petascale simulations

ABSTRACT

This paper presents some recent efforts carried out on the expansion of the scalability of TermoFluids multi-physics Computational Fluid Dynamics (CFD) code, aiming to achieve petascale capacity for a single simulation. We describe different aspects that we have improved in our code in order to efficiently run it on 131,072 CPU-cores. This work has been developed using the BlueGene/Q Mira supercomputer of the Argonne Leadership Computing Facility, where we have obtained feedback at the targeted scale. In summary, this is a practical paper showing our experience at reaching the petascale paradigm for a single simulation with TermoFluids.     

Contributions of the different water sources to the Elqui river runoff (northern Chile) evaluated by H/O isotopes

We present the results of an isotope (2H and 18O) and hydrogeochemical study in order to constrain the origin, recharge, and evolution of the surface and groundwater in the arid Andean realm of the Elqui watershed. The results of 2H and18O analyses of water samples obtained during our summer and winter campaigns indicate a generally meteoric origin of the river and spring waters of the watershed. The isotope signature of water of the Elqui river and its tributaries as well as that of groundwater in the coastal region fits the 2H-18O relation of delta2H =7.61delta18O+6.1. A relatively fast discharge and a quasi-closed catchment area can be asserted for water along the river flow path. The tributaries from the more arid coastal area, north of the Elqui river, differ in their isotopic signature due to evaporation and hydrochemically due to interactions with the strongly altered and fractured volcanic rocks of the basement. In the Andean zone, the18O-enriched hydrothermal spring of Ba?os del Toro exhibits the influence of water-rock interaction processes. The chemistry of the river water changes from sulphate- to chloride-rich along the river course from the high Andean mountains to the coast. The sulphate-rich character of these Andean waters reflects their passage through sulphide-rich rock massifs that were subjected to strong oxidation processes in the near superficial environment. This sulphate signature is enforced by past and present mining of precious metal epithermal deposits (e.g. those of El Indio-Tambo Au-Cu-As district), in which mineralised zones were developed during a series of Miocene magmatic-hydrothermal episodes in the Andean realm. Owing to the proximity of the lower Elqui river waters and its tributaries to the Pacific coast, the chloride character may be induced by agricultural and marine (sea spray, fog) sources. Generally, the main source of the Elqui river water is mainly attributed to surface runoff and less to contributions from the basement fractured aquifer.     

Portable implementation model for CFD simulations. Application to hybrid CPU/GPU supercomputers

Nowadays, high performance computing (HPC) systems experience a disruptive moment with a variety of novel architectures and frameworks, without any clarity of which one is going to prevail. In this context, the portability of codes across different architectures is of major importance. This paper presents a portable implementation model based on an algebraic operational approach for direct numerical simulation (DNS) and large eddy simulation (LES) of incompressible turbulent flows using unstructured hybrid meshes. The strategy proposed consists in representing the whole time-integration algorithm using only three basic algebraic operations: sparse matrix–vector product, a linear combination of vectors and dot product. The main idea is based on decomposing the nonlinear operators into a concatenation of two SpMV operations. This provides high modularity and portability. An exhaustive analysis of the proposed implementation for hybrid CPU/GPU supercomputers has been conducted with tests using up to 128 GPUs. The main objective consists in understanding the challenges of implementing CFD codes on new architectures.     

Facile Formation of Redox‐Active Totally Organic Nanoparticles in Water by In Situ Reduction of Organic Precursors Stabilized through Aromatic–Aromatic Interactions by Aromatic Polyelectrolytes

The formation of redox‐active, totally organic nanoparticles in water is achieved following a strategy similar to that used to form metal nanoparticles. It is based on two fundamental concepts: i) complexation through aromatic–aromatic interactions of a water‐soluble precursor aromatic molecule with polyelectrolytes bearing complementary charged aromatic rings, and ii) reduction of the precursor molecule to achieve stabilized nanoparticles. Thus, formazan nanoparticles are synthesized by reduction of a tetrazolium salt with ascorbic acid using polyelectrolytes bearing benzene sulfonate residues of high linear aromatic density, but cannot be formed in the presence of nonaromatic polyelectrolytes. The red colored nanoparticles are efficiently encapsulated in calcium alginate beads, showing macroscopic homogeneity. Bleaching kinetics with chlorine show linear rates on the order of tenths of milli­meters per minute. A linear behavior of the dependence of the rate of bleaching on the chlorine concentration is found, showing the potential of the nanoparticles for chlorine sensing.

     

Resistivity of thin gold films on mica induced by electron-surface scattering: Application of quantitative scanning tunneling microscopy

We report a comparison between the resistivity measured on thin gold films deposited on mica, with predictions based upon classical theories of size effects (Drude's, Sondheimer's and Calecki's), as well as predictions based upon quantum theories of electron-surface scattering (the modified theory of Sheng, Xing and Wang, the theory of Tesanovic, Jaric and Maekawa, and that of Trivedi and Aschroft). From topographic images of the surface recorded with a Scanning Tunneling Microscope, we determined the rms roughness amplitude, δ and the lateral correlation length, ξ corresponding to a Gaussian representation of the average height-height autocorrelation function, describing the roughness of each sample in the scale of length set by the Fermi wave length. Using (δ, ξ) as input data, we present a rigorous comparison between resistivity data and predictions based upon the theory of Calecki as well as quantum theoretical predictions without adjustable parameters. The resistivity was measured on gold films of different thickness evaporated onto mica substrates, between 4 K and 300 K. The resistivity data covers the range 0.1 < x(T) < 6.8, for 4 K < T < 300 K, where x(T) is the ratio between film thickness and electron mean free path in the bulk at temperature T. We experimentally identify electron-surface and electron-phonon scattering as the microscopic electron scattering mechanisms giving rise to the macroscopic resistivity. The different theories are all capable of estimating the thin film resistivity to an accuracy better than 10%; however the mean free path and the resistivity characterizing the bulk turn out to depend on film thickness. Surprisingly, only the Sondheimer theory and its quantum version, the modified theory of Sheng, Xing and Wang, predict and increase in resistivity induced by size effects that seems consistent with published galvanomagnetic phenomena also arising from electron-surface scattering measured at low temperatures.     

Polycaprolactone and poly-β-cyclodextrin polymer blend: a Biopolymers composite film for drug release

Nowadays, biomedical films containing drug carriers are preferred over conventional ones, since the protection of the injury and the therapy is joined within a single device. In the current work, we prepared polycaprolactone (PCL) composite films with β-cyclodextrin (βCD) or its epichlorohydrin crosslinked polymer (βCDP) as ibuprofen (Ibu) drug carrier. The composite films were prepared at different PCL/additive ratios (2, 5, 10 and 20 wt%). ATR-FTIR spectroscopy and water contact angle (WCA) measurements indicated a scarce presence of the additives on the surface. Cross-section scanning electron micrographs showed the presence of aggregates corresponding to βCD and βCDP in the inner regions of the films. The incorporation of βCD and βCDP into the PCL films did not affect their thermal properties as was determined from differential scanning calorimetry (DSC). PCL-films with 10 wt% of the inclusion complexes Ibu@βCD and Ibu@βCDP were prepared and the release studies were performed. At pH?=?7.2, PCL-Ibu@βCDP composite film released 55% of Ibu within the first six hours; eight times the amount released by PCL-Ibu@βCD within the same time interval. A plausible mechanism for ibuprofen release is discussed based on the cross-section SEM micrographs of composite films.