A fast semi-implicit difference method for a nonlinear two-sided space-fractional diffusion equation with variable diffusivity coefficients

In this paper, we derive a new nonlinear two-sided space-fractional diffusion equation with variable coefficients from the fractional Fick’s law. A semi-implicit difference method (SIDM) for this equation is proposed. The stability and convergence of the SIDM are discussed. For the implementation, we develop a fast accurate iterative method for the SIDM by decomposing the dense coefficient matrix into a combination of Toeplitz-like matrices. This fast iterative method significantly reduces the storage requirement of O(n2)O(n2) and computational cost of O(n3)O(n3) down to n and O(nlogn)O(nlogn), where n is the number of grid points. The method retains the same accuracy as the underlying SIDM solved with Gaussian elimination. Finally, some numerical results are shown to verify the accuracy and efficiency of the new method.

A Crank--Nicolson ADI spectral method for a two-dimensional riesz space fractional nonlinear reaction-diffusion equation

In this paper, a new alternating direction implicit Galerkin--Legendre spectral method for the two-dimensional Riesz space fractional nonlinear reaction-diffusion equation is developed. The temporal component is discretized by the Crank--Nicolson method. The detailed implementation of the method is presented. The stability and convergence analysis is strictly proven, which shows that the derived method is stable and convergent of order $2$ in time. An optimal error estimate in space is also obtained by introducing a new orthogonal projector. The present method is extended to solve the fractional FitzHugh--Nagumo model. Numerical results are provided to verify the theoretical analysis.

Stability and convergence of a new finite volume method for a two-sided space-fractional diffusion equation

In this paper, we consider a two-sided space-fractional diffusion equation with variable coefficients on a finite domain. Firstly, based on the nodal basis functions, we present a new fractional finite volume method for the two-sided space-fractional diffusion equation and derive the implicit scheme and solve it in matrix form. Secondly, we prove the stability and convergence of the implicit fractional finite volume method and conclude that the method is unconditionally stable and convergent. Finally, some numerical examples are given to show the effectiveness of the new numerical method, and the results are in excellent agreement with theoretical analysis.

A Novel High Order Space-Time Spectral Method for the Time Fractional Fokker--Planck Equation

The fractional Fokker-Planck equation is an important physical model for simulating anomalous diffusions with external forces. Because of the non-local property of the fractional derivative an interesting problem is to explore high accuracy numerical methods for fractional differential equations. In this paper, a space-time spectral method is presented for the numerical solution of the time fractional Fokker-Planck initial-boundary value problem. The proposed method employs the Jacobi polynomials for the temporal discretization and Fourier-like basis functions for the spatial discretization. Due to the diagonalizable trait of the Fourier-like basis functions, this leads to a reduced representation of the inner product in the Galerkin analysis. We prove that the time fractional Fokker-Planck equation attains the same approximation order as the time fractional diffusion equation developed in [23] by using the present method. That indicates an exponential decay may be achieved if the exact solution is sufficiently smooth. Finally, some numerical results are given to demonstrate the high order accuracy and efficiency of the new numerical scheme. The results show that the errors of the numerical solutions obtained by the space-time spectral method decay exponentially.

Numeracy counts in the statistical reasoning equation

The relationship between mathematics and statistical reasoning frequently receives comment (Vere-Jones 1995, Moore 1997); however most of the research into the area tends to focus on maths anxiety. Gnaldi (Gnaldi 2003) showed that in a statistics course for psychologists, the statistical understanding of students at the end of the course depended on students’ basic numeracy, rather than the number or level of previous mathematics courses the student had undertaken. As part of a study into the development of statistical thinking at the interface between secondary and tertiary education, students enrolled in an introductory data analysis subject were assessed regarding their statistical reasoning ability, basic numeracy skills and attitudes towards statistics. This work reports on the relationships between these factors and in particular the importance of numeracy to statistical reasoning.

Investigating the impact of project definition clarity on project performance: Structural equation modeling study

Having a clear project definition is crucial for successful construction projects. It affects design quality, project communication between stakeholders and final project performance in terms of cost, schedule and quality. This study examines the relationship between project definition and final project performance through a structural equation model comprising 4 latent constructs and 6 path hypotheses using data from a questionnaire survey of 120 general contractors in the Malaysian construction industry. The results show that in the study population, all three items impact the project performance, but the link between design quality and project performance is indirect. Instead, the clarity of project definition affects project performance indirectly through design quality and project communication and design quality affects project performance indirectly through project communication. The primary contribution is to provide quantitative confirmation of the more general statements made in the literature from around the world and therefore adds to and consolidates existing knowledge. Practical implications derived from the finding are also proposed for various project stakeholders. Furthermore, as lack of the clarity of project definition is a very common occurrence in construction projects globally, these findings have important ramifications for all construction projects in expanding and clarifying existing knowledge on what is needed for the successful delivery of construction projects.

Development of height-weight based equation for assessment of body composition in Sri Lankan children

Objective To develop a height and weight based equation to estimate total body water (TBW) in Sri Lankan children. Methods Cross sectional descriptive study done involving 5–15 year old healthy children. Height and weight were measured. TBW was assessed using isotope dilution method (D2O) and fat free mass (FFM) calculated. Multiple regression analysis was used to develop prediction equation and validated using PRESS statistical technique. Height, weight and sex code (male=1; female=0) were used as prediction variables. Results This study provides height and weight equation for the prediction of TBW in Sri Lankan children. To the best of our knowledge there are no published height weight prediction equations validated on South Asian populations. Conclusion Results of this study need to be affirmed by more studies on other closely related populations by using multicomponent body composition.

A stochastic exponential Euler scheme for simulation of stiff biochemical reaction systems

In order to simulate stiff biochemical reaction systems, an explicit exponential Euler scheme is derived for multidimensional, non-commutative stochastic differential equations with a semilinear drift term. The scheme is of strong order one half and A-stable in mean square. The combination with this and the projection method shows good performance in numerical experiments dealing with an alternative formulation of the chemical Langevin equation for a human ether a-go-go related gene ion channel mode

Numerical analysis of a new space-time variable fractional order advection-dispersion equation

Many physical processes appear to exhibit fractional order behavior that may vary with time and/or space. The continuum of order in the fractional calculus allows the order of the fractional operator to be considered as a variable. In this paper, we consider a new space–time variable fractional order advection–dispersion equation on a finite domain. The equation is obtained from the standard advection–dispersion equation by replacing the first-order time derivative by Coimbra’s variable fractional derivative of order α(x)∈(0,1]α(x)∈(0,1], and the first-order and second-order space derivatives by the Riemann–Liouville derivatives of order γ(x,t)∈(0,1]γ(x,t)∈(0,1] and β(x,t)∈(1,2]β(x,t)∈(1,2], respectively. We propose an implicit Euler approximation for the equation and investigate the stability and convergence of the approximation. Finally, numerical examples are provided to show that the implicit Euler approximation is computationally efficient.

Spatially explicit structural equation modeling

Structural equation modeling (SEM) is a powerful statistical approach for the testing of networks of direct and indirect theoretical causal relationships in complex data sets with intercorrelated dependent and independent variables. SEM is commonly applied in ecology, but the spatial information commonly found in ecological data remains difficult to model in a SEM framework. Here we propose a simple method for spatially explicit SEM (SE-SEM) based on the analysis of variance/covariance matrices calculated across a range of lag distances. This method provides readily interpretable plots of the change in path coefficients across scale and can be implemented using any standard SEM software package. We demonstrate the application of this method using three studies examining the relationships between environmental factors, plant community structure, nitrogen fixation, and plant competition. By design, these data sets had a spatial component, but were previously analyzed using standard SEM models. Using these data sets, we demonstrate the application of SE-SEM to regularly spaced, irregularly spaced, and ad hoc spatial sampling designs and discuss the increased inferential capability of this approach compared with standard SEM. We provide an R package, sesem, to easily implement spatial structural equation modeling.

Assessment of body composition in Sri Lankan children: Validation of a bioelectrical impedance prediction equation

Objective: To develop bioelectrical impedance analysis (BIA) equations to predict total body water (TBW) and fat-free mass (FFM) of Sri Lankan children. Subjects/Methods: Data were collected from 5- to 15-year-old healthy children. They were randomly assigned to validation (M/F: 105/83) and cross-validation (M/F: 53/41) groups. Height, weight and BIA were measured. TBW was assessed using isotope dilution method (D2 O). Multiple regression analysis was used to develop preliminary equations and cross-validated on an independent group. Final prediction equation was constructed combining the two groups and validated by PRESS (prediction of sum of squares) statistics. Impedance index (height2/impedance; cm2/Ω), weight and sex code (male = 1; female = 0) were used as variables. Results: Independent variables of the final prediction equation for TBW were able to predict 86.3% of variance with root means-squared error (RMSE) of 2.1l. PRESS statistics was 2.1l with press residuals of 1.2l. Independent variables were able to predict 86.9% of variance of FFM with RMSE of 2.7 kg. PRESS statistics was 2.8 kg with press residuals of 1.4 kg. Bland Altman technique showed that the majority of the residuals were within mean bias±1.96 s.d. Conclusions: Results of this study provide BIA equation for the prediction of TBW and FFM in Sri Lankan children. To the best of our knowledge there are no published BIA prediction equations validated on South Asian populations. Results of this study need to be affirmed by more studies on other closely related populations by using multi-component body composition assessment.

Relations between chronic regulatory focus and future time perspective: Results of a cross-lagged structural equation model

Future time perspective - the way individuals perceive their remaining time in life - importantly influences socio-emotional goals and motivational outcomes. Recently, researchers have called for studies that investigate relationships between personality and future time perspective. Using a cross-lagged panel design, this study investigated effects of chronic regulatory focus dimensions (promotion and prevention orientation) on future time perspective dimensions (focus on opportunities and limitations). Survey data were collected two times, separated by a 3. month time lag, from 85 participants. Results of structural equation modeling showed that promotion orientation had a positive lagged effect on focus on opportunities, and prevention orientation had a positive lagged effect on focus on limitations.

Characterization of anomalous relaxation using the time-fractional Bloch equation and multiple echo T2 *-weighted magnetic resonance imaging at 7T

PURPOSE To study the utility of fractional calculus in modeling gradient-recalled echo MRI signal decay in the normal human brain. METHODS We solved analytically the extended time-fractional Bloch equations resulting in five model parameters, namely, the amplitude, relaxation rate, order of the time-fractional derivative, frequency shift, and constant offset. Voxel-level temporal fitting of the MRI signal was performed using the classical monoexponential model, a previously developed anomalous relaxation model, and using our extended time-fractional relaxation model. Nine brain regions segmented from multiple echo gradient-recalled echo 7 Tesla MRI data acquired from five participants were then used to investigate the characteristics of the extended time-fractional model parameters. RESULTS We found that the extended time-fractional model is able to fit the experimental data with smaller mean squared error than the classical monoexponential relaxation model and the anomalous relaxation model, which do not account for frequency shift. CONCLUSIONS We were able to fit multiple echo time MRI data with high accuracy using the developed model. Parameters of the model likely capture information on microstructural and susceptibility-induced changes in the human brain.

Development of an efficiency equation for solar evaporator- collectors

Considering the growing energy needs and concern for environmental degradation, clean and inexhaustible energy sources, e.g solar energy are receiving greater attention for various applications. The use of solar energy systems for low temperature applications reduces the burden on conventional fossil fuels and has little or no harmful effects on the environment. The performance of a solar system depends to a great extent on the collector used for the conversion of solar radiant energy to thermal energy. A solar evaporatorcollector (SEC) is basically an unglazed flat plate collector where refrigerant, like R134a, is used as the working fluid. As the operating temperature of SEC is very low, it collects energy both from solar irradiation and ambient energy leading to a much higher efficiency than the conventional collectors. The capability of SEC to utilize ambient energy also enables the system to operate at night. Therefore it is not appropriate to use for the evaluation of performance of SEC by conventional efficiency equation where ambient energy and condensation is not considered as energy input in addition to irradiation. In the National University of Singapore, several Solar Assisted Heat Pump (SAHP) systems were built for the evaluation of performance under the metrological condition of Singapore for thermal applications of desalination and SEC was the main component to harness renewable energy. In this paper, the design and performance of SEC are explored. Furthermore, an attempt is made to develop an efficiency equation for SEC and maximum efficiency attained 98% under the meteorological condition of Singapore.

The Pandolf load carriage equation is a poor predictor of metabolic rate while wearing explosive ordnance disposal protective clothing

This investigation aimed to quantify metabolic rate when wearing an explosive ordnance disposal (EOD) ensemble (~33kg) during standing and locomotion; and determine whether the Pandolf load carriage equation accurately predicts metabolic rate when wearing an EOD ensemble during standing and locomotion. Ten males completed 8 trials with metabolic rate measured through indirect calorimetry. Walking in EOD at 2.5, 4.0 and 5.5km·h−1 was significantly (p < 0.05) greater than matched trials without the EOD ensemble by 49% (127W), 65% (213W) and 78% (345W), respectively. Mean bias (95% limits of agreement) between predicted and measured metabolism during standing, 2.5, 4 and 5.5km·h−1 were 47W (19 to 75W); −111W (−172 to −49W); −122W (−189 to −54W) and −158W (−245 to −72W), respectively. The Pandolf equation significantly underestimated measured metabolic rate during locomotion. These findings have practical implications for EOD technicians during training and operation and should be considered when developing maximum workload duration models and work-rest schedules.