Biochar type and pyrolysis temperature effects on soil quality indicators and structural stability.

Quality of soils of the arid zones with low organic matter can be improved through the application of natural amendments especially biochar from various available feedstocks. The objective of this study was to evaluate the impacts of corn residue and poultry manure and their biochars on soil organic carbon (SOC), hot-water soluble carbohydrates (HWSC), basal soil respiration (BSR),and structural stability determined by HEMC (high-energy moisture characteristic). A sandy loam soil in pots were thoroughly mixed with 1, 2 and 4% w/w of corn residues (CR) and poultry manure (PM) feedstock and their biochars prepared at 350 and 650 °C of slow pyrolysis. Maize seeds were planted in pots and grown until physiological maturity when soil characteristics were measured. Treatments considerably altered the means of studied soil quality indicators, and increased SOC (1.5-10 times) and HWSC (1-7 times), and HEMC indices: volume of drainable pores ratio (VDPR, 1.5- 3.5 times), and stability ratio (SR, 1-3 times). Increasing pyrolysis temperature, regardless of the type and rate of feedstock, significantly decreased the SOC, BSR and percent of water-stable aggregates, and consequently structural stability indices. Contribution of both PM feedstock and its biochars was less effective than the CR ones (particularly the biochar produced at higher pyrolysis temperature), due to elevated sodium adsorption ratio (SAR) associated with higher slaking, physico-chemical dispersion and lower aggregate and structural stability.

Modelling the effect of acoustic waves on the thermodynamics and kinetics of phase transformation in a solution: Including mass transportation.

Effects of acoustic waves on a phase transformation in a metastable phase were investigated in our previous work [S. R. Haqshenas, I. J. Ford, and N. Saffari, "Modelling the effect of acoustic waves on nucleation," J. Chem. Phys. 145, 024315 (2016)]. We developed a non-equimolar dividing surface cluster model and employed it to determine the thermodynamics and kinetics of crystallisation induced by an acoustic field in a mass-conserved system. In the present work, we developed a master equation based on a hybrid Szilard-Fokker-Planck model, which accounts for mass transportation due to acoustic waves. This model can determine the kinetics of nucleation and the early stage of growth of clusters including the Ostwald ripening phenomenon. It was solved numerically to calculate the kinetics of an isothermal sonocrystallisation process in a system with mass transportation. The simulation results show that the effect of mass transportation for different excitations depends on the waveform as well as the imposed boundary conditions and tends to be noticeable in the case of shock waves. The derivations are generic and can be used with any acoustic source and waveform.

The dynamic excitation of a granular chain: Contact mechanics finite element analysis and experimental validation.

There is currently interest in transmitting acoustic signals along granular chains to produce waveforms of relevance to biomedical ultrasound applications. The study of such a transduction mechanism is greatly aided by the use of validated theoretical models. In view of this, a finite element analysis is presented in this paper. The dynamics of a granular chain of six, 1 mm diameter chrome steel spherical beads, was excited at one end using a sinusoidal displacement signal at 73 kHz, and terminated by a rigid support. Output from this model was compared with the solution provided by the equivalent discrete dynamics model, and good agreement obtained. An experimental configuration involving the same chain, but terminated by an annular support made of a liquid photopolymer resin was also simulated and the velocity of the last sphere obtained through simulation was compared with laser vibrometer measurement, with good agreement. This model was then extended whereby the granular chain was coupled to an acoustic medium with the properties of water, via a thin vitreous carbon cylinder. Finite element predictions of the acoustic pressure indicate that, for a 73 kHz excitation frequency, harmonic rich acoustic pulses with harmonic content close to 1 MHz are predicted.

Modelling the effect of acoustic waves on nucleation.

A phase transformation in a metastable phase can be affected when it is subjected to a high intensity ultrasound wave. In this study we determined the effect of oscillation in pressure and temperature on a phase transformation using the Gibbs droplet model in a generic format. The developed model is valid for both equilibrium and non-equilibrium clusters formed through a stationary or non-stationary process. We validated the underlying model by comparing the predicted kinetics of water droplet formation from the gas phase against experimental data in the absence of ultrasound. Our results demonstrated better agreement with experimental data in comparison with classical nucleation theory. Then, we determined the thermodynamics and kinetics of nucleation and the early stage of growth of clusters in an isothermal sonocrystallisation process. This new contribution shows that the effect of pressure on the kinetics of nucleation is cluster size-dependent in contrast to classical nucleation theory.

A theoretical investigation of photoacoustic contrast agents.

Photoacoustic imaging offers significant potential as a biomedical imaging modality. For some applications, however, there is a need for contrast enhancement. In this paper, a theoretical comparison is presented of the efficacy of three different designs for photoacoustic contrast agents (PACAs), specifically, a droplet of dye, a bubble filled with gas coated by a dye loaded shell, and a droplet of volatile dye. For each case, the governing equations describing the dynamics of a single PACA in a homogenous incompressible fluid are derived. The coupled sets of equations describing the bubble oscillation and resulting radiated pressure, the photo-acoustic energy equation, and the equation of state are then solved numerically. The numerical results predict a stronger radiated acoustic signal for the same optical source energy density in the case of the volatile dye droplet by a factor of up to two orders of magnitude compared with the other two types of agent.

A first-order k-space model for elastic wave propagation in heterogeneous media.

A pseudospectral model of linear elastic wave propagation is described based on the first order stress-velocity equations of elastodynamics. k-space adjustments to the spectral gradient calculations are derived from the dyadic Green's function solution to the second-order elastic wave equation and used to (a) ensure the solution is exact for homogeneous wave propagation for timesteps of arbitrarily large size, and (b) also allows larger time steps without loss of accuracy in heterogeneous media. The formulation in k-space allows the wavefield to be split easily into compressional and shear parts. A perfectly matched layer (PML) absorbing boundary condition was developed to effectively impose a radiation condition on the wavefield. The staggered grid, which is essential for accurate simulations, is described, along with other practical details of the implementation. The model is verified through comparison with exact solutions for canonical examples and further examples are given to show the efficiency of the method for practical problems. The efficiency of the model is by virtue of the reduced point-per-wavelength requirement, the use of the fast Fourier transform (FFT) to calculate the gradients in k space, and larger time steps made possible by the k-space adjustments.

A fast full-wave solver for calculating ultrasound propagation in the body.

Therapeutic ultrasound is a promising non-invasive method for inducing various beneficial biological effects in the human body. In cancer treatment applications, high-power ultrasound is focused at a target tissue volume to ablate the malignant tumour. The success of the procedure depends on the ability to accurately focus ultrasound and destroy the target tissue volume through coagulative necrosis whilst preserving the surrounding healthy tissue. Patient-specific treatment planning strategies are therefore being developed to increase the efficacy of such therapies, while reducing any damage to healthy tissue. These strategies require to use high-performance computing methods to solve ultrasound wave propagation in the body quickly and accurately. For realistic clinical scenarios, all numerical methods which employ volumetric meshes require several hours or days to solve the full-wave propagation on a computer cluster. The boundary element method (BEM) is an efficient approach for modelling the wave field because only the boundaries of the hard and soft tissue regions require discretisation. This paper presents a multiple-domain BEM formulation with a novel preconditioner for solving the Helmholtz transmission problem (HTP). This new formulation is efficient at high-frequencies and where high-contrast materials are present. Numerical experiments are performed to solve the HTP in multiple domains comprising: (i) human ribs, an idealised abdominal fat layer and liver tissue, (ii) a human kidney with a perinephric fat layer, exposed to the acoustic field generated by a high-intensity focused ultrasound (HIFU) array transducer. The time required to solve the equations associated with these problems on a single workstation is of the order of minutes. These results demonstrate the great potential of this new BEM formulation for accurately and quickly solving ultrasound wave propagation problems in large anatomical domains which is essential for developing treatment planning strategies.

Preparation of microbubble suspensions by co-axial electrohydrodynamic atomization.

In this paper we report a novel method, based on co-axial electrohydrodynamic jetting, for the preparation of microbubble suspensions containing bubbles <10 microm in size and having a narrow size distribution. No selective filtration is necessary and the suspensions are produced directly by the process. To demonstrate the method, glycerol was used as the liquid medium, flowing in the outer needle of the co-axial twin needle arrangement and undergoing electrohydrodynamic atomization in the stable cone-jet mode while air flowed through the inner needle at the same time. At zero applied voltage a hollow stream of liquid flowed from the outer needle. When the applied voltage was increased, eventually the hollow stream became a stable cone-jet and emitted a microthread of bubbles, which were collected in a container of glycerol to obtain microbubble suspensions. The size of the microbubbles was measured via optical microscopy and laser diffractometry. Several microbubble suspensions were prepared and characterised and the size distribution was found to be critically dependent on the ratio (n) of flow rates of liquid/air and, in particular the flow rate of the air. At n=1.5, with the flow rate of air set at approximately 1.7 microl/s, a microbubble suspension containing bubbles in the size range 2-8 microm was obtained.

Analysis of solitary wave impulses in granular chains using ultrasonic excitation.

The propagation of broad bandwidth solitary wave impulses, generated within granular chains by narrow bandwidth ultrasonic excitation, is studied in detail. Theoretical predictions are compared to experimental results. It is demonstrated that the observed effects result from a sum of a solitary wave traveling out from the source with a wave that reflects from the far end of the chain. It is shown that this combination, when used with an excitation in the form of a long-duration tone burst, encourages the generation of multiple impulses with a characteristic periodicity. This study shows that the properties of the chain structure and the excitation can be adjusted so as to generate ultrasonic solitary wave impulses with a high amplitude and known frequency content, which are of interest in applications such as biomedical ultrasound.

Ultrasonic propagation in finite-length granular chains.

A narrowband ultrasound source has been used to generate solitary wave impulses in finite-length chains of spheres. Once the input signal is of sufficient amplitude, both harmonics and sub-harmonics of the input frequency can be generated as non-linear normal modes of the system, allowing a train of impulses to be established from a sinusoidal input. The characteristics of the response have been studied as a function of the physical properties of the chain, the input waveform and the level of static pre-compression. The results agree with the predictions of a theoretical model, based on a set of discrete dynamic equations for the spheres for finite-length chains. Impulses are only created for very small pre-compression forces of the order of 0.01N, where strongly non-linear behaviour is expected.

An assessment of the DORT method on simple scatterers using boundary element modelling.

The ability to focus through ribs overcomes an important limitation of a high-intensity focused ultrasound (HIFU) system for the treatment of liver tumours. Whilst it is important to generate high enough acoustic pressures at the treatment location for tissue lesioning, it is also paramount to ensure that the resulting ultrasonic dose on the ribs remains below a specified threshold, since ribs both strongly absorb and reflect ultrasound. The DORT (décomposition de l'opérateur de retournement temporel) method has the ability to focus on and through scatterers immersed in an acoustic medium selectively without requiring prior knowledge of their location or geometry. The method requires a multi-element transducer and is implemented via a singular value decomposition of the measured matrix of inter-element transfer functions. The efficacy of a method of focusing through scatterers is often assessed by comparing the specific absorption rate (SAR) at the surface of the scatterer, and at the focal region. The SAR can be obtained from a knowledge of the acoustic pressure magnitude and the acoustic properties of the medium and scatterer. It is well known that measuring acoustic pressures with a calibrated hydrophone at or near a hard surface presents experimental challenges, potentially resulting in increased measurement uncertainties. Hence, the DORT method is usually assessed experimentally by measuring the SAR at locations on the surface of the scatterer after the latter has been removed from the acoustic medium. This is also likely to generate uncertainties in the acoustic pressure measurement. There is therefore a strong case for assessing the efficacy of the DORT method through a validated theoretical model. The boundary element method (BEM) applied to exterior acoustic scattering problems is well-suited for such an assessment. In this study, BEM was used to implement the DORT method theoretically on locally reacting spherical scatterers, and to assess its focusing capability relative to the spherical focusing case, binarised apodisation based on geometric ray tracing and the phase conjugation method.

The potential for thermal damage posed by microbubble ultrasound contrast agents.

The development of coated microbubble ultrasound contrast agents for use in imaging applications and as carriers in drug and gene delivery applications has intensified the need for a clear understanding of their behaviour and potential bioeffects. Previous studies have focused on the risks posed by unencapsulated bubbles as representing the "worst case scenario". They have concluded that the risk of thermal damage should be minimal provided the threshold for inertial cavitation is not exceeded. However, these treatments have ignored the heating effects due to viscous dissipation in the coatings of contrast agent particles. Simulations indicate that the temperature rise due to this process may be sufficient to generate harmful bioeffects even under conventionally "safe" insonation conditions. The implications of these findings and strategies for addressing the risks posed by contrast agents are discussed.

Microbubble ultrasound contrast agents: a review.

The superior scattering properties of gas bubbles compared with blood cells have made microbubble ultrasound contrast agents important tools in ultrasound diagnosis. Over the past 2 years they have become the focus of a wide and rapidly expanding field of research, with their benefits being repeatedly demonstrated, both in ultrasound image enhancement, and more recently in drug and gene delivery applications. However, despite considerable investigation, their behaviour is by no means fully understood and, while no definite evidence of harmful effects has been obtained, there remain some concerns as to their safety. In this review the existing theoretical and experimental evidence is examined in order to clarify the extent to which contrast agents are currently understood and to identify areas for future research. In particular the disparity between the conditions considered in theoretical models and those encountered both in vitro, and more importantly in vivo is discussed, together with the controversy regarding the risk of harmful bio-effects.

Ultrasonic evaluation of ultra-high molecular weight polyethylene used in medical prostheses.

An ultrasonic non-destructive evaluation process is established to observe cracks and delaminations that occur below the surface of retrieved ultra-high molecular weight polyethylene joints. The result from the ultrasonic evaluation is compared with destructive (optical) test results of the actual cracks. Feasibility of characterizing different grades of bulk polyethylene is also established by means of ultrasonic attenuation measurements in the materials. This ultrasonic data can be used to give a better understanding of the failure mechanisms in the UHMWPE material.

A comparison of methods for focusing the field of a HIFU array transducer through human ribs.

A forward model, which predicts the scattering by human ribs of a multi-element high-intensity focused ultrasound transducer, was used to investigate the efficacy of a range of focusing approaches described in the literature. This forward model is based on the boundary element method and was described by Gélat et al (2011 Phys. Med. Biol. 56 5553-81; 2012 Phys. Med. Biol. 57 8471-97). The model has since been improved and features a complex surface impedance condition at the surface of the ribs. The inverse problem of focusing through the ribs was implemented on six transducer array-rib topologies and five methods of focusing were investigated, including spherical focusing, binarized apodization based on geometric ray tracing, phase conjugation and the decomposition of the time-reversal operator method. The excitation frequency was 1 MHz and the array was of spherical-section type. Both human and idealized rib topologies were considered. The merit of each method of focusing was examined. It was concluded that the constrained optimization approach offers greater potential than the other focusing methods in terms of maximizing the ratio of acoustic pressure magnitudes at the focus to those on the surface of the ribs whilst taking full advantage of the dynamic range of the phased array.

The optimization of acoustic fields for ablative therapies of tumours in the upper abdomen.

High intensity focused ultrasound (HIFU) enables highly localized, non-invasive tissue ablation and its efficacy has been demonstrated in the treatment of a range of cancers, including those of the kidney, prostate and breast. HIFU offers the ability to treat deep-seated tumours locally, and potentially bears fewer side effects than more invasive treatment modalities such as resection, chemotherapy and ionizing radiation. There remains however a number of significant challenges which currently hinder its widespread clinical application. One of these challenges is the need to transmit sufficient energy through the ribcage to ablate tissue at the required foci whilst minimizing the formation of side lobes and sparing healthy tissue. Ribs both absorb and reflect ultrasound strongly. This sometimes results in overheating of bone and overlying tissue during treatment, leading to skin burns. Successful treatment of a patient with tumours in the upper abdomen therefore requires a thorough understanding of the way acoustic and thermal energy is deposited. Previously, a boundary element approach based on a Generalized Minimal Residual (GMRES) implementation of the Burton-Miller formulation was developed to predict the field of a multi-element HIFU array scattered by human ribs, the topology of which was obtained from CT scan data (Gélat et al 2011 Phys. Med. Biol. 56 5553-81). The present paper describes the reformulation of the boundary element equations as a least-squares minimization problem with nonlinear constraints. The methodology has subsequently been tested at an excitation frequency of 1 MHz on a spherical multi-element array in the presence of ribs. A single array-rib geometry was investigated on which a 50% reduction in the maximum acoustic pressure magnitude on the surface of the ribs was achieved with only a 4% reduction in the peak focal pressure compared to the spherical focusing case. This method was then compared with a binarized apodization approach based on ray tracing and against the decomposition of the time-reversal operator (DORT). In both cases, the constrained optimization provided a superior ratio of focal peak pressure to maximum pressure magnitude on the surface of the ribs.

Nondestructive evaluation of bone cement and bone cement/metal interface failure.

To quantify the failure mechanisms related to the loosening of cemented hip joint replacements, novel techniques, capable of monitoring, nondestructively, the initiation and progression of failure during in vitro fatigue tests, were employed. Fatigue testing of model cement and cement-stem test pieces was monitored using acoustic emission (AE) sensors. Once damage was detected, an ultrasonic imaging system was used to obtain an image of the damage site and to measure the stiffness of the affected region. This method of examination provided a detailed insight into the internal crack propagation and delamination patterns. Initial work was conducted on bulk cement specimens subjected to bending and tension. The second stage of the work examined a model stem-cement interface under tensile opening loading conditions. A novel ultrasonic technique was used to measure the bond quality at the cement-metal interface. Progressive delamination was identified over time, and the AE technique was able to identify critical areas of delamination before they could be identified conclusively by ultrasonic imaging. The work has demonstrated the potential of the AE technique as a tool for the preclinical assessment of total hip replacements.

Crack initiation processes in acrylic bone cement.

A major constraint in improving the understanding of the micromechanics of the fatigue failure process and, hence, in optimizing bone cement performance is found in the uncertainties associated with monitoring the evolution of the internal defects that are believed to dominate in vivo failure. The present study aimed to synthesize high resolution imaging with complementary damage monitoring/detection techniques. As a result, evidence of the chronology of failure has been obtained. The earliest stages of crack initiation have been captured and it is proposed that, in the presence of a pore, crack initiation may occur away from the pore due to the combined influence of pore morphology and the presence of defects within regions of stress concentration. Furthermore, experimental evidence shows that large agglomerations of BaSO(4) are subject to microcracking during fatigue, although in the majority of cases, these are not the primary cause of failure. It is proposed that cracks may then remain contained within the agglomerations because of the clamping effect of the matrix during volumetric shrinkage upon curing.