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.

Description of body posture during Static load bearing exercises for individuals with transfemoral amputation fitted with bone-anchored prosthesis

The rehabilitation programs of bone-anchorage prostheses relying either on the OPRA (Integrum, Sweden) or the ILP (Orthodynamics, Germany) fixation involve some forms of static load bearing exercises (LBE). So far, most of biomechanical studies of these static LBEs focused on the direct measurements of the actual forces and moments applied on the OPRA fixation of individuals with transfemoral amputation (TFA). To date, the proof-of-concept of an apparatus to conduct these kinetic measurements has been presented, along with some preliminary data. The understanding of the kinetic data is essential to improve rehabilitation programs as well as the design of upcoming loading frames. However, kinetic information alone is difficult to interpret without concomitant kinematic data. The purpose of this preliminary study was to introduce a qualitative analysis describing the different body postures during LBE for a group of TFAs.

Numerical treatment of a two-dimensional variable-order fractional nonlinear reaction-diffusion model

A two-dimensional variable-order fractional nonlinear reaction-diffusion model is considered. A second-order spatial accurate semi-implicit alternating direction method for a two-dimensional variable-order fractional nonlinear reaction-diffusion model is proposed. Stability and convergence of the semi-implicit alternating direct method are established. Finally, some numerical examples are given to support our theoretical analysis. These numerical techniques can be used to simulate a two-dimensional variable order fractional FitzHugh-Nagumo model in a rectangular domain. This type of model can be used to describe how electrical currents flow through the heart, controlling its contractions, and are used to ascertain the effects of certain drugs designed to treat arrhythmia.

A kinematic analysis of the hand function

Background The evaluation of the hand function is an essential element within the clinical practice. The usual assessments are focus on the ability to perform activities of daily life. The inclusion of instruments to measure kinematic variables provides a new approach to the assessment. Inertial sensors adapted to the hand could be used as a complementary instrument to the traditional assessment. Material: clinimetric assessment (Upper Limb Functional Index, Quick Dash), antrophometric variables (eight and weight), dynamometry (palm preasure) was taken. Functional analysis was made with Acceleglove system for the right hand and computer system. The glove has six acceleration sensor, one on each finger and another one on the reverse palm. Method Analytic, transversal approach. Ten healthy subject made six task on evaluation table (tripod pinch, lateral pinch and tip pinch, extension grip, spherical grip and power grip). Each task was made and measure three times, the second one was analyze for the results section. A Matlab script was created for the analysis of each movement and detection phase based on module vector. Results The module acceleration vector offers useful information of the hand function. The data analysis obtained during the performance of functional gestures allows to identify five different phases within the movement, three static phase and tow dynamic, each module vector was allied to one task. Conclusion Module vector variables could be used for the analysis of the different task made by the hand. Inertial sensor could be use as a complement for the traditional assessment system.

Stochastic analysis of a radial-inflow turbine in the presence of parametric uncertainties

This paper presents an uncertainty quantification study of the performance analysis of the high pressure ratio single stage radial-inflow turbine used in the Sundstrand Power Systems T-100 Multi-purpose Small Power Unit. A deterministic 3D volume-averaged Computational Fluid Dynamics (CFD) solver is coupled with a non-statistical generalized Polynomial Chaos (gPC) representation based on a pseudo-spectral projection method. One of the advantages of this approach is that it does not require any modification of the CFD code for the propagation of random disturbances in the aerodynamic and geometric fields. The stochastic results highlight the importance of the blade thickness and trailing edge tip radius on the total-to-static efficiency of the turbine compared to the angular velocity and trailing edge tip length. From a theoretical point of view, the use of the gPC representation on an arbitrary grid also allows the investigation of the sensitivity of the blade thickness profiles on the turbine efficiency. The gPC approach is also applied to coupled random parameters. The results show that the most influential coupled random variables are trailing edge tip radius coupled with the angular velocity.

Information-theoretic analysis of brain white matter fiber orientation distribution functions

We propose a new information-theoretic metric, the symmetric Kullback-Leibler divergence (sKL-divergence), to measure the difference between two water diffusivity profiles in high angular resolution diffusion imaging (HARDI). Water diffusivity profiles are modeled as probability density functions on the unit sphere, and the sKL-divergence is computed from a spherical harmonic series, which greatly reduces computational complexity. Adjustment of the orientation of diffusivity functions is essential when the image is being warped, so we propose a fast algorithm to determine the principal direction of diffusivity functions using principal component analysis (PCA). We compare sKL-divergence with other inner-product based cost functions using synthetic samples and real HARDI data, and show that the sKL-divergence is highly sensitive in detecting small differences between two diffusivity profiles and therefore shows promise for applications in the nonlinear registration and multisubject statistical analysis of HARDI data.

The impact of nonlinear pedagogy on physical education teacher education students’ intrinsic motivation

Background: Providing motivationally supportive physical education experiences for learners is crucial since empirical evidence in sport and physical education research has associated intrinsic motivation with positive educational outcomes. Self-determination theory (SDT) provides a valuable framework for examining motivationally supportive physical education experiences through satisfaction of three basic psychological needs: autonomy, competence and relatedness. However, the capacity of the prescriptive teaching philosophy of the dominant traditional physical education teaching approach to effectively satisfy the psychological needs of students to engage in physical education has been questioned. The constraints-led approach (CLA) has been proposed as a viable alternative teaching approach that can effectively support students’ self-motivated engagement in physical education. Purpose: We sought to investigate whether adopting the learning design and delivery of the CLA, guided by key pedagogical principles of nonlinear pedagogy (NLP), would address basic psychological needs of learners, resulting in higher self-reported levels of intrinsic motivation. The claim was investigated using action research. The teacher/researcher delivered two lessons aimed at developing hurdling skills: one taught using the CLA and the other using the traditional approach. Participants and Setting: The main participant for this study was the primary researcher and lead author who is a PETE educator, with extensive physical education teaching experience. A sample of 54 pre-service PETE students undertaking a compulsory second year practical unit at an Australian university was recruited for the study, consisting of an equal number of volunteers from each of two practical classes. A repeated measures experimental design was adopted, with both practical class groups experiencing both teaching approaches in a counterbalanced order. Data collection and analysis: Immediately after participation in each lesson, participants completed a questionnaire consisting of 22 items chosen from validated motivation measures of basic psychological needs and indices of intrinsic motivation, enjoyment and effort. All questionnaire responses were indicated on a 7-point Likert scale. A two-tailed, paired-samples t-test was used to compare the groups’ motivation subscale mean scores for each teaching approach. The size of the effect for each group was calculated using Cohen’s d. To determine whether any significant differences between the subscale mean scores of the two groups was due to an order effect, a two-tailed, independent samples t test was used. Findings: Participants’ reported substantially higher levels of self-determination and intrinsic motivation during the CLA hurdles lesson compared to during the traditional hurdles lesson. Both groups reported significantly higher motivation subscale mean scores for competence, relatedness, autonomy, enjoyment and effort after experiencing the CLA than mean scores reported after experiencing the traditional approach. This significant difference was evident regardless of the order that each teaching approach was experienced. Conclusion: The theoretically based pedagogical principles of NLP that inform learning design and delivery of the CLA may provide teachers and coaches with tools to develop more functional pedagogical climates, which result in students exhibiting more intrinsically motivated behaviours during learning.

Analysis of Business Process Variants in Apromore

In this paper we illustrate a set of features of the Apromore process model repository for analyzing business process variants. Two types of analysis are provided: one is static and based on differences on the process control flow, the other is dynamic and based on differences in the process behavior between the variants. These features combine techniques for the management of large process model collections with those for mining process knowledge from process execution logs. The tool demonstration will be useful for researchers and practitioners working on large process model collections and process execution logs, and specifically for those with an interest in understanding, managing and consolidating business process variants both within and across organizational boundaries.

A preconditioned numerical solver for stiff nonlinear reaction-diffusion equations with fractional Laplacians that avoids dense matrices

The numerical solution of fractional partial differential equations poses significant computational challenges in regard to efficiency as a result of the spatial nonlocality of the fractional differential operators. The dense coefficient matrices that arise from spatial discretisation of these operators mean that even one-dimensional problems can be difficult to solve using standard methods on grids comprising thousands of nodes or more. In this work we address this issue of efficiency for one-dimensional, nonlinear space-fractional reaction–diffusion equations with fractional Laplacian operators. We apply variable-order, variable-stepsize backward differentiation formulas in a Jacobian-free Newton–Krylov framework to advance the solution in time. A key advantage of this approach is the elimination of any requirement to form the dense matrix representation of the fractional Laplacian operator. We show how a banded approximation to this matrix, which can be formed and factorised efficiently, can be used as part of an effective preconditioner that accelerates convergence of the Krylov subspace iterative solver. Our approach also captures the full contribution from the nonlinear reaction term in the preconditioner, which is crucial for problems that exhibit stiff reactions. Numerical examples are presented to illustrate the overall effectiveness of the solver.

A CCG virtual system for big data application communication costs analysis

Network topology and routing are two important factors in determining the communication costs of big data applications at large scale. As for a given Cluster, Cloud, or Grid system, the network topology is fixed and static or dynamic routing protocols are preinstalled to direct the network traffic. Users cannot change them once the system is deployed. Hence, it is hard for application developers to identify the optimal network topology and routing algorithm for their applications with distinct communication patterns. In this study, we design a CCG virtual system (CCGVS), which first uses container-based virtualization to allow users to create a farm of lightweight virtual machines on a single host. Then, it uses software-defined networking (SDN) technique to control the network traffic among these virtual machines. Users can change the network topology and control the network traffic programmingly, thereby enabling application developers to evaluate their applications on the same system with different network topologies and routing algorithms. The preliminary experimental results through both synthetic big data programs and NPB benchmarks have shown that CCGVS can represent application performance variations caused by network topology and routing algorithm.

Finite element analysis of deep transverse metatarsal ligaments mechanical response during landing

The deep transverse metatarsal ligaments (DTML) play an important role in stabilizing the metatarsal bones and manipulating foot transverse arch deformation. However, the biomechanical research about DTML in the foot maneuver is quite few. Due to the difficulties and lack of better measurement technology for these ligaments experimental monitor, the load transfer mechanism and internal stress state also hadn't been well addressed. The purpose of this study was to develop a detailing foot finite element model including DTML tissues, to investigate the mechanical response of DTML during the landing condition. The DTML was considered as hyperelastic material model was used to represent the nonlinear and nearly incompressible nature of the ligament tissue. From the simulation results, it is clearly to find that the peak maiximal principal stress of DTML was between the third and fourth metatarsals. Meanwhile, it seems the DTML in the middle position experienced higher tension than the sides DTML.

A revisit to Pope's cohort analysis

Gulland's [Gulland, J.A., 1965. Estimation of mortality rates. Annex to Arctic Fisheries Working Group Report (meeting in Hamburg, January 1965). ICES. C.M. 1965, Doc. No. 3 (mimeographed)] virtual population analysis (VPA) is commonly used for studying the dynamics of harvested fish populations. However, it necessitates the solving of a nonlinear equation for the instantaneous rate of fishing mortality of the fish in a population. Pope [Pope, J.G., 1972. An investigation of the accuracy of Virtual Population Analysis using cohort analysis. ICNAF Res. Bull. 9, 65-74. Also available in D.H. Cushing (ed.) (1983), Key Papers on Fish Populations, p. 291-301, IRL Press, Oxford, 405 p.] eliminated this necessity in his cohort analysis by approximating its underlying age- and time-dependent population model. His approximation has since become one of the most commonly used age- and time-dependent fish population models in fisheries science. However, some of its properties are not well understood. For example, many assert that it describes the dynamics of a fish population, from which the catch of fish is taken instantaneously in the middle of the year. Such an assertion has never been proven, nor has its implied instantaneous rate of fishing mortality of the fish of a particular age at a particular time been examined, nor has its implied catch equation been derived from a general catch equation. In this paper, we prove this assertion, examine its implied instantaneous rate of fishing mortality of the fish of a particular age at a particular time, derive its implied catch equation from a general catch equation, and comment on how to structure an age- and time-dependent population model to ensure its internal consistency. This work shows that Gulland's (1965) virtual population analysis and Pope's (1972) cohort analysis lie at the opposite end of a continuous spectrum as a general model for a seasonally occurring fishery; Pope's (1972) approximation implies an infinitely large instantaneous rate of fishing mortality of the fish of a particular age at a particular time in a fishing season of zero length; and its implied catch equation has an undefined instantaneous rate of fishing mortality of the fish in a population, but a well-defined cumulative instantaneous rate of fishing mortality of the fish in the population. This work also highlights a need for a more careful treatment of the times of start and end of a fishing season in fish population models.