Kinetic investigations of heterogeneously catalyzed reactions on the Ru(001) and Pt(110)-(1x2) surfaces

The initial probabilities of activated, dissociative chemisorption of methane and ethane on Pt(110)-(1 x 2) have been measured. The surface temperature was varied from 450 to 900 K with the reactant gas temperature constant at 300 K. Under these conditions, we probe the kinetics of dissociation via trapping-mediated (as opposed to 'direct') mechanism. It was found that the probabilities of dissociation of both methane and ethane were strong functions of the surface temperature with an apparent activation energies of 14.4 kcal/mol for methane and 2.8 kcal/mol for ethane, which implys that the methane and ethane molecules have fully accommodated to the surface temperature. Kinetic isotope effects were observed for both reactions, indicating that the C-H bond cleavage was involved in the rate-limiting step. A mechanistic model based on the trapping-mediated mechanism is used to explain the observed kinetic behavior. The activation energies for C-H bond dissociation of the thermally accommodated methane and ethane on the surface extracted from the model are 18.4 and 10.3 kcal/mol, respectively.

The studies of the catalytic decomposition of formic acid on the Ru(001) surface with thermal desorption mass spectrometry following the adsorption of DCOOH and HCOOH on the surface at 130 and 310 K are described. Formic acid (DCOOH) chemisorbs dissociatively on the surface via both the cleavage of its O-H bond to form a formate and a hydrogen adatom, and the cleavage of its C-O bond to form a carbon monoxide, a deuterium adatom and an hydroxyl (OH). The former is the predominant reaction. The rate of desorption of carbon dioxide is a direct measure of the kinetics of decomposition of the surface formate. It is characterized by a kinetic isotope effect, an increasingly narrow FWHM, and an upward shift in peak temperature with Ɵ_T, the coverage of the dissociatively adsorbed formic acid. The FWHM and the peak temperature change from 18 K and 326 K at Ɵ_T = 0.04 to 8 K and 395 K at Ɵ_T = 0.89. The increase in the apparent activation energy of the C-D bond cleavage is largely a result of self-poisoning by the formate, the presence of which on the surface alters the electronic properties of the surface such that the activation energy of the decomposition of formate is increased. The variation of the activation energy for carbon dioxide formation with Ɵ_T accounts for the observed sharp carbon dioxide peak. The coverage of surface formate can be adjusted over a relatively wide range so that the activation energy for C-D bond cleavage in the case of DCOOH can be adjusted to be below, approximately equal to, or well above the activation energy for the recombinative desorption of the deuterium adatoms. Accordingly, the desorption of deuterium was observed to be governed completely by the desorption kinetics of the deuterium adatoms at low Ɵ_T, jointly by the kinetics of deuterium desorption and C-D bond cleavage at intermediate Ɵ_T, and solely by the kinetics of C-D bond cleavage at high Ɵ_T. The overall branching ratio of the formate to carbon dioxide and carbon monoxide is approximately unity, regardless the initial coverage Ɵ_T, even though the activation energy for the production of carbon dioxide varies with Ɵ_T. The desorption of water, which implies C-O bond cleavage of the formate, appears at approximately the same temperature as that of carbon dioxide. These observations suggest that the cleavage of the C-D bond and that of the C-O bond of two surface formates are coupled, possibly via the formation of a short-lived surface complex that is the precursor to to the decomposition.

The measurement of steady-state rate is demonstrated here to be valuable in determining kinetics associated with short-lived, molecularly adsorbed precursor to further reactions on the surface, by determining the kinetic parameters of the molecular precursor of formaldehyde to its dissociation on the Pt(110)-(1 x 2) surface.

Overlayers of nitrogen adatoms on Ru(001) have been characterized both by thermal desorption mass spectrometry and low-energy electron diffraction, as well as chemically via the postadsorption and desorption of ammonia and carbon monoxide.

The nitrogen-adatom overlayer was prepared by decomposing ammonia thermally on the surface at a pressure of 2.8 x 10^(-6) Torr and a temperature of 480 K. The saturated overlayer prepared under these conditions has associated with it a (√247/10 x √247/10)R22.7° LEED pattern, has two peaks in its thermal desorption spectrum, and has a fractional surface coverage of 0.40. Annealing the overlayer to approximately 535 K results in a rather sharp (√3 x √3)R30° LEED pattern with an associated fractional surface coverage of one-third. Annealing the overlayer further to 620 K results in the disappearance of the low-temperature thermal desorption peak and the appearance of a rather fuzzy p(2x2) LEED pattern with an associated fractional surface coverage of approximately one-fourth. In the low coverage limit, the presence of the (√3 x √3)R30° N overlayer alters the surface in such a way that the binding energy of ammonia is increased by 20% relative to the clean surface, whereas that of carbon monoxide is reduced by 15%.

A general methodology for the indirect relative determination of the absolute fractional surface coverages has been developed and was utilized to determine the saturation fractional coverage of hydrogen on Ru(001). Formaldehyde was employed as a bridge to lead us from the known reference point of the saturation fractional coverage of carbon monoxide to unknown reference point of the fractional coverage of hydrogen on Ru(001), which is then used to determine accurately the saturation fractional coverage of hydrogen. We find that ƟSAT/H = 1.02 (±0.05), i.e., the surface stoichiometry is Ru : H = 1 : 1. The relative nature of the method, which cancels systematic errors, together with the utilization of a glass envelope around the mass spectrometer, which reduces spurious contributions in the thermal desorption spectra, results in high accuracy in the determination of absolute fractional coverages.

Particle fluid mechanics in shear flows, acoustic waves, and shock waves

Three different categories of flow problems of a fluid containing small particles are being considered here. They are: (i) a fluid containing small, non-reacting particles (Parts I and II); (ii) a fluid containing reacting particles (Parts III and IV); and (iii) a fluid containing particles of two distinct sizes with collisions between two groups of particles (Part V).

Part I

A numerical solution is obtained for a fluid containing small particles flowing over an infinite disc rotating at a constant angular velocity. It is a boundary layer type flow, and the boundary layer thickness for the mixture is estimated. For large Reynolds number, the solution suggests the boundary layer approximation of a fluid-particle mixture by assuming W = W_p. The error introduced is consistent with the Prandtl’s boundary layer approximation. Outside the boundary layer, the flow field has to satisfy the “inviscid equation” in which the viscous stress terms are absent while the drag force between the particle cloud and the fluid is still important. Increase of particle concentration reduces the boundary layer thickness and the amount of mixture being transported outwardly is reduced. A new parameter, β = 1/Ω τ_v, is introduced which is also proportional to μ. The secondary flow of the particle cloud depends very much on β. For small values of β, the particle cloud velocity attains its maximum value on the surface of the disc, and for infinitely large values of β, both the radial and axial particle velocity components vanish on the surface of the disc.

Part II

The “inviscid” equation for a gas-particle mixture is linearized to describe the flow over a wavy wall. Corresponding to the Prandtl-Glauert equation for pure gas, a fourth order partial differential equation in terms of the velocity potential ϕ is obtained for the mixture. The solution is obtained for the flow over a periodic wavy wall. For equilibrium flows where λ_v and λ_T approach zero and frozen flows in which λ_v and λ_T become infinitely large, the flow problem is basically similar to that obtained by Ackeret for a pure gas. For finite values of λ_v and λ_T, all quantities except v are not in phase with the wavy wall. Thus the drag coefficient C_D is present even in the subsonic case, and similarly, all quantities decay exponentially for supersonic flows. The phase shift and the attenuation factor increase for increasing particle concentration.

Part III

Using the boundary layer approximation, the initial development of the combustion zone between the laminar mixing of two parallel streams of oxidizing agent and small, solid, combustible particles suspended in an inert gas is investigated. For the special case when the two streams are moving at the same speed, a Green’s function exists for the differential equations describing first order gas temperature and oxidizer concentration. Solutions in terms of error functions and exponential integrals are obtained. Reactions occur within a relatively thin region of the order of λ_D. Thus, it seems advantageous in the general study of two-dimensional laminar flame problems to introduce a chemical boundary layer of thickness λ_D within which reactions take place. Outside this chemical boundary layer, the flow field corresponds to the ordinary fluid dynamics without chemical reaction.

Part IV

The shock wave structure in a condensing medium of small liquid droplets suspended in a homogeneous gas-vapor mixture consists of the conventional compressive wave followed by a relaxation region in which the particle cloud and gas mixture attain momentum and thermal equilibrium. Immediately following the compressive wave, the partial pressure corresponding to the vapor concentration in the gas mixture is higher than the vapor pressure of the liquid droplets and condensation sets in. Farther downstream of the shock, evaporation appears when the particle temperature is raised by the hot surrounding gas mixture. The thickness of the condensation region depends very much on the latent heat. For relatively high latent heat, the condensation zone is small compared with Ʌ_D.

For solid particles suspended initially in an inert gas, the relaxation zone immediately following the compression wave consists of a region where the particle temperature is first being raised to its melting point. When the particles are totally melted as the particle temperature is further increased, evaporation of the particles also plays a role.

The equilibrium condition downstream of the shock can be calculated and is independent of the model of the particle-gas mixture interaction.

Part V

For a gas containing particles of two distinct sizes and satisfying certain conditions, momentum transfer due to collisions between the two groups of particles can be taken into consideration using the classical elastic spherical ball model. Both in the relatively simple problem of normal shock wave and the perturbation solutions for the nozzle flow, the transfer of momentum due to collisions which decreases the velocity difference between the two groups of particles is clearly demonstrated. The difference in temperature as compared with the collisionless case is quite negligible.

Biocompatibilidade do polissacarídeo celulósico sintetizado por zoogloea sp.: um estudo em bexigas de coelho

Atualmente, o material utilizado para o tratamento endoscópico é o Deflux, porém este é um material não-biológico. Sabe-se que a substância ideal deve ser atóxica, biocompatível, não-migratória, não-antigênica e deve causar o mínimo possível de inflamação no local do implante. A bactéria Zoogloea sp. produz um exopolissacarídeo celulósico (CEP) com baixa citotoxicidade e alto biocompatibilidade. O objetivo deste estudo é investigar, na bexiga de coelho, a biocompatibilidade de implantes de exopolissacarídeo de celulose, produzidos pela Zooglea sp. Foram utilizados como modelo experimental, 20 coelhos adultos da raça Califórnia, com média de seis meses de idade. Os animais foram divididos em dois grupos, sendo o grupo G1, composto por animais mortos três dias após a aplicação do implante (n=9), e o grupo G2, composto por animais mortos três meses após a aplicação do implante (n=11). Cada animal recebeu, no total, quatro implantes, sendo dois de gel de biopolímero e dois de gel Deflux. Foram realizadas as técnicas imunohistoquímicas para marcação de colágeno tipos I e III, alfa-actina de músculo liso, PCNA e reação química TUNEL. Nas amostras de três dias, os implantes de CEP e deflux, eram estruturalmente homogêneos e livres de células inflamatórias ou vasos sanguíneos. Por outro lado, nas amostras de três meses, com exceção de algumas áreas, o CEP estava organizado como feixes curtos que eram sugestivos de um tecido fibroso. Apesar disso, o implante de CEP corou negativamente para colágenos tipos I e III, fibras elásticas, enquanto que o tricrômico de masson, não indicou a presença de colágeno. Em contraste as áreas de implante de deflux nas amostras de três meses estavam fragmentadas, mas ainda eram homogêneas, e ainda não havia nenhuma célula nem vaso sanguíneo em seu interior. As células positivas para PCNA podiam ser claramente percebidas dentro dessas ilhotas, dessa forma indicando um processo inflamatórioproliferativo, em curso. No grupo sacrificado aos três meses, os implantes de deflux ainda estavam negativos, mas em torno das áreas de CEP algumas células positivas para a técnica do TUNEL eram perceptíveis. Nos implantes de CEP de três meses, muitos vasos sanguíneos eram visualizados, e a sua densidade era de 23.865.48. A densidade de microvasos na lâmina própria (41.5111.19) foi significativamente diferente (p<0.001) daquela no implante de CEP. Nossos resultados mostraram que o CEP possui pouca imunogenicidade e se integra melhor no tecido hospedeiro quando comparado ao deflux. Portanto o CEP deve ser um material eficiente em casos em que a incorporação ao tecido é desejada como por exemplo em estruturas de suporte na cirurgia de reconstrução

Estimating material elasticity by spherical indentation load-relaxation tests on viscoelastic samples of finite thickness

A two-step viscoelastic spherical indentation method is proposed to compensate for 1) material relaxation and 2) sample thickness. In the first step, the indenter is moved at a constant speed and the reaction force is measured. In the second step, the indenter is held at a constant position and the relaxation response of the material is measured. Then the relaxation response is fit with a multi-exponential function which corresponds to a three-branch general Maxwell model. The relaxation modulus is derived by correcting the finite ramp time introduced in the first step. The proposed model takes into account the sample thickness, which is important for applications in which the sample thickness is less than ten times the indenter radius. The model is validated numerically by finite element simulations. Experiments are carried out on a 10% gelatin phantom and a chicken breast sample with the proposed method. The results for both the gelatin phantom and the chicken breast sample agree with the results obtained from a surface wave method. Both the finite element simulations and experimental results show improved elasticity estimations by incorporating the sample thickness into the model. The measured shear elasticities of the 10% gelatin sample are 6.79 and 6.93 kPa by the proposed finite indentation method at sample thickness of 40 and 20 mm, respectively. The elasticity of the same sample is estimated to be 6.53 kPa by the surface wave method. For the chicken breast sample, the shear elasticity is measured to be 4.51 and 5.17 kPa by the proposed indentation method at sample thickness of 40 and 20 mm, respectively. Its elasticity is measured by the surface wave method to be 4.14 kPa. © 2011 IEEE.

Nanodiamond formation by hot-filament chemical vapor deposition on carbon ions bombarded Si

Nanocrystalline diamond films were grown by a two-step process on Si(1 0 0) substrate, which was first pretreated by pure carbon ions bombardment. The bombarded Si substrate was then transformed into a hot-filament chemical vapor deposition (HFCVD) system for further growth. Using the usual CH4/H-3 feed gas ratio for micro crystalline diamond growth, nanodiamond crystallites were obtained. The diamond nucleation density is comparable to that obtained by biasing the substrate. The uniformly distributed lattice damage is proposed to be responsible for the formation of the nanodiamond. (C) 2002 Elsevier Science B.V. All rights reserved.

Preparation of highly active Pt/C cathode electrocatalysts for DMFCs by an improved aqueous impregnation method

An improved aqueous impregnation method was used to prepare 40 wt% Pt/C electrocatalysts. TEM analysis of the samples showed that the Pt particles impregnated for a short time have a very narrow size distribution in the range of 1-4 nm with an average size of 2.6 nm. UV-vis spectroscopy measurements verified that the redox reaction between PtCl62- and formaldehyde took place with a slow rate at ambient temperature via a two-step reaction path, where PtCl42- serves as an intermediate. The use of the short-time-impregnated 40 wt% Pt/C as cathode electrocatalysts in direct methanol fuel cells yields better performance than that of commercial 40 wt% Pt/C electrocatalyst. Experimental evidence provides clues for the fundamental understanding of elementary steps of the redox reactions, which helps in guiding the design and preparation of highly dispersed Pt catalyst for fuel cells.

Facile synthesis of highly uniform octahedral LuVO4 microcrystals by a facile chemical conversion method

Uniform octahedral LuVO4 microcrystals have been successfully prepared through a designed two-step hydrothermal method. One-dimensional lutetium precursor was first prepared through a simple hydrothermal route. Subsequently, a well-shaped octahedral LuVO4 sample was synthesized at the expense of the wirelike precursors during the hydrothermal process. The whole process in this method was carried out in aqueous conditions without the use of any organic solvents, surfactant, or catalyst. The conversion process from nanowire precursor to octahedral product has been investigated in detail. The LuVO4 : Ln(3+) (Ln Eu, Dy, Sm, and Er) phosphors show strong light emissions with different colors coming from different activator ions under ultraviolet light excitation or low-voltage electron beam excitation. Furthermore, this general and facile method may be of much significance in the synthesis of many other lanthanide compounds with polyhedral morphology.

High-efficiency and low-cost hybrid nanomaterial as enhancing electrocatalyst: Spongelike AWN core/shell nanomaterial with hollow cavity

A high-efficiency and low-cost spongelike Au/Pt core/shell electrocatalyst with hollow cavity has been facilely obtained via a simple two-step wet chemical process. Hollow gold nanospheres were first synthesized via a modified galvanic replacement reaction between Co nanoparticles in situ produced and HAUCl(4). The as-prepared gold hollow spheres were employed as seeds to further grow spongelike Pt shell. It is found that the surface of this hybrid nanomaterial owns many Pt nanospikes, which form a spongelike nanostructure. All experimental data including scanning electron microscopy, transmission electron microscopy, X-ray photoelectron spectroscopy, and UV-vis-near-infrared spectroscopy have been employed to characterize the obtained Au/Pt hybrid nanomaterial. The rapid development of fuel cell has inspired us to investigate the electrocatalytic properties for dioxygen and methanol of this novel hybrid nanomaterial. Spongelike hybrid nanomaterial mentioned here exhibits much higher catalytic activity for dioxygen reduction and methanol oxidation than the common Pt electrode.

Analysis of chemical calorific effect during reactive extrusion processes for free radical polymerization

In the reactive extrusion process for polymerization, the chemical calorific effect has a great influence on the temperature. In order to quantitatively analyze the polymerization trend and optimize the processing conditions, the phenomena of the chemical calorific effect during reactive extrusion processes for free radical polymerization were analyzed. Numerical computation expressions of the heat of chemical reaction and the reactive calorific intensity were deduced, and then a numerical simulation of the reactive extrusion process for the polymerization of n-butyl methacrylate was carried out. The evolutions of the heat of chemical reaction and the reactive calorific intensity along the! axial direction of the extruder are presented, on the basis of which reactive processing conditions can be optimized.

Shape-memory and biocompatibility properties of segmented polyurethanes based on poly(L-lactide)

A series of segmented poly (L-lactide)-polyurethanes (PLA-PU) were synthesized by a two-step method, with oligo-poly(L-lactide) (PLA) as the soft segments and the reaction product of 2,4-toluene diisocyanate(TDI) and ethylene glycol(EG) as the hard segments. The shape memory properties of PLA-PUs were examined. The processed PLA-PUs could recover almost 100% to their original shape within 10 degrees C from the lowest recovery temperature. In the recovery process, the PLA-PUs showed a maximum contracting stress of shape change in the range of 1.5-4 MPa depending on the PLA segmental length and the hard-segmental content and higher than that of poly (e-caprolactone polyurethane) (PCL-PU). Besides, the influence of deforming and fixing temperatures on shape memory properties of PLA-PU was studied in detail. They could affect not only the recovery temperature but also the maximum contracting stress. The experiments of cell incubation were used to evaluate the biocompatibility of PLA-PU. The results show that the biocompatibility of PLA-PU is comparable to that of the pure PLA. This kind of polyurethane can be used as implanted medical devices with a shape memory property.

Synthesis and properties of novel polyimides derived from 2,2 ',3,3 '-benzophenonetetracarboxylic dianhydride

A new synthetic route to 2,2',3,3'-BTDA (where BTDA is benzophenonetetracarboxylic dianhydride), an isomer of 2,3',3',4'-BTDA and 3,3',4,4'-BTDA, is described. Single-crystal X-ray diffraction analysis of 2,2',3,3'-BTDA has shown that this dianhydride has a bent and noncoplanar structure. The polymerizations of 2,2',3,3'-BTDA with 4,4'-oxydianiline (ODA) and 4,4'-bis(4-aminophenoxy)benzene (TPEQ) have been investigated with a conventional two-step process. A trend of cyclic oligomers forming in the reaction of 2,2',3,3'-BTDA and ODA has been found and characterized with IR, NMR, matrix-assisted laser desorption/ionization time-of-flight mass spectrometry, and elemental analyses. Films based on 2,2',3,3'-BTDA/TPEQ can only be obtained from corresponding polyimide (PI) solutions prepared by chemical imidization because those from their polyamic acids by thermal imidization are brittle. PIS from 2,2',3,3'-BTDA have lower inherent viscosities and worse thermal and mechanical properties than the corresponding 2,3',3',4'-BTDA- and 3,3',4,4'-BTDA-based PIs. PIs from 2,2',3,3'-BTDA and 2,3',3',4'-BTDA are amorphous, whereas those from 3,3',4,4'-BTDA have some crystallinity, according to wide-angle X-ray diffraction.

Syntheses and properties of novel polyimides derived from 2-(4-aminophenyl)-5-aminopyrimidine

2-(4-Aminophenyl)-5-aminopyrimidine (4) is synthesized via a condensation reaction of vinamidium salts and amidine chloride salts, followed by hydrazine palladium catalyzed reduction. A series of novel homo- and copolyimides containing pyrimidine unit are prepared from the diamine and 1,4-phenylenediamine (PDA) with pyromellitic dianhydride (PMDA) or 3,3',4,4'-biphenyl tertracarboxylic dianhydride (BPDA) via a conventional two-step thermal imidization method. The poly(amic acid) precursors had inherent viscosities of 0.97-4.38 dL/g (c = 0.5 g/dL, in DMAc, 30 degrees C) and all of them could be cast and thermally converted into flexible and tough polyimide films. All of the polyimides showed excellent thermal stability and mechanical properties. The glass transition temperatures of the resulting polyimides are in the range of 307-434 degrees C and the 10% weight loss temperature is in the range of 556-609 degrees C under air. The polyimide films possess strength at break in the range of 185-271 MPa, elongations at break in the range of 6.8-51%, and tensile modulus in the range of 3.5-6.46 GPa. The polymer films are insoluble in common organic solvents, exhibiting high chemical resistance.

Kinetics of Y(III) stripping for the system Y/HCl/Cyanex 272-P507/heptane with constant interfacial area cell with laminar flow

Kinetics and mechanism of stripping of yttrium(III) previously extracted by mixtures of bis(2,4,4-trimethylpentyl)phosphinic acid (Cyanex 272, HA), and 2-ethylhexyl phosphonic acid mono-2-ethylhexl ester (P507, HB) dissolved in heptane have been investigated by constant interfacial-area cell by laminar flow. The corresponding equilibrium stripping equation and equilibrium constant were obtained. The studies of effects of the stirring rate and temperature on the stripping rate show that the stripping regime is dependent on the stripping conditions. The plot of interfacial area on the rate has shown a linear relationship. This fact together with the strong surface activity of mixtures of Cyanex 272 and P507 at heptane-water interfaces makes the interface the most probable locale for the chemical reactions. The stripping rate constant is obtained, and the value is compared with that of the system with Cyanex 272 and P507 alone. It is concluded that the stripping ability with the mixtures is easier than that of P507 due to lower the activation energy of the mixtures. The stripping rate equation has also been obtained, and the rate-determining steps are the two-step interfacial chemical reactions as predicted from interfacial reaction models.

Extraction kinetics of rare earth elements with sec-octylphenoxy acetic acid

The kinetics of RE (La, Gd, Er, Yb and Y) extraction with sec-octylphenoxy acetic acid was investigated using a constant interfacial area cell with laminar flow at 303 K. The natures of the extracted complexes have some effect on the extraction rate which is controlled by the reaction rate of M(III) and extractant molecules at two-phase interface for Er(III), Yb(III) and Y(III), by a mixed chemical reaction-diffusion for Gd(III) and a diffusion for La( III). The extractant molecules tend to adsorb at the interface. So an interfacial extraction reaction model was derived.

中低纬电离层暴的统计分析与经验建模

As a key issue of ionospheric weather study, systemic studies on ionospheric storms can not only further improve our understanding of the response of the ionosphere to solar and geomagnetic disturbances, but also help us to reveal the chemical, dynamic and electro-dynamic mechanisms during storms. Empirical modelling for regional ionospheric storm is also very useful, because it can provide us with tools and references for the forecasting and further practical application of ionospheric activity. In this thesis, we focus on describing and forecasting of ionospheric storms at middle and low latitudes. The main points of my investigations are listed as follows. (1) By using magnetic storms during the period over 50 years, the dependence of the type, onset time and time delay of the ionospheric storms on magnetic latitude, season and local time at middle and low latitudes in the East-Asian sector are studied. The results show that the occurrences of the types of ionospheric disturbances differ in latitude and season. The onset of the ionospheric storms depends on local time. At middle latitudes, most negative phase onsets are within the local time interval from night to early morning, and they rarely occurred in the local noon and afternoon sectors. At low latitudes, positive phases commence most frequently in the daytime sector as well as pre-midnight sector. The average time delays for both the positive and negative ionospheric storms increase with descending latitudes. The time delay has significant dependence on the local time of main phase onset (MPO). The time delay of positive response is shorter for daytime MPO and longer for night-time MPO, whereas the opposite applies for negative response. (2) Based on some previous researches, a primary empirical model for mid-latitude ionospheric disturbance is set up. By fitting to the observed data, we get a high accuracy with a mean RMSE of only 12-14% in summer and equinox. The model output has been compared with the output of STORM model, and the results show that, our model is much better than STORM in summer and a little better for some mid-latitude stations at equinox. Especially, for the type of two-step geomagnetic storm, our model can present twice descending of foF2 very well. In addition, our model can forecast positive ionospheric storms.