A variable-stepsize Jacobian-free exponential integrator for simulating transport in heterogeneous porous media : application to wood drying

A Jacobian-free variable-stepsize method is developed for the numerical integration of the large, stiff systems of differential equations encountered when simulating transport in heterogeneous porous media. Our method utilises the exponential Rosenbrock-Euler method, which is explicit in nature and requires a matrix-vector product involving the exponential of the Jacobian matrix at each step of the integration process. These products can be approximated using Krylov subspace methods, which permit a large integration stepsize to be utilised without having to precondition the iterations. This means that our method is truly "Jacobian-free" - the Jacobian need never be formed or factored during the simulation. We assess the performance of the new algorithm for simulating the drying of softwood. Numerical experiments conducted for both low and high temperature drying demonstrates that the new approach outperforms (in terms of accuracy and efficiency) existing simulation codes that utilise the backward Euler method via a preconditioned Newton-Krylov 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.

Caustic potash pulping study of water efficient Arudno donax and Sorghum bicolor; a comparison to sugarcane bagasse

Numerous crops grow in sugar regions that have the potential to increase the amount of biomass available to a small bagasse-based pulp factory. Arundo donax and Sorghum offer unique advantages to farmers compared to other agricultural crops. Sorghum bicolour requires only 1/3 of the water of sugarcane. Arundo donax is a very high yield crop, it can also grow with little water but it has the further advantage in that it is also highly stress tolerant, making it suitable for land which is unsuited to other crops. Pulps produced from these crops were benchmarked against sugarcane bagasse pulp. Arundo, sorghum and bagasse were pulped using KOH and anthraquinone to 20 Kappa number so as to produce a bleachable pulp. The unbleached sorghum pulp has better tensile strength properties than the unbleached Arundo pulp (43.8 Nm/g compared to 21.4 Nm/g) and the bleached sorghum pulp tensile strength was similar to bagasse (28.4 Nm/g). At 20 Kappa number, sorghum pulp had acceptable yield for a non-wood fibre (45% c.f. 55% for bagasse), Arundo donax pulp had low tensile strength, and relatively low yield (38.7%), even for an agricultural fibre and required severe cooking conditions to achieve similar delignification to sugarcane bagasse or sorghum. Sorghum and Arundo donax produced thicker handsheets than bagasse (>160 μm c.f. 122 μm for bagasse). In preliminary experiments sorghum and bagasse responded slightly better to Totally Chlorine Free bleaching (QPP), although none achieved a satisfactory brightness level and more optimisation is needed.

A comprehensive dual-scale wood torrefaction model : application to the analysis of thermal run-away in industrial heat treatment processes

A dual-scale model of the torrefaction of wood was developed and used to study industrial configurations. At the local scale, the computational code solves the coupled heat and mass transfer and the thermal degradation mechanisms of the wood components. At the global scale, the two-way coupling between the boards and the stack channels is treated as an integral component of the process. This model is used to investigate the effect of the stack configuration on the heat treatment of the boards. The simulations highlight that the exothermic reactions occurring in each single board can be accumulated along the stack. This phenomenon may result in a dramatic eterogeneity of the process and poses a serious risk of thermal runaway, which is often observed in industrial plants. The model is used to explain how thermal runaway can be lowered by increasing the airflow velocity, the sticker thickness or by gas flow reversal.

Breaking the code of silence : the Wood Royal Commission into New South Wales police : a brief overview

The three-volume Final Report of the Wood inquiry into NSW Police (Royal Commission Into the New South Wales Police Service, 'Final Report, Vol I: Corruption; Vol II: Reform; Vol III: Appendices', May 1997) was publicly released on 15 May 1997, to much media fanfare. The Sydney Morning Herald (SMH) devoted an 8-page special report on I May to the pending release of the Inquiry Report, headed The Police Purge. On the day of the public release of the Report, the SMH five-page 'Special Report' under the banner The Police Verdict was headlined Wood, Carr Split on Drugs. The Australian led with Call for Drug Law Revamp, Force Overhaul to Fight Corruption, Wood Attacks Culture of Greed, and the Daily Telegraph front page 'Final Verdict' was True Blue Strategy for an Honest Police Force...

A comparison between bagasse and water efficient alternatives using KOH/AQ pulping

Numerous crops grow in sugar regions that have the potential to increase the amount of biomass available to a small bagasse-based pulp factory. Arundo donax and Sorghum offer unique advantages to farmers compared to other agricultural crops. Sorghum bicolour requires only 1/3 of the water of sugarcane. Arundo donax is a very high yield crop, it can also grow with little water but it has the further advantage in that it is also highly stress tolerant, making it suitable for land which is unsuited to other crops. Pulps produced from these crops were benchmarked against sugarcane bagasse pulp. Arundo, sorghum and bagasse were pulped using KOH and anthraquinone to 20 Kappa number so as to produce a bleachable pulp which is suitable for making photocopier paper and tissue products. The unbleached sorghum pulp has better tensile strength properties than the unbleached Arundo pulp (43.8 Nm/g compared to 21.4 Nm/g) and the bleached sorghum pulp tensile strength was similar to bagasse (28.4 Nm/g). At 20 Kappa number, sorghum pulp had acceptable yield for a non-wood fibre (45% c.f. 55% for bagasse), Arundo donax pulp had low tensile strength, and relatively low yield (38.7%), even for an agricultural fibre and required severe cooking conditions to achieve similar delignification to sugarcane bagasse or sorghum. Sorghum and Arundo donax produced thicker handsheets than bagasse (>160 µm c.f. 122 µm for bagasse). In preliminary experiments sorghum and bagasse responded slightly better to Totally Chlorine Free peroxide bleaching (QPP), although none achieved a satisfactory brightness level and further improvement would be required to produce a bleached pulp.

An investigation into the production of chloromethly furfual from bagasse and bagasse pulp

There are many attractive alternatives to produce chemicals similar to those currently produced from fossil fuel resources. The most viable renewable resource of fixed carbon is biomass. This paper examines processing conditions for the production and recovery of furanics from bagasse as well as bagasse pulp. It is shown that bio-oil consisting mainly of furanics (~84% chloromethly furfural) may be obtained in yields of ~78% and ~87% by weight from bagasse and bagasse pulp respectively using a biphasic acid hydrolysis system. The biphasic system consists of an organic layer of dichloroethane and an aqueous phase of concentrated hydrochloric acid. Generally the lower the impurity content and the higher the cellulose content, the higher the furanics yield.

Survival of rat functional dental pulp cells in vascularized tissue engineering chambers

Regenerative endodontics aims to preserve, repair or regenerate the dental pulp tissue. Dental pulp stem cells, have a potential use in dental tissue generation. However, specific requirements to drive the dental tissue generation are still obscured. We established an in vivo model for studying the survival of dental pulp cells (DPC) and their potential to generate dental pulp tissue. DPC were mixed with collagen scaffold with or without slow release bone morphogenic protein 4 (BMP-4) and fibroblast growth factor 2 (FGF2). The cell suspension was transplanted into a vascularized tissue engineering chamber in the rat groin. Tissue constructs were harvested after 2, 4, 6, and 8 weeks and processed for histomorphological and immunohistochemical analysis. After 2 weeks newly formed tissue with new blood vessel formation were observed inside the chamber. DPC were found around dentin, particularly around the vascular pedicle and also close to the gelatin microspheres. Cell survival, was confirmed up to 8 weeks after transplantation. Dentin Sialophosphoprotein (DSPP) positive matrix production was detected in the chamber, indicating functionality of dental pulp progenitor cells. This study demonstrates the potential of our tissue engineering model to study rat dental pulp cells and their behavior in dental pulp regeneration, for future development of an alternative treatment using these techniques.

Which industry is greener? An empirical study of nine industries in OECD countries

This study analyzed the relationship between the CO2 emissions of different industries and economic growth in OECD countries from 1970 to 2005. We tested an environmental Kuznets curve (EKC) hypothesis and found that total CO2 emissions from nine industries show an N-shaped trend instead of an inverted U or monotonic increasing trend with increasing income. The EKC hypothesis for sector-level CO2 emissions was supported in the (1) paper, pulp, and printing industry; (2) wood and wood products industry; and (3) construction industry. We also found that emissions from coal and oil increase with economic growth in the steel and construction industries. In addition, the non-metallic minerals, machinery, and transport equipment industries tend to have increased emissions from oil and electricity with economic growth. Finally, the EKC turning point and the relationship between GDP per capita and sectoral CO2 emissions differ among industries according to the fuel type used. Therefore, environmental policies for CO2 reduction must consider these differences in industrial characteristics. © 2013 Elsevier Ltd.

Wood density : a tool to find complementary species for the design of mixed species plantations

The successful establishment and growth of mixed-species forest plantations requires that complementary or facilitatory species be identified. This can be difficult in many tropical areas because the growth characteristics of endemic species are often unknown, particularly when grown at potentially higher densities in plantations than in natural forests. Here, we investigate whether wood density is a useful and readily accessible trait for choosing complementary species for mixed species plantations. Wood density represents the carbon investment per unit volume of stem with a trade-off generally found between fast (low wood density) and slow (high wood density) growing species. To do this, we use data collected from 18 highly diverse mixed species plantations (4–23 mostly native species) aged from 6 to 11 years at the time of data collection located on Leyte Island, Philippines. We found significant negative correlations between wood densities and the height of the most abundant species, as well as with measures of overall stand growth and tree diameter size distribution. Not only do species with denser woods have slower growth rates, but also mixed-species plantations with higher average wood density and higher stem density were also less productive, at least in these young plantations. Similarly, stands with a high diversity in wood densities were less productive. There is growing interest in making greater use of native multi-species mixtures in smallholder and community planting programs in the tropics, and our results show databases of wood density values may help improve their design. In the early development stages of plantations, canopy closure and rapid height growth are usually key silvicultural targets, and wood density values can predict the rapid height development of species. If plantations are being grown for the livelihood of small landholders then the best target is to choose some species with different wood densities. This allows an early harvest of low-wood density species for early income, and will also reduce competition for slower growing trees with higher wood densities for later income generation.

Comparing the effectiveness of education and technology in reducing wood smoke pollution: A field experiment

This study describes a field experiment assessing the effectiveness of education and technological innovation in reducing air pollution generated by domestic wood heaters. Two-hundred and twenty four households from a small regional center in Australia were randomly assigned to one of four experimental conditions: (1) Education only – households received a wood smoke reduction education pack containing information about the negative health impacts of wood smoke pollution, and advice about wood heater operation and firewood management; (2) SmartBurn only – households received a SmartBurn canister designed to improve combustion and help wood fires burn more efficiently, (3) Education and SmartBurn, and (4) neither Education nor SmartBurn (control). Analysis of covariance, controlling for pre-intervention household wood smoke emissions, wood moisture content, and wood heater age, revealed that education and SmartBurn were both associated with significant reduction in wood smoke emissions during the post-intervention period. Follow-up mediation analyses indicated that education reduced emissions by improving wood heater operation practices, but not by increasing health risk perceptions. As predicted, SmartBurn exerted a direct effect on emission levels, unmediated by wood heater operation practices or health risk perceptions.

“This is not a burning issue for me”: How citizens justify their use of wood heaters in a city with a severe air pollution problem

Although wood smoke pollution has been linked to health problems, wood burning remains a popular form of domestic heating in many countries across the world. In this paper, we describe the rhetoric of resistance to wood heater regulation amongst citizens in the regional Australian town of Armidale, where wood smoke levels regularly exceed national health advisory limits. We discuss how this is related to particular sources of resistance, such as affective attachment to wood heating and socio-cultural norms. The research draws on six focus groups with participants from households with and without wood heating. With reference to practice theory, we argue that citizen discourses favouring wood burning draw upon a rich suite of justifications and present this activity as a natural and traditional activity promoting comfort and cohesion. Such discourses also emphasise the identity of the town as a rural community and the supposed gemeinschaft qualities of such places. We show that, in this domain of energy policy, it is not enough to present ‘facts’ which have little emotional association or meaning for the populace. Rather, we need understand how social scripts, often localised, inform identity and practice.

The affect heuristic and public support for three types of wood smoke mitigation policies

This study applied the affect heuristic model to investigate key psychological factors (affective associations, perceived benefits, and costs of wood heating) contributing to public support for three distinct types of wood smoke mitigation policies: education, incentives, and regulation. The sample comprised 265 residents of Armidale, an Australian regional community adversely affected by winter wood smoke pollution. Our results indicate that residents with stronger positive affective associations with wood heating expressed less support for wood smoke mitigation policies involving regulation. This relationship was fully mediated by expected benefits and costs associated with wood heating. Affective associations were unrelated to public support for policies involving education and incentives, which were broadly endorsed by all segments of the community, and were more strongly associated with rational considerations. Latent profile analysis revealed no evidence to support the proposition that some community members experience internal “heart versus head” conflicts in which their positive affective associations with wood heating would be at odds with their risk judgments about the dangers of wood smoke pollution. Affective associations and cost/benefit judgments were very consistent with each other.