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 method of lines approach for modelling saltwater intrusion in coastal aquifiers

The equations governing saltwater intrusion in coastal aquifers are complex. Backward Euler time stepping approaches are often used to advance the solution to these equations in time, which typically requires that small time steps be taken in order to ensure that an accurate solution is obtained. We show that a method of lines approach incorporating variable order backward differentiation formulas can greatly improve the efficiency of the time stepping process.

Efficient simulation of unsaturated flow using exponential time integration

We assess the performance of an exponential integrator for advancing stiff, semidiscrete formulations of the unsaturated Richards equation in time. The scheme is of second order and explicit in nature but requires the action of the matrix function φ(A) where φ(z) = [exp(z) - 1]/z on a suitability defined vector v at each time step. When the matrix A is large and sparse, φ(A)v can be approximated by Krylov subspace methods that require only matrix-vector products with A. We prove that despite the use of this approximation the scheme remains second order. Furthermore, we provide a practical variable-stepsize implementation of the integrator by deriving an estimate of the local error that requires only a single additional function evaluation. Numerical experiments performed on two-dimensional test problems demonstrate that this implementation outperforms second-order, variable-stepsize implementations of the backward differentiation formulae.

Performance assessment of exponential Rosenbrock method for large systems of ODE

This paper studies time integration methods for large stiff systems of ordinary differential equations (ODEs) of the form u'(t) = g(u(t)). For such problems, implicit methods generally outperform explicit methods, since the time step is usually less restricted by stability constraints. Recently, however, explicit so-called exponential integrators have become popular for stiff problems due to their favourable stability properties. These methods use matrix-vector products involving exponential-like functions of the Jacobian matrix, which can be approximated using Krylov subspace methods that require only matrix-vector products with the Jacobian. In this paper, we implement exponential integrators of second, third and fourth order and demonstrate that they are competitive with well-established approaches based on the backward differentiation formulas and a preconditioned Newton-Krylov solution strategy.

Efficient solution of two-sided nonlinear space-fractional diffusion equations using fast Poisson preconditioners

We develop a fast Poisson preconditioner for the efficient numerical solution of a class of two-sided nonlinear space fractional diffusion equations in one and two dimensions using the method of lines. Using the shifted Gr¨unwald finite difference formulas to approximate the two-sided(i.e. the left and right Riemann-Liouville) fractional derivatives, the resulting semi-discrete nonlinear systems have dense Jacobian matrices owing to the non-local property of fractional derivatives. We employ a modern initial value problem solver utilising backward differentiation formulas and Jacobian-free Newton-Krylov methods to solve these systems. For efficient performance of the Jacobianfree Newton-Krylov method it is essential to apply an effective preconditioner to accelerate the convergence of the linear iterative solver. The key contribution of our work is to generalise the fast Poisson preconditioner, widely used for integer-order diffusion equations, so that it applies to the two-sided space fractional diffusion equation. A number of numerical experiments are presented to demonstrate the effectiveness of the preconditioner and the overall solution strategy.

Exponential integrators and a dual-scale model for wood drying

For the timber industry, the ability to simulate the drying of wood is invaluable for manufacturing high quality wood products. Mathematically, however, modelling the drying of a wet porous material, such as wood, is a diffcult task due to its heterogeneous and anisotropic nature, and the complex geometry of the underlying pore structure. The well{ developed macroscopic modelling approach involves writing down classical conservation equations at a length scale where physical quantities (e.g., porosity) can be interpreted as averaged values over a small volume (typically containing hundreds or thousands of pores). This averaging procedure produces balance equations that resemble those of a continuum with the exception that effective coeffcients appear in their deffnitions. Exponential integrators are numerical schemes for initial value problems involving a system of ordinary differential equations. These methods differ from popular Newton{Krylov implicit methods (i.e., those based on the backward differentiation formulae (BDF)) in that they do not require the solution of a system of nonlinear equations at each time step but rather they require computation of matrix{vector products involving the exponential of the Jacobian matrix. Although originally appearing in the 1960s, exponential integrators have recently experienced a resurgence in interest due to a greater undertaking of research in Krylov subspace methods for matrix function approximation. One of the simplest examples of an exponential integrator is the exponential Euler method (EEM), which requires, at each time step, approximation of φ(A)b, where φ(z) = (ez - 1)/z, A E Rnxn and b E Rn. For drying in porous media, the most comprehensive macroscopic formulation is TransPore [Perre and Turner, Chem. Eng. J., 86: 117-131, 2002], which features three coupled, nonlinear partial differential equations. The focus of the first part of this thesis is the use of the exponential Euler method (EEM) for performing the time integration of the macroscopic set of equations featured in TransPore. In particular, a new variable{ stepsize algorithm for EEM is presented within a Krylov subspace framework, which allows control of the error during the integration process. The performance of the new algorithm highlights the great potential of exponential integrators not only for drying applications but across all disciplines of transport phenomena. For example, when applied to well{ known benchmark problems involving single{phase liquid ow in heterogeneous soils, the proposed algorithm requires half the number of function evaluations than that required for an equivalent (sophisticated) Newton{Krylov BDF implementation. Furthermore for all drying configurations tested, the new algorithm always produces, in less computational time, a solution of higher accuracy than the existing backward Euler module featured in TransPore. Some new results relating to Krylov subspace approximation of '(A)b are also developed in this thesis. Most notably, an alternative derivation of the approximation error estimate of Hochbruck, Lubich and Selhofer [SIAM J. Sci. Comput., 19(5): 1552{1574, 1998] is provided, which reveals why it performs well in the error control procedure. Two of the main drawbacks of the macroscopic approach outlined above include the effective coefficients must be supplied to the model, and it fails for some drying configurations, where typical dual{scale mechanisms occur. In the second part of this thesis, a new dual{scale approach for simulating wood drying is proposed that couples the porous medium (macroscale) with the underlying pore structure (microscale). The proposed model is applied to the convective drying of softwood at low temperatures and is valid in the so{called hygroscopic range, where hygroscopically held liquid water is present in the solid phase and water exits only as vapour in the pores. Coupling between scales is achieved by imposing the macroscopic gradient on the microscopic field using suitably defined periodic boundary conditions, which allows the macroscopic ux to be defined as an average of the microscopic ux over the unit cell. This formulation provides a first step for moving from the macroscopic formulation featured in TransPore to a comprehensive dual{scale formulation capable of addressing any drying configuration. Simulation results reported for a sample of spruce highlight the potential and flexibility of the new dual{scale approach. In particular, for a given unit cell configuration it is not necessary to supply the effective coefficients prior to each simulation.

Efficient simulation of unsaturated flow in heterogeneous soils using an exponential integrator

- Provided a practical variable-stepsize implementation of the exponential Euler method (EEM). - Introduced a new second-order variant of the scheme that enables the local error to be estimated at the cost of a single additional function evaluation. - New EEM implementation outperformed sophisticated implementations of the backward differentiation formulae (BDF) of order 2 and was competitive with BDF of order 5 for moderate to high tolerances.

Elevated concentration of TNF-a induces trophoblast differentiation in mouse blastocyst outgrowths

Elevated expression of tumour necrosis factora (TNF-a) is associated with adverse pregnancy outcome. This study has examined the expression of TNF-a and its receptors (TNF-Rs) by mouse blastocysts and blastocyst outgrowths from day 4 to 9.5 of pregnancy and investigated the effects of elevated TNF-a on the inner cell mass (ICM) and trophoblast cells of blastocyst outgrowths. RTPCR demonstrated TNF-a mRNA expression from day 7.5 to 9.5, TNF-R1 from day 6.5 to 9.5 and TNF-R2 from day 5.5 to 7.5 of pregnancy, and in situ hybridisation revealed the trophoblast giant cells (TGCs) of the early placenta as the site of TNF-a expression. Day 4 blastocysts were cultured in a physiologically high concentration of TNF-a (100 ng/ml) for 72 h to the outgrowth stage and then compared to blastocysts cultured in media alone. TNF-a-treated blastocyst outgrowths exhibited a significant reduction in ICM cells (mean € SD 23.90€10.42 vs 9.37€7.45, t-test, P<0.0001) with no significant change in the numbers of trophoblast cells (19.97€8.14 vs 21.73€7.79, t-test, P=0.39). Within the trophoblast cell population, the TNF-a-treated outgrowths exhibited a significant increase in multinucleated cells (14.10€5.53 vs 6.37€5.80, t-test, P<0.0001) and a corresponding significant decrease in mononucleated cells (5.87€3.60 vs 15.37€5.87, t-test, P<0.0001). In summary, this study describes the expression of TNF-a and its receptors during the peri-implantation period in the mouse. It also reports that elevated TNF-a restricts ICM proliferation in the blastocyst and changes the ratio of mononucleated to multinucleated trophoblast cells. These findings suggest a mechanism by which increased

Structural differentiation and corporate venturing : the moderating role of formal and informal integration mechanisms

Research has suggested that corporate venturing is crucial to strategic renewal and firm performance, yet scholars still debate the appropriate organizational configurations to facilitate the creation of new businesses in existing organizations. Our study investigates the effectiveness of combining structural differentiation with formal and informal organizational as well as top management team integration mechanisms in establishing an appropriate context for venturing activities. Our findings suggest that structural differentiation has a positive effect on corporate venturing. In addition, our study indicates that a shared vision has a positive effect on venturing in a structurally differentiated context. Socially integrated senior teams and cross-functional interfaces, however, are ineffective integration mechanisms for establishing linkages across differentiated units and for successfully pursuing corporate venturing.

A forward and backward secure key management in wireless sensor networks for PCS/SCADA

Process Control Systems (PCSs) or Supervisory Control and Data Acquisition (SCADA) systems have recently been added to the already wide collection of wireless sensor networks applications. The PCS/SCADA environment is somewhat more amenable to the use of heavy cryptographic mechanisms such as public key cryptography than other sensor application environments. The sensor nodes in the environment, however, are still open to devastating attacks such as node capture, which makes designing a secure key management challenging. In this paper, a key management scheme is proposed to defeat node capture attack by offering both forward and backward secrecies. Our scheme overcomes the pitfalls which Nilsson et al.'s scheme suffers from, and is not more expensive than their scheme.

Multilineage differentiation potential of bone and cartilage cells derived from explant culture

To date, mesenchymal stem cells (MSCs) from various tissues have been reported, but the yield and differentiation potential of different tissue-derived MSCs is still not clear. This study was undertaken in an attempt to investigate the multilineage stem cell potential of bone and cartilage explant cultures in comparison with bone marrow derived mesenchymal stem cells (BMSCs). The results showed that the surface antigen expression of tissue-derived cells was consistent with that of mesenchymal stem cells, such as lacking the haematopoietic and common leukocyte markers (CD34, CD45) while expressing markers related to adhesion (CD29, CD166) and stem cells (CD90, CD105). The tissue-derived cells were able to differentiate into osteoblast, chondrocyte and adipocyte lineage pathways when stimulated in the appropriate differentiating conditions. However, compared with BMSCs, tissue-derived cells showed less capacity for multilineage differentiation when the level of differentiation was assessed in monolayer culture by analysing the expression of tissue-specific genes by reverse transcription polymerase chain reaction (RT-PCR) and histology. In high density pellet cultures, tissue-derived cells were able to differentiate into chondrocytes, expressing chondrocyte markers such as proteoglycans, type II collagen and aggrecan. Taken together, these results indicate that cells derived from tissue explant cultures reserved certain degree of differentiation properties of MSCs in vitro.

Stem cell regulatory gene expression in human adult dental pulp and periodontal ligament cells undergoing odontogenic/osteogenic differentiation

Introduction During development and regeneration, odontogenesis and osteogenesis are initiated by a cascade of signals driven by several master regulatory genes. Methods In this study, we investigated the differential expression of 84 stem cell–related genes in dental pulp cells (DPCs) and periodontal ligament cells (PDLCs) undergoing odontogenic/osteogenic differentiation. Results Our results showed that, although there was considerable overlap, certain genes had more differential expression in PDLCs than in DPCs. CCND2, DLL1, and MME were the major upregulated genes in both PDLCs and DPCs, whereas KRT15 was the only gene significantly downregulated in PDLCs and DPCs in both odontogenic and osteogenic differentiation. Interestingly, a large number of regulatory genes in odontogenic and osteogenic differentiation interact or crosstalk via Notch, Wnt, transforming growth factor β (TGF-β)/bone morphogenic protein (BMP), and cadherin signaling pathways, such as the regulation of APC, DLL1, CCND2, BMP2, and CDH1. Using a rat dental pulp and periodontal defect model, the expression and distribution of both BMP2 and CDH1 have been verified for their spatial localization in dental pulp and periodontal tissue regeneration. Conclusions This study has generated an overview of stem cell–related gene expression in DPCs and PDLCs during odontogenic/osteogenic differentiation and revealed that these genes may interact through the Notch, Wnt, TGF-β/BMP, and cadherin signalling pathways to play a crucial role in determining the fate of dental derived cell and dental tissue regeneration. These findings provided a new insight into the molecular mechanisms of the dental tissue mineralization and regeneration

'Anorexia Nervosa' : asceticism, differentiation, government

This paper looks at the severe fasting practices most commonly found among young women. Almost all explanations for this behaviour centre around the notion of the pathological condition 'anorexia nervosa'. However, food asceticism has a well-documented history, particularly when it concerns religious fasting. In ancient Greece, dietary asceticism constituted an important part of the means by which individuals constructed an acceptable 'self'. Ascetic fasting then later resurfaced at various historical moments and in various different places — such as amongst medieval religious women and, in a broader way, amongst contemporary young women. It is argued that these practices have traditionally formed part of the mechanisms by which differentiation by age and sex occurs. Overall, it is hoped that this analysis will permit not only a different focus on 'anorexia nervosa', but also on some of the ways in which young people become gendered.

Asceticism, differentiation, government : 'anorexia nervosa' as an achievement

As with the ancient Greeks, the `self' is not something to be discovered, it is something to be created. Practices such as ascetic fasting are not expressions of the struggle between the authentic self and the external world, they are the very practices by which a `self' is formed.