The analytical solution and numerical solution of the fractional diffusion-wave equation with damping

Fractional partial differential equations have been applied to many problems in physics, finance, and engineering. Numerical methods and error estimates of these equations are currently a very active area of research. In this paper we consider a fractional diffusionwave equation with damping. We derive the analytical solution for the equation using the method of separation of variables. An implicit difference approximation is constructed. Stability and convergence are proved by the energy method. Finally, two numerical examples are presented to show the effectiveness of this approximation.

Numerical solution of an optimal control model of dendritic cell treatment of a growing tumour

A new optimal control model of the interactions between a growing tumour and the host immune system along with an immunotherapy treatment strategy is presented. The model is based on an ordinary differential equation model of interactions between the growing tu- mour and the natural killer, cytotoxic T lymphocyte and dendritic cells of the host immune system, extended through the addition of a control function representing the application of a dendritic cell treat- ment to the system. The numerical solution of this model, obtained from a multi species Runge–Kutta forward-backward sweep scheme, is described. We investigate the effects of varying the maximum al- lowed amount of dendritic cell vaccine administered to the system and find that control of the tumour cell population is best effected via a high initial vaccine level, followed by reduced treatment and finally cessation of treatment. We also found that increasing the strength of the dendritic cell vaccine causes an increase in the number of natural killer cells and lymphocytes, which in turn reduces the growth of the tumour.

Lateral impact response and parametric studies of axially loaded square concrete filled steel tube columns

This paper presents a numerical study on the response of axially loaded slender square concrete filled steel tube (CFST) columns under low velocity lateral impact loading. A finite element analysis (FEA) model was developed using the explicit dynamic nonlinear finite element code LS -DYNA in which the strain rate effects of both steel and concrete, contact between steel tube and concrete and confinement effect provided by the steel tube for the concrete were considered. The model also benefited from a relatively recent feature of LS-DYNA for applying a pre-loading in the explicit solver. The developed numerical model was verified for its accuracy and adequacy by comparing the results with experimental results available in the literature. The verified model was then employed to conduct a parametric study to investigate the influence of axial load level, impact location, support conditions, and slenderness ratio on the response of the CFST columns. A good agreement between the numerical and experimental results was achieved. The model could reasonably predict the impact load-deflection history and deformed shape of the column at the end of the impact event. The results of the parametric study showed that whilst impact location, axial load level and slenderness ratio can have a significant effect on the peak impact force, residual lateral deflection and maximum lateral deflection, the influence of support fixity is minimal. With an increase of axial load to up to a certain level, the peak force increases; however, a further increase in the axial load causes a decrease in the peak force. Both residual lateral deflection and maximum lateral deflection increase as axial load level increases. Shifting the impact location towards the supports increases the peak force and reduces both residual and maximum lateral deflections. A rise in slenderness ratio decreases the peak force and increases the residual and maximum lateral deflections.

Study of the underlying electrochemistry of polycrystalline gold electrodes in aqueous solution and electrocatalysis by large amplitude Fourier transformed alternating current voltammetry

Polycrystalline gold electrodes of the kind that are routinely used in analysis and catalysis in aqueous media are often regarded as exhibiting relatively simple double-layer charging/discharging and monolayer oxide formation/ removal in the positive potential region. Application of the large amplitude Fourier transformed alternating current (FT-ac) voltammetric technique that allows the faradaic current contribution of fast electron-transfer processes to be emphasized in the higher harmonic components has revealed the presence of well-defined faradaic (premonolayer oxidation) processes at positive potentials in the double-layer region in acidic and basic media which are enhanced by electrochemical activation. These underlying quasi-reversible interfacial electron-transfer processes may mediate the course of electrocatalytic oxidation reactions of hydrazine, ethylene glycol, and glucose on gold electrodes in aqueous media. The observed responses support key assumptions associated with the incipient hydrous oxide adatom mediator (IHOAM) model of electrocatalysis.

Auto-inhibition of hydrogen gas evolution on gold in aqueous acid solution

It was demonstrated recently that dramatic changes in the redox behaviour of gold/aqueous solution interfaces may be observed following either cathodic or thermal electrode pretreatment. Further work on the cathodic pretreatment of gold in acid solution revealed that as the activity of the gold surface was increased, its performance as a substrate for hydrogen gas evolution under constant potential conditions deteriorated. The change in activity of the gold atoms at the interface, which was attributed to a hydrogen embrittlement process (the occurrence of the latter was subsequently checked by surface microscopy), was confirmed, as in earlier work, by the appearance of a substantial anodic peak at ca. 0.5 V (RHE) in a post-activation positive sweep. Changes in the catalytic activity of a metal surface reflect the fact that the structure (or topography), thermodynamic activity and electronic properties of a surface are dependent not only on pretreatment but also, in the case of the hydrogen evolution reaction, vary with time during the course of reaction. As will be reported shortly, similar (and often more dramatic) time-dependent behaviour was observed for hydrogen gas evolution on other metal electrodes.

The role of defects, or active states, in surface electrochemistry with particular reference to gold in neutral solution

Metastable, active, or nonequilibrium states due to the presence of abnormal structures and various types of defects are well known in metallurgy. The role of such states at gold surfaces in neutral aqueous media (an important electrode system in the microsensor area) was explored using cyclic voltammetry. It was demonstrated that, as postulated in earlier work from this laboratory, there is a close relationship between premonolayer oxidation, multilayer hydrous oxide reduction and electrocatalytic behaviour in the case of this and other metal electrode systems. Some of the most active, and therefore most important, entities at surfaces (e.g., metal adatoms) are not readily imageable or detectable by high resolution surface microscopy techniques. Cyclic voltammetry, however, provides significant, though not highly specific, information about such species. The main conclusion is that further practical and theoretical work on active states of metal surfaces is highly desirable as their behaviour is not simple and is of major importance in many electrocatalytic processes.

The formation of gold nanoparticles using hydroquinone as a reducing agent through a localized pH change upon addition of NaOH to a solution of HAuCl4

We demonstrate a rapid synthesis of gold nanoparticles using hydroquinone as a reducing agent under acidic conditions without the need for precursor seed particles. The nanoparticle formation process is facilitated by the addition of NaOH to a solution containing HAuCl4 and hydroquinone to locally change the pH; this enhances the reducing capability of hydroquinone to form gold nucleation centres, after which further growth of gold can take place through an autocatalytic mechanism. The stability of the nanoparticles is highly dependent on the initial solution pH, and both the concentration of added NaOH and hydroquinone present in solution. The gold nanoparticles were characterized by UV–visible spectroscopy, transmission electron microscopy, Fourier transform infrared spectroscopy, atomic force microscopy, dynamic light scattering, and zeta potential measurements. It was found that under optimal conditions that stable aqueous suspensions of 20 nm diameter nanoparticles can be achieved where benzoquinone, the oxidized product of hydroquinone, acts as a capping agent preventing nanoparticles aggregation.

The electrochemical corrosion behaviour of quaternary gold alloys when exposed to 3.5% NaCl solution

Lower carat gold alloys, specifically 9 carat gold alloys, containing less than 40 % gold, and alloying additions of silver, copper and zinc, are commonly used in many jewellery applications, to offset high costs and poor mechanical properties associated with pure gold. While gold is considered to be chemically inert, the presence of active alloying additions raises concerns about certain forms of corrosion, particularly selective dissolution of these alloys. The purpose of this study was to systematically study the corrosion behaviour of a series of quaternary gold–silver–copper–zinc alloys using dc potentiodynamic scanning in saline (3.5 % NaCl) environment. Full anodic/cathodic scans were conducted to determine the overall corrosion characteristics of the alloy, followed by selective anodic scans and subsequent morphological and compositional analysis of the alloy surface and corroding media to determine the extent of selective dissolution. Varying degrees of selective dissolution and associated corrosion rates were observed after anodic polarisation in 3.5 % NaCl, depending on the alloy composition. The corrosion behaviour of the alloys was determined by the extent of anodic reactions which induce (1) formation of oxide scales on the alloy surface and or (2) dissolution of Zn and Cu species. In general, the improved corrosion characteristics of alloy #3 was attributed to the composition of Zn/Cu in the alloy and thus favourable microstructure promoting the formation of protective oxide/chloride scales and reducing the extent of Cu and Zn dissolution.

A review of the removal of anions and oxyanions of the halogen elements from aqueous solution by layered double hydroxides

The application of layered double hydroxides (LDHs) and thermally activated LDHs for the removal of various fluorine (F-, BF-4), chlorine (Cl-,ClO-4), bromine (Br-, BrO-3) and iodine (I-, IO-3) species from aqueous solutions has been reviewed in this article. LDHs and thermally activated LDHs were able to significantly reduce the concentration of selected anions in laboratory scale experiments. The M2+:M3+ cation ratio of the LDH adsorbent was an important factor which influenced anion uptake. Though LDHs were able to remove some target anion species through anion exchange and surface adsorption thermal activation and reformation generally produced better results. The presence of competing anions including carbonate, phosphate and sulphate had a significant impact on uptake of the target anion as LDHs typically exhibit lower affinity towards monovalent anions compared to anions with multiple charges. The removal of fluoride and perchlorate from aqueous solution by a continuous flow system utilising fixed bed columns packed with LDH adsorbents has also been investigated. The adsorption capacity of the columns at breakpoint was heavily dependent on the flow rate and lower than result reported for the corresponding batch methods. There is still considerable scope for future research on numerous topics summarised in this article.

Numerical simulations of particle deposition in metal foam heat exchangers

Australia is a high-potential country for geothermal power with reserves currently estimated in the tens of millions of petajoules, enough to power the nation for at least 1000 years at current usage. However, these resources are mainly located in isolated arid regions where water is scarce. Therefore, wet cooling systems for geothermal plants in Australia are the least attractive solution and thus air-cooled heat exchangers are preferred. In order to increase the efficiency of such heat exchangers, metal foams have been used. One issue raised by this solution is the fouling caused by dust deposition. In this case, the heat transfer characteristics of the metal foam heat exchanger can dramatically deteriorate. Exploring the particle deposition property in the metal foam exchanger becomes crucial. This paper is a numerical investigation aimed to address this issue. Two dimensional (2D) numerical simulations of a standard one-row tube bundle wrapped with metal foam in cross-flow are performed and highlight preferential particle deposition areas.

Numerical simulations of particle deposition in metal foam heat exchangers

Australia is a high potential country for geothermal power with reserves currently estimated in the tens of millions of petajoules, enough to power the nation for at least 1000 years at current usage.However, these resources are mainly located in isolated arid regions where water is scarce. Therefore, wet cooling systems for geothermal plants in Australia are the least attractive solution and thus air-cooled heat exchangers are preferred. In order to increase the efficiency of such heat exchangers, metal foams have been used. One issue raised by this solution is the fouling caused by dust deposition. In this case, the heat transfer characteristics of the metal foam heat exchanger can dramatically deteriorate. Exploring the particle deposition property in the metal foam exchanger becomes crucial. This paper is a numerical investigation aimed to address this issue. Two-dimensional(2D numerical simulations of a standard one-row tube bundle wrapped with metal foam in cross-flow are performed and highlight preferential particle deposition areas.

Observation of solution-induced corneal staining with fluorescein, rose bengal and lissamine green

Ever since sodium fluorescein (‘fluorescein’ [FL]) was first used to investigate the ocular surface over a century ago, the term ‘staining’ has been taken to mean the presence of ocular surface fluorescence [1]. This term has not been necessarily taken to infer any particular mechanism of causation, and indeed, can be attributed to a variety of possible aetiologies [2]. In recent times, there has been considerable interest in a form of ocular surface fluorescence seen in association with the use of certain combinations of soft contact lenses and multipurpose solutions. The first clinical account of this phenomenon was reported by Jones et al. [3], which was followed by a more formal investigation by the same author in 2002 [4]. Jones et al described this appearance as a ‘classic solution-based toxicity reaction’. Subsequently, this appearance has come to be known as ‘solution-induced corneal staining’ or more recently by the acronym ‘SICS’ [5]. The term SICS is potentially problematic in that from a cell biology point of view, there is an inference that ‘staining’ means the entry of a dye into corneal epithelial cells. Morgan and Maldonado-Codina [2] noted there was no foundation of solid scientific literature underpinning our understanding of the true basic causative mechanisms of this phenomenon; since that time, further work has been published in this field [6] and [7] but questions still remain about the precise aetiology of this phenomenon...

Strong and bioactive tri-calcium phosphate scaffolds with tube-like macropores

Calcium phosphate ceramic scaffolds have been widely investigated for bone tissue engineering due to their excellent biocompatibility and biodegradation. Unfortunately, they have the shortcoming of low mechanical properties. In order to provide strong, bioactive, and biodegradable scaffolds, a new approach of infiltrating the macro-tube ABS (acrylontrile butadiene styrene) templates with a hydroxyapatite/bioactive glass mixed slurry was developed to fabricate porous Si-doped TCP (tri-calcium phosphate) scaffolds. The porous Si-doped TCP ceramics with a high porosity (~65%) and with interconnected macrotubes (~0.8mm in diameter) and micropores (5-100 m) had a high compressive strength (up to 14.68+0.2MPa), which was comparable to that of a trabecular bone and was much higher than those of pure TCP scaffolds. Additional cell attachment study and MTT cytotoxicity assay proved the bioactivity and biocompatibility of the new scaffolds. Thus a potential bioceramic material and a new approach to make the potential scaffolds were developed for bone tissue engineering.

A finite volume method with linearisation in time for the solution of advection-reaction-diffusion systems

The numerical solution in one space dimension of advection--reaction--diffusion systems with nonlinear source terms may invoke a high computational cost when the presently available methods are used. Numerous examples of finite volume schemes with high order spatial discretisations together with various techniques for the approximation of the advection term can be found in the literature. Almost all such techniques result in a nonlinear system of equations as a consequence of the finite volume discretisation especially when there are nonlinear source terms in the associated partial differential equation models. This work introduces a new technique that avoids having such nonlinear systems of equations generated by the spatial discretisation process when nonlinear source terms in the model equations can be expanded in positive powers of the dependent function of interest. The basis of this method is a new linearisation technique for the temporal integration of the nonlinear source terms as a supplementation of a more typical finite volume method. The resulting linear system of equations is shown to be both accurate and significantly faster than methods that necessitate the use of solvers for nonlinear system of equations.

Measurement of the first neutron diffraction peak of D2O flowing in a 2 mm diameter aluminium tube as a function of temperature

A pilot experiment was performed using the WOMBAT powder diffraction instrument at ANSTO in which the first neutron diffraction peak (Q0) was measured for D2O flowing in a 2 mm internal diameter aluminium tube. Measurements of Q0 were made at -9, 4.3, 6.9, 12, 18.2 and 21.5 °C. The D2O was circulated using a siphon with water in the lower reservoir returned to the upper reservoir using a small pump. This enabled stable flow to be maintained for several hours. For example, if the pump flow increased slightly, the upper reservoir level rose, increasing the siphon flow until it matched the return flow. A neutron wavelength of 2.4 Å was used and data integrated over 60 minutes for each temperature. A jet of nitrogen from a liquid N2 Dewar was directed over the aluminium tube to vary water temperature. After collection of the data, the d spacing of the aluminium peaks was used to calculate the temperature of the aluminium within the neutron beam and therefore was considered to be an accurate measure of water temperature within the beam. Sigmaplot version 12.3 was used to fit a Weibull five parameter peak fit to the first neutron diffraction peak. The values of Q0 obtained in this experiment showed an increase with temperature consistent with data in the literature [1] but were consistently higher than published values for bulk D20. For example at 21.5 °C we obtained a value of 2.008 Å-1 for Q0 compared to a literature value of 1.988 Å-1 for bulk D2O at 20 °C, a difference of 1%. Further experiments are required to see if this difference is real or artifactual.