Voronoi, Delaunay, and Block-Structured Mesh Refinement for Solution of the Shallow-Water Equations on the Sphere

Alternative meshes of the sphere and adaptive mesh refinement could be immensely beneficial for weather and climate forecasts, but it is not clear how mesh refinement should be achieved. A finite-volume model that solves the shallow-water equations on any mesh of the surface of the sphere is presented. The accuracy and cost effectiveness of four quasi-uniform meshes of the sphere are compared: a cubed sphere, reduced latitude–longitude, hexagonal–icosahedral, and triangular–icosahedral. On some standard shallow-water tests, the hexagonal–icosahedral mesh performs best and the reduced latitude–longitude mesh performs well only when the flow is aligned with the mesh. The inclusion of a refined mesh over a disc-shaped region is achieved using either gradual Delaunay, gradual Voronoi, or abrupt 2:1 block-structured refinement. These refined regions can actually degrade global accuracy, presumably because of changes in wave dispersion where the mesh is highly nonuniform. However, using gradual refinement to resolve a mountain in an otherwise coarse mesh can improve accuracy for the same cost. The model prognostic variables are height and momentum collocated at cell centers, and (to remove grid-scale oscillations of the A grid) the mass flux between cells is advanced from the old momentum using the momentum equation. Quadratic and upwind biased cubic differencing methods are used as explicit corrections to a fast implicit solution that uses linear differencing.

A fast two-grid and finite section method for a class of integral equations on the real line with application to an acoustic scattering problem in the half-plane

We consider the numerical treatment of second kind integral equations on the real line of the form ∅(s) = ∫_(-∞)^(+∞)▒〖κ(s-t)z(t)ϕ(t)dt,s=R〗 (abbreviated ϕ= ψ+K_z ϕ) in which K ϵ L_1 (R), z ϵ L_∞ (R) and ψ ϵ BC(R), the space of bounded continuous functions on R, are assumed known and ϕ ϵ BC(R) is to be determined. We first derive sharp error estimates for the finite section approximation (reducing the range of integration to [-A, A]) via bounds on (1-K_z )^(-1)as an operator on spaces of weighted continuous functions. Numerical solution by a simple discrete collocation method on a uniform grid on R is then analysed: in the case when z is compactly supported this leads to a coefficient matrix which allows a rapid matrix-vector multiply via the FFT. To utilise this possibility we propose a modified two-grid iteration, a feature of which is that the coarse grid matrix is approximated by a banded matrix, and analyse convergence and computational cost. In cases where z is not compactly supported a combined finite section and two-grid algorithm can be applied and we extend the analysis to this case. As an application we consider acoustic scattering in the half-plane with a Robin or impedance boundary condition which we formulate as a boundary integral equation of the class studied. Our final result is that if z (related to the boundary impedance in the application) takes values in an appropriate compact subset Q of the complex plane, then the difference between ϕ(s)and its finite section approximation computed numerically using the iterative scheme proposed is ≤C_1 [kh log⁡〖(1⁄kh)+(1-Θ)^((-1)⁄2) (kA)^((-1)⁄2) 〗 ] in the interval [-ΘA,ΘA](Θ<1) for kh sufficiently small, where k is the wavenumber and h the grid spacing. Moreover this numerical approximation can be computed in ≤C_2 N log⁡N operations, where N = 2A/h is the number of degrees of freedom. The values of the constants C1 and C2 depend only on the set Q and not on the wavenumber k or the support of z.

Does reactive surface area depend on grain size? Results from pH 3, 25 circle C far-from-equilibrium flow-through dissolution experiments on anorthite and biotite

Laboratory determined mineral weathering rates need to be normalised to allow their extrapolation to natural systems. The principle normalisation terms used in the literature are mass, and geometric- and BET specific surface area (SSA). The purpose of this study was to determine how dissolution rates normalised to these terms vary with grain size. Different size fractions of anorthite and biotite ranging from 180-150 to 20-10 mu m were dissolved in pH 3, HCl at 25 degrees C in flow through reactors under far from equilibrium conditions. Steady state dissolution rates after 5376 h (anorthite) and 4992 h (biotite) were calculated from Si concentrations and were normalised to initial- and final- mass and geometric-, geometric edge- (biotite), and BET SSA. For anorthite, rates normalised to initial- and final-BET SSA ranged from 0.33 to 2.77 X 10(-10) mol(feldspar) m(-2) s(-1), rates normalised to initial- and final-geometric SSA ranged from 5.74 to 8.88 X 10(-10) mol(feldspar) m(-2) s(-1) and rates normalised to initial- and final-mass ranged from 0.11 to 1.65 mol(feldspar) g(-1) s(-1). For biotite, rates normalised to initial- and final-BET SSA ranged from 1.02 to 2.03 X 10(-12) mol(biotite) m(-2) s(-1), rates normalised to initial- and final-geometric SSA ranged from 3.26 to 16.21 X 10(-12) mol(biotite) m(-2) s(-1), rates normalised to initial- and final-geometric edge SSA ranged from 59.46 to 111.32 x 10(-12) mol(biotite) m(-2) s(-1) and rates normalised to initial- and final-mass ranged from 0.81 to 6.93 X 10(-12) mol(biotite) g(-1) s(-1). For all normalising terms rates varied significantly (p <= 0.05) with grain size. The normalising terms which gave least variation in dissolution rate between grain sizes for anorthite were initial BET SSA and initial- and final-geometric SSA. This is consistent with: (1) dissolution being dominated by the slower dissolving but area dominant non-etched surfaces of the grains and, (2) the walls of etch pits and other dissolution features being relatively unreactive. These steady state normalised dissolution rates are likely to be constant with time. Normalisation to final BET SSA did not give constant ratios across grain size due to a non-uniform distribution of dissolution features. After dissolution coarser grains had a greater density of dissolution features with BET-measurable but unreactive wall surface area than the finer grains. The normalising term which gave the least variation in dissolution rates between grain sizes for biotite was initial BET SSA. Initial- and final-geometric edge SSA and final BET SSA gave the next least varied rates. The basal surfaces dissolved sufficiently rapidly to influence bulk dissolution rate and prevent geometric edge SSA normalised dissolution rates showing the least variation. Simple modelling indicated that biotite grain edges dissolved 71-132 times faster than basal surfaces. In this experiment, initial BET SSA best integrated the different areas and reactivities of the edge and basal surfaces of biotite. Steady state dissolution rates are likely to vary with time as dissolution alters the ratio of edge to basal surface area. Therefore they would be more properly termed pseudo-steady state rates, only appearing constant because the time period over which they were measured (1512 h) was less than the time period over wich they would change significantly. (c) 2006 Elsevier Inc. All rights reserved.

On the numerical integration of a class of singular perturbation problems

A three-point difference scheme recently proposed in Ref. 1 for the numerical solution of a class of linear, singularly perturbed, two-point boundary-value problems is investigated. The scheme is derived from a first-order approximation to the original problem with a small deviating argument. It is shown here that, in the limit, as the deviating argument tends to zero, the difference scheme converges to a one-sided approximation to the original singularly perturbed equation in conservation form. The limiting scheme is shown to be stable on any uniform grid. Therefore, no advantage arises from using the deviating argument, and the most accurate and efficient results are obtained with the deviation at its zero limit.

Comparison of pollinators and natural enemies: a meta-analysis of landscape and local effects on abundance and richness in crops

To manage agroecosystems for multiple ecosystem services, we need to know whether the management of one service has positive, negative, or no effects on other services. We do not yet have data on the interactions between pollination and pest-control services. However, we do have data on the distributions of pollinators and natural enemies in agroecosystems. Therefore, we compared these two groups of ecosystem service providers, to see if the management of farms and agricultural landscapes might have similar effects on the abundance and richness of both. In a meta-analysis, we compared 46 studies that sampled bees, predatory beetles, parasitic wasps, and spiders in fields, orchards, or vineyards of food crops. These studies used the proximity or proportion of non-crop or natural habitats in the landscapes surrounding these crops (a measure of landscape complexity), or the proximity or diversity of non-crop plants in the margins of these crops (a measure of local complexity), to explain the abundance or richness of these beneficial arthropods. Compositional complexity at both landscape and local scales had positive effects on both pollinators and natural enemies, but different effects on different taxa. Effects on bees and spiders were significantly positive, but effects on parasitoids and predatory beetles (mostly Carabidae and Staphylinidae) were inconclusive. Landscape complexity had significantly stronger effects on bees than it did on predatory beetles and significantly stronger effects in non-woody rather than in woody crops. Effects on richness were significantly stronger than effects on abundance, but possibly only for spiders. This abundance-richness difference might be caused by differences between generalists and specialists, or between arthropods that depend on non-crop habitats (ecotone species and dispersers) and those that do not (cultural species). We call this the ‘specialist-generalist’ or ‘cultural difference’ mechanism. If complexity has stronger effects on richness than abundance, it might have stronger effects on the stability than the magnitude of these arthropod-mediated ecosystem services. We conclude that some pollinators and natural enemies seem to have compatible responses to complexity, and it might be possible to manage agroecosystems for the benefit of both. However, too few studies have compared the two, and so we cannot yet conclude that there are no negative interactions between pollinators and natural enemies, and no trade-offs between pollination and pest-control services. Therefore, we suggest a framework for future research to bridge these gaps in our knowledge.

A benchmark-driven modelling approach for evaluating deployment choices on a multi-core architecture

The complexity of current and emerging architectures provides users with options about how best to use the available resources, but makes predicting performance challenging. In this work a benchmark-driven model is developed for a simple shallow water code on a Cray XE6 system, to explore how deployment choices such as domain decomposition and core affinity affect performance. The resource sharing present in modern multi-core architectures adds various levels of heterogeneity to the system. Shared resources often includes cache, memory, network controllers and in some cases floating point units (as in the AMD Bulldozer), which mean that the access time depends on the mapping of application tasks, and the core's location within the system. Heterogeneity further increases with the use of hardware-accelerators such as GPUs and the Intel Xeon Phi, where many specialist cores are attached to general-purpose cores. This trend for shared resources and non-uniform cores is expected to continue into the exascale era. The complexity of these systems means that various runtime scenarios are possible, and it has been found that under-populating nodes, altering the domain decomposition and non-standard task to core mappings can dramatically alter performance. To find this out, however, is often a process of trial and error. To better inform this process, a performance model was developed for a simple regular grid-based kernel code, shallow. The code comprises two distinct types of work, loop-based array updates and nearest-neighbour halo-exchanges. Separate performance models were developed for each part, both based on a similar methodology. Application specific benchmarks were run to measure performance for different problem sizes under different execution scenarios. These results were then fed into a performance model that derives resource usage for a given deployment scenario, with interpolation between results as necessary.

Wind-tunnel modelling of dispersion from a scalar area source in urban-like roughness

A wind-tunnel study was conducted to investigate ventilation of scalars from urban-like geometries at neighbourhood scale by exploring two different geometries a uniform height roughness and a non-uniform height roughness, both with an equal plan and frontal density of λ p = λ f = 25%. In both configurations a sub-unit of the idealized urban surface was coated with a thin layer of naphthalene to represent area sources. The naphthalene sublimation method was used to measure directly total area-averaged transport of scalars out of the complex geometries. At the same time, naphthalene vapour concentrations controlled by the turbulent fluxes were detected using a fast Flame Ionisation Detection (FID) technique. This paper describes the novel use of a naphthalene coated surface as an area source in dispersion studies. Particular emphasis was also given to testing whether the concentration measurements were independent of Reynolds number. For low wind speeds, transfer from the naphthalene surface is determined by a combination of forced and natural convection. Compared with a propane point source release, a 25% higher free stream velocity was needed for the naphthalene area source to yield Reynolds-number-independent concentration fields. Ventilation transfer coefficients w T /U derived from the naphthalene sublimation method showed that, whilst there was enhanced vertical momentum exchange due to obstacle height variability, advection was reduced and dispersion from the source area was not enhanced. Thus, the height variability of a canopy is an important parameter when generalising urban dispersion. Fine resolution concentration measurements in the canopy showed the effect of height variability on dispersion at street scale. Rapid vertical transport in the wake of individual high-rise obstacles was found to generate elevated point-like sources. A Gaussian plume model was used to analyse differences in the downstream plumes. Intensified lateral and vertical plume spread and plume dilution with height was found for the non-uniform height roughness

A high-order arbitrarily unstructured finite-volume model of the global atmosphere: tests solving the shallow-water equations

Simulations of the global atmosphere for weather and climate forecasting require fast and accurate solutions and so operational models use high-order finite differences on regular structured grids. This precludes the use of local refinement; techniques allowing local refinement are either expensive (eg. high-order finite element techniques) or have reduced accuracy at changes in resolution (eg. unstructured finite-volume with linear differencing). We present solutions of the shallow-water equations for westerly flow over a mid-latitude mountain from a finite-volume model written using OpenFOAM. A second/third-order accurate differencing scheme is applied on arbitrarily unstructured meshes made up of various shapes and refinement patterns. The results are as accurate as equivalent resolution spectral methods. Using lower order differencing reduces accuracy at a refinement pattern which allows errors from refinement of the mountain to accumulate and reduces the global accuracy over a 15 day simulation. We have therefore introduced a scheme which fits a 2D cubic polynomial approximately on a stencil around each cell. Using this scheme means that refinement of the mountain improves the accuracy after a 15 day simulation. This is a more severe test of local mesh refinement for global simulations than has been presented but a realistic test if these techniques are to be used operationally. These efficient, high-order schemes may make it possible for local mesh refinement to be used by weather and climate forecast models.

Effects of Delta9-tetrahydrocannabivarin on [35S]GTPgammaS binding in mouse brain cerebellum and piriform cortex membranes

BACKGROUND AND PURPOSE: We have recently shown that the phytocannabinoid Delta9-tetrahydrocannabivarin (Delta9-THCV) and the CB1 receptor antagonist AM251 increase inhibitory neurotransmission in mouse cerebellum and also exhibit anticonvulsant activity in a rat piriform cortical (PC) model of epilepsy. Possible mechanisms underlying cannabinoid actions in the CNS include CB1 receptor antagonism (by displacing endocannabinergic tone) or inverse agonism at constitutively active CB1 receptors. Here, we investigate the mode of cannabinoid action in [35S]GTPgammaS binding assays. EXPERIMENTAL APPROACH: Effects of Delta9-THCV and AM251 were tested either alone or against WIN55,212-2-induced increases in [35S]GTPgammaS binding in mouse cerebellar and PC membranes. Effects on non-CB receptor expressing CHO-D2 cell membranes were also investigated. KEY RESULTS :Delta9-THCV and AM251 both acted as potent antagonists of WIN55,212-2-induced increases in [35S]GTPgammaS binding in cerebellar and PC membranes (Delta9-THCV: pA2=7.62 and 7.44 respectively; AM251: pA2=9.93 and 9.88 respectively). At micromolar concentrations, Delta9-THCV or AM251 alone caused significant decreases in [35S]GTPgammaS binding; Delta9-THCV caused larger decreases than AM251. When applied alone in CHO-D2 membranes, Delta9-THCV and AM251 also caused concentration-related decreases in G protein activity. CONCLUSIONS AND IMPLICATIONS: Delta9-THCV and AM251 act as CB1 receptors antagonists in the cerebellum and PC, with AM251 being more potent than Delta9-THCV in both brain regions. Individually, Delta9-THCV or AM251 exhibited similar potency at CB1 receptors in the cerebellum and the PC. At micromolar concentrations, Delta9-THCV and AM251 caused a non-CB receptor-mediated depression of basal [35S]GTPgammaS binding.

Driving places - Marc Augé, non-places, and the geographies of England's M1 motorway

In this article I provide a critical account of the 'placing' of England's M1 motor-way. I start by critiquing Marc Auge's anthropological writings on 'non-places' which have provided a common point of reference for academics discussing spaces of travel, consumption and exchange in the contemporary world. I argue that Auge's ethnology of supermodernity results in a rather partial account of these sites, that he overstates the novelty of contemporary experiences of these spaces, and that he fails to acknowledge the heterogeneity and materiality of the social networks bound up with the production of non-places/places. I suggest that, rather than focusing on the presences and absences associated with the polarities of place and non-place, academics should examine the multiple, partial, dynamic and relational 'placings' which arise through the diverse performances and movements associated with travel, consumption and exchange. I then trace the topologies of England's M1 motorway, examining some of the different ways in which the motorway has been assembled, performed and placed over the past 45 years.

Does gym use impact upon nutritional knowledge?

Purpose – The purpose of this paper is to investigate gym and non-gym users' use and understanding of nutrition labels. Design/methodology/approach – A consumer survey in the form of a questionnaire conducted in the Greater London area in February/March 2005. Subject recruitment process took place in both a gym and university setting. Frequency tables and ?2-test are used to assess relationships between variables (p=0.05). Findings – The resulting sample consisted of 187 subjects, with predominance of females and gym users. Of the subjects, 88 per cent reported to at least occasionally read nutrition labels, with higher reading rates amongst women, irrespective of gym user status. Total and saturated fats are the most often information viewed on labels, however the overall knowledge of the calorie content of fat is low, with 53 per cent of subjects responding saturated fat contains more calories per gram when compared with other types of fats. This paper does not find significant differences in the use and understanding of nutrition labels between gym and non-gym users, but highlights the publics' continued lack of understanding of nutrition labels. Originality/value – This paper is unique as it investigates whether there is any difference between gym/non-gym users' use and interpretation of use of nutrition labels. It finds gender impacted more on nutritional labels knowledge than gym user's status. This points to a gender issue and questions the quality of information available to the general public. This paper is valuable as it highlights and identifies an area that requires further research and assessment, and is therefore useful to key stakeholders responsible for public health nutrition.

The EU's Export Refunds on Processed Foods: Legitimate in the WTO?

Export subsidies on processed foods are an important trade policy instrument for the European Union. GATT Article XVI legitimised the use of export subsidies on primary agricultural products, under certain circumstances, but forbade the use of export subsidies on non-primary products. However it was never satisfactorily resolved whether export subsidies could be paid on the primary agricultural products incorporated into processed products, such as pasta. The Uruguay Round Agreements, and particularly the Agreement on Agriculture (the URAA), apparently legitimised the EU’s practice of paying export subsidies on incorporated agricultural products, at least while the Peace Clause was in force. With the demise of the Peace Clause the question arises whether GATT Article XVI has any residual force, given that the range of primary agricultural products exempted by Article XVI from the ban on export subsidies is narrower than the list of agricultural products covered by the URAA.

Rate-limiting steps in the development of atherosclerosis: the response-to-influx theory

A large number of processes are involved in the pathogenesis of atherosclerosis but it is unclear which of them play a rate-limiting role. One way of resolving this problem is to investigate the highly non-uniform distribution of disease within the arterial system; critical steps in lesion development should be revealed by identifying arterial properties that differ between susceptible and protected sites. Although the localisation of atherosclerotic lesions has been investigated intensively over much of the 20th century, this review argues that the factor determining the distribution of human disease has only recently been identified. Recognition that the distribution changes with age has, for the first time, allowed it to be explained by variation in transport properties of the arterial wall; hitherto, this view could only be applied to experimental atherosclerosis in animals. The newly discovered transport variations which appear to play a critical role in the development of adult disease have underlying mechanisms that differ from those elucidated for the transport variations relevant to experimental atherosclerosis: they depend on endogenous NO synthesis and on blood flow. Manipulation of transport properties might have therapeutic potential. Copyright (C) 2004 S. Karger AG, Basel.

Does gym use impact upon nutritional knowledge?

Purpose - The purpose of this paper is to investigate gym and non-gym users' use and understanding of nutrition labels. Design/methodology/approach - A consumer survey in the form of a questionnaire conducted in the Greater London area in February/March 2005. Subject recruitment process took place in both a gym and university setting. Frequency tables and chi(2)-test are used to assess relationships between variables (p = 0.05). Findings - The resulting sample consisted of 187 subjects, with predominance of females and gym users. Of the subjects, 88 per cent reported to at least occasionally read nutrition labels, with higher reading rates amongst women, irrespective of gym user status. Total and saturated fats are the most often information viewed on labels, however the overall knowledge of the calorie content of fat is low, with 53 per cent of subjects responding saturated fat contains more calories per gram when compared with other types of fats. This paper does not find significant differences in the use and understanding of nutrition labels between gym and non-gym users, but highlights the publics' continued lack of understanding of nutrition labels. Originality/value - This paper is unique as it investigates whether there is any difference between gym/non-gym users' use and interpretation of use of nutrition labels. It finds gender impacted more on nutritional labels knowledge than gym user's status. This points to a gender issue and questions the quality of information available to the general public. This paper is valuable as it highlights and identifies an area that requires further research and assessment, and is therefore useful to key stakeholders responsible for public health nutrition.

Finding the point on Bezier Curves with the normal vector passing an external point

Two algorithms for finding the point on non-rational/rational Bezier curves of which the normal vector passes through a given external point are presented. The algorithms are based on Bezier curves generation algorithms of de Casteljau's algorithm for non-rational Bezier curve or Farin's recursion for rational Bezier curve, respectively. Orthogonal projections from the external point are used to guide the directional search used in the proposed iterative algorithms. Using Lyapunov's method, it is shown that each algorithm is able to converge to a local minimum for each case of non-rational/rational Bezier curves. It is also shown that on convergence the distance between the point on curves to the external point reaches a local minimum for both approaches. Illustrative examples are included to demonstrate the effectiveness of the proposed approaches.