Two novel numerical methods for solving the two-dimensional fractional Fitzhugh-Nagumo-monodomain model

Fractional mathematical models represent a new approach to modelling complex spatial problems in which there is heterogeneity at many spatial and temporal scales. In this paper, a two-dimensional fractional Fitzhugh-Nagumo-monodomain model with zero Dirichlet boundary conditions is considered. The model consists of a coupled space fractional diffusion equation (SFDE) and an ordinary differential equation. For the SFDE, we first consider the numerical solution of the Riesz fractional nonlinear reaction-diffusion model and compare it to the solution of a fractional in space nonlinear reaction-diffusion model. We present two novel numerical methods for the two-dimensional fractional Fitzhugh-Nagumo-monodomain model using the shifted Grunwald-Letnikov method and the matrix transform method, respectively. Finally, some numerical examples are given to exhibit the consistency of our computational solution methodologies. The numerical results demonstrate the effectiveness of the methods.

Computational analysis of laminated glass panels under blast loads: A comparison of two dimensional and three dimensional modelling approaches

This paper illustrates the use of finite element (FE) technique to investigate the behaviour of laminated glass (LG) panels under blast loads. Two and three dimensional (2D and 3D) modelling approaches available in LS-DYNA FE code to model LG panels are presented. Results from the FE analysis for mid-span deflection and principal stresses compared well with those from large deflection plate theory. The FE models are further validated using the results from a free field blast test on a LG panel. It is evident that both 2D and 3D LG models predict the experimental results with reasonable accuracy. The 3D LG models give slightly more accurate results but require considerably more computational time compared to the 2D LG models.

Two-dimensional graphene heterojunctions: The tunable mechanical properties

We report the mechanical properties of different two-dimensional carbon heterojunctions (HJs) made from graphene and various stable graphene allotropes, including α-, β-, γ- and 6612-graphyne (GY), and graphdiyne (GDY). It is found that all HJs exhibit a brittle behaviour except the one with α-GY, which however shows a hardening process due to the formation of triple carbon rings. Such hardening process has greatly deferred the failure of the structure. The yielding of the HJs is usually initiated at the interface between graphene and graphene allotropes, and monoatomic carbon rings are normally formed after yielding. By varying the locations of graphene (either in the middle or at the two ends of the HJs), similar mechanical properties have been obtained, suggesting insignificant impacts from location of graphene allotropes. Whereas, changing the types and percentages of the graphene allotropes, the HJs exhibit vastly different mechanical properties. In general, with the increasing graphene percentage, the yield strain decreases and the effective Young’s modulus increases. Meanwhile, the yield stress appears irrelevant with the graphene percentage. This study provides a fundamental understanding of the tensile properties of the heterojunctions that are crucial for the design and engineering of their mechanical properties, in order to facilitate their emerging future applications in nanoscale devices, such as flexible/stretchable electronics.

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.

Resource allocation in two species systems: Is it worth acknowledging species interactions?

It is becoming increasingly popular to consider species interactions when managing ecological foodwebs. Such an approach is useful in determining how management can affect multiple species, with either beneficial or detrimental consequences. Identifying such actions is particularly valuable in the context of conservation decision making as funding is severely limited. This paper outlines a new approach that simplifies the resource allocation problem in a two species system for a range of species interactions: independent, mutualism, predator-prey, and competitive exclusion. We assume that both species are endangered and we do not account for decisions over time. We find that optimal funding allocation is to the conservation of the species with the highest marginal gain in expected probability of survival and that, across all except mutualist interaction types, optimal conservation funding allocation differs between species. Loss in efficiency from ignoring species interactions was most severe in predator-prey systems. The funding problem we address, where an ecosystem includes multiple threatened species, will only become more commonplace as increasing numbers of species worldwide become threatened. © 2011 Elsevier B.V.

Two-dimensional hydrogen-bonded polymers in the crystal structures of the ammonium salts of phenoxyacetic acid, (4-fluorophenoxy)acetic acid and (4-chloro-2-methylphenoxy)acetic acid

The structures of the ammonium salts of phenoxyacetic acid, NH4+ C8H6O3- (I), (4-fluorophenoxy)acetic acid NH4+ C8H5FO3- (II) and the herbicidally active (4-chloro-2-methylphenoxy)acetic acid (MCPA), NH4+ C9H8ClO3-. 0.5(H2O) (III) have been determined. All have two-dimensional layered structures based on inter-species ammonium N-H...O hydrogen-bonding associations which give core substructures consisting primarily of conjoined cyclic motifs. Crystals of (I) and (II) are isomorphous with the core comprising R2/1(5), R2/1(4) and centrosymmetric R2/4(8) ring motifs, giving two-dimensional layers lying parallel to (100). In (III), the water molecule of solvation lies on a crystallographic twofold rotation axis and bridges two carboxyl O-atoms in an R4/4(12) hydrogen-bonded motif, creating two R3/4(10) rings which together with a conjoined centrosymmetric R2/4(8) ring incorporating both ammonium cations, generate two-dimensional layers lying parallel to (100). No pi-pi ring associations are present in any of the structures.