A finite volume scheme with preconditioned Lanczos method for two-dimensional space-fractional reaction–diffusion equations

Fractional differential equations have been increasingly used as a powerful tool to model the non-locality and spatial heterogeneity inherent in many real-world problems. However, a constant challenge faced by researchers in this area is the high computational expense of obtaining numerical solutions of these fractional models, owing to the non-local nature of fractional derivatives. In this paper, we introduce a finite volume scheme with preconditioned Lanczos method as an attractive and high-efficiency approach for solving two-dimensional space-fractional reaction–diffusion equations. The computational heart of this approach is the efficient computation of a matrix-function-vector product f(A)bf(A)b, where A A is the matrix representation of the Laplacian obtained from the finite volume method and is non-symmetric. A key aspect of our proposed approach is that the popular Lanczos method for symmetric matrices is applied to this non-symmetric problem, after a suitable transformation. Furthermore, the convergence of the Lanczos method is greatly improved by incorporating a preconditioner. Our approach is show-cased by solving the fractional Fisher equation including a validation of the solution and an analysis of the behaviour of the model.

Diagnostics and two-dimensional simulation of low-frequency inductively coupled plasmas with neutral gas heating and electron heat fluxes

This article presents the results on the diagnostics and numerical modeling of low-frequency (∼460 KHz) inductively coupled plasmas generated in a cylindrical metal chamber by an external flat spiral coil. Experimental data on the electron number densities and temperatures, electron energy distribution functions, and optical emission intensities of the abundant plasma species in low/intermediate pressure argon discharges are included. The spatial profiles of the plasma density, electron temperature, and excited argon species are computed, for different rf powers and working gas pressures, using the two-dimensional fluid approach. The model allows one to achieve a reasonable agreement between the computed and experimental data. The effect of the neutral gas temperature on the plasma parameters is also investigated. It is shown that neutral gas heating (at rf powers≥0.55kW) is one of the key factors that control the electron number density and temperature. The dependence of the average rf power loss, per electron-ion pair created, on the working gas pressure shows that the electron heat flux to the walls appears to be a critical factor in the total power loss in the discharge.

Reactive species in Ar+H2 plasma-aided nanofabrication : two-dimensional discharge modelling

Understanding the generation of reactive species in a plasma is an important step towards creating reliable and robust plasma-aided nanofabrication processes. A two-dimensional fluid simulation of the number densities of surface preparation species in a low-temperature, low-pressure, non-equilibrium Ar+H2 plasma is conducted. The operating pressure and H2 partial pressure have been varied between 70-200 mTorr and 0.1-50%, respectively. An emphasis is placed on the application of these results to nanofabrication. A reasonable balance between operating pressures and H 2 partial pressures that would optimize the number densities of the two working units largely responsible for activation and passivation of surface dangling bonds (Ar+ and H respectively) in order to achieve acceptable rates of surface activation and passivation is obtained. It is found that higher operating pressures (150-200 mTorr) and lower H2 partial pressures (∼5%) are required in order to ensure high number densities of Ar+ and H species. This paper contributes to the improvement of the controllability and predictability of plasma-based nanoassembly processes.

Two-dimensional simulation of nanoassembly precursor species in Ar+H2+C2H2 reactive plasmas

The results of two-dimensional fluid simulation of number densities and fluxes of the main building blocks and surface preparation species involved in nanoassembly of carbon-based nanopatterns in Ar+H2+C2H2 reactive plasmas are reported. It is shown that the process parameters and non-uniformity of surface fluxes of each particular species may affect the targeted nanopattern quality. The results can be used to improve predictability of plasma-aided nanofabrication processes and optimize the parameters of plasma nanotools.KGaA, Weinheim.

Two-dimensional simulation of nanoparticle deposition from high-density plasmas on microstructured surfaces

Selective and controlled deposition of plasma-grown nanoparticles is one of the pressing problems of plasma-aided nanofabrication. The results of advanced numerical simulations of motion of charge-variable nanoparticles in the plasma presheath and sheath areas and in localized microscopic electric fields created by surface microstructures are reported. Conditions for site-selective deposition of such nanoparticles onto individual microstructures and open surface areas within a periodic micropattern are formulated. The effects of plasma parameters, surface potential, and micropattern features on nanoparticle deposition are investigated and explained using particle charging and plasma force arguments. The results are generic and applicable to a broad range of nanoparticle-generating plasmas and practical problems ranging from management of nanoparticle contamination in microelectronics to site-selective nanoparticle deposition into specified device locations, and synthesis of advanced microporous materials and nanoparticle superlattices. © 2007 American Institute of Physics.

Thermal conductivity of configurable two-dimensional carbon nanotube architecture and strain modulation

We reported the thermal conductivity of the two-dimensional carbon nanotube (CNT)-based architecture, which can be constructed through welding of single-wall CNTs by electron beam. Using large-scale nonequilibrium molecular dynamics simulations, the thermal conductivity is found to vary with different junction types due to their different phonon scatterings at the junction. The strong length and strain dependence of the thermal conductivity suggests an effective avenue to tune the thermal transport properties of the CNT-based architecture, benefiting the design of nanoscale thermal rectifiers or phonon engineering.

A Crank--Nicolson ADI spectral method for a two-dimensional riesz space fractional nonlinear reaction-diffusion equation

In this paper, a new alternating direction implicit Galerkin--Legendre spectral method for the two-dimensional Riesz space fractional nonlinear reaction-diffusion equation is developed. The temporal component is discretized by the Crank--Nicolson method. The detailed implementation of the method is presented. The stability and convergence analysis is strictly proven, which shows that the derived method is stable and convergent of order $2$ in time. An optimal error estimate in space is also obtained by introducing a new orthogonal projector. The present method is extended to solve the fractional FitzHugh--Nagumo model. Numerical results are provided to verify the theoretical analysis.

Two-dimensional bow and stern flows in a finite-depth fluid with constant vorticity

This thesis is concerned with two-dimensional free surface flows past semi-infinite surface-piercing bodies in a fluid of finite-depth. Throughout the study, it is assumed that the fluid in question is incompressible, and that the effects of viscosity and surface tension are negligible. The problems considered are physically important, since they can be used to model the flow of water near the bow or stern of a wide, blunt ship. Alternatively, the solutions can be interpreted as describing the flow into, or out of, a horizontal slot. In the past, all research conducted on this topic has been dedicated to the situation where the flow is irrotational. The results from such studies are extended here, by allowing the fluid to have constant vorticity throughout the flow domain. In addition, new results for irrotational flow are also presented. When studying the flow of a fluid past a surface-piercing body, it is important to stipulate in advance the nature of the free surface as it intersects the body. Three different possibilities are considered in this thesis. In the first of these possibilities, it is assumed that the free surface rises up and meets the body at a stagnation point. For this configuration, the nonlinear problem is solved numerically with the use of a boundary integral method in the physical plane. Here the semi-infinite body is assumed to be rectangular in shape, with a rounded corner. Supercritical solutions which satisfy the radiation condition are found for various values of the Froude number and the dimensionless vorticity. Subcritical solutions are also found; however these solutions violate the radiation condition and are characterised by a train of waves upstream. In the limit that the height of the body above the horizontal bottom vanishes, the flow approaches that due to a submerged line sink in a $90^\circ$ corner. This limiting problem is also examined as a special case. The second configuration considered in this thesis involves the free surface attaching smoothly to the front face of the rectangular shaped body. For this configuration, nonlinear solutions are computed using a similar numerical scheme to that used in the stagnant attachment case. It is found that these solution exist for all supercritical Froude numbers. The related problem of the cusp-like flow due to a submerged sink in a corner is also considered. Finally, the flow of a fluid emerging from beneath a semi-infinite flat plate is examined. Here the free surface is assumed to detach from the trailing edge of the plate horizontally. A linear problem is formulated under the assumption that the elevation of the plate is close to the undisturbed free surface level. This problem is solved exactly using the Wiener-Hopf technique, and subcritical solutions are found which are characterised by a train of sinusoidal waves in the far field. The nonlinear problem is also considered. Exact relations between certain parameters for supercritical flow are derived using conservation of mass and momentum arguments, and these are confirmed numerically. Nonlinear subcritical solutions are computed, and the results are compared to those predicted by the linear theory.

Numerical treatment of a two-dimensional variable-order fractional nonlinear reaction-diffusion model

A two-dimensional variable-order fractional nonlinear reaction-diffusion model is considered. A second-order spatial accurate semi-implicit alternating direction method for a two-dimensional variable-order fractional nonlinear reaction-diffusion model is proposed. Stability and convergence of the semi-implicit alternating direct method are established. Finally, some numerical examples are given to support our theoretical analysis. These numerical techniques can be used to simulate a two-dimensional variable order fractional FitzHugh-Nagumo model in a rectangular domain. This type of model can be used to describe how electrical currents flow through the heart, controlling its contractions, and are used to ascertain the effects of certain drugs designed to treat arrhythmia.

The role of tunnel junction resistances and defects on electron transport mechanism in networks of two-dimensional disordered conductors

The effect of tunnel junction resistances on the electronic property and the magneto-resistance of few-layer graphene sheet networks is investigated. By decreasing the tunnel junction resistances, transition from strong localization to weak localization occurs and magneto-resistance changes from positive to negative. It is shown that the positive magneto-resistance is due to Zeeman splitting of the electronic states at the Fermi level as it changes with the bias voltage. As the tunnel junction resistances decrease, the network resistance is well described by 2D weak localization model. Sensitivity of the magneto-resistance to the bias voltage becomes negligible and diminishes with increasing temperature. It is shown 2D weak localization effect mainly occurs inside of the few-layer graphene sheets and the minimum temperature of 5 K in our experiments is not sufficiently low to allow us to observe 2D weak localization effect of the networks as it occurs in 2D disordered metal films. Furthermore, defects inside the few-layer graphene sheets have negligible effect on the resistance of the networks which have small tunnel junction resistances between few-layer graphene sheets

Two-dimensional coordination polymeric structures in caesium complexes with ring-substituted phenoxyacetic acids

The two-dimensional polymeric structures of the caesium complexes with the phenoxyacetic acid analogues (4-fluorophenoxy)acetic acid, (3-chloro-2-methylphenoxy)acetic acid and the herbicidally active (2,4-dichlorophen­oxy)acetic acid (2,4-D), namely poly[[5-(4-fluorophenoxy)acetato][4-(4-fluorophenoxy)acetato]dicaesium], [Cs2(C8H6FO3)2]n, (I), poly[aqua[5-(3-chloro-2-methylphenoxy)acetato]caesium], [Cs(C9H8ClO3)(H2O)]n, (II), and poly[[7-(2,4-di­chlorophenoxy)acetato][(2,4-dichlorphenoxy)acetic acid]caesium], [Cs(C8H5Cl2O3)(C8H6Cl2O3)]n, (III), are described. In (I), the Cs+ cations of the two individual irregular coordination polyhedra in the asymmetric unit (one CsO7 and the other CsO8) are linked by bridging carboxylate O-atom donors from the two ligand molecules, both of which are involved in bidentate chelate Ocarboxy,Ophenoxy interactions, while only one has a bidentate carboxylate O,O'-chelate inter­action. Polymeric extension is achieved through a number of carboxylate O-atom bridges, with a minimum CsCs separation of 4.3231 (9) Å, giving layers which lie parallel to (001). In hydrated complex (II), the irregular nine-coordination about the Cs+ cation comprises a single monodentate water molecule, a bidentate Ocarboxy,Ophenoxy chelate interaction and six bridging carboxylate O-atom bonding interactions, giving a CsCs separation of 4.2473 (3) Å. The water mol­ecule forms intralayer hydrogen bonds within the two-dimensional layers, which lie parallel to (100). In complex (III), the irregular centrosymmetric CsO6Cl2 coordination environment comprises two O-atom donors and two ring-substituted Cl-atom donors from two hydrogen bis[(2,4-dichlorophenoxy)acetate] ligand species in a bidentate chelate mode, and four O-atom donors from bridging carboxyl groups. The duplex ligand species lie across crystallographic inversion centres, linked through a short O-HO hydrogen bond involving the single acid H atom. Structure extension gives layers which lie parallel to (001). The present set of structures of Cs salts of phenoxyacetic acids show previously demonstrated trends among the alkali metal salts of simple benzoic acids with no stereochemically favourable interactive substituent groups for formation of two-dimensional coordination polymers.