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.

Multicomponent charge transport in electrolyte solutions

The work presented in this thesis investigates the mathematical modelling of charge transport in electrolyte solutions, within the nanoporous structures of electrochemical devices. We compare two approaches found in the literature, by developing onedimensional transport models based on the Nernst-Planck and Maxwell-Stefan equations. The development of the Nernst-Planck equations relies on the assumption that the solution is infinitely dilute. However, this is typically not the case for the electrolyte solutions found within electrochemical devices. Furthermore, ionic concentrations much higher than those of the bulk concentrations can be obtained near the electrode/electrolyte interfaces due to the development of an electric double layer. Hence, multicomponent interactions which are neglected by the Nernst-Planck equations may become important. The Maxwell-Stefan equations account for these multicomponent interactions, and thus they should provide a more accurate representation of transport in electrolyte solutions. To allow for the effects of the electric double layer in both the Nernst-Planck and Maxwell-Stefan equations, we do not assume local electroneutrality in the solution. Instead, we model the electrostatic potential as a continuously varying function, by way of Poisson’s equation. Importantly, we show that for a ternary electrolyte solution at high interfacial concentrations, the Maxwell-Stefan equations predict behaviour that is not recovered from the Nernst-Planck equations. The main difficulty in the application of the Maxwell-Stefan equations to charge transport in electrolyte solutions is knowledge of the transport parameters. In this work, we apply molecular dynamics simulations to obtain the required diffusivities, and thus we are able to incorporate microscopic behaviour into a continuum scale model. This is important due to the small size scales we are concerned with, as we are still able to retain the computational efficiency of continuum modelling. This approach provides an avenue by which the microscopic behaviour may ultimately be incorporated into a full device-scale model. The one-dimensional Maxwell-Stefan model is extended to two dimensions, representing an important first step for developing a fully-coupled interfacial charge transport model for electrochemical devices. It allows us to begin investigation into ambipolar diffusion effects, where the motion of the ions in the electrolyte is affected by the transport of electrons in the electrode. As we do not consider modelling in the solid phase in this work, this is simulated by applying a time-varying potential to one interface of our two-dimensional computational domain, thus allowing a flow field to develop in the electrolyte. Our model facilitates the observation of the transport of ions near the electrode/electrolyte interface. For the simulations considered in this work, we show that while there is some motion in the direction parallel to the interface, the interfacial coupling is not sufficient for the ions in solution to be "dragged" along the interface for long distances.

Simultaneous determination of 6,7-dimethylesculetin and geniposide in rat plasma and its application to pharmacokinetic studies of Yin Chen Hao Tang preparation

A method for the rapid and simultaneous determination of 6,7-dimethylesculetin (CAS 120-08-1) and geniposide (CAS 24512-63-8) in rat plasma has been developed, using validated high performance liquid chromatography (HPLC) with solid phase extraction (SPE). The HPLC analysis was performed on a commercially available column (200 mm x 4.6 mm, 5 microm) with acetonitrile-methanol-0.1% aqueous formic acid as mobile phase and the UV detection at 343 nm and 238 nm for 6,7-dimethylesculetin and geniposide, respectively. The calibration curves for 6,7-dimethylesculetin and geniposide were linear over the range 0.4-25.6 microg/mL and 1.12-71.68 microg/mL, respectively. The lower limits of quantitation were 0.40 microg/ mL and 1.12 microg/mL, and the lower limits of detection were 0.06 microg/mL and 0.09 microg/ mL, respectively. The intra-day and inter-day precision for 6,7-dimethylesculetin and geniposide were < 5%, whereas the absolute recovery percentages were > 74%. A successful application of the developed HPLC analysis was demonstrated for the pharmacokinetic study of a Traditional Chinese Medicine formula of Yin Chen Hao Tang preparation.

Voltammetric, spectroscopic, and microscopic investigations of electrocrystallized forms of semiconducting AgTCNQ (TCNQ=7,7,8,8-tetracyanoquinodimethane) exhibiting different morphologies and colors

Chemically synthesized AgTCNQ exists in two forms that differ in their morphologies (needles and microcrystals) and colors (red and blue). It is now shown that both forms exhibit essentially indistinguishable X-ray diffraction, spectroscopic, and thermochemical data, implying that they are not separate phases, as implied in some literature. Electrochemical reduction of TCNQ((MeCN)) in the presence of Ag+((MeCN)) generates both red and blue AgTCNQ. On glassy carbon, platinum, or indium tin oxide electrodes and at relatively positive deposition potentials, slow growth of high aspect ratio, red needle AgTCNQ crystals occurs. After longer times and at more negative deposition potentials, blue microcrystalline AgTCNQ thin films are favored. Blue AgTCNQ is postulated to be generated via reduction of a Ag+\[(TCNQ(center dot-))(TCNQ)]((MeCN)) intermediate. At even more negative potentials, Ag-(metal) formation inhibits further growth of AgTCNQ. On a gold electrode, Ag-(metal)) deposition occurs at more positive potentials than on the other electrode materials examined. However, surface plasmon resonance data indicate (hat a small potential region is available between the stripping of Ag-(metal)) and the oxidation of TCNQ(center dot-)(MeCN) back to TCNQ(MeCN) where AgTCNQ may form. AgTCNQ in both the red and blue forms also can be prepared electrochemically on a TCNQ((s)) modified electrode in -0.1 M AgNO3(aq) where deposition of Ag(m,,,I) onto the TCNQ((s)) crystals allows a charge transfer process to occur. However, the morphology formed in this solid-solid phase transformation is more difficult to control.

Short term variability in urinary bisphenol A in Australian children

Used frequently in food contact materials, bisphenol A (BPA) has been studied extensively in recent years, and ubiquitous exposure in the general population has been demonstrated worldwide. Characterising within- and between-individual variability of BPA concentrations is important for characterising exposure in biomonitoring studies, and this has been investigated previously in adults, but not in children. The aim of this study was to characterise the short-term variability of BPA in spot urine samples in young children. Children aged ≥2-<4 years (n = 25) were recruited from an existing cohort in Queensland Australia, and donated four spot urine samples each over a two day period. Samples were analysed for total BPA using isotope dilution online solid phase extraction-liquid chromatography-tandem mass spectrometry, and concentrations ranged from 0.53–74.5 ng/ml, with geometric mean and standard deviation of 2.70 ng/ml and 2.94 ng/ml, respectively. Sex and time of sample collection were not significant predictors of BPA concentration. The between-individual variability was approximately equal to the within-individual variability (ICC = 0.51), and this ICC is somewhat higher than previously reported literature values. This may be the result of physiological or behavioural differences between children and adults or of the relatively short exposure window assessed. Using a bootstrapping methodology, a single sample resulted in correct tertile classification approximately 70% of the time. This study suggests that single spot samples obtained from young children provide a reliable characterization of absolute and relative exposure over the short time window studied, but this may not hold true over longer timeframes.

Analysis of albumin-associated peptides and proteins from ovarian cancer patients

Albumin binds low–molecular-weight molecules, including proteins and peptides, which then acquire its longer half-life, thereby protecting the bound species from kidney clearance. We developed an experimental method to isolate albumin in its native state and to then identify [mass spectrometry (MS) sequencing] the corresponding bound low–molecular-weight molecules. We used this method to analyze pooled sera from a human disease study set (high-risk persons without cancer, n= 40; stage I ovarian cancer, n = 30; stage III ovarian cancer, n = 40) to demonstrate the feasibility of this approach as a discovery method. Methods Albumin was isolated by solid-phase affinity capture under native binding and washing conditions. Captured albumin-associated proteins and peptides were separated by gel electrophoresis and subjected to iterative MS sequencing by microcapillary reversed-phase tandem MS. Selected albumin-bound protein fragments were confirmed in human sera by Western blotting and immunocompetition. Results In total, 1208 individual protein sequences were predicted from all 3 pools. The predicted sequences were largely fragments derived from proteins with diverse biological functions. More than one third of these fragments were identified by multiple peptide sequences, and more than one half of the identified species were in vivo cleavage products of parent proteins. An estimated 700 serum peptides or proteins were predicted that had not been reported in previous serum databases. Several proteolytic fragments of larger molecules that may be cancer-related were confirmed immunologically in blood by Western blotting and peptide immunocompetition. BRCA2, a 390-kDa low-abundance nuclear protein linked to cancer susceptibility, was represented in sera as a series of specific fragments bound to albumin. Conclusion Carrier-protein harvesting provides a rich source of candidate peptides and proteins with potential diverse tissue and cellular origins that may reflect important disease-related information.

Methods for extracting genomic DNA from whole blood samples : current perspectives

Deoxyribonucleic acid (DNA) extraction has considerably evolved since it was initially performed back in 1869. It is the first step required for many of the available downstream applications used in the field of molecular biology. Whole blood samples are one of the main sources used to obtain DNA, and there are many different protocols available to perform nucleic acid extraction on such samples. These methods vary from very basic manual protocols to more sophisticated methods included in automated DNA extraction protocols. Based on the wide range of available options, it would be ideal to determine the ones that perform best in terms of cost-effectiveness and time efficiency. We have reviewed DNA extraction history and the most commonly used methods for DNA extraction from whole blood samples, highlighting their individual advantages and disadvantages. We also searched current scientific literature to find studies comparing different nucleic acid extraction methods, to determine the best available choice. Based on our research, we have determined that there is not enough scientific evidence to support one particular DNA extraction method from whole blood samples. Choosing a suitable method is still a process that requires consideration of many different factors, and more research is needed to validate choices made at facilities around the world.

Bisphenol A exposure is not associated with area-level socioeconomic index in Australian children using pooled urine samples

Bisphenol A (BPA) is used extensively in food-contact materials and has been detected routinely in populations worldwide, and this exposure has been linked to a range of negative health outcomes in humans. There is some evidence of an association between BPA and different socioeconomic variables which may be the result of different dietary patterns. The aim of this study was to conduct a preliminary investigation of the association between BPA and socioeconomic status in Australian children using pooled urine specimens and an area level socioeconomic index. Surplus pathology urine specimens collected from children aged 0-15 years in Queensland, Australia as samples of convenience (n = 469) were pooled by age, sex and area level socioeconomic index (n = 67 pools), and analysed for total BPA using online solid phase extraction LC-MS/MS. Concentration ranged from 1.08-27.4 ng/ml with geometric mean 2.57 ng/ml, and geometric mean exposure was estimated as 70.3 ng/kg d-1. Neither BPA concentration nor excretion was associated with age or sex, and the authors found no evidence of an association with socioeconomic status. These results suggest that BPA exposure is not associated with socioeconomic status in the Australian population due to relatively homogenous exposures in Australia, or that the socioeconomic gradient is relatively slight in Australia compared with other OECD countries.

Decline in perfluorooctane sulfonate and perfluorooctanoate serum concentrations in an Australian population from 2002 to 2011

Some perfluoroalkyl and polyfluoroalkyl substances (PFASs) have become widespread pollutants detected in human and wildlife samples worldwide. The main objective of this study was to assess temporal trends of PFAS concentrations in human blood in Australia over the last decade (2002–2011), taking into consideration age and sex trends. Pooled human sera from 2002/03 (n=26); 2008/09 (n=24) and 2010/11 (n=24) from South East Queensland, Australia were obtained from de-identified surplus pathology samples and compared with samples collected previously from 2006/07 (n=84). A total of 9775 samples in 158 pools were available for assessment of PFASs. Stratification criteria included sex and age: <16 years (2002/03 only); 0–4 (2006/07, 2008/09, 2010/11); 5–15 (2006/07, 2008/09, 2010/11); 16–30; 31–45; 46–60; and >60 years (all collection periods). Sera were analyzed using on-line solid-phase extraction coupled to high-performance liquid chromatography-isotope dilution-tandem mass spectrometry. Perfluorooctane sulfonate (PFOS) was detected in the highest concentrations ranging from 5.3–19.2 ng/ml (2008/09) to 4.4–17.4 ng/ml (2010/11). Perfluorooctanoate (PFOA) was detected in the next highest concentration ranging from 2.8–7.3 ng/ml (2008/09) to 3.1–6.5 ng/ml (2010/11). All other measured PFASs were detected at concentrations <1 ng/ml with the exception of perfluorohexane sulfonate which ranged from 1.2–5.7 ng/ml (08/09) and 1.4–5.4 ng/ml (10/11). The mean concentrations of both PFOS and PFOA in the 2010/11 period compared to 2002/03 were lower for all adult age groups by 56%. For 5-15 year olds, the decrease was 66% (PFOS) and 63% (PFOA) from 2002/03 to 2010/11. For 0-4 year olds the decrease from 2006/07 (when data were first available for this age group) was 50% (PFOS) and 22% (PFOA). This study provides strong evidence for decreasing serum PFOS and PFOA concentrations in an Australian population from 2002 through 2011. Age trends were variable and concentrations were higher in males than females. Global use has been in decline since around 2002 and hence primary exposure levels are expected to be decreasing. Further biomonitoring will allow assessment of PFAS exposures to confirm trends in exposure as primary and eventually secondary sources are depleted.

Breaking the DNA-vaccine bottleneck

Monash University in Australia has developed a new approach towards DNA vaccine development that has the potential to cut the time it takes to produce a vaccine from up to nine months to four weeks or less. The university has designed and filed a patent on a commercially viable, single-stage technology for manufacturing DNA molecules. The technology was used to produce malaria and measles DNA vaccines, which were tested to be homogeneous supercoiled DNA, free from RNA and protein contaminations and meeting FDA regulatory standards for DNA vaccines. The technique is based on customized, smart, polymeric, monolithic adsorbents that can purify DNA very rapidly. The design criteria of solid-phase adsorbent include rapid adsorption and desorption kinetics, physical composition, and adequate selectivity , capacity and recovery. The new show technology significantly improved binding capacities, higher recovery, drastically reduced use of buffers and processing time, less clogging, and higher yields of DNA.

Preparation, Analysis and Use of an Affinity Adsorbent for the Purification of GST Fusion Protein

Methods are presented for the preparation, ligand density analysis and use of an affinity adsorbent for the purification of a glutathione S-transferase (GST) fusion protein in packed and expanded bed chromatographic processes. The protein is composed of GST fused to a zinc finger transcription factor (ZnF). Glutathione, the affinity ligand for GST purification, is covalently immobilized to a solid-phase adsorbent (Streamline™). The GST–ZnF fusion protein displays a dissociation constant of 0.6 x10-6 M to glutathione immobilized to Streamline™. Ligand density optimization, fusion protein elution conditions (pH and glutathione concentration) and ligand orientation are briefly discussed.

Plasmid DNA purification via the use of a dual affinity protein

Methods are presented for the production, affinity purification and analysis of plasmid DNA (pDNA). Batch fermentation is used for the production of the pDNA, and expanded bed chromatography, via the use of a dual affinity glutathione S-transferase (GST) fusion protein, is used for the capture and purification of the pDNA. The protein is composed of GST, which displays affinity for glutathione immobilized to a solid-phase adsorbent, fused to a zinc finger transcription factor, which displays affinity for a target 9-base pair sequence contained within the target pDNA. A Picogreen™ fluorescence assay and/or anx ethidium bromide agarose gel electrophoresis assay can be used to analyze the eluted pDNA.

Microscale consolidation analysis of relaxation behavior of single living chondrocytes subjected to varying strain-rates

Besides the elastic stiffness, the relaxation behavior of single living cells is also of interest of various researchers when studying cell mechanics. It is hypothesized that the relaxation response of the cells is governed by both intrinsic viscoelasticity of the solid phase and fluid-solid interactions mechanisms. There are a number of mechanical models have been developed to investigate the relaxation behavior of single cells. However, there is lack of model enable to accurately capture both of the mechanisms. Therefore, in this study, the porohyperelastic (PHE) model, which is an extension of the consolidation theory, combined with inverse Finite Element Analysis (FEA) technique was used at the first time to investigate the relaxation response of living chondrocytes. This model was also utilized to study the dependence of relaxation behavior of the cells on strain-rates. The stress-relaxation experiments under the various strain-rates were conducted with the Atomic Force Microscopy (AFM). The results have demonstrated that the PHE model could effectively capture the stress-relaxation behavior of the living chondrocytes, especially at intermediate to high strain-rates. Although this model gave some errors at lower strain-rates, its performance was acceptable. Therefore, the PHE model is properly a promising model for single cell mechanics studies. Moreover, it has been found that the hydraulic permeability of living chondrocytes reduced with decreasing of strain-rates. It might be due to the intracellular fluid volume fraction and the fluid pore pressure gradients of chondrocytes were higher when higher strain-rates applied.

Characterization of epitopes for virus-neutralizing monoclonal antibodies to Ross River virus E2 using phage-displayed random peptide libraries

Ross River virus (RRV) is the predominant cause of epidemic polyarthritis in Australia, yet the antigenic determinants are not well defined. We aimed to characterize epitope(s) on RRV-E2 for a panel of monoclonal antibodies (MAbs) that recognize overlapping conformational epitopes on the E2 envelope protein of RRV and that neutralize virus infection of cells in vitro. Phage-displayed random peptide libraries were probed with the MAbs T1E7, NB3C4, and T10C9 using solution-phase and solid-phase biopanning methods. The peptides VSIFPPA and KTAISPT were selected 15 and 6 times, respectively, by all three of the MAbs using solution-phase biopanning. The peptide LRLPPAP was selected 8 times by NB3C4 using solid-phase biopanning; this peptide shares a trio of amino acids with the peptide VSIFPPA. Phage that expressed the peptides VSIFPPA and LRLPPAP were reactive with T1E7 and/or NB3C4, and phage that expressed the peptides VSIFPPA, LRLPPAP, and KTAISPT partially inhibited the reactivity of T1E7 with RRV. The selected peptides resemble regions of RRV-E2 adjacent to sites mutated in neutralization escape variants of RRV derived by culture in the presence of these MAbs (E2 210-219 and 238-245) and an additional region of E2 172-182. Together these sites represent a conformational epitope of E2 that is informative of cellular contact sites on RRV.

Immunochemical approach to the mapping of an assembled epitope of human chorionic gonadotropin: Proximity of CTP-alpha to the receptor binding region of the beta-subunit

A single-step solid-phase RIA (SS-SPRIA) developed in our laboratory using hybridoma culture supernatants has been utilised for the quantitation of epitope-paratope interactions. Using SS-SPRIA as a quantitative tool for the assessment of epitope stability, it was found that several assembled epitopes of human chorionic gonadotropin (hCG) are differentially stable to proteolysis and chemical modification. Based on these observations an approach has been developed for identifying the amino acid residues constituting an epitopic region. This approach has now been used to map an assembled epitope at/near the receptor binding region of the hormone. The mapped site forms a part of the seat belt region and the cystine knot region (C34-C38-C88-C90-H106). The carboxy terminal region of the alpha-subunit forms a part of the epitope indicating its proximity to the receptor binding region. These results are in agreement with the reported receptor binding region identified through other approaches and the X-ray crystal structure of hCG.