A multilevel multiscale finite volume method

The Multiscale Finite Volume (MsFV) method has been developed to efficiently solve reservoir-scale problems while conserving fine-scale details. The method employs two grid levels: a fine grid and a coarse grid. The latter is used to calculate a coarse solution to the original problem, which is interpolated to the fine mesh. The coarse system is constructed from the fine-scale problem using restriction and prolongation operators that are obtained by introducing appropriate localization assumptions. Through a successive reconstruction step, the MsFV method is able to provide an approximate, but fully conservative fine-scale velocity field. For very large problems (e.g. one billion cell model), a two-level algorithm can remain computational expensive. Depending on the upscaling factor, the computational expense comes either from the costs associated with the solution of the coarse problem or from the construction of the local interpolators (basis functions). To ensure numerical efficiency in the former case, the MsFV concept can be reapplied to the coarse problem, leading to a new, coarser level of discretization. One challenge in the use of a multilevel MsFV technique is to find an efficient reconstruction step to obtain a conservative fine-scale velocity field. In this work, we introduce a three-level Multiscale Finite Volume method (MlMsFV) and give a detailed description of the reconstruction step. Complexity analyses of the original MsFV method and the new MlMsFV method are discussed, and their performances in terms of accuracy and efficiency are compared.

Geochemical study of vertebrate fossils from the Upper Cretaceous (Santonian) Csehbanya Formation (Hungary): Evidence for a freshwater habitat of mosasaurs and pycnodont fish

The diverse vertebrate remains from the Upper Cretaceous freshwater settings at Iharkut, Hungary, contain two fossil groups, Pycnodontiformes fish and Mosasauridae that are almost exclusively known from marine palaeo-environments. Hence, their appearance in alluvial sediments is very unusual. Trace element and isotope compositions of the remains have been analyzed to investigate the taphonomy and the ecological differences among the different fossil groups present at Iharkut. All examined fossils have undergone post-depositional diagenetic alteration, which resulted in high concentrations of REE, U, and Fe, together with almost complete homogenization of delta(18)O(CO3) values. Similar REE patterns in different fossils suggest a common origin for all remains, hence the discovered species most likely lived in the same local ecosystem. Despite partial diagenetic overprinting, the delta(18)O(PO4) values of the fossils indicate sufficient taxon-specific isotopic diversity to permit some broad conclusions on the palaeo-environment of the fossils. In particular, it is apparent that the isotopic composition of the Pycnodontiformes fish and Mosasauridae remains is most compatible with a freshwater palaeo-habitat and incompatible with a marine palaeo-environment. In addition, the Sr concentration and isotope data indicate that the Pycnodontiformes and Mosasauridae likely lived predominantly in a freshwater environment and were not simply occasional visitors to the Iharkut river ecosystem. Regarding other fossil groups, high delta(18)O(PO4) values of Alligatoroidea and Iharkutosuchus teeth suggest that these small crocodile species might have inhabited swamps and ponds where the water was relatively rich in (18)O due to evaporation. (c) 2009 Elsevier B.V. All rights reserved.

Model development for capercaillie (Tetrao urogallus) habitat in the Jura Mountains (Western Switzerland)

Capercaillie, Tetrao urogallus, is a threatened species in central Europe, and Swiss populations declined 40 to 50 % between 1970 and 1985. Capercaillie are sensitive to forest structure, and loss of habitat is a major cause of their decline. Knowledge of habitat characteristics is therefore essential for capercaillie conservation. Here, we present models predicting capercaillie probability of occurrence, based on relevant structural habitat variables. Models were built using multiple logistic regression analyses on capercaillie presence/absence data. Vegetation survey was carried out in July 1999 in a 170-km2 forested area (Jura mountains, canton de Vaud, western Switzerland) inhabited by capercaillie and presence/absence of the species was assessed according to dropping presence/absence. The survey was based on 10-m-radius sample plots each in a 1-km2 forest patch (n = 76 with capercaillie droppings, n = 80 without). A first model included seven out of 27 measured habitat variables and a second model only four. The latter model best represents practical needs. It includes three variables which had a negative impact on capercaillie presence: tree and shrub covers and spruce, Picea excelsa, shrub cover, and one which had a positive effect: bilberry, Vaccinium myrtillus, cover, highlighting that capaercaillie selected open forest with high bilberry abundance. The model can be used to map potential capercaillie habitat distribution and to manage the habitat in favour of capercaillie (protection and adapted forestry practices) in the Swiss Jura mountains.

Importance of polyethylene thickness in total shoulder arthroplasty: a finite element analysis.

BACKGROUND: Articular surfaces reconstruction is essential in total shoulder arthroplasty. Because of the limited glenoid bone support, thin glenoid component could improve anatomical reconstruction, but adverse mechanical effects might appear. METHODS: With a numerical musculoskeletal shoulder model, we analysed and compared three values of thickness of a typical all-polyethylene glenoid component: 2, 4 (reference) and 6mm. A loaded movement of abduction in the scapular plane was simulated. We evaluated the humeral head translation, the muscle moment arms, the joint force, the articular contact pattern, and the polyethylene and cement stress. Findings Decreasing polyethylene thickness from 6 to 2mm slightly increased humeral head translation and muscle moment arms. This induced a small decreased of the joint reaction force, but important increase of stress within the polyethylene and the cement mantel. Interpretation The reference thickness of 4mm seems a good compromise to avoid stress concentration and joint stuffing.

Nonlinear Finite Element Study of Piles in Integral Abutment Bridges, HR-227, Part 2, 1982

The highway departments of the states which use integral abutments in bridge design were contacted in order to study the extent of integral abutment use in skewed bridges and to survey the different guidelines used for analysis and design of integral abutments in skewed bridges. The variation in design assumptions and pile orientations among the various states in their approach to the use of integral abutments on skewed bridges is discussed. The problems associated with the treatment of the approach slab, backfill, and pile cap, and the reason for using different pile orientations are summarized in the report. An algorithm based on a state-of-the-art nonlinear finite element procedure previously developed by the authors was modified and used to study the influence of different factors on behavior of piles in integral abutment bridges. An idealized integral abutment was introduced by assuming that the pile is rigidly cast into the pile cap and that the approach slab offers no resistance to lateral thermal expansion. Passive soil and shear resistance of the cap are neglected in design. A 40-foot H pile (HP 10 X 42) in six typical Iowa soils was analyzed for fully restrained pile head and pinned pile head. According to numerical results, the maximum safe length for fully restrained pile head is one-half the maximum safe length for pinned pile head. If the pile head is partially restrained, the maximum safe length will lie between the two limits. The numerical results from an investigation of the effect of predrilled oversized holes indicate that if the length of the predrilled oversized hole is at least 4 feet below the ground, the vertical load-carrying capacity of the H pile is only reduced by 10 percent for 4 inches of lateral displacement in very stiff clay. With no predrilled oversized hole, the pile failed before the 4-inch lateral displacement was reached. Thus, the maximum safe lengths for integral abutment bridges may be increased by predrilling. Four different typical Iowa layered soils were selected and used in this investigation. In certain situations, compacted soil (> 50 blow count in standard penetration tests) is used as fill on top of natural soil. The numerical results showed that the critical conditions will depend on the length of the compacted soil. If the length of the compacted soil exceeds 4 feet, the failure mechanism for the pile is similar to one in a layer of very stiff clay. That is, the vertical load-carrying capacity of the H pile will be greatly reduced as the specified lateral displacement increases.

Woodland Wildlife Habitat Improvement Guide, February 4, 2010

Most non-grazed Iowa woodlands and surrounding openings are excellent habitat for a variety of birds and animals if there is a diversity of over story tree species and understory vegetation. As vegetative structures of woodlands change over time, the diversity of the woodland will change, and some species of birds and animals will benefit more then others. To optimize habitat development for the wildest range of bird and animal species, concentrate on maintaining as much vegetative diversity in the woodland as possible. To make improvement for individual species, the special needs of those species will have to be met by targeting precise woodland activities in specific areas.

Habitat selection in zones of parapatric contact between the common shrew Sorex araneaus and Millet's shrew S. coronatus

(1) The common shrew Sorex araneus and Millet's shrew S. coronatusare sibling species.They are morphologically and genetically very similar but do not hybridize.Their parapatric distribution throughout south-western Europe, with a few narrow zones of distributional overlap, suggests that they are in competitive parapatry. (2) Two of these contact zones were studied; there was evidence of coexistence over periods of 2 years as well as habitat segregation. In both zones, the species segregated on litter thickness and humidity variables. (3) A simple analysis of spatial distribution showed that habitats visible in the field corresponded to the habitats selected by the species. Habitat selection was found throughout the annual life-cycle of the shrews. (4) In one contact zone, a removal experiment was performed to test whether habitat segregation is induced by interspecific interactions. The experiment showed that the species select habitats differentially when both are present and abandon habitat selection when their competitor removed. (5)These results confirm the role of resource partitioning in promoting narrow ranges of distributional overlap between such parapatric species and qualitatively support the prediction of habitat selection theory that, in a two-species system, coexistence may be achieved by differential habitat selection to avoid competition. The results also support the view that the common shrew and Millet's shrew are in competitive parapatry.

An operator formulation of the multiscale finite-volume method with correction function

The multiscale finite-volume (MSFV) method has been derived to efficiently solve large problems with spatially varying coefficients. The fine-scale problem is subdivided into local problems that can be solved separately and are coupled by a global problem. This algorithm, in consequence, shares some characteristics with two-level domain decomposition (DD) methods. However, the MSFV algorithm is different in that it incorporates a flux reconstruction step, which delivers a fine-scale mass conservative flux field without the need for iterating. This is achieved by the use of two overlapping coarse grids. The recently introduced correction function allows for a consistent handling of source terms, which makes the MSFV method a flexible algorithm that is applicable to a wide spectrum of problems. It is demonstrated that the MSFV operator, used to compute an approximate pressure solution, can be equivalently constructed by writing the Schur complement with a tangential approximation of a single-cell overlapping grid and incorporation of appropriate coarse-scale mass-balance equations.

Interpolating and sampling sequences in finite Riemann surfaces

We provide a description of the interpolating and sampling sequences on a space of holomorphic functions on a finite Riemann surface, where a uniform growth restriction is imposed on the holomorphic functions.

Global habitat suitability models of terrestrial mammals.

Detailed large-scale information on mammal distribution has often been lacking, hindering conservation efforts. We used the information from the 2009 IUCN Red List of Threatened Species as a baseline for developing habitat suitability models for 5027 out of 5330 known terrestrial mammal species, based on their habitat relationships. We focused on the following environmental variables: land cover, elevation and hydrological features. Models were developed at 300 m resolution and limited to within species' known geographical ranges. A subset of the models was validated using points of known species occurrence. We conducted a global, fine-scale analysis of patterns of species richness. The richness of mammal species estimated by the overlap of their suitable habitat is on average one-third less than that estimated by the overlap of their geographical ranges. The highest absolute difference is found in tropical and subtropical regions in South America, Africa and Southeast Asia that are not covered by dense forest. The proportion of suitable habitat within mammal geographical ranges correlates with the IUCN Red List category to which they have been assigned, decreasing monotonically from Least Concern to Endangered. These results demonstrate the importance of fine-resolution distribution data for the development of global conservation strategies for mammals.

The number of planar central configurations for the 4-body problem is finite when 3 mass positions are fixed

In the n{body problem a central con guration is formed when the position vector of each particle with respect to the center of mass is a common scalar multiple of its acceleration vector. Lindstrom showed for n = 3 and for n > 4 that if n ? 1 masses are located at xed points in the plane, then there are only a nite number of ways to position the remaining nth mass in such a way that they de ne a central con guration. Lindstrom leaves open the case n = 4. In this paper we prove the case n = 4 using as variables the mutual distances between the particles.