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.

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.

Three-dimensional off-design numerical analysis of an organic Rankine cycle radial-inflow turbine

Optimisation of organic Rankine cycles(ORCs for binary cycle applications could play a major role in determining the competitiveness of low to moderate renewable sources. An important aspect of the optimisation is to maximise the turbine output power for a given resource. This requires careful attention to the turbine design notably through numerical simulations. Challenges in the numerical modelling of radial-inflow turbines using high-density working fluids still need to be addressed in order to improve the turbine design and better optimise ORCs. Thispaper presents preliminary 3D numerical simulations of a high-density radial-inflow ORC turbine in sensible geothermal conditions. Following extensive investigation of the operating conditions and thermodynamic cycle analysis, therefrigerant R143a is chosen as the high-density working fluid. The 1D design of the candidate radial-inflow turbine is presented in details. Furthermore, commercially-available software Ansys-CFX is used to perform preliminary steady-state 3D CFD simulations of the candidate R143a radial-inflow turbine for a number of operating conditions including off-design conditions. The real-gas properties are obtained using the Peng–Robinson equations of state.The thermodynamic ORC cycle is presented. The preliminary design created using dedicated radial-inflow turbine software Concepts-Rital is discussed and the 3D CFD results are presented and compared against the meanline analysis.

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 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.

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