Numerical methods for solving a two-dimensional variable-order anomalous subdiffusion equation

Anomalous subdiffusion equations have in recent years received much attention. In this paper, we consider a two-dimensional variable-order anomalous subdiffusion equation. Two numerical methods (the implicit and explicit methods) are developed to solve the equation. Their stability, convergence and solvability are investigated by the Fourier method. Moreover, the effectiveness of our theoretical analysis is demonstrated by some numerical examples. © 2011 American Mathematical Society.

Analysis of acoustic emission data for accurate damage assessment for structural health monitoring applications

Structural health monitoring (SHM) refers to the procedure used to assess the condition of structures so that their performance can be monitored and any damage can be detected early. Early detection of damage and appropriate retrofitting will aid in preventing failure of the structure and save money spent on maintenance or replacement and ensure the structure operates safely and efficiently during its whole intended life. Though visual inspection and other techniques such as vibration based ones are available for SHM of structures such as bridges, the use of acoustic emission (AE) technique is an attractive option and is increasing in use. AE waves are high frequency stress waves generated by rapid release of energy from localised sources within a material, such as crack initiation and growth. AE technique involves recording these waves by means of sensors attached on the surface and then analysing the signals to extract information about the nature of the source. High sensitivity to crack growth, ability to locate source, passive nature (no need to supply energy from outside, but energy from damage source itself is utilised) and possibility to perform real time monitoring (detecting crack as it occurs or grows) are some of the attractive features of AE technique. In spite of these advantages, challenges still exist in using AE technique for monitoring applications, especially in the area of analysis of recorded AE data, as large volumes of data are usually generated during monitoring. The need for effective data analysis can be linked with three main aims of monitoring: (a) accurately locating the source of damage; (b) identifying and discriminating signals from different sources of acoustic emission and (c) quantifying the level of damage of AE source for severity assessment. In AE technique, the location of the emission source is usually calculated using the times of arrival and velocities of the AE signals recorded by a number of sensors. But complications arise as AE waves can travel in a structure in a number of different modes that have different velocities and frequencies. Hence, to accurately locate a source it is necessary to identify the modes recorded by the sensors. This study has proposed and tested the use of time-frequency analysis tools such as short time Fourier transform to identify the modes and the use of the velocities of these modes to achieve very accurate results. Further, this study has explored the possibility of reducing the number of sensors needed for data capture by using the velocities of modes captured by a single sensor for source localization. A major problem in practical use of AE technique is the presence of sources of AE other than crack related, such as rubbing and impacts between different components of a structure. These spurious AE signals often mask the signals from the crack activity; hence discrimination of signals to identify the sources is very important. This work developed a model that uses different signal processing tools such as cross-correlation, magnitude squared coherence and energy distribution in different frequency bands as well as modal analysis (comparing amplitudes of identified modes) for accurately differentiating signals from different simulated AE sources. Quantification tools to assess the severity of the damage sources are highly desirable in practical applications. Though different damage quantification methods have been proposed in AE technique, not all have achieved universal approval or have been approved as suitable for all situations. The b-value analysis, which involves the study of distribution of amplitudes of AE signals, and its modified form (known as improved b-value analysis), was investigated for suitability for damage quantification purposes in ductile materials such as steel. This was found to give encouraging results for analysis of data from laboratory, thereby extending the possibility of its use for real life structures. By addressing these primary issues, it is believed that this thesis has helped improve the effectiveness of AE technique for structural health monitoring of civil infrastructures such as bridges.

Analysis of the damping properties of high-transmission-ratio Cycloidal drives

Cycloidal drives are widely used in today’s industries for drives where large reduction ratios are required. Drive-train dynamics plays an important role in their design. This paper presents a new methodology for assessing damping characteristics of Cycloidal drives and compares the natural frequencies obtained from experiments and theoretical/numerical calculations using Fast-Fourier-Transforms.

Automatic region-of-interest detection and prioritisation for visually optimised coding of low bit rate videos

The increasing popularity of video consumption from mobile devices requires an effective video coding strategy. To overcome diverse communication networks, video services often need to maintain sustainable quality when the available bandwidth is limited. One of the strategy for a visually-optimised video adaptation is by implementing a region-of-interest (ROI) based scalability, whereby important regions can be encoded at a higher quality while maintaining sufficient quality for the rest of the frame. The result is an improved perceived quality at the same bit rate as normal encoding, which is particularly obvious at the range of lower bit rate. However, because of the difficulties of predicting region-of-interest (ROI) accurately, there is a limited research and development of ROI-based video coding for general videos. In this paper, the phase spectrum quaternion of Fourier Transform (PQFT) method is adopted to determine the ROI. To improve the results of ROI detection, the saliency map from the PQFT is augmented with maps created from high level knowledge of factors that are known to attract human attention. Hence, maps that locate faces and emphasise the centre of the screen are used in combination with the saliency map to determine the ROI. The contribution of this paper lies on the automatic ROI detection technique for coding a low bit rate videos which include the ROI prioritisation technique to give different level of encoding qualities for multiple ROIs, and the evaluation of the proposed automatic ROI detection that is shown to have a close performance to human ROI, based on the eye fixation data.

Infrared spectra of Mg-SiO smokes : comparison with analytical electron microscopy studies

An important component of current models for interstellar and circumstellar evolution is the infrared (IR)spectral data collected from stellar outflows around oxygen-rich stars and from the general interstellar medium [1]. IR spectra from these celestial bodies are usually interpreted as showing the general properties of sub-micron sized silicate grains [2]. Two major features at 10 and 20 microns are reasonably attributed to amorphous olivine or pyroxene (e.g. Mg2Si04 or MgSi03) on the basis of comparisons with natural standards and vapor condensed silicates [3-6]. In an attempt to define crystallisation rates for spectrally amorphous condensates, Nuth and Donn [5] annealed experimentally produced amorphous magnesium silicate smokes at 1000K. On analysing these smokes at various annealing times, Nuth and Donn [5] showed that changes in crystallinity measured by bulk X-ray diffraction occured at longer annealing times (days) than changes measured by IR spectra (a few hours). To better define the onset of crystallinity in these magnesium silicates, we have examined each annealed product using a JEOL 1OOCX analytical electron microscope (AEM). In addition, the development of chemical diversity with annealing has been monitored using energy dispersive spectroscopy of individual grains from areas <20nm in diameter. Furthermore, the crystallisation kinetics of these smokes under ambient, room temperature conditions have been examined using bulk and fourier transform infrared (FTIR)spectra.

Average structure of an An (sub 48) plagioclase from the Hogarth Ranges

The average structure (CI) of a volcanic plagioclase megacryst with composition Ano, from the Hogarth Ranges, Australia, has been determined using three-dimensional, singlecrystal neutron and X-ray diffraction data. Least squaresr efinements, incorporating anisotropic thermal motion of all atoms and an extinction correction, resulted in weighted R factors (based on intensities) of 0.076 and 0.056, respectively, for the neutron and X-ray data. Very weak e reflections could be detected in long-exposure X-ray and electron diffraction photographs of this crystal, but the refined average structure is believed to be unaffected by the presence of such a weak superstructure. The ratio of the scattering power of Na to that of Ca is different for X ray and neutron radiation, and this radiation-dependence of scattering power has been used to determine the distribution of Na and Ca over a split-atom M site (two sites designated M' and M") in this Ano, plagioclase. Relative peak-height ratios M'/M", revealed in difference Fourier sections calculated from neutron and X-ray data, formed the basis for the cation-distribution analysis. As neutron and X-ray data sets were directly compared in this analysis, it was important that systematic bias between refined neutron and X-ray positional parameters could be demonstrated to be absent. In summary, with an M-site model constrained only by the electron-microprobedetermined bulk composition of the crystal, the following values were obtained for the M-site occupanciesN: ar, : 0.29(7),N ar. : 0.23(7),C ar, : 0.15(4),a nd Car" : 0.33(4). These results indicate that restrictive assumptions about M sites, on which previous plagioclase refinements have been based, are not applicable to this Ano, and possibly not to the entire compositional range. T-site ordering determined by (T-O) bond-length variation-t,o : 0.51(l), trm = t2o = t2m = 0.32(l)-is weak, as might be expectedf rom the volcanic origin of this megacryst.

Identification of copper species present in Cu-ZSM-5 catalysts for NOx reduction

The structure of Cu-ZSM-5 catalysts that show activity for direct NO decomposition and selective catalytic reduction of NOx by hydrocarbons has been investigated by a multitude of modern surface analysis and spectroscopy techniques including X-ray photoelectron spectroscopy, thermogravimetric analysis, and in situ Fourier transform infrared spectroscopy. A series of four catalysts were prepared by exchange of Na-ZSM-5 with dilute copper acetate, and the copper loading was controlled by variation of the solution pH. Underexchanged catalysts contained isolated Cu2+OH-(H2O) species and as the copper loading was increased Cu2+ ions incorporated into the zeolite lattice appeared. The sites at which the latter two copper species were located were fundamentally different. The Cu2+OH-(H2O) moieties were bound to two lattice oxygen ions and associated with one aluminum framework species. In contrast, the Cu2+ ions were probably bound to four lattice oxygen ions and associated with two framework aluminum ions. Once the Cu-ZSM-5 samples attained high levels of exchange, the development of [Cu(μ-OH)2Cu]n2+OH-(H2O) species along with a small concentration of Cu(OH)2 was observed. On activation in helium to 500°C the Cu2+OH-(H2O) species transformed into Cu2+O- and Cu+ moieties, whereas the Cu2+ ions were apparently unaffected by this treatment (apart from the loss of ligated water molecules). Calcination of the precursors resulted in the formation of Cu2+O2- and a one-dimensional CuO species. Temperature-programmed desorption studies revealed that oxygen was removed from the latter two species at 407 and 575°C, respectively. © 1999 Academic Press.

FT Raman spectroscopic characterization of oxalate precursors to YBCO superconductors

FT Raman spectroscopy has been used to characterise the composition of the oxalate precursor to YBCO superconductors. By comparison to spectra of barium, copper and yttrium oxalate it is concluded that the co-precipitate incorporates not only the individual oxalate species but also a species ascribed to a mixed oxalate system. Significantly, Raman spectroscopy demonstrated that the precursor was not amorphous as previously deduced from XRD studies. In contrast, it is hypothesised that the sample consists of very small crystalline particles.

Swept source OCT with air puff chamber for corneal dynamics measurements

None of currently used tonometers produce estimated IOP values that are free of errors. Measurement incredibility arises from indirect measurement of corneal deformation and the fact that pressure calculations are based on population averaged parameters of anterior segment. Reliable IOP values are crucial for understanding and monitoring of number of eye pathologies e.g. glaucoma. We have combined high speed swept source OCT with air-puff chamber. System provides direct measurement of deformation of cornea and anterior surface of the lens. This paper describes in details the performance of air-puff ssOCT instrument. We present different approaches of data presentation and analysis. Changes in deformation amplitude appears to be good indicator of IOP changes. However, it seems that in order to provide accurate intraocular pressure values an additional information on corneal biomechanics is necessary. We believe that such information could be extracted from data provided by air-puff ssOCT.

Time-varying bispectral analysis of visually evoked multi-channel EEG

Theoretical foundations of higher order spectral analysis are revisited to examine the use of time-varying bicoherence on non-stationary signals using a classical short-time Fourier approach. A methodology is developed to apply this to evoked EEG responses where a stimulus-locked time reference is available. Short-time windowed ensembles of the response at the same offset from the reference are considered as ergodic cyclostationary processes within a non-stationary random process. Bicoherence can be estimated reliably with known levels at which it is significantly different from zero and can be tracked as a function of offset from the stimulus. When this methodology is applied to multi-channel EEG, it is possible to obtain information about phase synchronization at different regions of the brain as the neural response develops. The methodology is applied to analyze evoked EEG response to flash visual stimulii to the left and right eye separately. The EEG electrode array is segmented based on bicoherence evolution with time using the mean absolute difference as a measure of dissimilarity. Segment maps confirm the importance of the occipital region in visual processing and demonstrate a link between the frontal and occipital regions during the response. Maps are constructed using bicoherence at bifrequencies that include the alpha band frequency of 8Hz as well as 4 and 20Hz. Differences are observed between responses from the left eye and the right eye, and also between subjects. The methodology shows potential as a neurological functional imaging technique that can be further developed for diagnosis and monitoring using scalp EEG which is less invasive and less expensive than magnetic resonance imaging.

Application of TG-FTIR to study SO2 evolved during the thermal decomposition of coal-derived pyrite

The thermal decomposition of the coal-derived pyrite was studied using thermogravimetry combining with Fourier-transform infrared spectroscopy (TG-FTIR) techniques to gain knowledge on the SO2 gas evolution process and formation mechanism during the thermal decomposition of the coal-derived pyrite. The results showed that the thermal decomposition of the coal-derived pyrite which started at about 400 ◦C was complete at 600 ◦C; the gas evolved can be established by combining the DTG peak, the Gram–Schmidt curve and in situ FTIR spectroscopic evolved gas analysis. It can be observed from the spectra that the pyrolysis products for the sample mainly vary in quantity, but not in species. It was proposed that the oxidation of the coal-derived pyrite started at about 400 ◦C and that pyrrhotite and hematite were formed as primary products. The SO2 released by the thermal decomposition of the coal-derived pyrite mainly occurred in the first pyrolysis stage between 410 and 470 ◦C with the maximum rate at 444 ◦C. Furthermore, the SO2 gas evolution and formation mechanism during the thermal decomposition of the coal-derived pyrite has been proposed.

Preparation and thermal energy storage properties of paraffin/calcined diatomite composites as form-stable phase change materials

A composite paraffin-based phase change material (PCM) was prepared by blending composite paraffin and calcined diatomite through the fusion adsorption method. In this study, raw diatomite was purified by thermal treatment in order to improve the adsorption capacity of diatomite, which acted as a carrier material to prepare shape-stabilized PCMs. Two forms of paraffin (paraffin waxes and liquid paraffin) with different melting points were blended together by the fusion method, and the optimum mixed proportion with a suitable phase-transition temperature was obtained through differential scanning calorimetry (DSC) analysis. Then the prepared composite paraffin was adsorbed in calcined diatomite. The prepared paraffin/calcined diatomite composites were characterized by the scanning electron microscope (SEM) and Fourier transformation infrared (FT-IR) analysis techniques. Thermal energy storage properties of the composite PCMs were determined by DSC method. DSC results showed that there was an optimum adsorption ratio between composite paraffin and calcined diatomite and the phase-transition temperature and the latent heat of the composite PCMs were 33.04 ◦C and 89.54 J/g, respectively. Thermal cycling test of composite PCMs showed that the prepared material is thermally reliable and chemically stable. The obtained paraffin/calcined diatomite composites have proper latent heat and melting temperatures, and show practical significance and good potential application value.

Numerical methods of the variable-order Rayleigh-Stokes problem for a heated generalized second grade fluid with fractional derivative

Rayleigh–Stokes problems have in recent years received much attention due to their importance in physics. In this article, we focus on the variable-order Rayleigh–Stokes problem for a heated generalized second grade fluid with fractional derivative. Implicit and explicit numerical methods are developed to solve the problem. The convergence, stability of the numerical methods and solvability of the implicit numerical method are discussed via Fourier analysis. Moreover, a numerical example is given and the results support the effectiveness of the theoretical analysis.

Numerical approximation for a variable-order nonlinear reaction–subdiffusion equation

Fractional reaction–subdiffusion equations are widely used in recent years to simulate physical phenomena. In this paper, we consider a variable-order nonlinear reaction–subdiffusion equation. A numerical approximation method is proposed to solve the equation. Its convergence and stability are analyzed by Fourier analysis. By means of the technique for improving temporal accuracy, we also propose an improved numerical approximation. Finally, the effectiveness of the theoretical results is demonstrated by numerical examples.