Calculation of volume from crank angle using reconfigurable hardware

Fast calculation of quantities such as in-cylinder volume and indicated power is important in internal combustion engine research. Multiple channels of data including crank angle and pressure were collected for this purpose using a fully instrumented diesel engine research facility. Currently, existing methods use software to post-process the data, first calculating volume from crank angle, then calculating the indicated work and indicated power from the area enclosed by the pressure-volume indicator diagram. Instead, this work investigates the feasibility of achieving real-time calculation of volume and power via hardware implementation on Field Programmable Gate Arrays (FPGAs). Alternative hardware implementations were investigated using lookup tables, Taylor series methods or the CORDIC (CoOrdinate Rotation DIgital Computer) algorithm to compute the trigonometric operations in the crank angle to volume calculation, and the CORDIC algorithm was found to use the least amount of resources. Simulation of the hardware based implementation showed that the error in the volume and indicated power is less than 0.1%.

Calculation of engine parameters using reconfigurable hardware

The feasibility of real-time calculation of parameters for an internal combustion engine via reconfigurable hardware implementation is investigated as an alternative to software computation. A detailed in-hardware field programmable gate array (FPGA)-based design is developed and evaluated using input crank angle and in-cylinder pressure data from fully instrumented diesel engines in the QUT Biofuel Engine Research Facility (BERF). Results indicate the feasibility of employing a hardware-based implementation for real-time processing for speeds comparable to the data sampling rate currently used in the facility, with acceptably low level of discrepancies between hardware and software-based calculation of key engine parameters.

Hydrogen spillover mechanism on a pd-doped mg surface as revealed by ab initio density functional calculation

The hydrogenation kinetics of Mg is slow, impeding its application for mobile hydrogen storage. We demonstrate by ab initio density functional theory (DFT) calculations that the reaction path can be greatly modified by adding transition metal catalysts. Contrasting with Ti doping, a Pd dopant will result in a very small activation barrier for both dissociation of molecular hydrogen and diffusion of atomic H on the Mg surface. This new computational finding supports for the first time by ab initio simulationthe proposed hydrogen spillover mechanism for rationalizing experimentally observed fast hydrogenation kinetics for Pd-capped Mg materials.

Managing memory and reducing I/O cost for correlation matrix calculation in bioinformatics

The generation of a correlation matrix from a large set of long gene sequences is a common requirement in many bioinformatics problems such as phylogenetic analysis. The generation is not only computationally intensive but also requires significant memory resources as, typically, few gene sequences can be simultaneously stored in primary memory. The standard practice in such computation is to use frequent input/output (I/O) operations. Therefore, minimizing the number of these operations will yield much faster run-times. This paper develops an approach for the faster and scalable computing of large-size correlation matrices through the full use of available memory and a reduced number of I/O operations. The approach is scalable in the sense that the same algorithms can be executed on different computing platforms with different amounts of memory and can be applied to different problems with different correlation matrix sizes. The significant performance improvement of the approach over the existing approaches is demonstrated through benchmark examples.

Creative city index 2012 - Data sources, calculation methods and tools [Report spreadsheet]

Spreadsheet for Creative City Index 2012

Calculation of dose kernels for radiotherapy applications using the GEANT 4 Monte Carlo toolkit

Dose kernels may be used to calculate dose distributions in radiotherapy (as described by Ahnesjo et al., 1999). Their calculation requires use of Monte Carlo methods, usually by forcing interactions to occur at a point. The Geant4 Monte Carlo toolkit provides a capability to force interactions to occur in a particular volume. We have modified this capability and created a Geant4 application to calculate dose kernels in cartesian, cylindrical, and spherical scoring systems. The simulation considers monoenergetic photons incident at the origin of a 3 m x 3 x 9 3 m water volume. Photons interact via compton, photo-electric, pair production, and rayleigh scattering. By default, Geant4 models photon interactions by sampling a physical interaction length (PIL) for each process. The process returning the smallest PIL is then considered to occur. In order to force the interaction to occur within a given length, L_FIL, we scale each PIL according to the formula: PIL_forced = L_FIL 9 (1 - exp(-PIL/PILo)) where PILo is a constant. This ensures that the process occurs within L_FIL, whilst correctly modelling the relative probability of each process. Dose kernels were produced for an incident photon energy of 0.1, 1.0, and 10.0 MeV. In order to benchmark the code, dose kernels were also calculated using the EGSnrc Edknrc user code. Identical scoring systems were used; namely, the collapsed cone approach of the Edknrc code. Relative dose difference images were then produced. Preliminary results demonstrate the ability of the Geant4 application to reproduce the shape of the dose kernels; median relative dose differences of 12.6, 5.75, and 12.6 % were found for an incident photon energy of 0.1, 1.0, and 10.0 MeV respectively.

The catalytic role of an isolated-Ti atom in the hydrogenation of Ti-doped Al(001) surface : An ab initio density functional theory calculation

Recent work [S. Chaudhuri, J.T. Muckerman, J. Phys. Chem. B 109 (2005) 6952] reported that two Ti-substituted atoms on an Al(0 0 1) surface can form a catalytically active site for the dissociation of H2, but the diffusion barrier of atomic H away from Ti site is as high as 1.57 eV. By using ab initio density functional calculations, we found that two hydrogen molecules can dissociate on isolated-Ti atom doped Al(0 0 1) surface with small activation barriers (0.21 and 0.235 eV for first and second H2, respectively). Additionally, the diffusion barrier of atomic H away from Ti site is also moderate (0.47 eV). These results contribute further towards understanding the improved kinetics observed in recycling of hydrogen with Ti-doped NaAlH4.

Effects of small field output factors on IMRT optimisation and dose calculation

Introduction Given the known challenges of obtaining accurate measurements of small radiation fields, and the increasing use of small field segments in IMRT beams, this study examined the possible effects of referencing inaccurate field output factors in the planning of IMRT treatments. Methods This study used the Brainlab iPlan treatment planning system to devise IMRT treatment plans for delivery using the Brainlab m3 microMLC (Brainlab, Feldkirchen, Germany). Four pairs of sample IMRT treatments were planned using volumes, beams and prescriptions that were based on a set of test plans described in AAPM TG 119’s recommendations for the commissioning of IMRT treatment planning systems [1]: • C1, a set of three 4 cm volumes with different prescription doses, was modified to reduce the size of the PTV to 2 cm across and to include an OAR dose constraint for one of the other volumes. • C2, a prostate treatment, was planned as described by the TG 119 report [1]. • C3, a head-and-neck treatment with a PTV larger than 10 cm across, was excluded from the study. • C4, an 8 cm long C-shaped PTV surrounding a cylindrical OAR, was planned as described in the TG 119 report [1] and then replanned with the length of the PTV reduced to 4 cm. Both plans in each pair used the same beam angles, collimator angles, dose reference points, prescriptions and constraints. However, one of each pair of plans had its beam modulation optimisation and dose calculation completed with reference to existing iPlan beam data and the other had its beam modulation optimisation and dose calculation completed with reference to revised beam data. The beam data revisions consisted of increasing the field output factor for a 0.6 9 0.6 cm2 field by 17 % and increasing the field output factor for a 1.2 9 1.2 cm2 field by 3 %. Results The use of different beam data resulted in different optimisation results with different microMLC apertures and segment weightings between the two plans for each treatment, which led to large differences (up to 30 % with an average of 5 %) between reference point doses in each pair of plans. These point dose differences are more indicative of the modulation of the plans than of any clinically relevant changes to the overall PTV or OAR doses. By contrast, the maximum, minimum and mean doses to the PTVs and OARs were smaller (less than 1 %, for all beams in three out of four pairs of treatment plans) but are more clinically important. Of the four test cases, only the shortened (4 cm) version of TG 119’s C4 plan showed substantial differences between the overall doses calculated in the volumes of interest using the different sets of beam data and thereby suggested that treatment doses could be affected by changes to small field output factors. An analysis of the complexity of this pair of plans, using Crowe et al.’s TADA code [2], indicated that iPlan’s optimiser had produced IMRT segments comprised of larger numbers of small microMLC leaf separations than in the other three test cases. Conclusion: The use of altered small field output factors can result in substantially altered doses when large numbers of small leaf apertures are used to modulate the beams, even when treating relatively large volumes.

Optimizing I/O cost and managing memory for composition vector method based on correlation matrix calculation in bioinformatics

The generation of a correlation matrix for set of genomic sequences is a common requirement in many bioinformatics problems such as phylogenetic analysis. Each sequence may be millions of bases long and there may be thousands of such sequences which we wish to compare, so not all sequences may fit into main memory at the same time. Each sequence needs to be compared with every other sequence, so we will generally need to page some sequences in and out more than once. In order to minimize execution time we need to minimize this I/O. This paper develops an approach for faster and scalable computing of large-size correlation matrices through the maximal exploitation of available memory and reducing the number of I/O operations. The approach is scalable in the sense that the same algorithms can be executed on different computing platforms with different amounts of memory and can be applied to different bioinformatics problems with different correlation matrix sizes. The significant performance improvement of the approach over previous work is demonstrated through benchmark examples.

Calculation of the oxygen isotope discrimination factor for studying plant respiration

Measurement of discrimination against 18O during dark respiration in plants is currently accepted as the only reliable method of estimating the partitioning of electrons between the cytochrome and alternative pathways. In this paper, we review the theory of the technique and its application to a gas-phase system. We extend it to include sampling effects and show that the isotope discrimination factor, D, is calculated as –dln(1 + δ)/dlnO*, where δ is isotopic composition of the substrate oxygen and O*=[O2]/[N2] in a closed chamber containing tissue respiring in the dark. It is not necessary to integrate the expression but, if the integrated form is used, the resultant regression should not be constrained through the origin. This is important since any error in D will have significant effects on the estimation of the flux of electrons through the two pathways.

Hydrolysis of genotoxic methyl-substituted oxiranes : Experimental kinetic and semiempirical studies

The kinetics of acid-catalyzed hydrolysis of seven methylated aliphatic epoxides - R1R2C(O)CR3R4 (A: R1=R2=R3=R4=H; B: R1=R2=R3=H, R4=Me; C: R1=R2=H, R3=R4=Me; D: R1=R3=H, R2=R4=Me(trans); E: R1=R3=H, R2=R4=Me(cis); F: R1=R3=R4=Me, R2=H; G: R1=R2=R3=R4=Me) - has been studied at 36 ± 1.5°C. Compounds with two methyl groups at the same carbon atom of the oxirane ring exhibit highest rate constants (k(eff) in reciprocal molar concentration per second: 11.0 ± 1.3 for C, 10.7 ± 2.1 for F, and 8.7 ± 0.7 for G as opposed to 0.124 ± 0.003 for B, 0.305 ± 0.003 for D, and 0.635 ± 0.036 for E). Ethylene oxide (A) displays the lowest rate of hydrolysis (0.027 M-1 s-1). The results are consistent with literature data available for compounds A, B, and C. To model the reactivities we have employed quantum chemical calculations (MNDO, AM1, PM3, and MINDO/3) of the main reaction species. There is a correlation of the logarithm k(eff) with the total energy of epoxide ring opening. The best correlation coefficients (r) were obtained using the AM1 and MNDO methods (0.966 and 0.957, respectively). However, unlike MNDO, AM1 predicts approximately zero energy barriers for the oxirane ring opening of compounds B, C, E and G, which is not consistent with published kinetic data. Thus, the MNDO method provides a preferential means of modeling the acidic hydrolysis of the series of methylated oxiranes. The general ranking of mutagenicity in vitro, A > B > C, is in line with the concept that this sequence also gradually leaves the expoxide reactivity optimal for genotoxicity toward reactivities leading to higher biological detoxifications.

Further development of discrete computational techniques for calculation of restricted diffusion propagators in porous media

Magnetic resonance is a well-established tool for structural characterisation of porous media. Features of pore-space morphology can be inferred from NMR diffusion-diffraction plots or the time-dependence of the apparent diffusion coefficient. Diffusion NMR signal attenuation can be computed from the restricted diffusion propagator, which describes the distribution of diffusing particles for a given starting position and diffusion time. We present two techniques for efficient evaluation of restricted diffusion propagators for use in NMR porous-media characterisation. The first is the Lattice Path Count (LPC). Its physical essence is that the restricted diffusion propagator connecting points A and B in time t is proportional to the number of distinct length-t paths from A to B. By using a discrete lattice, the number of such paths can be counted exactly. The second technique is the Markov transition matrix (MTM). The matrix represents the probabilities of jumps between every pair of lattice nodes within a single timestep. The propagator for an arbitrary diffusion time can be calculated as the appropriate matrix power. For periodic geometries, the transition matrix needs to be defined only for a single unit cell. This makes MTM ideally suited for periodic systems. Both LPC and MTM are closely related to existing computational techniques: LPC, to combinatorial techniques; and MTM, to the Fokker-Planck master equation. The relationship between LPC, MTM and other computational techniques is briefly discussed in the paper. Both LPC and MTM perform favourably compared to Monte Carlo sampling, yielding highly accurate and almost noiseless restricted diffusion propagators. Initial tests indicate that their computational performance is comparable to that of finite element methods. Both LPC and MTM can be applied to complicated pore-space geometries with no analytic solution. We discuss the new methods in the context of diffusion propagator calculation in porous materials and model biological tissues.

Calculation of the speed of walking of individuals with transfemoral amputation fitted with bone-anchored prosthesis

Typically, the walking ability of individuals with a transfemoral amputation (TFA) can be represented by the speed of walking (SofW) obtained in experimental settings. Recent developments in portable kinetic systems allow assessing the level of activity of TFA during actual daily living outside the confined space of a gait lab. Unfortunately, only minimal spatio-temporal characteristics could be extracted from the kinetic data including the cadence and the duration on gait cycles. Therefore, there is a need for a way to use some of these characteristics to assess the instantaneous speed of walking during daily living. The purpose of the study was to compare several methods to determine SofW using minimal spatial gait characteristics.

trans-Dichloro-bis-(arylazoimidazole)palladium(II): Synthesis, Structure, Photoisomerization, and DFT Calculation

Reaction between PdCl2 and 1-alkyl-2-(arylazo)imidazole (RaaiR') or 1-alkyl-2-(naphthyl-alpha/beta-azo)imidazole (alpha/beta-NaiR') under reflux in ethanol has isolated complexes of compositions Pd(RaaiR')(2)Cl-2 (5, 6) and Pd(alpha/beta-NaiR')(2)Cl-2 (7, 8). The X-ray structure determination of one of the molecules, Pd(alpha-NaiBz)(2)Cl-2 (7c), has reported a trans-PdCl2 configuration, and alpha-NaiBz acts as monodentate N(imidazole) donor ligand. The spectral (IR, UV-vis, H-1 NMR) data support the structure. UV light irradiation (light source: Perkin-Elmer LS 55 spectrofluorimeter, Xenon discharge lamp, lambda = 360-396 nm) in a MeCN solution of the complexes shows E-to-Z isomerization of the coordinated azoimidazole unit. The reverse transformation, Z-to-E, is very slow with visible light irradiation. Quantum yields (phi(E-Z)) of E-to-Z isomerization are calculated, and phi is lower than that of the free ligand but comparable with those of Cd(II) and Hg(II) complexes of the same ligand. The Z-to-E isomerization is a thermally induced process. The activation energy (E-a) of Z-to-E isomerization is calculated by controlled-temperature experimentation. cis-Pd(azoimidazole)Cl-2 complexes (azomidazole acts as N(imidazole) and N(azo) Chelating ligand) do not respond upon light irradiation, which supports the idea that the presence of noncoordinated azo-N to make free azo (-N=N-) function is important to reveal photochromic activity. DFT calculation of Pd(alpha-NaiBz)(2)Cl-2 (7c) has suggested that the HOMO of the molecule is constituted of Pd (32%) and Cl (66%), and hence photo excitation may use the energy of Pd and Cl instead of that of the photofunctional -N=N-Ar motif; thus, the rate of photoisomerization and quantum yield decrease versus the free ligand values.

Statistical power calculation and sample size determination for environmental studies with data below detection limits

Power calculation and sample size determination are critical in designing environmental monitoring programs. The traditional approach based on comparing the mean values may become statistically inappropriate and even invalid when substantial proportions of the response values are below the detection limits or censored because strong distributional assumptions have to be made on the censored observations when implementing the traditional procedures. In this paper, we propose a quantile methodology that is robust to outliers and can also handle data with a substantial proportion of below-detection-limit observations without the need of imputing the censored values. As a demonstration, we applied the methods to a nutrient monitoring project, which is a part of the Perth Long-Term Ocean Outlet Monitoring Program. In this example, the sample size required by our quantile methodology is, in fact, smaller than that by the traditional t-test, illustrating the merit of our method.