ترغب بنشر مسار تعليمي؟ اضغط هنا

Simulation of reaction-diffusion processes in three dimensions using CUDA

109   0   0.0 ( 0 )
 نشر من قبل Istvan Lagzi
 تاريخ النشر 2010
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

Numerical solution of reaction-diffusion equations in three dimensions is one of the most challenging applied mathematical problems. Since these simulations are very time consuming, any ideas and strategies aiming at the reduction of CPU time are important topics of research. A general and robust idea is the parallelization of source codes/programs. Recently, the technological development of graphics hardware created a possibility to use desktop video cards to solve numerically intensive problems. We present a powerful parallel computing framework to solve reaction-diffusion equations numerically using the Graphics Processing Units (GPUs) with CUDA. Four different reaction-diffusion problems, (i) diffusion of chemically inert compound, (ii) Turing pattern formation, (iii) phase separation in the wake of a moving diffusion front and (iv) air pollution dispersion were solved, and additionally both the Shared method and the Moving Tiles method were tested. Our results show that parallel implementation achieves typical acceleration values in the order of 5-40 times compared to CPU using a single-threaded implementation on a 2.8 GHz desktop computer.


قيم البحث

اقرأ أيضاً

Snowflake growth provides us with a fascinating example of spontaneous pattern formation in nature. Attempts to understand this phenomenon have led to important insights in non-equilibrium dynamics observed in various active scientific fields, rangin g from pattern formation in physical and chemical systems, to self-assembly problems in biology. Yet, very few models currently succeed in reproducing the diversity of snowflake forms in three dimensions, and the link between model parameters and thermodynamic quantities is not established. Here, we report a modified phase field model that describes the subtlety of the ice vapour phase transition, through anisotropic water molecules attachment and condensation, surface diffusion, and strong anisotropic surface tension, that guarantee the anisotropy, faceting and dendritic growth of snowflakes. We demonstrate that this model reproduces the growth dynamics of the most challenging morphologies of snowflakes from the Nakaya diagram. We find that the growth dynamics of snow crystals matches the selection theory, consistently with previous experimental observations.
We study the decay process for the reaction-diffusion process of three species on the small-world network. The decay process is manipulated from the deterministic rate equation of three species in the reaction-diffusion system. The particle density a nd the global reaction rate on a two dimensional small-world network adding new random links is discussed numerically, and the global reaction rate before and after the crossover is also found by means of the Monte Carlo simulation. The time-dependent global reaction rate scales as a power law with the scaling exponent 0.66 at early time regime while it scales with -0.50 at long time regime, in all four cases of the added probability $p=0.2-0.8$. Especially, our result presented is compared with the numerical calculation of regular networks.
CELES is a freely available MATLAB toolbox to simulate light scattering by many spherical particles. Aiming at high computational performance, CELES leverages block-diagonal preconditioning, a lookup-table approach to evaluate costly functions and ma ssively parallel execution on NVIDIA graphics processing units using the CUDA computing platform. The combination of these techniques allows to efficiently address large electrodynamic problems ($>10^4$ scatterers) on inexpensive consumer hardware. In this paper, we validate near- and far-field distributions against the well-established multi-sphere $T$-matrix (MSTM) code and discuss the convergence behavior for ensembles of different sizes, including an exemplary system comprising $10^5$ particles.
We formulate, solve computationally and study experimentally the problem of collecting solar energy in three dimensions(1-5). We demonstrate that absorbers and reflectors can be combined in the absence of sun tracking to build three-dimensional photo voltaic (3DPV) structures that can generate measured energy densities (energy per base area, kWh/m2) higher by a factor of 2-20 than stationary flat PV panels, versus an increase by a factor of 1.3-1.8 achieved with a flat panel using dual-axis sun tracking(6). The increased energy density is countered by a higher solar cell area per generated energy for 3DPV compared to flat panel design (by a factor of 1.5-4 in our conditions), but accompanied by a vast range of improvements. 3DPV structures are steadier sources of solar energy generation at all latitudes: they can double the number of peak power generation hours and dramatically reduce the seasonal, latitude and weather variations of solar energy generation compared to a flat panel design. Self-supporting 3D shapes can create new schemes for PV installation and the increased energy density can facilitate the use of cheaper thin film materials in area-limited applications. Our findings suggest that harnessing solar energy in three dimensions can open new avenues towards Terawatt-scale generation.
The Graphics Processing Unit (GPU) is a powerful tool for parallel computing. In the past years the performance and capabilities of GPUs have increased, and the Compute Unified Device Architecture (CUDA) - a parallel computing architecture - has been developed by NVIDIA to utilize this performance in general purpose computations. Here we show for the first time a possible application of GPU for environmental studies serving as a basement for decision making strategies. A stochastic Lagrangian particle model has been developed on CUDA to estimate the transport and the transformation of the radionuclides from a single point source during an accidental release. Our results show that parallel implementation achieves typical acceleration values in the order of 80-120 times compared to CPU using a single-threaded implementation on a 2.33 GHz desktop computer. Only very small differences have been found between the results obtained from GPU and CPU simulations, which are comparable with the effect of stochastic transport phenomena in atmosphere. The relatively high speedup with no additional costs to maintain this parallel architecture could result in a wide usage of GPU for diversified environmental applications in the near future.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا