Comparison of SRIM, MCNPX and GEANT simulations with experimental data for thick Al absorbers

Proton computerized tomography deals with relatively thick targets like the human head or trunk. In this case precise analytical calculation of the proton final energy is a rather complicated task, thus the Monte Carlo simulation stands out as a solution. We used the GEANT4.8.2 code to calculate the proton final energy spectra after passing a thick Al absorber and compared it with the same conditions of the experimental data. The ICRU49, Ziegler85 and Ziegler2000 models from the low energy extension pack were used. The results were also compared with the SRIM2008 and MCNPX2.4 simulations, and with solutions of the Boltzmann transport equation in the Fokker-Planck approximation. (C) 2009 Elsevier Ltd. All rights reserved.

Integration Methods Used in Numerical Simulations of Electromagnetic Transients

This paper made an analysis of some numerical integration methods that can be used in electromagnetic transient simulations. Among the existing methods, we analyzed the trapezoidal integration method (or Heun formula), Simpson's Rule and Runge-Kutta. These methods were used in simulations of electromagnetic transients in power systems, resulting from switching operations and maneuvers that occur in transmission lines. Analyzed the characteristics such as accuracy, computation time and robustness of the methods of integration.

A case study of the 9 August 1988 South Atlantic storm: Numerical simulations of the wave activity

During 9-11 August 1988, a cyclone developed over Uruguay in the lee of the Andes Mountains and moved over the South Atlantic Ocean, where it redeveloped into an intense storm. This storm was responsible for unusual wave activity along the Brazilian shoreline from 22° to 32°S. The Brazilian news media reported the loss of at least one life, waves of 3 m and higher, and the disappearance of a drainage pipe, which weighed 8000 kg, off the shores of Rio de Janeiro. In this paper, the evolution of this intense storm and the associated ocean wave response is studied through European Centre for Medium-Range Weather Forecasts analyses, a hydrostatic limited-area meteorological model, and a second-generation prognostic wave model. The atmospheric model results indicated the presence of a long-lived and large fetch with surface wind velocities higher than 12 m s -1 directed toward the coast. Some areas with velocities of 20 m s -1 were embedded in the fetch. The wave model forced by this wind field was able to simulate waves with a significant height of 8 m far from the coast and about 4 m in regions very close to the Brazilian coast in agreement with the occurrence reported at Rio de Janeiro. The swell propagation toward the coast of Rio de Janeiro was obstructed by a northeastward 10-m wind during the first 24-h period of the model's integration. During the second 24-h period, the fetch was still large and strong, but the obstacle was removed by a counterclockwise rotation of wind direction favoring the swell and windsea propagation toward the Rio de Janeiro coast.

Coil to extended coil transition in polygalacturonic acid: Conductometric titration and monte carlo simulations

The conformational transition from coil to extended coil for polygalacturonic acid has been studied by conductometric titrations and Monte Carlo simulations. The results of conductometric titrations at different polymer concentrations have been analyzed using the model proposed by Manning,1 which describes the conductivity of polyelectrolitic solutions. This experimental approach provides the transport factor and the average distance between charged groups at different degrees of ionization (α). The mean distances between charged groups have been compared with the values obtained by Monte Carlo simulations. In these simulations the polymer chain is modeled as a self-avoiding random walk in a cubic lattice. The monomers interact through the unscreened Coulombic potential. The ratio between the end-to-end distance and the number of ionized beads provides the average distance between charged monomers. The experimental and theoretical values are in good agreement for the whole range of ionization degrees accessed by conductometric titrations. These results suggest that the electrostatic interactions seem to be the major contribution for the coil to extended coil conformational change. The small deviations for α ≤ 0.5 suggests that the stiffness of the chain, associated with local interactions, becomes increasingly significant as the fraction of charged groups is decreased. © 2000 American Chemical Society.

Simulations of incompressible fluid flows by a least squares finite element method

In this work simulations of incompressible fluid flows have been done by a Least Squares Finite Element Method (LSFEM) using velocity-pressure-vorticity and velocity-pressure-stress formulations, named u-p-ω) and u-p-τ formulations respectively. These formulations are preferred because the resulting equations are partial differential equations of first order, which is convenient for implementation by LSFEM. The main purposes of this work are the numerical computation of laminar, transitional and turbulent fluid flows through the application of large eddy simulation (LES) methodology using the LSFEM. The Navier-Stokes equations in u-p-ω and u-p-τ formulations are filtered and the eddy viscosity model of Smagorinsky is used for modeling the sub-grid-scale stresses. Some benchmark problems are solved for validate the numerical code and the preliminary results are presented and compared with available results from the literature. Copyright © 2005 by ABCM.

On numerical simulations of a nonlinear self-excited system with two non-ideal sources

In this work, the dynamic behavior of self-synchronization and synchronization through mechanical interactions between the nonlinear self-excited oscillating system and two non-ideal sources are examined by numerical simulations. The physical model of the system vibrating consists of a non-linear spring of Duffing type and a nonlinear damping described by Rayleigh's term. This system is additional forced by two unbalanced identical direct current motors with limited power (non-ideal excitations). The present work mathematically implements the parametric excitation described by two periodically changing stiffness of Mathieu type that are switched on/off. Copyright © 2005 by ASME.

A study of the effect of overshooting deep convection on the water content of the TTL and lower stratosphere from Cloud Resolving Model simulations

Simulations of overshooting, tropical deep convection using a Cloud Resolving Model with bulk microphysics are presented in order to examine the effect on the water content of the TTL (Tropical Tropopause Layer) and lower stratosphere. This case study is a subproject of the HIBISCUS (Impact of tropical convection on the upper troposphere and lower stratosphere at global scale) campaign, which took place in Bauru, Brazil (22° S, 49° W), from the end of January to early March 2004. Comparisons between 2-D and 3-D simulations suggest that the use of 3-D dynamics is vital in order to capture the mixing between the overshoot and the stratospheric air, which caused evaporation of ice and resulted in an overall moistening of the lower stratosphere. In contrast, a dehydrating effect was predicted by the 2-D simulation due to the extra time, allowed by the lack of mixing, for the ice transported to the region to precipitate out of the overshoot air. Three different strengths of convection are simulated in 3-D by applying successively lower heating rates (used to initiate the convection) in the boundary layer. Moistening is produced in all cases, indicating that convective vigour is not a factor in whether moistening or dehydration is produced by clouds that penetrate the tropopause, since the weakest case only just did so. An estimate of the moistening effect of these clouds on an air parcel traversing a convective region is made based on the domain mean simulated moistening and the frequency of convective events observed by the IPMet (Instituto de Pesquisas Meteorológicas, Universidade Estadual Paulista) radar (S-band type at 2.8 Ghz) to have the same 10 dBZ echo top height as those simulated. These suggest a fairly significant mean moistening of 0.26, 0.13 and 0.05 ppmv in the strongest, medium and weakest cases, respectively, for heights between 16 and 17 km. Since the cold point and WMO (World Meteorological Organization) tropopause in this region lies at ∼ 15.9 km, this is likely to represent direct stratospheric moistening. Much more moistening is predicted for the 15-16 km height range with increases of 0.85-2.8 ppmv predicted. However, it would be required that this air is lofted through the tropopause via the Brewer Dobson circulation in order for it to have a stratospheric effect. Whether this is likely is uncertain and, in addition, the dehydration of air as it passes through the cold trap and the number of times that trajectories sample convective regions needs to be taken into account to gauge the overall stratospheric effect. Nevertheless, the results suggest a potentially significant role for convection in determining the stratospheric water content. Sensitivity tests exploring the impact of increased aerosol numbers in the boundary layer suggest that a corresponding rise in cloud droplet numbers at cloud base would increase the number concentrations of the ice crystals transported to the TTL, which had the effect of reducing the fall speeds of the ice and causing a ∼13% rise in the mean vapour increase in both the 15-16 and 16-17 km height ranges, respectively, when compared to the control case. Increases in the total water were much larger, being 34% and 132% higher for the same height ranges, but it is unclear whether the extra ice will be able to evaporate before precipitating from the region. These results suggest a possible impact of natural and anthropogenic aerosols on how convective clouds affect stratospheric moisture levels.

Comparison of GEANT4 simulations with experimental data for thick al absorbers

Proton beams in medical applications deal with relatively thick targets like the human head or trunk. Therefore, relatively small differences in the total proton stopping power given, for example, by the different models provided by GEANT4 can lead to significant disagreements in the final proton energy spectra when integrated along lengthy proton trajectories. This work presents proton energy spectra obtained by GEANT4.8.2 simulations using ICRU49, Ziegler1985 and Ziegler2000 models for 19.68MeV protons passing through a number of Al absorbers with various thicknesses. The spectra were compared with the experimental data, with TRIM/SRIM2008 and MCNPX2.4.0 simulations, and with the Payne analytical solution for the transport equation in the Fokker-Plank approximation. It is shown that the MCNPX simulations reasonably reproduce well all experimental spectra. For the relatively thin targets all the methods give practically identical results but this is not the same for the thick absorbers. It should be noted that all the spectra were measured at the proton energies significantly above 2MeV, i.e., in the so-called Bethe-Bloch region. Therefore the observed disagreements in GEANT4 results, simulated with different models, are somewhat unexpected. Further studies are necessary for better understanding and definitive conclusions. © 2009 American Institute of Physics.

Determination of minimum fluidization velocity for gas-solid beds by experimental data and numerical simulations

The purpose of this work is to predict the minimum fluidization velocity Umf in a gas-solid fluidized bed. The study was carried out with an experimental apparatus for sand particles with diameters between 310μm and 590μm, and density of 2,590kg/m3. The experimental results were compared with numerical simulations developed in MFIX (Multiphase Flow with Interphase eXchange) open source code [1], for three different sizes of particles: 310mum, 450μm and 590μm. A homogeneous mixture with the three kinds of particles was also studied. The influence of the particle diameter was presented and discussed. The Ergun equation was also used to describe the minimum fluidization velocity. The experimental data presented a good agreement with Ergun equation and numerical simulations. Copyright © 2011 by ASME.

The fractional-nonlinear robotic manipulator: Modeling and dynamic simulations

In this paper, we applied the Riemann-Liouville approach and the fractional Euler-Lagrange equations in order to obtain the fractional-order nonlinear dynamics equations of a two link robotic manipulator. The aformentioned equations have been simulated for several cases involving: integer and non-integer order analysis, with and without external forcing acting and some different initial conditions. The fractional nonlinear governing equations of motion are coupled and the time evolution of the angular positions and the phase diagrams have been plotted to visualize the effect of fractional order approach. The new contribution of this work arises from the fact that the dynamics equations of a two link robotic manipulator have been modeled with the fractional Euler-Lagrange dynamics approach. The results reveal that the fractional-nonlinear robotic manipulator can exhibit different and curious behavior from those obtained with the standard dynamical system and can be useful for a better understanding and control of such nonlinear systems. © 2012 American Institute of Physics.

Molecular Dynamics simulations of track formation at different ensembles

A series of Molecular Dynamics simulations of thermal spikes has been run in zircon. For two different ensembles: microcanonical one and a combination of microcanonical one acting on the simulation core with Langevin one on the side walls of simulation. Depending on the used ensemble, different track-formation threshold energies were found. When the combined ensemble is carried out, the total energy of the simulations varies with the temperature which can influence how annealing fission-track models should deal with the lattice recovery. A fission-track annealing model is tested with the simulation results. © 2012 Elsevier Ltd. All rights reserved.

On the unzipping mechanisms of carbon nanotubes: Insights from reactive molecular dynamics simulations

Unzipping carbon nanotubes (CNTs) is considered one of the most promising approaches for the controlled and large-scale production of graphene nanoribbons (GNR). These structures are considered of great importance for the development of nanoelectronics because of its dimensions and intrinsic nonzero band gap value. Despite many years of investigations some details on the dynamics of the CNT fracture/unzipping processes remain unclear. In this work we have investigated some of these process through molecular dynamics simulations using reactive force fields (ReaxFF), as implemented in the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code. We considered multi-walled CNTs of different dimensions and chiralities and under induced mechanical stretching. Our preliminary results show that the unzipping mechanisms are highly dependent on CNT chirality. Well-defined and distinct fracture patterns were observed for the different chiralities. Armchair CNTs favor the creation of GNRs with well-defined armchair edges, while zigzag and chiral ones produce GNRs with less defined and defective edges. © 2012 Materials Research Society.

Mathematical analysis and simulations involving chemotherapy and surgery on large human tumours under a suitable cell-kill functional response

Dosage and frequency of treatment schedules are important for successful chemotherapy. However, in this work we argue that cell-kill response and tumoral growth should not be seen as separate and therefore are essential in a mathematical cancer model. This paper presents a mathematical model for sequencing of cancer chemotherapy and surgery. Our purpose is to investigate treatments for large human tumours considering a suitable cell-kill dynamics. We use some biological and pharmacological data in a numerical approach, where drug administration occurs in cycles (periodic infusion) and surgery is performed instantaneously. Moreover, we also present an analysis of stability for a chemotherapeutic model with continuous drug administration. According to Norton & Simon [22], our results indicate that chemotherapy is less eficient in treating tumours that have reached a plateau level of growing and that a combination with surgical treatment can provide better outcomes.

Assembling a xylanase-lichenase chimera through all-atom molecular dynamics simulations

Multifunctional enzyme engineering can improve enzyme cocktails for emerging biofuel technology. Molecular dynamics through structure-based models (SB) is an effective tool for assessing the tridimensional arrangement of chimeric enzymes as well as for inferring the functional practicability before experimental validation. This study describes the computational design of a bifunctional xylanase-lichenase chimera (XylLich) using the xynA and bglS genes from Bacillus subtilis. In silico analysis of the average solvent accessible surface area (SAS) and the root mean square fluctuation (RMSF) predicted a fully functional chimera, with minor fluctuations and variations along the polypeptide chains. Afterwards, the chimeric enzyme was built by fusing the xynA and bglS genes. XylLich was evaluated through small-angle X-ray scattering (SAXS) experiments, resulting in scattering curves with a very accurate fit to the theoretical protein model. The chimera preserved the biochemical characteristics of the parental enzymes, with the exception of a slight variation in the temperature of operation and the catalytic efficiency (k cat/Km). The absence of substantial shifts in the catalytic mode of operation was also verified. Furthermore, the production of chimeric enzymes could be more profitable than producing a single enzyme separately, based on comparing the recombinant protein production yield and the hydrolytic activity achieved for XylLich with that of the parental enzymes. © 2013 Elsevier B.V. All rights reserved.

Ab initio protein folding simulations using atomic burials as informational intermediates between sequence and structure

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)