A Generalised RBF Finite Difference Approach to Solve Nonlinear Heat Conduction Problems on Unstructured Datasets

Radial Basis Functions have traditionally been used to provide a continuous interpolation of scattered data sets. However, this interpolation also allows for the reconstruction of partial derivatives throughout the solution field, which can then be used to drive the solution of a partial differential equation. Since the interpolation takes place on a scattered dataset with no local connectivity, the solution is essentially meshless. RBF-based methods have been successfully used to solve a wide variety of PDEs in this fashion. Such full-domain RBF methods are highly flexible and can exhibit spectral convergence rates Madych & Nelson (1990). However, in their traditional implementation the fully-populated matrix systems which are produced lead to computational complexities of at least order-N2 with datasets of size N. In addition, they suffer fromincreasingly poor numerical conditioning as the size of the dataset grows, and also with increasingly flat interpolating functions. This is a consequence of ill-conditioning in the determination of RBF weighting coefficients (as demonstrated in Driscoll & Fornberg (2002)), and is described by Robert Schaback Schaback (1995) as the uncertainty relation; better conditioning is associated with worse accuracy, and worse conditioning is associated with improved accuracy. Many techniques have been developed to reduce the effect of the uncertainty relation in the traditional RBF formulation, such as RBF-specific preconditioners Baxter (2002); Beatson et al. (1999); Brown (2005); Ling & Kansa (2005), or adaptive selection of data centres Ling et al. (2006); Ling & Schaback (2004). However, at present the only reliable methods of controlling numerical ill-conditioning and computational cost as problem size increases are domain decomposition Hernandez Rosales & Power (2007); Wong et al. (1999); Zhang (2007); Zhou et al. (2003), or the use of locally supported basis functions Fasshauer (1999); Schaback (1997); Wendland (1995); Wu (1995). In this work the domain decomposition principle is applied, forming a large number of heavily overlapping systems that cover the solution domain. A small RBF collocation system is formed around each global data centre, with each collocation system used to approximate the governing PDE at its centrepoint, in terms of the solution value at surrounding collocation points. This leads to a sparse global linear system which may be solved using a variety of standard solvers. In this way, the proposed formulation emulates a finite difference method, with the RBF collocation systems replacing the polynomial interpolation functions used in traditional finite difference methods. However, unlike such polynomial functions RBF collocation is well suited to scattered data, and the method may be applied to both structured and unstructured datasets without modification. The method is applied here to solve the nonlinear heat conduction equation. In order to reduce the nonlinearity in the governing equation the Kirchhoff integral transformation is applied, and the transformed equation is solved using a Picard iterative process. The application of the Kirchhoff transform necessitates that the thermal property functions be transformed to Kirchhoff space also. If the thermal properties are a known and integrable function of temperature then the transformation may be performed analytically. Otherwise, an integration-interpolation procedure can be performed using 1D radial basis functions, as described in Stevens & Power (2010). In recent years a number of local RBF collocation techniques have been proposed, and applied a wide variety of problems (for example; Divo & Kassab (2007); Lee et al. (2003); Sarler & Vertnik (2006); Wright & Fornberg (2006)). A more comprehensive review of such methods is given in Stevens et al. (2009). Unlike most local RBF collocation methods that are used in the literature, the technique described here utilises the Hermitian RBF collocation formulation (see section 2 for more details), and allows both the PDE-boundary and PDE-governing operators to be included within in the local collocation systems. This inclusion of the governing PDE within the basis functions is shown in Stevens et al. (2009) to significantly improve the accuracy and stability of solutions obtained for linear transport problems. Additionally, the incorporation of information about the convective velocity field into the basis functionswas shown to have a stabilising effect, similar to traditional upwinding methods but without the requirement to alter the stencil configuration based on the local convective field. The standard approach to the solution of linear and nonlinear heat conduction problems is the use of finite difference and finite volume methods with simple polynomial interpolants Bejan (1993); Holman (2002); Kreith & Bohn (2000). Due to the dominance of diffusion in most cases, central differencing techniques are commonly used to compute the heat fluxes. However, limiter methods (such as the unconditionally stable TVD schemes) may be used for nonlinear heat conduction problems where the effective convection term, which results from the non-zero variation of thermal conductivity with temperature, can be expected to approach the magnitude of the diffusive term (see, for example, Shen & Han (2002)). Full-domain RBF methods have also been examined for use with nonlinear heat conduction problems (see Chantasiriwan (2007)), however such methods are restricted to small dataset sizes, due to the computational cost and numerical conditioning experienced by full-domain RBF techniques on large datasets. The present work demonstrates how local RBF collocation may be used as an alternative to traditional finite difference and finite volume methods, for nonlinear heat conduction problems. The described method retains freedom from a volumetric mesh, while allowing solution over unstructured datasets. A central stencil configuration is used in each case, and the solution is stabilised via the inclusion of the governing and boundary PDEs within the local collocation systems (\u201cimplicit upwinding\u201d), rather than by adjusting the stencil configuration based on the local solution field (\u201ctraditional upwinding\u201d). The method is validated using a transient numerical example with a known analytical solution (see section 4), and the ability of the formulation to handle strongly nonlinear problems is demonstrated in the solution of a food freezing problem (see section 5)

Alteration effects of volcanic ash in seawater: Anomalous Y/Ho ratios in coastal waters of the Central Mediterranean sea

This paper presents the results of a study based on data collected during the oceanographic cruise ANSIC 2001 carried out in the Ionian Sea during the explosive activity of Mount Etna in the summer of 2001. Anomalous low values of Y/Ho ratios in seawater suggest extensive scavenging processes on the surfaces of smectitic alteration products, with Y and Ho fractionation controlled by the differences in their electronic configurations and behaviour during solution/surface complexation equilibria. These processes can also be traced through the presence of significant tetrad effects recorded in the chondrite-normalised Rare Earth Elements and Yttrium (YREEs) patterns of suspended particulate matter. This suggests that the preferential Y scavenging from seawater is due to the formation of inner-sphere complexes with OH- groups on montmorillonite crystal surfaces. The preliminary results of kinetic experiments of YREE released from volcanic ash to coexisting seawater, and the related effects on Y/Ho ratios and Ce anomalies, are consistent with the fractionation of Light Rare Earth Elements (LREEs) with respect to Heavy Rare Earth Elements (HREEs) observed in dissolved phase. They suggest a behaviour of Y similar to that reported for LREEs, particularly for Ce and Pr. © 2007 Elsevier Ltd. All rights reserved

Earth conductivity structures and their effects on geomagnetic induction in pipelines

Anomalous, large pipe-to-soil potentials (PSP) have been observed along a natural gas pipeline in eastern Ontario, Canada, where there is a major geological contact between the highly resistive rocks of the Precambrian Shield to the west and the more conductive Paleozoic sediments to the east. This study tested the hypothesis that large variations of PSP are related to lateral changes of Earth conductivity under the pipeline. Concurrent and co-located PSP and magnetotelluric (MT) geophysical data were acquired in the study area. Results from the MT survey were used to model PSP variations based on distributed-source transmission line theory, using a spatially-variant surface geoelectric field. Different models were built to investigate the impact of different subsurface features. Good agreement between modelled and observed PSP was reached when impedance peaks related to major changes of subsurface geological conditions were included. The large PSP could therefore be attributed to the presence of resistive intrusive bodies in the upper crust and/or boundaries between tectonic terranes. This study demonstrated that combined PSP-MT investigations are a useful tool in the identification of potential hazards caused by geomagnetically induced currents in pipelines

Wellness and Multiple Sclerosis: The National MS Society Establishes a Wellness Research Working Group and Research Priorities

Background: People with multiple sclerosis (MS) have identified “wellness” and associated behaviors as a high priority based on “social media listening” undertaken by the National MS Society (i.e. the Society). Objective: The Society recently convened a group that consisted of researchers with experience in MS and wellness-related research, Society staff members, and an individual with MS for developing recommendations regarding a wellness research agenda. Method: The members of the group engaged in focal reviews and discussions involving the state of science within three approaches for promoting wellness in MS, namely diet, exercise, and emotional wellness. Results: That process informed a group-mediated activity for developing and prioritizing research goals for wellness in MS. This served as a background for articulating the mission and objectives of the Society’s Wellness Research Working Group. Conclusion: The primary mission of the Wellness Research Working Group is the provision of scientific evidence supporting the application of lifestyle, behavioral, and psychosocial approaches for promoting optimal health of mind, body, and spirit (i.e. wellness) in people with MS as well as managing the disease and its consequences