This study aims to develop a chloride diffusion simulation method that considers the hydration microstructure and pore solution properties during the hydration of tricalcium silicate(C3S).The method combines the hydra...This study aims to develop a chloride diffusion simulation method that considers the hydration microstructure and pore solution properties during the hydration of tricalcium silicate(C3S).The method combines the hydration simulation,thermodynamic calculation,and finite element analysis to examine the effects of pore solution,including effect of electrochemical potential,effect of chemical activity,and effect of mechanical interactions between ions,on the chloride effective diffusion coefficient of hydrated C3S paste.The results indicate that the effect of electrochemical potential on chloride diffusion becomes stronger with increasing hydration age due to the increase in the content of hydrated calcium silicate;as the hydration age increases,the effect of chemical activity on chloride diffusion weakens when the number of diffusible elements decreases;the effect of mechanical interactions between ions on chloride diffusion decreases with the increase of hydration age.展开更多
For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geomet...For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.展开更多
The interfacial reactions in partial transient liquid-phase bonding of Si3N4 ceramics with Ti/Ni/Ti interlayers were studied by means of scanning electron microscopy （SEM）, energy dispersive spectrometry （EDS） and...The interfacial reactions in partial transient liquid-phase bonding of Si3N4 ceramics with Ti/Ni/Ti interlayers were studied by means of scanning electron microscopy （SEM）, energy dispersive spectrometry （EDS） and X-ray diffractometry （XRD）. It was shown that the interfacial structure of Si3N4/TiN/Ti5Si3+Ti5Si4 + Ni3Si/ （NiTi ） /Ni3Ti/ Ni was formed after bonding. The activation energies for TiN layer and the mixed reaction layer of Ti5Si3 + Ti5Si4 + Ni3Si are 546. 8 kJ/mol and 543. 9 kJ/mol, respectively. The formation and transition processes of interface layer sequence in the joint were clarified by diffusion path. An important characteristic, which is different from the conventional brazing and soid-state diffusion bonding, has been found, i. e., during the partial transient liquid-phase bonding, not only the reaction layers which have formed grow, but also the diffusion path in the subsequent reaction changes because of the remarkable variation of the concentration on the metal side.展开更多
In this paper, two finite difference streamline diffusion (FDSD) schemes for solving two-dimensional time-dependent convection-diffusion equations are constructed. Stability and optimal order error estimati-ions for c...In this paper, two finite difference streamline diffusion (FDSD) schemes for solving two-dimensional time-dependent convection-diffusion equations are constructed. Stability and optimal order error estimati-ions for considered schemes are derived in the norm stronger than L^2-norm.展开更多
In this article the travelling wave solution for a class of nonlinear reaction diffusion problems are considered. Using the homotopic method and the theory of travelling wave transform, the approximate solution for th...In this article the travelling wave solution for a class of nonlinear reaction diffusion problems are considered. Using the homotopic method and the theory of travelling wave transform, the approximate solution for the corresponding problem is obtained.展开更多
The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte ...The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.展开更多
Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the movi...Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.展开更多
Coal seam gas content is frequently measured in quantity during underground coal mining operation and coalbed methane(CBM)exploration as a significant basic parameter.Due to the calculation error of lost gas and resid...Coal seam gas content is frequently measured in quantity during underground coal mining operation and coalbed methane(CBM)exploration as a significant basic parameter.Due to the calculation error of lost gas and residual gas in the direct method,the efficiency and accuracy of the current methods are not inadequate to the large area multi-point measurement of coal seam gas content.This paper firstly deduces a simplified theoretical dynamic model for calculating lost gas based on gas dynamic diffusion theory.Secondly,the effects of various factors on gas dynamic diffusion from coal particle are experimentally studied.And sampling procedure of representative coal particle is improved.Thirdly,a new estimation method of residual gas content based on excess adsorption and competitive adsorption theory is proposed.The results showed that the maximum error of calculating the losing gas content by using the new simplified model is only 4%.Considering the influence of particle size on gas diffusion law,the particle size of the collected coal sample is below 0.25 mm,which improves the measurement speed and reflects the safety representativeness of the sample.The determination time of gas content reduced from 36 to 3 h/piece.Moreover,the absolute error is 0.15–0.50 m^3/t,and the relative error is within 5%.A new engineering method for determining the coal seam gas content is developed according to the above research.展开更多
We consider the functional separation of variables to the nonlinear diffusion equation with source and convection term： ut = （A（x）D（u）ux）x ＋ B（x）Q（u）, Ax ≠ 0. The functional separation of variables to thi...We consider the functional separation of variables to the nonlinear diffusion equation with source and convection term： ut = （A（x）D（u）ux）x ＋ B（x）Q（u）, Ax ≠ 0. The functional separation of variables to this equation is studied by using the group foliation method. A classification is carried out for the equations which admit the function separable solutions. As a consequence, some solutions to the resulting equations are obtained.展开更多
In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existenc...In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results.展开更多
In this paper,we consider a Cauchy problem of the time fractional diffusion equation(TFDE)in x∈[0,L].This problem is ubiquitous in science and engineering applications.The illposedness of the Cauchy problem is explai...In this paper,we consider a Cauchy problem of the time fractional diffusion equation(TFDE)in x∈[0,L].This problem is ubiquitous in science and engineering applications.The illposedness of the Cauchy problem is explained by its solution in frequency domain.Furthermore,the problem is formulated into a minimization problem with a modified Tikhonov regularization method.The gradient of the regularization functional based on an adjoint problem is deduced and the standard conjugate gradient method is presented for solving the minimization problem.The error estimates for the regularized solutions are obtained under Hp norm priori bound assumptions.Finally,numerical examples illustrate the effectiveness of the proposed method.展开更多
In this article,we consider to solve the inverse initial value problem for an inhomogeneous space-time fractional diffusion equation.This problem is ill-posed and the quasi-boundary value method is proposed to deal wi...In this article,we consider to solve the inverse initial value problem for an inhomogeneous space-time fractional diffusion equation.This problem is ill-posed and the quasi-boundary value method is proposed to deal with this inverse problem and obtain the series expression of the regularized solution for the inverse initial value problem.We prove the error estimates between the regularization solution and the exact solution by using an a priori regularization parameter and an a posteriori regularization parameter choice rule.Some numerical results in one-dimensional case and two-dimensional case show that our method is efficient and stable.展开更多
A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretizatio...A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming （P1）2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.展开更多
Aim The purpose of this study was to develop a mathe-matical model to quantitatively describe the passive trans-port of macromolecules within dental biofilms. Methodology Fluorescently labeled dextrans with different ...Aim The purpose of this study was to develop a mathe-matical model to quantitatively describe the passive trans-port of macromolecules within dental biofilms. Methodology Fluorescently labeled dextrans with different molecular mass （3 kD,10 kD,40 kD,70 kD,2 000 kD） were used as a series of diffusion probes. Streptococcus mutans,Streptococcus sanguinis,Actinomyces naeslundii and Fusobacterium nucleatum were used as inocula for biofilm formation. The diffusion processes of different probes through the in vitro biofilm were recorded with a confocal laser microscope. Results Mathematical function of biofilm penetration was constructed on the basis of the inverse problem method. Based on this function,not only the relationship between average concentration of steady-state and molecule weights can be analyzed,but also that between penetrative time and molecule weights. Conclusion This can be used to predict the effective concentration and the penetrative time of anti-biofilm medicines that can diffuse through oral biofilm. Further-more,an improved model for large molecule is proposed by considering the exchange time at the upper boundary of the dental biofilm.展开更多
In this paper, the cubic and quintic diffusion equation under stochastic non homogeneity is solved using Wiener- Hermite expansion and perturbation (WHEP) technique, Homotopy perturbation method (HPM) and Pickard appr...In this paper, the cubic and quintic diffusion equation under stochastic non homogeneity is solved using Wiener- Hermite expansion and perturbation (WHEP) technique, Homotopy perturbation method (HPM) and Pickard approximation technique. The analytic solution of the linear case is obtained using Eigenfunction expansion .The Picard approximation method is used to introduce the first and second order approximate solution for the non linear case. The WHEP technique is also used to obtain approximate solution under different orders and different corrections. The Homotopy perturbation method (HPM) is also used to obtain some approximation orders for mean and variance. Using mathematica-5, the methods of solution are illustrated through figures, comparisons among different methods and some parametric studies.展开更多
A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order...A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.展开更多
The formation mechanism of calcium vanadate and manganese vanadate and the difference between calcium and manganese in the reaction with vanadium are basic issues in the calcification roasting and manganese roasting p...The formation mechanism of calcium vanadate and manganese vanadate and the difference between calcium and manganese in the reaction with vanadium are basic issues in the calcification roasting and manganese roasting process with vanadium slag.In this work,CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) diffusion couples were prepared and roasted for different time periods to illustrate and compare the diffusion reaction mechanisms.Then,the changes in the diffusion product and diffusion coefficient were investigated and calculated based on scanning electron microscopy (SEM) with energy dispersive X-ray spectroscopy (EDS) analysis.Results show that with the extension of the roasting time,the diffusion reaction gradually proceeds among the CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) diffusion couples.The regional boundaries of calcium and vanadium are easily identifiable for the CaO–V_(2)O_(5) diffusion couple.Meanwhile,for the MnO_(2)–V_(2)O_(5) diffusion couple,MnO_(2) gradually decomposes to form Mn_(2)O_(3),and vanadium diffuses into the interior of Mn_(2)O_(3).Only a part of vanadium combines with manganese to form the diffusion production layer.CaV_(2)O_(6) and MnV_(2)O_(6) are the interfacial reaction products of the CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) diffusion couples,respectively,whose thicknesses are 39.85 and 32.13μm when roasted for 16 h.After 16 h,both diffusion couples reach the reaction equilibrium due to the limitation of diffusion.The diffusion coefficient of the CaO–V_(2)O_(5) diffusion couple is higher than that of the MnO_(2)–V_(2)O_(5) diffusion couple for the same roasting time,and the diffusion reaction between vanadium and calcium is easier than that between vanadium and manganese.展开更多
This paper aims to present a new streamline diffusion method with low order rectangular Bernardi-Raugel elements to solve the generalized Oseen equations. With the help of the Bramble-Hilbert lemma, the optimal errors...This paper aims to present a new streamline diffusion method with low order rectangular Bernardi-Raugel elements to solve the generalized Oseen equations. With the help of the Bramble-Hilbert lemma, the optimal errors of the velocity and pressure are estimated, which are independent of the considered parameter e. With an interpolation postprocessing approach, the superconvergent error of the pressure is obtained. Finally, a numerical experiment is carried out to confirm the theoretical results.展开更多
With the rapid development of Mg alloys,deeper understanding to the thermodynamic and diffusional kinetic behavior of intermetallic compounds(IMCs)is important for studying the effect of alloying elements to the micro...With the rapid development of Mg alloys,deeper understanding to the thermodynamic and diffusional kinetic behavior of intermetallic compounds(IMCs)is important for studying the effect of alloying elements to the microstructure evolution.Specially,a systematic quantitative investigation on the diffusional growth of IMCs is of great necessity.However,the works studying the elemental diffusion behaviors of multiple-element IMCs are rare in magnesium alloy systems.The current work takes the ternary Mg-Al-Zn system as research target,and combines the diffusion couple technique,phase stability diagrams,in-situ observation technique and numerical inverse method to investigate the temperature-dependent kinetic coefficients.The parabolic growth constant(PGC)and interdiffusion coefficients for Mg solid-solution phase andγ-Mg_(17)Al_(12),β-Mg_(2)Al_(3),ε-Mg_(23)Al_(30),MgZn_(2),Mg_(2)Zn_(3),τ-Mg_(32)(Zn,Al)49 andφ-Mg_(5)Zn_(2)Al_(2) IMCs in the Mg-Al-Zn alloy system are determined.By comparing the current experimental with calculation results,the rate-controlling factor of the temperature-dependent diffusion growth ofφ,τandεternary IMCs in the Mg-Al-Zn system is further discussed in detail.展开更多
A localized space-time method of fundamental solutions(LSTMFS)is extended for solving three-dimensional transient diffusion problems in this paper.The interval segmentation in temporal direction is developed for the a...A localized space-time method of fundamental solutions(LSTMFS)is extended for solving three-dimensional transient diffusion problems in this paper.The interval segmentation in temporal direction is developed for the accurate simulation of long-time-dependent diffusion problems.In the LSTMFS,the whole space-time domain with nodes arranged i divided into a series of overlapping subdomains with a simple geometry.In each subdomain,the conventional method of fundamental solutions is utilized and the coefficients associated with the considered domain can be explicitly determined.By calculating a combined sparse matrix system,the value at any node inside the space-time domain can be obtained.Numerical experi-ments demonstrate that high accuracy and efficiency can be simultaneously achieved via the LSTMFS,even for the problems defined on a long-time and quite complex computational domain.展开更多
基金Funded by the Natural Science Foundation of Jiangsu Province(No.BK20241529)China Postdoctoral Science Foundation(No.2024M750736)。
文摘This study aims to develop a chloride diffusion simulation method that considers the hydration microstructure and pore solution properties during the hydration of tricalcium silicate(C3S).The method combines the hydration simulation,thermodynamic calculation,and finite element analysis to examine the effects of pore solution,including effect of electrochemical potential,effect of chemical activity,and effect of mechanical interactions between ions,on the chloride effective diffusion coefficient of hydrated C3S paste.The results indicate that the effect of electrochemical potential on chloride diffusion becomes stronger with increasing hydration age due to the increase in the content of hydrated calcium silicate;as the hydration age increases,the effect of chemical activity on chloride diffusion weakens when the number of diffusible elements decreases;the effect of mechanical interactions between ions on chloride diffusion decreases with the increase of hydration age.
文摘For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.
文摘The interfacial reactions in partial transient liquid-phase bonding of Si3N4 ceramics with Ti/Ni/Ti interlayers were studied by means of scanning electron microscopy （SEM）, energy dispersive spectrometry （EDS） and X-ray diffractometry （XRD）. It was shown that the interfacial structure of Si3N4/TiN/Ti5Si3+Ti5Si4 + Ni3Si/ （NiTi ） /Ni3Ti/ Ni was formed after bonding. The activation energies for TiN layer and the mixed reaction layer of Ti5Si3 + Ti5Si4 + Ni3Si are 546. 8 kJ/mol and 543. 9 kJ/mol, respectively. The formation and transition processes of interface layer sequence in the joint were clarified by diffusion path. An important characteristic, which is different from the conventional brazing and soid-state diffusion bonding, has been found, i. e., during the partial transient liquid-phase bonding, not only the reaction layers which have formed grow, but also the diffusion path in the subsequent reaction changes because of the remarkable variation of the concentration on the metal side.
基金Project supported by National Natural Science Foundation of China and China State Key project for Basic Researchcs.
文摘In this paper, two finite difference streamline diffusion (FDSD) schemes for solving two-dimensional time-dependent convection-diffusion equations are constructed. Stability and optimal order error estimati-ions for considered schemes are derived in the norm stronger than L^2-norm.
基金Supported by the National Natural Sciences Foundation of China(40676016 and 10471039)the National Key Project for Basic Research(2003CB415101-03 and 2004CB418304)+2 种基金the Key Project of the Chinese Academy of Sciences(KZCX3-SW-221)in part by E-Institutes of Shanghai Municipal Education Commission(N.E03004)the Natural Science Foundation of Zeijiang,China(Y606268).
文摘In this article the travelling wave solution for a class of nonlinear reaction diffusion problems are considered. Using the homotopic method and the theory of travelling wave transform, the approximate solution for the corresponding problem is obtained.
文摘The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.
基金Project supported by the National Natural Science Foundation of China(Grant No.11072117)the Natural Science Foundation of Ningbo City,China(GrantNo.2013A610103)+2 种基金the Natural Science Foundation of Zhejiang Province,China(Grant No.Y6090131)the Disciplinary Project of Ningbo City,China(GrantNo.SZXL1067)the K.C.Wong Magna Fund in Ningbo University,China
文摘Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.
基金the National Natural Science Foundation of China(51774119,51374095,and 51604092)the primary research projects of critical scientific research in Henan Colleges and Universities(19zx003)+1 种基金Program for Innovative Research Team in University of Ministry of Education of China(IRT_16R22)State Key Laboratory Cultivation Base for Gas Geology and Gas Control(Henan Polytechnic University)(WS2018A02)。
文摘Coal seam gas content is frequently measured in quantity during underground coal mining operation and coalbed methane(CBM)exploration as a significant basic parameter.Due to the calculation error of lost gas and residual gas in the direct method,the efficiency and accuracy of the current methods are not inadequate to the large area multi-point measurement of coal seam gas content.This paper firstly deduces a simplified theoretical dynamic model for calculating lost gas based on gas dynamic diffusion theory.Secondly,the effects of various factors on gas dynamic diffusion from coal particle are experimentally studied.And sampling procedure of representative coal particle is improved.Thirdly,a new estimation method of residual gas content based on excess adsorption and competitive adsorption theory is proposed.The results showed that the maximum error of calculating the losing gas content by using the new simplified model is only 4%.Considering the influence of particle size on gas diffusion law,the particle size of the collected coal sample is below 0.25 mm,which improves the measurement speed and reflects the safety representativeness of the sample.The determination time of gas content reduced from 36 to 3 h/piece.Moreover,the absolute error is 0.15–0.50 m^3/t,and the relative error is within 5%.A new engineering method for determining the coal seam gas content is developed according to the above research.
基金The project supported by National Natural Science Foundation of China under Grant No. 10371098 and the Program for New Century Excellent Talents in Universities under Grant No. NCET-04-0968
文摘We consider the functional separation of variables to the nonlinear diffusion equation with source and convection term： ut = （A（x）D（u）ux）x ＋ B（x）Q（u）, Ax ≠ 0. The functional separation of variables to this equation is studied by using the group foliation method. A classification is carried out for the equations which admit the function separable solutions. As a consequence, some solutions to the resulting equations are obtained.
基金supported by the National Basic Research Program under the Grant 2005CB321701the National Natural Science Foundation of China under the Grants 60474027 and 10771211.
文摘In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results.
基金Supported by the National Natural Science Foundation of China(Grant No.11471253 and No.11571311)
文摘In this paper,we consider a Cauchy problem of the time fractional diffusion equation(TFDE)in x∈[0,L].This problem is ubiquitous in science and engineering applications.The illposedness of the Cauchy problem is explained by its solution in frequency domain.Furthermore,the problem is formulated into a minimization problem with a modified Tikhonov regularization method.The gradient of the regularization functional based on an adjoint problem is deduced and the standard conjugate gradient method is presented for solving the minimization problem.The error estimates for the regularized solutions are obtained under Hp norm priori bound assumptions.Finally,numerical examples illustrate the effectiveness of the proposed method.
基金The project is supported by the National Natural Science Foundation of China(11561045,11961044)the Doctor Fund of Lan Zhou University of Technology.
文摘In this article,we consider to solve the inverse initial value problem for an inhomogeneous space-time fractional diffusion equation.This problem is ill-posed and the quasi-boundary value method is proposed to deal with this inverse problem and obtain the series expression of the regularized solution for the inverse initial value problem.We prove the error estimates between the regularization solution and the exact solution by using an a priori regularization parameter and an a posteriori regularization parameter choice rule.Some numerical results in one-dimensional case and two-dimensional case show that our method is efficient and stable.
基金supported by the National Natural Science Foundation of China(No.10771150)the National Basic Research Program of China(No.2005CB321701)+1 种基金the Program for New Century Excellent Talents in University(No.NCET-07-0584)the Natural Science Foundation of Sichuan Province(No.07ZB087)
文摘A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming （P1）2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.
基金supported by a grant from the National Natural Science Foundation of China (NSFC) No. 81070826/30872886/30400497Sponsored by Shanghai Rising-Star Program No. 09QA1403700+1 种基金funded by Shanghai Leading Academic Discipline Project (Project Number: S30206)the Science and Technology Commission of Shanghai (08DZ2271100)
文摘Aim The purpose of this study was to develop a mathe-matical model to quantitatively describe the passive trans-port of macromolecules within dental biofilms. Methodology Fluorescently labeled dextrans with different molecular mass （3 kD,10 kD,40 kD,70 kD,2 000 kD） were used as a series of diffusion probes. Streptococcus mutans,Streptococcus sanguinis,Actinomyces naeslundii and Fusobacterium nucleatum were used as inocula for biofilm formation. The diffusion processes of different probes through the in vitro biofilm were recorded with a confocal laser microscope. Results Mathematical function of biofilm penetration was constructed on the basis of the inverse problem method. Based on this function,not only the relationship between average concentration of steady-state and molecule weights can be analyzed,but also that between penetrative time and molecule weights. Conclusion This can be used to predict the effective concentration and the penetrative time of anti-biofilm medicines that can diffuse through oral biofilm. Further-more,an improved model for large molecule is proposed by considering the exchange time at the upper boundary of the dental biofilm.
文摘In this paper, the cubic and quintic diffusion equation under stochastic non homogeneity is solved using Wiener- Hermite expansion and perturbation (WHEP) technique, Homotopy perturbation method (HPM) and Pickard approximation technique. The analytic solution of the linear case is obtained using Eigenfunction expansion .The Picard approximation method is used to introduce the first and second order approximate solution for the non linear case. The WHEP technique is also used to obtain approximate solution under different orders and different corrections. The Homotopy perturbation method (HPM) is also used to obtain some approximation orders for mean and variance. Using mathematica-5, the methods of solution are illustrated through figures, comparisons among different methods and some parametric studies.
基金supported by the National Natural Science Foundation of China (No. 10601022)NSF ofInner Mongolia Autonomous Region of China (No. 200607010106)513 and Science Fund of InnerMongolia University for Distinguished Young Scholars (No. ND0702)
文摘A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.
基金supported by the National Natural Science Foundation of China(Nos.52174277 and 51874077)the Fundamental Funds for the Central Universities,China(No.N2225032)+1 种基金the China Postdoctoral Science Foundation(No.2022M720683)the Postdoctoral Fund of Northeastern University,China。
文摘The formation mechanism of calcium vanadate and manganese vanadate and the difference between calcium and manganese in the reaction with vanadium are basic issues in the calcification roasting and manganese roasting process with vanadium slag.In this work,CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) diffusion couples were prepared and roasted for different time periods to illustrate and compare the diffusion reaction mechanisms.Then,the changes in the diffusion product and diffusion coefficient were investigated and calculated based on scanning electron microscopy (SEM) with energy dispersive X-ray spectroscopy (EDS) analysis.Results show that with the extension of the roasting time,the diffusion reaction gradually proceeds among the CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) diffusion couples.The regional boundaries of calcium and vanadium are easily identifiable for the CaO–V_(2)O_(5) diffusion couple.Meanwhile,for the MnO_(2)–V_(2)O_(5) diffusion couple,MnO_(2) gradually decomposes to form Mn_(2)O_(3),and vanadium diffuses into the interior of Mn_(2)O_(3).Only a part of vanadium combines with manganese to form the diffusion production layer.CaV_(2)O_(6) and MnV_(2)O_(6) are the interfacial reaction products of the CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) diffusion couples,respectively,whose thicknesses are 39.85 and 32.13μm when roasted for 16 h.After 16 h,both diffusion couples reach the reaction equilibrium due to the limitation of diffusion.The diffusion coefficient of the CaO–V_(2)O_(5) diffusion couple is higher than that of the MnO_(2)–V_(2)O_(5) diffusion couple for the same roasting time,and the diffusion reaction between vanadium and calcium is easier than that between vanadium and manganese.
基金supported by the National Natural Science Foundation of China(Nos.11271340 and11671369)
文摘This paper aims to present a new streamline diffusion method with low order rectangular Bernardi-Raugel elements to solve the generalized Oseen equations. With the help of the Bramble-Hilbert lemma, the optimal errors of the velocity and pressure are estimated, which are independent of the considered parameter e. With an interpolation postprocessing approach, the superconvergent error of the pressure is obtained. Finally, a numerical experiment is carried out to confirm the theoretical results.
基金funded by the National Natural Science Foundation of China(No.51801116 and 52001176)the Shandong Province Key Research and Development Plan(No.2019GHZ019 and 2021SFGC1001)the Youth Innovation and Technology Support Program of Shandong Provincial Colleges and Universities(No.2020KJA002).
文摘With the rapid development of Mg alloys,deeper understanding to the thermodynamic and diffusional kinetic behavior of intermetallic compounds(IMCs)is important for studying the effect of alloying elements to the microstructure evolution.Specially,a systematic quantitative investigation on the diffusional growth of IMCs is of great necessity.However,the works studying the elemental diffusion behaviors of multiple-element IMCs are rare in magnesium alloy systems.The current work takes the ternary Mg-Al-Zn system as research target,and combines the diffusion couple technique,phase stability diagrams,in-situ observation technique and numerical inverse method to investigate the temperature-dependent kinetic coefficients.The parabolic growth constant(PGC)and interdiffusion coefficients for Mg solid-solution phase andγ-Mg_(17)Al_(12),β-Mg_(2)Al_(3),ε-Mg_(23)Al_(30),MgZn_(2),Mg_(2)Zn_(3),τ-Mg_(32)(Zn,Al)49 andφ-Mg_(5)Zn_(2)Al_(2) IMCs in the Mg-Al-Zn alloy system are determined.By comparing the current experimental with calculation results,the rate-controlling factor of the temperature-dependent diffusion growth ofφ,τandεternary IMCs in the Mg-Al-Zn system is further discussed in detail.
基金the Fundamental Research Funds for the Central Universities(Grants B200203009 and B200202126)the Natural Science Foundation of Jiangsu Province(Grant BK20190073)+2 种基金the State Key Laboratory of Acoustics,Chinese Academy of Sciences(Grant SKLA202001)the State Key Laboratory of Mechanical Behavior and System Safety of Traffic Engineering Structures,Shijiazhuang Tiedao University(Grant KF2020-22)the China Postdoctoral Science Foundation(Grants 2017M611669 and 2018T110430).
文摘A localized space-time method of fundamental solutions(LSTMFS)is extended for solving three-dimensional transient diffusion problems in this paper.The interval segmentation in temporal direction is developed for the accurate simulation of long-time-dependent diffusion problems.In the LSTMFS,the whole space-time domain with nodes arranged i divided into a series of overlapping subdomains with a simple geometry.In each subdomain,the conventional method of fundamental solutions is utilized and the coefficients associated with the considered domain can be explicitly determined.By calculating a combined sparse matrix system,the value at any node inside the space-time domain can be obtained.Numerical experi-ments demonstrate that high accuracy and efficiency can be simultaneously achieved via the LSTMFS,even for the problems defined on a long-time and quite complex computational domain.
