We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure coreplatform@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
Computer-aided design (CAD) plays an essential role in creative idea generation on 2D screens during the design process. In most CAD scenarios, virtual object translation is an essential operation, and it is commonly used when designers simulate their innovative solutions. The degrees of freedom (DoF) of virtual object translation modes have been found to directly impact users’ task performance and psychological aspects in simulated environments. Little is known in the existing literature about the sense of agency (SoA), which is a critical psychological aspect emphasizing the feeling of control, in translation modes on 2D screens during the design process. Hence, this study aims to assess users’ SoA in virtual object translation modes on mouse-based, touch-based, and handheld augmented reality (AR) interfaces through subjective and objective measures, such as self-report, task performance, and electroencephalogram (EEG) data. Based on our findings in this study, users perceived a greater feeling of control in 1DoF translation mode, which may help them come up with more creative ideas, than in 3DoF translation mode in the design process; additionally, the handheld AR interface offers less control feel, which may have a negative impact on design quality and creativity, as compared with mouse- and touch-based interfaces. This research contributes to the current literature by analyzing the association between virtual object translation modes and SoA, as well as the relationship between different 2D interfaces and SoA in CAD. As a result of these findings, we propose several design considerations for virtual object translation on 2D screens, which may enable designers to perceive a desirable feeling of control during the design process.
Complex fault patterns associated with rift development in the Hailar Basin were largely influenced by the Mongolia–Okhotsk Ocean and Palaeo-Pacific tectonic regimes during the Late Jurassic to Early Cretaceous periods. Based on 3D seismic data from the Beier Depression in the Hailar Basin, we characterized the reactivation history of multi-trend major faults and examined the evolution of the Beier Depression during the Early Cretaceous period. NE–SW-, NW–SE- and ENE–WSW-oriented major faults originated from strike-slip-associated structures that were pre-existing fabrics and then were reactivated and propagated upward under extensional regimes in the Late Jurassic. During the syn-rift stage (K1t–K1n), the Hailar Basin was in a NNW–SSE- to NW–SE-oriented extensional setting, and major faults of all orientations were active. There was tectonic quiescence (K1n1L) between the syn-rift stages (a rifting transition stage). The short compression stage after the syn-rift stage caused regional compressional deformation. During the post-rift stage (K1d–K1y), the extension direction rotated to an E–W orientation, and a new population of N–S-trending faults formed together with the reactivation of NE–SW- and ENE–WSW-trending major faults. Structural analysis shows that the major ENE–WSW-trending major faults were polycyclic growth faults reactivated via an upward propagation mode and that the NE–SW-trending faults were dip linkage faults reactivated via a dip linkage mode. The reactivation intensity of the NE–SW-trending major faults was stronger than that of the ENE–WSW-trending major faults. These results demonstrate the differences in the evolution of the different trending faults in the same tectonic regime, and the complexity of the final fault patterns in the Beier Depression was produced by differences in the reactivation of major faults. The originate interpretation of the multi-trend major faults in the Hailar Basin provides new insights into fault generation, and the classification of fault growth also has useful implications for future research on multiphase rifts.
Dynamics of two-dimensional flow past a rigid flat plate with a trailing closed flexible filament acting as a deformable afterbody are investigated numerically by an immersed boundary-lattice Boltzmann method for the fluid flow and a finite element method for the filament motion. The effects of Reynolds number ($Re$) and length ratio ($Lr$) on the flow patterns and dynamics of the rigid-flexible coupling system are studied. Based on our numerical results, five typical state modes have been identified in $Lr\unicode{x2013}Re$ plane in terms of the filament shape and corresponding dynamics, i.e. static deformation, micro-vibration, multi-frequency flapping, periodic flapping and chaotic flapping modes, respectively. Benefiting from the passive flow control by using the flexible filament as a deformable afterbody, the coupled system may enjoy a significant drag reduction (up to $22\,\%$) compared with bare plate scenarios ($Lr=1$). Maximum drag reduction achieved at $L_{c,{min}} \in [1.8, 2]$ is often accompanied by the onset of the system state transition. The flow characteristic and its relation to the change in hydrodynamic drag are further explored in order to reveal the underlying mechanisms of the counterintuitive dynamical behaviour of the coupled system. The scaling laws for the form drag and the friction drag, which arise from the pressure and viscous effects, respectively, are proposed to estimate the overall drag acting on the system. The results obtained in the present study may shed some light on understanding the dynamical behaviour of rigid-flexible coupling systems.
The rapid and accurate taxonomic identification of fossils is of great significance in paleontology, biostratigraphy, and other fields. However, taxonomic identification is often labor-intensive and tedious, and the requisition of extensive prior knowledge about a taxonomic group also requires long-term training. Moreover, identification results are often inconsistent across researchers and communities. Accordingly, in this study, we used deep learning to support taxonomic identification. We used web crawlers to collect the Fossil Image Dataset (FID) via the Internet, obtaining 415,339 images belonging to 50 fossil clades. Then we trained three powerful convolutional neural networks on a high-performance workstation. The Inception-ResNet-v2 architecture achieved an average accuracy of 0.90 in the test dataset when transfer learning was applied. The clades of microfossils and vertebrate fossils exhibited the highest identification accuracies of 0.95 and 0.90, respectively. In contrast, clades of sponges, bryozoans, and trace fossils with various morphologies or with few samples in the dataset exhibited a performance below 0.80. Visual explanation methods further highlighted the discrepancies among different fossil clades and suggested similarities between the identifications made by machine classifiers and taxonomists. Collecting large paleontological datasets from various sources, such as the literature, digitization of dark data, citizen-science data, and public data from the Internet may further enhance deep learning methods and their adoption. Such developments will also possibly lead to image-based systematic taxonomy to be replaced by machine-aided classification in the future. Pioneering studies can include microfossils and some invertebrate fossils. To contribute to this development, we deployed our model on a server for public access at www.ai-fossil.com.
The propulsion of a pitching flexible plate in a uniform flow is investigated numerically. The effects of bending stiffness ($K$), pitching amplitude ($A_L$) and frequency ($St$) on the wake patterns, thrust generations and propulsive performances of the fluid–plate system are analysed. Four typical wake patterns, i.e. von Kármán, reversed von Kármán, deflected and chaotic wakes, emerge from various kinematics, and the $St-A_L$ wake maps are given for various $K$. The drag-to-thrust transitions (DTT) and the wake transitions (WT) between the von Kármán and reversed von Kármán wakes are examined. Results indicate that the WT and DTT boundaries can be scaled by the chord-averaged distance of travel, $\mathcal {L}$, which leads to $\mathcal {L}\times St \approx 1$ and $\mathcal {L}\times St \approx 1.2$, respectively. Further, the resonance mechanism for the performance enhancement is revealed and confirmed in a wide range of parameters. The dimensionless average speed of plate, $\mathcal {U^*}\left (=\mathcal {L}\times St\right )$, is adopted merely to characterize the propulsive performances. For the first time, the $\mathcal {U^*}$-based scaling laws for the thrust and power are revealed in pitching rigid and flexible plates for various $A_L$ and $St$. This study may deepen our understanding of biological swimming and flying, and provide a guide for bionic design.
The association between blood transfusion and ventilator-associated events (VAEs) has not been fully understood. We sought to determine whether blood transfusion increases the risk of a VAE.
Design:
Nested case-control study.
Setting:
This study was based on a registry of healthcare-associated infections in intensive care units at West China Hospital system.
Patients:
1,657 VAE cases and 3,293 matched controls were identified.
Methods:
For each case, 2 controls were randomly selected using incidence density sampling. We defined blood transfusion as a time-dependent variable, and we used weighted Cox models to calculate hazard ratios (HRs) for all 3 tiers of VAEs.
Results:
Blood transfusion was associated with increased risk of ventilator-associated complication-plus (VAC-plus; HR, 1.47; 95% CI, 1.22–1.77; P <.001), VAC-only (HR, 1.29; 95% CI, 1.01–1.65; P = .038), infection-related VAC-plus (IVAC-plus; HR, 1.78; 95% CI, 1.33–2.39; P < .001), and possible ventilator-associated pneumonia (PVAP; HR, 2.10; 95% CI, 1.10–3.99; P = .024). Red blood cell (RBC) transfusion was also associated with increased risk of VAC-plus (HR, 1.34; 95% CI, 1.08–1.65; P = .007), IVAC-plus (HR, 1.70; 95% CI, 1.22–2.36; P = .002), and PVAP (HR, 2.49; 95% CI, 1.17–5.28; P = .018). Compared to patients without transfusion, the risk of VAE was significantly higher in patients with RBC transfusions of >3 units (HR, 1.73; 95% CI, 1.25–2.40; P = .001) but not in those with RBC transfusions of 0–3 units.
Conclusion:
Blood transfusions were associated with increased risk of all tiers of VAE. The risk was significantly higher among patients who were transfused with >3 units of RBCs.
To investigate the clinical impact of ventilator-associated events (VAEs) on adverse prognoses and risk factors for mortality among intensive care unit (ICU) patients receiving invasive mechanical ventilation (IMV) based on an ICU healthcare-associated infection (ICU-HAI) registry.
Design:
A cohort study was conducted based on an ICU-HAI registry including 30,830 patients between 2015 and 2018.
Setting:
The study was conducted using data from 5 adult ICUs of a referral hospital.
Patients:
Adult patients in the ICU-HAI registry who received ≥4 consecutive IMV days.
Methods:
Clinical outcomes and mortality risk factors for VAEs were analyzed using propensity score matching (PSM), multivariate regression models, and sensitivity analyses.
Results:
Of 6,426 included patients, 1,803 developed 1,899 VAEs. After PSM, patients with VAEs did have prolonged length of stay in the ICU and in the hospital, increased hospitalization costs, longer days on mechanical ventilation, higher proportion of ≥9 days on mechanical ventilation, higher rate of failure in extubating mechanical ventilation, and excess all-cause mortality in the ICU. Older age (adjusted OR [aOR], 1.02), higher APACHE II score on ICU admission (aOR, 1.06), pneumonia (aOR, 1.49), blood transfusion (aOR 1.43), immunosuppressive drugs (aOR, 1.69), central-line catheter (aOR, 2.06), and ≥2 VAEs in the ICU (aOR, 1.99) were associated with higher risks for all-cause mortality in an ICU.
Conclusions:
Patients with VAEs indeed had poorer clinical outcomes. Older age, higher APACHE II score on ICU admission, pneumonia, blood transfusion, immunosuppressive drugs, central-line catheter, and ≥2 VAEs in the ICU were risk factors for all-cause mortality of VAE patients in the ICU.
A fever clinic within a hospital plays a vital role in pandemic control because it serves as an outpost for pandemic discovery, monitoring and handling. As the outbreak of coronavirus disease 2019 (COVID-19) in Wuhan was gradually brought under control, the fever clinic in the West Campus of Wuhan Union Hospital introduced a new model for construction and management of temporary mobile isolation wards. A traditional battlefield hospital model was combined with pandemic control regulations, to build a complex of mobile isolation wards that used adaptive design and construction for medical operational, medical waste management and water drainage systems. The mobile isolation wards allowed for the sharing of medical resources with the fever clinic. This increased the capacity and efficiency of receiving, screening, triaging and isolation and observation of patients with fever. The innovative mobile isolation wards also controlled new sudden outbreaks of COVID-19. We document the adaptive design and construction model of the novel complex of mobile isolation wards and explain its characteristics, functions and use.
The convergent Richtmyer–Meshkov instability (RMI) of an $\textrm {SF}_6$ layer with a uniform outer surface and a sinusoidal inner surface surrounded by air (generated by a novel soap film technique) is studied in a semiannular convergent shock tube using high-speed schlieren photography. The outer interface initially suffers only a slight deformation over a long period of time, but distorts quickly at late stages when the inner interface is close to it and produces strong coupling effects. The development of the inner interface can be divided into three stages. At stage I, the interface amplitude first drops suddenly to a lower value due to shock compression, then decreases gradually to zero (phase reversal) and later increases sustainedly in the negative direction. After the reshock (stage II), the perturbation amplitude exhibits long-term quasi-linear growth with time. The quasi-linear growth rate depends weakly on the pre-reshock amplitude and wavelength, but strongly on the pre-reshock growth rate. An empirical model for the growth of convergent RMI under reshock is proposed, which reasonably predicts the present results and those in the literature. At stage III, the perturbation growth is promoted by the Rayleigh–Taylor instability caused by a rarefaction wave reflected from the outer interface. It is found that both the Rayleigh–Taylor effect and the interface coupling depend heavily on the layer thickness. Therefore, controlling the layer thickness is an effective way to modulate the late-stage instability growth, which may be useful for the target design.
The evolution of an $\text{SF}_{6}$ layer surrounded by air is experimentally studied in a semi-annular convergent shock tube by high-speed schlieren photography. The gas layer with a sinusoidal outer interface and a circular inner interface is realized by the soap-film technique such that the initial condition is well controlled. Results show that the thicker the gas layer, the weaker the interface–coupling effect and the slower the evolution of the outer interface. Induced by the distorted transmitted shock and the interface coupling, the inner interface exhibits a slow perturbation growth which can be largely suppressed by reducing the layer thickness. After the reshock, the inner perturbation increases linearly at a growth rate independent of the initial layer thickness as well as of the outer perturbation amplitude and wavelength, and the growth rate can be well predicted by the model of Mikaelian (Physica D, vol. 36, 1989, pp. 343–357) with an empirical coefficient of 0.31. After the linear stage, the growth rate decreases continuously, and finally the perturbation freezes at a constant amplitude caused by the successive stagnation of spikes and bubbles. The convergent geometry constraint as well as the very weak compressibility at late stages are responsible for this instability freeze-out.
The onset of thermal convection in a rapidly rotating spherical shell is studied by linear stability analysis based on the fully compressible Navier–Stokes equations. Compressibility is quantified by the number of density scale heights $N_{\unicode[STIX]{x1D70C}}$, which measures the intensity of density stratification of the motionless, polytropic base state. The nearly adiabatic flow with polytropic index $n=1.499<n_{a}=1.5$ is considered, where $n_{a}$ is the adiabatic polytropic index. By investigating the stability of the base state with respect to the disturbance of specified wavenumber, the instability process is found to be sensitive to the Prandtl number $Pr$ and to $N_{\unicode[STIX]{x1D70C}}$. For large $Pr$ and small $N_{\unicode[STIX]{x1D70C}}$, the quasi-geostrophic columnar mode loses stability first; while for relatively small $Pr$ a new quasi-geostrophic compressible mode is identified, which becomes unstable first under strong density stratification. The inertial mode can also occur first for relatively small $Pr$ and a certain intensity of density stratification in the parameter range considered. Although the Rayleigh numbers $Ra$ for the onsets of the quasi-geostrophic compressible mode and columnar mode are different by several orders of magnitude, we find that they follow very similar scaling laws with the Taylor number. The critical $Ra$ for convection onset is found to be always positive, in contrast with previous results based on the widely used anelastic model that convection can occur at negative $Ra$. By evaluating the relative magnitude of the time derivative of density perturbation in the continuity equation, we show that the anelastic approximation in the present system cannot be applied in the small-$Ra$ and large-$N_{\unicode[STIX]{x1D70C}}$ regime.
Reynolds-averaged Navier–Stokes (RANS) simulations with turbulence closure models continue to play important roles in industrial flow simulations. However, the commonly used linear eddy-viscosity models are intrinsically unable to handle flows with non-equilibrium turbulence (e.g. flows with massive separation). Reynolds stress models, on the other hand, are plagued by their lack of robustness. Recent studies in plane channel flows found that even substituting Reynolds stresses with errors below 0.5 % from direct numerical simulation databases into RANS equations leads to velocities with large errors (up to 35 %). While such an observation may have only marginal relevance to traditional Reynolds stress models, it is disturbing for the recently emerging data-driven models that treat the Reynolds stress as an explicit source term in the RANS equations, as it suggests that the RANS equations with such models can be ill-conditioned. So far, a rigorous analysis of the condition of such models is still lacking. As such, in this work we propose a metric based on local condition number function for a priori evaluation of the conditioning of the RANS equations. We further show that the ill-conditioning cannot be explained by the global matrix condition number of the discretized RANS equations. Comprehensive numerical tests are performed on turbulent channel flows at various Reynolds numbers and additionally on two complex flows, i.e. flow over periodic hills, and flow in a square duct. Results suggest that the proposed metric can adequately explain observations in previous studies, i.e. deteriorated model conditioning with increasing Reynolds number and better conditioning of the implicit treatment of the Reynolds stress compared to the explicit treatment. This metric can play critical roles in the future development of data-driven turbulence models by enforcing the conditioning as a requirement on these models.
Direct numerical simulation is conducted to uncover the response of a supersonic turbulent boundary layer to streamwise concave curvature and the related physical mechanisms at a Mach number of 2.95. Streamwise variations of mean flow properties, turbulence statistics and turbulent structures are analysed. A method to define the boundary layer thickness based on the principal strain rate is proposed, which is applicable for boundary layers subjected to wall-normal pressure and velocity gradients. While the wall friction grows with the wall turning, the friction velocity decreases. A logarithmic region with constant slope exists in the concave boundary layer. However, with smaller slope, it is located lower than that of the flat boundary layer. Streamwise varying trends of the velocity and the principal strain rate within different wall-normal regions are different. The turbulence level is promoted by the concave curvature. Due to the increased turbulence generation in the outer layer, secondary bumps are noted in the profiles of streamwise and spanwise turbulence intensity. Peak positions in profiles of wall-normal turbulence intensity and Reynolds shear stress are pushed outward because of the same reason. Attributed to the Görtler instability, the streamwise extended vortices within the hairpin packets are intensified and more vortices are generated. Through accumulations of these vortices with a similar sense of rotation, large-scale streamwise roll cells are formed. Originated from the very large-scale motions and by promoting the ejection, sweep and spanwise events, the formation of large-scale streamwise roll cells is the physical cause of the alterations of the mean properties and turbulence statistics. The roll cells further give rise to the vortex generation. The large number of hairpin vortices formed in the near-wall region lead to the improved wall-normal correlation of turbulence in the concave boundary layer.
The multipath effect and Non-Line-Of-Sight (NLOS) reception of Global Positioning System (GPS) signals both serve to degrade performance, particularly in urban areas. Although receiver design continues to evolve, residual multipath errors and NLOS signals remain a challenge in built-up areas. It is therefore desirable to identify direct, multipath-affected and NLOS GPS measurements in order improve ranging-based position solutions. The traditional signal strength-based methods to achieve this, however, use a single variable (for example, Signal to Noise Ratio (C/N0)) as the classifier. As this single variable does not completely represent the multipath and NLOS characteristics of the signals, the traditional methods are not robust in the classification of signals received. This paper uses a set of variables derived from the raw GPS measurements together with an algorithm based on an Adaptive Neuro Fuzzy Inference System (ANFIS) to classify direct, multipath-affected and NLOS measurements from GPS. Results from real data show that the proposed method could achieve rates of correct classification of 100%, 91% and 84%, respectively, for LOS, Multipath and NLOS based on a static test with special conditions. These results are superior to the other three state-of-the-art signal reception classification methods.
Heavy metal contamination in the paddy soils of China is a serious concern because of its health risk through transfer in food chains. A field experiment was conducted in 2014–2015 to investigate the long-term effects of different biochar amendments on cadmium (Cd) and arsenic (As) immobilisation in a contaminated paddy field in southern China. Two types of biochar, a rice-straw-derived biochar (RB) and a coconut-by-product-derived biochar (CB), were amended separately to determine their impacts on rice yield and their efficacy in reducing Cd and As in rice. The two-year field experiment showed that biochar amendments significantly improved the rice yields and that CB is superior to RB, especially in the first growth season. Using a large amount of biochar amendment (22.5tha–1) significantly increased soil pH and total organic carbon, and concomitantly decreased the Cd content in rice grains over the four growth seasons, regardless of biochar type and application rate. Arsenic levels in rice were similar to the control, and results from this study suggest that there was a sustainable effect of biochar on Cd sequestration in soil and reduction of Cd accumulation in rice for at least two years. Biochar amendment in soil could be considered as a sustainable, reliable and cost-effective option to remediate heavy metal contamination in paddy fields for long periods.
Sterol regulatory element binding protein 1 (SREBP1) has a central regulatory effect on milk fat synthesis. Lipopolysaccharides (LPS) can induce mastitis and cause milk fat depression in cows. SREBP1 is also known to be associated with inflammatory regulation. Thus, in the current study, we hypothesized that LPS-induced milk fat depression in dairy cow mammary epithelial cells (DCMECs) operates via decreased SREBP1 expression and activity. To examine the hypothesis, DCMECs were isolated and purified from dairy cow mammary tissue and treated with LPS (10 µg/ml). LPS treatment of DCMECs suppressed lipid-metabolism-related transcription factor SREBP1 mRNA expression, nuclear translocation and protein expression, leading to reduced triglyceride content. The transcription levels of acetyl-CoA carboxylase-1 and fatty acid synthetase were significantly down-regulated in DCMECs after LPS treatment, suggesting that acetyl-CoA carboxylase-1 and fatty acid synthetase involved in de novo milk fat synthesis was regulated by SREBP1. In summary, these results suggest that LPS induces milk fat depression in dairy cow mammary epithelial cells via decreased expression of SREBP1 in a time-dependent manner.
The influences of non-Oberbeck–Boussinesq (NOB) effects on flow instabilities and bifurcation characteristics of Rayleigh–Bénard convection are examined. The working fluid is air with reference Prandtl number $Pr=0.71$ and contained in two-dimensional rigid cavities of finite aspect ratios. The fluid flow is governed by the low-Mach-number equations, accounting for the NOB effects due to large temperature difference involving flow compressibility and variations of fluid viscosity and thermal conductivity with temperature. The intensity of NOB effects is measured by the dimensionless temperature differential $\unicode[STIX]{x1D716}$. Linear stability analysis of the thermal conduction state is performed. An $\unicode[STIX]{x1D716}^{2}$ scaling of the leading-order corrections of critical Rayleigh number $Ra_{cr}$ and disturbance growth rate $\unicode[STIX]{x1D70E}$ due to NOB effects is identified, which is a consequence of an intrinsic symmetry of the system. The influences of weak NOB effects on flow instabilities are further studied by perturbation expansion of linear stability equations with regard to $\unicode[STIX]{x1D716}$, and then the influence of aspect ratio $A$ is investigated in detail. NOB effects are found to enhance (weaken) flow stability in large (narrow) cavities. Detailed contributions of compressibility, viscosity and buoyancy actions on disturbance kinetic energy growth are identified quantitatively by energy analysis. Besides, a weakly nonlinear theory is developed based on centre-manifold reduction to investigate the NOB influences on bifurcation characteristics near convection onset, and amplitude equations are constructed for both codimension-one and -two cases. Rich bifurcation regimes are observed based on amplitude equations and also confirmed by direct numerical simulation. Weakly nonlinear analysis is useful for organizing and understanding these simulation results.
Somatic cell nuclear transfer (SCNT) is an important technique for life science research. However, most SCNT embryos fail to develop to term due to undefined reprogramming defects. Here, we show that abnormal Xi occurs in somatic cell NT blastocysts, whereas in female blastocysts derived from cumulus cell nuclear transfer, both X chromosomes were inactive. H3K27me3 removal by Kdm6a mRNA overexpression could significantly improve preimplantation development of NT embryos, and even reached a 70.2% blastocyst rate of cleaved embryos compared with the 38.5% rate of the control. H3K27me3 levels were significantly reduced in blastomeres from cloned blastocysts after overexpression of Kdm6a. qPCR indicated that rDNA transcription increased in both NT embryos and 293T cells after overexpression of Kdm6a. Our findings demonstrate that overexpression of Kdm6a improved the development of cloned mouse embryos by reducing H3K27me3 and increasing rDNA transcription.
This article is devoted to the study of some high-order difference schemes for the distributed-order time-fractional equations in both one and two space dimensions. Based on the composite Simpson formula and Lubich second-order operator, a difference scheme is constructed with convergence in the L1(L∞)-norm for the one-dimensional case, where τ,h and σ are the respective step sizes in time, space and distributed-order. Unconditional stability and convergence are proven. An ADI difference scheme is also derived for the two-dimensional case, and proven to be unconditionally stable and convergent in the L1(L∞)-norm, where h1 and h2 are the spatial step sizes. Some numerical examples are also given to demonstrate our theoretical results.