Towards effective indexing for very large video sequence database

With rapid advances in video processing technologies and ever fast increments in network bandwidth, the popularity of video content publishing and sharing has made similarity search an indispensable operation to retrieve videos of user interests. The video similarity is usually measured by the percentage of similar frames shared by two video sequences, and each frame is typically represented as a high-dimensional feature vector. Unfortunately, high complexity of video content has posed the following major challenges for fast retrieval: (a) effective and compact video representations, (b) efficient similarity measurements, and (c) efficient indexing on the compact representations. In this paper, we propose a number of methods to achieve fast similarity search for very large video database. First, each video sequence is summarized into a small number of clusters, each of which contains similar frames and is represented by a novel compact model called Video Triplet (ViTri). ViTri models a cluster as a tightly bounded hypersphere described by its position, radius, and density. The ViTri similarity is measured by the volume of intersection between two hyperspheres multiplying the minimal density, i.e., the estimated number of similar frames shared by two clusters. The total number of similar frames is then estimated to derive the overall similarity between two video sequences. Hence the time complexity of video similarity measure can be reduced greatly. To further reduce the number of similarity computations on ViTris, we introduce a new one dimensional transformation technique which rotates and shifts the original axis system using PCA in such a way that the original inter-distance between two high-dimensional vectors can be maximally retained after mapping. An efficient B+-tree is then built on the transformed one dimensional values of ViTris' positions. Such a transformation enables B+-tree to achieve its optimal performance by quickly filtering a large portion of non-similar ViTris. Our extensive experiments on real large video datasets prove the effectiveness of our proposals that outperform existing methods significantly.

Characterization of a volume-sensitive chloride current in rat osteoblast-like (ROS 17/2.8) cells

1. During osmotic swelling, cultured osteoblastic cells (ROS 17/2.8) exhibited activation of large amplitude Cl- currents in the whole-cell configuration of the patch-clamp technique. Effects of hypotonic shock on cell volume and membrane conductance were rapidly reversed on return to isotonic conditions. 2. Voltage command pulses in the range -80 to +50 mV produce instantaneous activation of Cl- currents. At potentials more positive than +50 mV the current exhibited time-dependent inactivation. The instantaneous current-voltage relationship was outwardly rectifying. 3. The anion permeability sequence of the induced current was SCN- (2.2) > I- (1.9) > Br- (1.5) > Cl- (1.0) > F- (0.8) > gluconate- (0.2). This corresponds to Eisenman's sequence I. 4. The volume-sensitive Cl- current was effectively inhibited by the Cl- channel blockers 4,4'-diisothiocyanatostilbene-2,2-disulphonic acid (DIDS) and 5-nitro-2-(3-phenylpropylamino) benzoic acid (NPPB). Outward currents were more effectively suppressed by DIDS than inward currents. The concentrations for 50% inhibition (IC50) of outward and inward currents were 81 and 298 μM, respectively. NPPB was equally effective at inhibiting outward and inward currents (IC50 of 64 μM). The current was relatively insensitive to diphenylamine-2-carboxylate (DPC), 500 μM producing only 22.5 ± 4.0% inhibition. 5. Inhibitors of protein kinase A (H-89, 1 μM) and tyrosine kinase (tyrphostin A25, 200 μM) were without effect upon activation of Cl- currents in response to hypotonic shock. Under isotonic conditions, elevation of intracellular Ca2+ by ionomycin (1 μM) or activation of protein kinase C by 12-O-tetradecanoylphorbol 13-acetate (TPA, 0.1 μM) failed to evoke increases in basal Cl- conductance levels. 6. It is concluded that an outwardly rectifying Cl- conductance is activated upon osmotic swelling and may be involved in cell volume regulation of ROS 17/2.8 cells.

The prediction of vibration in large electric machines

This thesis reports the development of a reliable method for the prediction of response to electromagnetically induced vibration in large electric machines. The machines of primary interest are DC ship-propulsion motors but much of the work reported has broader significance. The investigation has involved work in five principal areas. (1) The development and use of dynamic substructuring methods. (2) The development of special elements to represent individual machine components. (3) Laboratory scale investigations to establish empirical values for properties which affect machine vibration levels. (4) Experiments on machines on the factory test-bed to provide data for correlation with prediction. (5) Reasoning with regard to the effect of various design features. The limiting factor in producing good models for machines in vibration is the time required for an analysis to take place. Dynamic substructuring methods were adopted early in the project to maximise the efficiency of the analysis. A review of existing substructure- representation and composite-structure assembly methods includes comments on which are most suitable for this application. In three appendices to the main volume methods are presented which were developed by the author to accelerate analyses. Despite significant advances in this area, the limiting factor in machine analyses is still time. The representation of individual machine components was addressed as another means by which the time required for an analysis could be reduced. This has resulted in the development of special elements which are more efficient than their finite-element counterparts. The laboratory scale experiments reported were undertaken to establish empirical values for the properties of three distinct features - lamination stacks, bolted-flange joints in rings and cylinders and the shimmed pole-yoke joint. These are central to the preparation of an accurate machine model. The theoretical methods are tested numerically and correlated with tests on two machines (running and static). A system has been devised with which the general electromagnetic forcing may be split into its most fundamental components. This is used to draw some conclusions about the probable effects of various design features.

Assembly accuracy analysis for small components with a planar surface in large-scale metrology

Large-scale mechanical products, such as aircraft and rockets, consist of large numbers of small components, which introduce additional difficulty for assembly accuracy and error estimation. Planar surfaces as key product characteristics are usually utilised for positioning small components in the assembly process. This paper focuses on assembly accuracy analysis of small components with planar surfaces in large-scale volume products. To evaluate the accuracy of the assembly system, an error propagation model for measurement error and fixture error is proposed, based on the assumption that all errors are normally distributed. In this model, the general coordinate vector is adopted to represent the position of the components. The error transmission functions are simplified into a linear model, and the coordinates of the reference points are composed by theoretical value and random error. The installation of a Head-Up Display is taken as an example to analyse the assembly error of small components based on the propagation model. The result shows that the final coordination accuracy is mainly determined by measurement error of the planar surface in small components. To reduce the uncertainty of the plane measurement, an evaluation index of measurement strategy is presented. This index reflects the distribution of the sampling point set and can be calculated by an inertia moment matrix. Finally, a practical application is introduced for validating the evaluation index.

Spatial scale and abundance patterns of large fish communities in freshwater marshes of the Florida Everglades

Anthropogenic habitat alterations and water-management practices have imposed an artificial spatial scale onto the once contiguous freshwater marshes of the Florida Everglades. To gain insight into how these changes may affect biotic communities, we examined whether variation in the abundance and community structure of large fishes (SL . 8 cm) in Everglades marshes varied more at regional or intraregional scales, and whether this variation was related to hydroperiod, water depth, floating mat volume, and vegetation density. From October 1997 to October 2002, we used an airboat electrofisher to sample large fishes at sites within three regions of the Everglades. Each of these regions is subject to unique watermanagement schedules. Dry-down events (water depth , 10 cm) occurred at several sites during spring in 1999, 2000, 2001, and 2002. The 2001 dry-down event was the most severe and widespread. Abundance of several fishes decreased significantly through time, and the number of days post-dry-down covaried significantly with abundance for several species. Processes operating at the regional scale appear to play important roles in regulating large fishes. The most pronounced patterns in abundance and community structure occurred at the regional scale, and the effect size for region was greater than the effect size for sites nested within region for abundance of all species combined, all predators combined, and each of the seven most abundant species. Non-metric multi-dimensional scaling revealed distinct groupings of sites corresponding to the three regions. We also found significant variation in community structure through time that correlated with the number of days post-dry-down. Our results suggest that hydroperiod and water management at the regional scale influence large fish communities of Everglades marshes.

A Hydrological Budget (2002–2008) for a Large Subtropical Wetland Ecosystem Indicates Marine Groundwater Discharge Accompanies Diminished Freshwater Flow

Water budget parameters are estimated for Shark River Slough (SRS), the main drainage within Everglades National Park (ENP) from 2002 to 2008. Inputs to the water budget include surface water inflows and precipitation while outputs consist of evapotranspiration, discharge to the Gulf of Mexico and seepage losses due to municipal wellfield extraction. The daily change in volume of SRS is equated to the difference between input and outputs yielding a residual term consisting of component errors and net groundwater exchange. Results predict significant net groundwater discharge to the SRS peaking in June and positively correlated with surface water salinity at the mangrove ecotone, lagging by 1 month. Precipitation, the largest input to the SRS, is offset by ET (the largest output); thereby highlighting the importance of increasing fresh water inflows into ENP for maintaining conditions in terrestrial, estuarine, and marine ecosystems of South Florida.

Fig. 2. Lithology, water and total carbon content, cumulative grain-size distribution, and quartz, feldspar, and mica contents of core sequence Lz1013 from south-western Beaver Lake

A 100 cm long sediment sequence was recovered from Beaver Lake in Amery Oasis, East Antarctica, using gravity and piston corers. Sedimentological and mineralogical analyses and the absence of micro and macrofossils indicate that the sediments at the base of the sequence formed under glacial conditions, probably prior to c. 12 500 cal. yr BP. The sediments between c. 81 and 31 cm depth probably formed under subaerial conditions, indicating that isostatic uplift since deglaciation has been substantially less than eustatic sea-level rise and that large areas of the present-day floor of Beaver Lake must have been subaerially exposed following deglaciation. The upper 31 cm of the sediment sequence were deposited under glaciomarine conditions similar to those of today, supporting geomorphic observations that the Holocene was a period of relative sea-level highstand in Amery Oasis.

Oceanographic investiations from the Sea of Ochotsk

The Sea of Okhotsk is a marginal sea of the Pacific Ocean, which is characterized by strong variations in the productivity and sediment supply due to sea ice transport and river input. Furthermore the variations in the hydrological cycle determine the formation of the SOIW (Sea of Okhotsk Intermediate Water) which plays an important role in the ventilation processes in the intermediate water of the N-Pacific. Isotope data measured on planktonic and benthic foraminifera, sedimentological and geochemical studies of sediment cores and surface samples from the Sea of Okhotsk are used to reconstruct the paleoceanography during the past 350.000 years. The dating and correlation of the sediments are based on oxygen isotope stratigraphy, absolute ages, magnetic susceptibility as well as a detailled tephrachronology of the entire basin. The sedimentation rates are characterized by temporal and spatial variations. The maximum sedimentation rate takes place at the continental slope off Sakhalin due to the input of the Amur River, the sea ice drift and the high productivity. The sedimentation rate in the eastern part of the Sea of Okhotsk is generelly high because of the influence of the nutrient-rich Kamchatka Current. In the central and northern parts of the Sea of Okhotsk, areas with low productivity and reduced terrestrial supply, the sedimentation rate is the lowest. The analyses of the surface sediment samples make it possible to characterize the (sub)- recent sediment supply and transportation processes. The bulk sediment measurements, isotope data and the accumulation rate of ice-rafted debris (IRD) show a dominant sea ice cover and a region with a high productivity as well as a high Amur River input in the western part of the sea. The eastern part of the Sea of Okhotsk, however, is marked by the predominance of warm and nutrient-rich water masses coming from the Kamchatka Current which restricts the sea ice cover. This is reflected in low content of ice-rafted debris and high productivity proxies as well as in isotope data. The deposits of the Sea of Okhotsk are characterized by terrestrial, biogenic and volcanogenic sediment input which varies temporally and spatially. Here, the sedimentation pattern is dominated by the terrestrial input. Bulk sediment measurements and sample analyses of the > 63 micron particle input make it possible to distinguish glacial and interglacial fluctuations. The sedimentation processes during glacial times are determined by a high content of ice-rafted debris, whereas the primary production is higher during interglacial periods. During the last glacial/interglacial cycle the IRD-distribution pattern indicates a strong sea ice transport in the western part and in large areas of the open sea in the eastern part of the Sea of Okhotsk with a relatively constant ice-drift system. The IRD flux in sediments of the oxygen isotope Stage 6 reflects a new sedimentation pattern in the eastern part of the sea. This high IRD accumulation rate indicates ice advances beyond the shelf margin and an iceberg transport from NE-E direction into the Sea of Okhotsk. The several large, brief, negative anomalies in d13C values of Neogloboquadrina pachyderma (s) show releases of methane from basin sediments which correspond to periods of relative sea level falls. The high sedimentation rates on the Sakhalin slope allow insights into the climatic history in Holocene and indicate shorter-scale variations oscillation in Stage 3, which correlate with the global climatic changes. These variations are described as Dansgaard-Oeschger cycles in Greenland ice cores and as Heinrich-Events in several marine sediment cores from the N-Atlantic.

Occurrence of microtextures recorded from quartz sand grains of ODP Hole 178-1101A (Table 1)

Sediment drifts on the continental rise are located proximal to the western side of the Antarctic Peninsula and recorded changes in glacial volume and thermal regime over the last ca. 15 m.y. At Ocean Drilling Program (ODP) Site 1101 (Leg 178), which recovered sediments back to 3.1 Ma, glacial-interglacial cyclicity was identified based on the biogenic component and sedimentary structures observed in X-radiographs, magnetic susceptibility and lithofacies descriptions. Glacial intervals are dominated by fine-grained laminated mud and interglacial units consist of bioturbated muds enriched in biogenic components. From 2.2 to 0.76 Ma, planktonic foraminifera and calcareous nannofossils dominate in the interglacials suggesting a shift of the Antarctic Polar Front (APF) to the south near the drifts. Prior to 2.2 Ma, cyclicity cannot be identified and diatoms dominate the biogenic component and high percent opal suggests warmer conditions south of the APF and reduced sea ice over the drifts. Analyses of the coarse-grained terrigenous fraction (pebbles and coarse sand) from Sites 1096 and 1101 record glaciers at sea-level releasing iceberg-rafted debris (IRD) throughout the last 3.1 m.y. Analyses of quartz sand grains in IRD with the scanning electron microscope (SEM) show an abrupt change in the frequency of occurrence of microtextures at ~1.35 Ma. During the Late Pliocene to Early Pleistocene, the population of quartz grains included completely weathered grains and a low frequency of crushing and abrasion, suggesting that glaciers were small and did not inundate the topography. Debris shed from mountain peaks was transported supraglacially or englacially allowing weathered grains to pass through the glacier unmodified. During glacial periods from 1.35-0.76 Ma, glaciers expanded in size. The IRD flux was very high and dropstones have diverse lithologies. Conditions resembling those at the Last Glacial Maximum (LGM) have been episodically present on the Antarctic Peninsula since ~0.76 Ma. Quartz sand grains show high relief, fracture and abrasion common under thick ice and the IRD flux is low with a more restricted range of dropstone lithologies.

Feasibility of MR-Based Body Composition Analysis in Large Scale Population Studies

Introduction Quantitative and accurate measurements of fat and muscle in the body are important for prevention and diagnosis of diseases related to obesity and muscle degeneration. Manually segmenting muscle and fat compartments in MR body-images is laborious and time-consuming, hindering implementation in large cohorts. In the present study, the feasibility and success-rate of a Dixon-based MR scan followed by an intensity-normalised, non-rigid, multi-atlas based segmentation was investigated in a cohort of 3,000 subjects. Materials and Methods 3,000 participants in the in-depth phenotyping arm of the UK Biobank imaging study underwent a comprehensive MR examination. All subjects were scanned using a 1.5 T MR-scanner with the dual-echo Dixon Vibe protocol, covering neck to knees. Subjects were scanned with six slabs in supine position, without localizer. Automated body composition analysis was performed using the AMRA Profiler™ system, to segment and quantify visceral adipose tissue (VAT), abdominal subcutaneous adipose tissue (ASAT) and thigh muscles. Technical quality assurance was performed and a standard set of acceptance/rejection criteria was established. Descriptive statistics were calculated for all volume measurements and quality assurance metrics. Results Of the 3,000 subjects, 2,995 (99.83%) were analysable for body fat, 2,828 (94.27%) were analysable when body fat and one thigh was included, and 2,775 (92.50%) were fully analysable for body fat and both thigh muscles. Reasons for not being able to analyse datasets were mainly due to missing slabs in the acquisition, or patient positioned so that large parts of the volume was outside of the field-of-view. Discussion and Conclusions In conclusion, this study showed that the rapid UK Biobank MR-protocol was well tolerated by most subjects and sufficiently robust to achieve very high success-rate for body composition analysis. This research has been conducted using the UK Biobank Resource.

Second Quarter 2016: Slowdown for Large Hotels Continues; Small Hotels Have Now Slowed as Well

Our Standardized Unexpected Price (SUP) metric continues to show a decline in the price of large hotels, and now also the price of small hotels has eased—even though hotel transaction volume has increased. Although debt and equity financing for hotels remain relatively inexpensive, we are concerned that the total volatility of hotel returns is greater relative to the return volatility for other commercial real estate. If this trend continues, lenders will eventually start to tighten hotel lending standards. Our early warning indicators all continue to suggest that the downward trend in hotel prices should continue into the next quarter. This is report number 19 of the index series.

A semi-Lagrangian finite volume method for solving viscoelastic flow problems

Recently, there has been considerable interest in solving viscoelastic problems in 3D particularly with the improvement in modern computing power. In many applications the emphasis has been on economical algorithms which can cope with the extra complexity that the third dimension brings. Storage and computer time are of the essence. The advantage of the finite volume formulation is that a large amount of memory space is not required. Iterative methods rather than direct methods can be used to solve the resulting linear systems efficiently.

Development of ultra-high-performance concrete (UHPC) using waste glass materials ─ towards innovative eco-friendly concrete

Le béton conventionnel (BC) a de nombreux problèmes tels que la corrosion de l’acier d'armature et les faibles résistances des constructions en béton. Par conséquent, la plupart des structures fabriquées avec du BC exigent une maintenance fréquent. Le béton fibré à ultra-hautes performances (BFUP) peut être conçu pour éliminer certaines des faiblesses caractéristiques du BC. Le BFUP est défini à travers le monde comme un béton ayant des propriétés mécaniques, de ductilité et de durabilité supérieures. Le BFUP classique comprend entre 800 kg/m³ et 1000 kg/m³ de ciment, de 25 à 35% massique (%m) de fumée de silice (FS), de 0 à 40%m de poudre de quartz (PQ) et 110-140%m de sable de quartz (SQ) (les pourcentages massiques sont basés sur la masse totale en ciment des mélanges). Le BFUP contient des fibres d'acier pour améliorer sa ductilité et sa résistance aux efforts de traction. Les quantités importantes de ciment utilisées pour produire un BFUP affectent non seulement les coûts de production et la consommation de ressources naturelles comme le calcaire, l'argile, le charbon et l'énergie électrique, mais affectent également négativement les dommages sur l'environnement en raison de la production substantielle de gaz à effet de serre dont le gas carbonique (CO[indice inférieur 2]). Par ailleurs, la distribution granulométrique du ciment présente des vides microscopiques qui peuvent être remplis avec des matières plus fines telles que la FS. Par contre, une grande quantité de FS est nécessaire pour combler ces vides uniquement avec de la FS (25 à 30%m du ciment) ce qui engendre des coûts élevés puisqu’il s’agit d’une ressource limitée. Aussi, la FS diminue de manière significative l’ouvrabilité des BFUP en raison de sa surface spécifique Blaine élevée. L’utilisation du PQ et du SQ est également coûteuse et consomme des ressources naturelles importantes. D’ailleurs, les PQ et SQ sont considérés comme des obstacles pour l’utilisation des BFUP à grande échelle dans le marché du béton, car ils ne parviennent pas à satisfaire les exigences environnementales. D’ailleurs, un rapport d'Environnement Canada stipule que le quartz provoque des dommages environnementaux immédiats et à long terme en raison de son effet biologique. Le BFUP est généralement vendu sur le marché comme un produit préemballé, ce qui limite les modifications de conception par l'utilisateur. Il est normalement transporté sur de longues distances, contrairement aux composantes des BC. Ceci contribue également à la génération de gaz à effet de serre et conduit à un coût plus élevé du produit final. Par conséquent, il existe le besoin de développer d’autres matériaux disponibles localement ayant des fonctions similaires pour remplacer partiellement ou totalement la fumée de silice, le sable de quartz ou la poudre de quartz, et donc de réduire la teneur en ciment dans BFUP, tout en ayant des propriétés comparables ou meilleures. De grandes quantités de déchets verre ne peuvent pas être recyclées en raison de leur fragilité, de leur couleur, ou des coûts élevés de recyclage. La plupart des déchets de verre vont dans les sites d'enfouissement, ce qui est indésirable puisqu’il s’agit d’un matériau non biodégradable et donc moins respectueux de l'environnement. Au cours des dernières années, des études ont été réalisées afin d’utiliser des déchets de verre comme ajout cimentaire alternatif (ACA) ou comme granulats ultrafins dans le béton, en fonction de la distribution granulométrique et de la composition chimique de ceux-ci. Cette thèse présente un nouveau type de béton écologique à base de déchets de verre à ultra-hautes performances (BEVUP) développé à l'Université de Sherbrooke. Les bétons ont été conçus à l’aide de déchets verre de particules de tailles variées et de l’optimisation granulaire de la des matrices granulaires et cimentaires. Les BEVUP peuvent être conçus avec une quantité réduite de ciment (400 à 800 kg/m³), de FS (50 à 220 kg/m³), de PQ (0 à 400 kg/m³), et de SQ (0-1200 kg/m³), tout en intégrant divers produits de déchets de verre: du sable de verre (SV) (0-1200 kg/m³) ayant un diamètre moyen (d[indice inférieur 50]) de 275 µm, une grande quantité de poudre de verre (PV) (200-700 kg/m³) ayant un d50 de 11 µm, une teneur modérée de poudre de verre fine (PVF) (50-200 kg/m³) avec d[indice inférieur] 50 de 3,8 µm. Le BEVUP contient également des fibres d'acier (pour augmenter la résistance à la traction et améliorer la ductilité), du superplastifiants (10-60 kg/m³) ainsi qu’un rapport eau-liant (E/L) aussi bas que celui de BFUP. Le remplacement du ciment et des particules de FS avec des particules de verre non-absorbantes et lisse améliore la rhéologie des BEVUP. De plus, l’utilisation de la PVF en remplacement de la FS réduit la surface spécifique totale nette d’un mélange de FS et de PVF. Puisque la surface spécifique nette des particules diminue, la quantité d’eau nécessaire pour lubrifier les surfaces des particules est moindre, ce qui permet d’obtenir un affaissement supérieur pour un même E/L. Aussi, l'utilisation de déchets de verre dans le béton abaisse la chaleur cumulative d'hydratation, ce qui contribue à minimiser le retrait de fissuration potentiel. En fonction de la composition des BEVUP et de la température de cure, ce type de béton peut atteindre des résistances à la compression allant de 130 à 230 MPa, des résistances à la flexion supérieures à 20 MPa, des résistances à la traction supérieure à 10 MPa et un module d'élasticité supérieur à 40 GPa. Les performances mécaniques de BEVUP sont améliorées grâce à la réactivité du verre amorphe, à l'optimisation granulométrique et la densification des mélanges. Les produits de déchets de verre dans les BEVUP ont un comportement pouzzolanique et réagissent avec la portlandite générée par l'hydratation du ciment. Cependant, ceci n’est pas le cas avec le sable de quartz ni la poudre de quartz dans le BFUP classique, qui réagissent à la température élevée de 400 °C. L'addition des déchets de verre améliore la densification de l'interface entre les particules. Les particules de déchets de verre ont une grande rigidité, ce qui augmente le module d'élasticité du béton. Le BEVUP a également une très bonne durabilité. Sa porosité capillaire est très faible, et le matériau est extrêmement résistant à la pénétration d’ions chlorure (≈ 8 coulombs). Sa résistance à l'abrasion (indice de pertes volumiques) est inférieure à 1,3. Le BEVUP ne subit pratiquement aucune détérioration aux cycles de gel-dégel, même après 1000 cycles. Après une évaluation des BEVUP en laboratoire, une mise à l'échelle a été réalisée avec un malaxeur de béton industriel et une validation en chantier avec de la construction de deux passerelles. Les propriétés mécaniques supérieures des BEVUP a permis de concevoir les passerelles avec des sections réduites d’environ de 60% par rapport aux sections faites de BC. Le BEVUP offre plusieurs avantages économiques et environnementaux. Il réduit le coût de production et l’empreinte carbone des structures construites de béton fibré à ultra-hautes performances (BFUP) classique, en utilisant des matériaux disponibles localement. Il réduit les émissions de CO[indice inférieur 2] associées à la production de clinkers de ciment (50% de remplacement du ciment) et utilise efficacement les ressources naturelles. De plus, la production de BEVUP permet de réduire les quantités de déchets de verre stockés ou mis en décharge qui causent des problèmes environnementaux et pourrait permettre de sauver des millions de dollars qui pourraient être dépensés dans le traitement de ces déchets. Enfin, il offre une solution alternative aux entreprises de construction dans la production de BFUP à moindre coût.

Finite difference calculations of permeability in large domains in a wide porosity range.

Determining effective hydraulic, thermal, mechanical and electrical properties of porous materials by means of classical physical experiments is often time-consuming and expensive. Thus, accurate numerical calculations of material properties are of increasing interest in geophysical, manufacturing, bio-mechanical and environmental applications, among other fields. Characteristic material properties (e.g. intrinsic permeability, thermal conductivity and elastic moduli) depend on morphological details on the porescale such as shape and size of pores and pore throats or cracks. To obtain reliable predictions of these properties it is necessary to perform numerical analyses of sufficiently large unit cells. Such representative volume elements require optimized numerical simulation techniques. Current state-of-the-art simulation tools to calculate effective permeabilities of porous materials are based on various methods, e.g. lattice Boltzmann, finite volumes or explicit jump Stokes methods. All approaches still have limitations in the maximum size of the simulation domain. In response to these deficits of the well-established methods we propose an efficient and reliable numerical method which allows to calculate intrinsic permeabilities directly from voxel-based data obtained from 3D imaging techniques like X-ray microtomography. We present a modelling framework based on a parallel finite differences solver, allowing the calculation of large domains with relative low computing requirements (i.e. desktop computers). The presented method is validated in a diverse selection of materials, obtaining accurate results for a large range of porosities, wider than the ranges previously reported. Ongoing work includes the estimation of other effective properties of porous media.