Tomographic reconstruction using heuristic Monte Carlo methods

A tomographic reconstruction method based on Monte Carlo random searching guided by the information contained in the projections of radiographed objects is presented. In order to solve the optimization problem, a multiscale algorithm is proposed to reduce computation. The reconstruction is performed in a coarse-to-fine multigrid scale that initializes each resolution level with the reconstruction of the previous coarser level, which substantially improves the performance. The method was applied to a real case reconstructing the internal structure of a small metallic object with internal components, showing excellent results.     

可压稠密气固两相流动过程的模拟计算

基于稠密气体分子运动论和颗粒动理学,建立可压稠密气固两相流动模型。采用梯度模拟来考虑气相可压缩性对气相湍流的影响。模拟计算表明气固两相射流速度沿轴向和径向减小,颗粒浓度下降。气固两相射流具有高的颗粒温度,呈现强烈的气固两相湍流流动特性。     

General expressions of peaked traveling wave solutions of CH-γ and CH equations

We use qualitative analysis and numerical simulation to study peaked traveling wave solutions of CH-γ and CH equations. General expressions of peakon and periodic cusp wave solutions are obtained. Some previous results become our special cases.     

激光诱导间质肿瘤热疗的数值模拟和实验研究

本文在考虑生物组织物性动态变化的情况下建立了激光诱导间质肿瘤热疗(LITT)的物理数学模型,采用MonteCarlo方法数值模拟了LITT中激光能量在生物组织内的传输过程,基于Pennes生物传热方程和Arrhenius方程数值求解了组织内的温度分布和热损伤体积的变化,分析了热物性及血液灌注率的动态变化对LITT过程的影响,并与相应的离体实验结果进行了对比。数值模拟结果表明,组织的热物性及血液灌注率的动态变化对于热损伤体积的变化具有重要的影响。因此在激光诱导间质肿瘤热疗的数值模拟中应该考虑热物性及血液灌注率的动态变化以期为临床治疗方案的制定提供更为准确的依据。     

Anisotropic friction model of DNA electrophoresis in polymer solutions: Comparison with direct observations

The electrophoresis of DNA chains in uncrosslinked polymer solutions with a Brownian dynamics simulation with an anisotropic friction tensor was analyzed. According to the degree of anisotropy, three types of migration behavior are obtained: fluctuation without or with periodicity between U‐shaped and compact conformations, or migration with linear conformation. We found good agreement between our simulation results and the direct observations of DNA by fluorescence microscopy. © 2003 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 41: 1316–1322, 2003     

Transport properties of acetone aqueous solutions: molecular dynamics simulation and NMR studies

Two quantities η_rel and are applied to study the nonideal acetone–water association mixture. An all-atom acetone model and a TIP5P water model have been adopted for molecular dynamics simulation. We study the transport properties of the system comparing the 's of strong hydrogen bond and weak contact based on transport properties, MD simulations together with NMR experimental data and find good agreement of concentration dependence, which exhibits the cooperation effect.     

Methods for implementing the stream boundary condition in DSMC computations

In the direct simulation Monte‐Carlo (DSMC) method for simulating rarefied gas flows, the velocities of simulator particles that cross a simulation boundary and enter the simulation space are typically generated using the acceptance–rejection procedure that samples the velocities from a truncated theoretical velocity distribution that excludes low and high velocities. This paper analyses an alternative technique, where the velocities of entering particles are obtained by extending the simulation procedures to a region adjacent to the simulation space, and considering the movement of particles generated within that region during the simulation time step. The alternative method may be considered as a form of acceptance–rejection procedure, and permits the generation of all possible velocities, although the population of high velocities is depleted with respect to the theoretical distribution. Nevertheless, this is an improvement over the standard acceptance–rejection method. Previous implementations of the alternative method gave a number flux lower than the theoretical number required. Two methods for obtaining the correct number flux are presented. For upstream boundaries in high‐speed flows, the alternative method is more computationally efficient than the acceptance–rejection method. However, for downstream boundaries, the alternative method is extremely inefficient. The alternative method, with the correct theoretical number flux, should therefore be used in DSMC computations in favour of the acceptance–rejection method for upstream boundaries in high‐speed flows. Copyright © 2003 John Wiley & Sons, Ltd.     

A numerical and experimental investigation of the interactions between a non‐uniform planar array of incompressible free jets

The interaction between multiple incompressible air jets has been studied numerically and experimentally. The numerical predictions have been first validated using experimental data for a single jet configuration. The spreading features of five unequal jets in the configuration of one larger central jet surrounded by four smaller equi‐distant jets, have been studied, for different lateral spacing ratios of 1.5, 2.0 and 2.5 and a central jet Reynolds number of 1.24×105 (corresponding to a Mach number of 0.16). Flow of five equal jets has also been simulated, for the sake of comparison. The jet interactions commence at an axial distance of about 3–4 diameters and complete by an axial distance of about 10 diameters for the lowest spacing ratio of 1.5. For larger spacing ratios, the length required for the start and completion of jet interaction increase. Peripheral jets bend more towards the central jet and merge at a smaller distance, when their sizes are smaller than that of the central jet. The entrainment ratio for multiple jets is higher than that for a single jet. Excellent agreement is observed between the experimental data and theoretical predictions for both mean flow field and turbulent quantities, at regions away from the jet inlet. The potential core length and initial jet development, however, are not predicted very accurately due to differences in the assumed and actual velocity profiles at the jet inlet. Copyright © 2004 John Wiley & Sons, Ltd.     

Dependence of the kinetics of the anionic polymerization of methyl methacrylate on the concentration in active centers

Studies on the anionic polymerization of methyl methacrylate in tetrahydrofuran and in the presence of sparteine have revealed a beneficial effect due to this additive, resulting in a decrease in the extent of termination. Better control of the definition of the polymers formed can thus be achieved in the presence of this additive. On the other hand, macromolecular engineering requires a range of active species concentrations lower than 10^?3 mol L^?1 and particularly the synthesis of polymers of high molar masses. For a better understanding of the mechanism of chain growth under such concentration conditions, the kinetics of polymerization have been investigated with a technique based on adiabatic calorimetry. Sparteine has been found to lack sufficient cation‐binding power to prevent the propagating enolate ion pairs from aggregating. The rate constant of propagation of nonaggregated species has been estimated, as well as the aggregation constant of equilibrium. For very low initiator concentrations, termination reactions have been shown to profoundly alter the control of the polymerization and to prevent a quantitative monomer conversion. Theoretical maximal conversions have been calculated from kinetic data and compare well with the experimental values. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 4964–4975, 2004     

Stereological simulations of asymmetric polymerization in channels of bile acid inclusion compounds: A detailed analysis of the monomer arrangements in the channels

We studied simulations by computer graphics to estimate the steric mechanism of the asymmetric polymerization of prochiral diene monomers in channels of inclusion compounds of steroidal bile acids, such as deoxycholic acid (DCA) and cholic acid. We applied a hierarchization method to interpret the crystal structures of bile acids, clarifying that the chiral host molecules associated to form characteristic 2₁-helical assemblies with uneven surfaces. A detailed analysis of the uneven channels in a close-packing state indicated that there were many possible arrangements of the monomers in the channels. The plausible arrangements in the channel could explain a previous study, which showed that the polymerization in the DCA channel yielded chiral polymers with a predominant configuration from prochiral diene monomers, such as 2-methyl-trans-1,3-pentadiene. On the basis of such simulation studies of the arrangements of guest monomers in the channel, we examined a plausible steric mechanism for asymmetric inclusion polymerization. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 4648–4655, 2004