In this study,a straightforward one-step hydrothermal method was successfully utilized to synthesize the solid solution Na_(0.9)Mg_(0.45)Ti_(3.55)O_(8)-Na_(2)Ni_(2)Ti_(6)O_(16)(NNMTO-x),where x denotes the molar perce...In this study,a straightforward one-step hydrothermal method was successfully utilized to synthesize the solid solution Na_(0.9)Mg_(0.45)Ti_(3.55)O_(8)-Na_(2)Ni_(2)Ti_(6)O_(16)(NNMTO-x),where x denotes the molar percentage of Na_(2)Ni_(2)Ti_(6)O_(16)(NNTO)within Na_(0.9)Mg_(0.45)Ti_(3.55)O_(8)(NMTO),with x values of 10,20,30,40,and 50.Both XPS(X-ray Photoelectron Spectroscopy)and EDX(Energy Dispersive X-ray Spectroscopy)analyses unequivocally validated the formation of the NNMTO-x solid solutions.It was observed that when x is below 40,the NNMTO-x solid solution retains the structural characteristics of the original NMTO.However,beyond this threshold,significant alterations in crystal morphology were noted,accompanied by a noticeable decline in photocatalytic activity.Notably,the absorption edge of NNMTO-x(x<40)exhibited a shift towards the visible-light spectrum,thereby substantially broadening the absorption range.The findings highlight that NNMTO-30 possesses the most pronounced photocatalytic activity for the reduction of CO_(2).Specifically,after a 6 h irradiation period,the production rates of CO and CH_(4)were recorded at 42.38 and 1.47μmol/g,respectively.This investigation provides pivotal insights that are instrumental in the advancement of highly efficient and stable photocatalysts tailored for CO_(2)reduction processes.展开更多
In this investigation,a hybrid approach integrating the IDDES turbulence model and FW-H is employed to forecast the hydroacoustic of the rim driven thruster(RDT)under non-cavitation and uniform flow conditions at vary...In this investigation,a hybrid approach integrating the IDDES turbulence model and FW-H is employed to forecast the hydroacoustic of the rim driven thruster(RDT)under non-cavitation and uniform flow conditions at varying loading conditions(J=0.3 and J=0.6).It is revealed that the quadrupole term contribution in the P-FWH method significantly affects the monopole term in the low-frequency region,while it mainly affects the dipole term in the high-frequency region.Specifically,the overall sound pressure levels(SPL)of the RDT using the P-FWH method are 2.27 dB,10.03 dB,and 16.73 dB at the receiving points from R1 to R3 under the heavy-loaded condition,while they increase by 0.67 dB at R1,and decrease by 14.93 dB at R2,and 22.20 dB at R3,for the light-loaded condition.The study also utilizes the pressure-time derivatives to visualize the numerical noise and to pinpoint the dynamics of the vortex cores,and the optimization of the grid design can significantly reduce the numerical noise.The computational accuracy of the P-FWH method can meet the noise requirements for the preliminary design of rim driven thrusters.展开更多
Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The t...Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.展开更多
Rare earth oxide was prepared via direct pyrolysis of rare earth chloride solution. Based on this technique, a new-type jet-flow pyrolysis reactor was designed, and then the fluid dynamics (pressure and velocity) insi...Rare earth oxide was prepared via direct pyrolysis of rare earth chloride solution. Based on this technique, a new-type jet-flow pyrolysis reactor was designed, and then the fluid dynamics (pressure and velocity) inside the reactor was numerically simulated using a computational fluid dynamics method. The self-produced pressure (p) and the fuel inlet velocity (v) satisfied a quadratic function,p=0.06v2+0.23v?4.49. To fully utilize the combustion-generated heat in pyrolysis of rare earth chloride, an appropriate external pressure p=v2+3v?4.27 should be imposed at the feed inlet. The 1.25- and 1.5-fold increase of feed inlet diameter resulted in decline of adsorption dynamic pressure, but the intake of rare earth chloride increased by more than 30% and 60%, respectively. The fluid flow in the reactor was affected by the feeding rate; the fluid flow peaked near the throat of venturi and gradually smoothed down at the jet-flow reactor’s terminal along with the sharp decline of feeding rate.展开更多
A numerical technique of the target-region locating （TRL） solver in conjunction with the wave-front method is presented for the application of the finite element method （FEM） for 3-D electromagnetic computation. F...A numerical technique of the target-region locating （TRL） solver in conjunction with the wave-front method is presented for the application of the finite element method （FEM） for 3-D electromagnetic computation. First, the principle of TRL technique is described. Then, the availability of TRL solver for nonlinear application is particularly discussed demonstrating that this solver can be easily used while still remaining great efficiency. The implementation on how to apply this technique in FEM based on magnetic vector potential （MVP） is also introduced. Finally, a numerical example of 3-D magnetostatic modeling using the TRL solver and FEMLAB is given. It shows that a huge computer resource can be saved by employing the new solver.展开更多
An efficient numerical model for wave refraction, diffraction and reflection is presented in this paper. In the model the modified time-dependent mild-slope equation is transformed into an evolution equation and an im...An efficient numerical model for wave refraction, diffraction and reflection is presented in this paper. In the model the modified time-dependent mild-slope equation is transformed into an evolution equation and an improved ADI method involving a relaxation factor is adopted to solve it. The method has the advantage of improving the numerical stability and convergence rate by properly determining the relaxation factor. The range of the relaxation factor making the differential scheme unconditionally stable is determined by stability analysis. Several verifications are performed to examine the accuracy of the present model. The numerical results coincide with the analytic solutions or experimental data very well and the computer time is reduced.展开更多
The numerical analysis of heat transfer of laminar nanofluid flow over a fiat stretching sheet is presented. Two sets of boundary conditions （BCs） axe analyzed, i.e., a constant （Case 1） and a linear streamwise va...The numerical analysis of heat transfer of laminar nanofluid flow over a fiat stretching sheet is presented. Two sets of boundary conditions （BCs） axe analyzed, i.e., a constant （Case 1） and a linear streamwise variation of nanopaxticle volume fraction and wall temperature （Case 2）. The governing equations and BCs axe reduced to a set of nonlinear ordinary differential equations （ODEs） and the corresponding BCs, respectively. The dependencies of solutions on Prandtl number Pr, Lewis number Le, Brownian motion number Nb, and thermophoresis number Nt are studied in detail. The results show that the reduced Nusselt number and the reduced Sherwood number increase for the BCs of Case 2 compared with Case 1. The increases of Nb, Nt, and Le numbers cause a decrease of the reduced Nusselt number, while the reduced Sherwood number increases with the increase of Nb and Le numbers. For low Prandtl numbers, an increase of Nt number can cause to decrease in the reduced Sherwood number, while it increases for high Prandtl numbers.展开更多
In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarant...In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarantee of computations with a given precision. The equations of programmed constraints and those of constraint perturbations are defined. The stability of the programmed manifold for numerical solutions of the kinematical and dynamical equations is obtained by corresponding construction of the constraint perturbation equations. The dynamical equations of system with programmed constraints are set up in the form of Lagrange’s equations in generalized coordinates. Certain inverse problems of rigid body dynamics are examined.展开更多
The numerical solution for a type of quasilinear wave equation is studied.The three-level difference scheme for quasi-linear waver equation with strong dissipative term is constructed and the convergence is proved.The...The numerical solution for a type of quasilinear wave equation is studied.The three-level difference scheme for quasi-linear waver equation with strong dissipative term is constructed and the convergence is proved.The error of the difference solution is estimated.The theoretical results are controlled on a numerical example.展开更多
This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solutio...This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.展开更多
In this paper, authors discuss the numerical methods of general discontinuous boundary value problems for elliptic complex equations of first order, They first give the well posedness of general discontinuous boundary...In this paper, authors discuss the numerical methods of general discontinuous boundary value problems for elliptic complex equations of first order, They first give the well posedness of general discontinuous boundary value problems, reduce the discontinuous boundary value problems to a variation problem, and then find the numerical solutions of above problem by the finite element method. Finally authors give some error-estimates of the foregoing numerical solutions.展开更多
In this paper, the Adomian decomposition method is developed for the numerical solutions of a class of nonlinear evolution equations with nonlinear term of any order, utt＋auxx ＋ bu ＋ cu^p＋ du^2p-1=0, which contain...In this paper, the Adomian decomposition method is developed for the numerical solutions of a class of nonlinear evolution equations with nonlinear term of any order, utt＋auxx ＋ bu ＋ cu^p＋ du^2p-1=0, which contains some important famous equations. When setting the initial conditions in different forms, some new generalized numerical solutions： numerical hyperbolic solutions, numerical doubly periodic solutions are obtained. The numerical solutions are compared with exact solutions. The scheme is tested by choosing different values of p, positive and negative, integer and fraction, to illustrate the efficiency of the ADM method and the generalization of the solutions.展开更多
In this paper, T-stability of the Euler-Maruyama method is taken into account for linear stochastic delay differential equations with multiplicative noise and constant time lag in the Under a certain condition for coe...In this paper, T-stability of the Euler-Maruyama method is taken into account for linear stochastic delay differential equations with multiplicative noise and constant time lag in the Under a certain condition for coefficients, T-stability of the numerical scheme is researched. Moreover, some numerical examples will be presented to support the theoretical results.展开更多
An initial value problem was considered for a coupled differential system with multi-term Caputo type fractional derivatives. By means of nonlinear alternative of Leray-Schauder and Banach contraction principle,the ex...An initial value problem was considered for a coupled differential system with multi-term Caputo type fractional derivatives. By means of nonlinear alternative of Leray-Schauder and Banach contraction principle,the existence and uniqueness of solutions for the system were derived. Using a fractional predictorcorrector method, a numerical method was presented for the specified system. An example was given to illustrate the obtained results.展开更多
The exothermic efficiency of microwave heating an electrolyte/water solution is remarkably high due to the dielectric heating by orientation polarization of water and resistance heating by the Joule process occurred s...The exothermic efficiency of microwave heating an electrolyte/water solution is remarkably high due to the dielectric heating by orientation polarization of water and resistance heating by the Joule process occurred simultaneously compared with pure water.A three-dimensional finite element numerical model of multi-feed microwave heating industrial liquids continuously flowing in a meter-scale circular tube is presented.The temperature field inside the applicator tube in the cavity is solved by COMSOL Multiphysics and professional programming to describe the momentum,energy and Maxwell's equations.The evaluations of the electromagnetic field,the temperature distribution and the velocity field are simulated for the fluids dynamically heated by singleand multi-feed microwave system,respectively.Both the pilot experimental investigations and numerical results of microwave with single-feed heating for fluids with different effective permittivity and flow rates show that the presented numerical modeling makes it possible to analyze dynamic process of multi-feed microwave heating the industrial liquid.The study aids in enhancing the understanding and optimizing of dynamic process in the use of multi-feed microwave heating industrial continuous flow for a variety of material properties and technical parameters.展开更多
The control volume formulation with the QUICK finite difference scheme is used to solveincompressible liquid flow past a solid sphere in terms of stream function and vorticity.Several tech-nical points are addressed o...The control volume formulation with the QUICK finite difference scheme is used to solveincompressible liquid flow past a solid sphere in terms of stream function and vorticity.Several tech-nical points are addressed on improving the accuracy and efficiency of numerical simulation of similarproblems of fluid flow.In particular,the importance of suitable specification of the distortion func-tion to enforcing the far field boundarv conditions is emphasized.展开更多
The boundary knot method(BKM) is a boundary-type meshfree method. Only non-singular general solutions are used during the whole solution procedures. The effective condition number(ECN), which depends on the right-hand...The boundary knot method(BKM) is a boundary-type meshfree method. Only non-singular general solutions are used during the whole solution procedures. The effective condition number(ECN), which depends on the right-hand side vector of a linear system, is considered as an alternative criterion to the traditional condition number. In this paper, the effective condition number is used to help determine the position and distribution of the collocation points as well as the quasi-optimal collocation point numbers. During the solution process, we propose an NMN-search algorithm. Numerical examples show that the ECN is reliable to measure the feasibility of the BKM.展开更多
Photoinduced birefringence is investigated in a new amorphous copolymer containing azobenzene groups. The levels of birefringence signal are found to depend an the polarization angle between the pump beam and the prob...Photoinduced birefringence is investigated in a new amorphous copolymer containing azobenzene groups. The levels of birefringence signal are found to depend an the polarization angle between the pump beam and the probe beam, and on the ellipticity of the pump beam. Both the growth and decay processes of the birefringence signal can be described by known biexponential equations. The rate constants and the amplitudes associated with the growth process of the photoinduced birefringence are observed to display a linear dependence with the pump beam intensity. A new dynamic model of the photoinduced birefringence is presented taking into account the contributions of both the bans and cis isomers of azobenzene groups and the local polymer segments. The numerical treatment of this model shows good agreement with the experimental data in the whole writing-erasing processes of the photoinduced birefringence conducted in our polymer samples.展开更多
Practical resolution of consolidation problems that we often face requires an extensive and solid knowledge of the different parameters highlighted by the Terzaghi one-dimensional consolidation theory. This theory, wi...Practical resolution of consolidation problems that we often face requires an extensive and solid knowledge of the different parameters highlighted by the Terzaghi one-dimensional consolidation theory. This theory, with its assumptions, leads to a partial differential equation of second order in space and first order in time of pore water pressure. Analytical and numerical resolutions of this equation allow determining the water pressure variation before and after the application of a charge. Numerical modeling has enabled the simulation of the whole results obtained by the two methods of resolution (pressure, degree of consolidation, time factor, among others) to have a physical analysis and a lawful observation that lead to a suitable understanding of the phenomenon of Terzaghi one-dimensional consolidation.展开更多
A numerical solution of the weight function for a two-electrode electromagnetic flowmeter was proposed. The solution was obtained by using the finite element method based on the basic equation of a traditional two-ele...A numerical solution of the weight function for a two-electrode electromagnetic flowmeter was proposed. The solution was obtained by using the finite element method based on the basic equation of a traditional two-electrode electromagnetic flowmeter. The two-dimensional distribution of the weight function of the electromagnetic flowmeter obtained was verified by the analytical solution. Three-dimensional distribution of the weight function was also presented in the paper. It can be employed to analyze the sensitivity and linearity of the electromagnetic flowmeter with non-uniform magnetic field, and even to assist the design of the excitation coil pair.展开更多
基金Supported by the Doctoral Research Start-up Project of Yuncheng University(YQ-2023067)Project of Shanxi Natural Science Foundation(202303021211189)+1 种基金Fund Program for the Scientific Activities of Selected Returned Overseas Professionals in Shanxi Provinces(20220036)Shanxi ProvinceIntelligent Optoelectronic Sensing Application Technology Innovation Center and Shanxi Province Optoelectronic Information Science and TechnologyLaboratory,Yuncheng University.
文摘In this study,a straightforward one-step hydrothermal method was successfully utilized to synthesize the solid solution Na_(0.9)Mg_(0.45)Ti_(3.55)O_(8)-Na_(2)Ni_(2)Ti_(6)O_(16)(NNMTO-x),where x denotes the molar percentage of Na_(2)Ni_(2)Ti_(6)O_(16)(NNTO)within Na_(0.9)Mg_(0.45)Ti_(3.55)O_(8)(NMTO),with x values of 10,20,30,40,and 50.Both XPS(X-ray Photoelectron Spectroscopy)and EDX(Energy Dispersive X-ray Spectroscopy)analyses unequivocally validated the formation of the NNMTO-x solid solutions.It was observed that when x is below 40,the NNMTO-x solid solution retains the structural characteristics of the original NMTO.However,beyond this threshold,significant alterations in crystal morphology were noted,accompanied by a noticeable decline in photocatalytic activity.Notably,the absorption edge of NNMTO-x(x<40)exhibited a shift towards the visible-light spectrum,thereby substantially broadening the absorption range.The findings highlight that NNMTO-30 possesses the most pronounced photocatalytic activity for the reduction of CO_(2).Specifically,after a 6 h irradiation period,the production rates of CO and CH_(4)were recorded at 42.38 and 1.47μmol/g,respectively.This investigation provides pivotal insights that are instrumental in the advancement of highly efficient and stable photocatalysts tailored for CO_(2)reduction processes.
基金The National Natural Science Foundation of China(Grant No.52201376)the Natural Science Foundation of Hubei Province,China(Grant No.2023AFB683).
文摘In this investigation,a hybrid approach integrating the IDDES turbulence model and FW-H is employed to forecast the hydroacoustic of the rim driven thruster(RDT)under non-cavitation and uniform flow conditions at varying loading conditions(J=0.3 and J=0.6).It is revealed that the quadrupole term contribution in the P-FWH method significantly affects the monopole term in the low-frequency region,while it mainly affects the dipole term in the high-frequency region.Specifically,the overall sound pressure levels(SPL)of the RDT using the P-FWH method are 2.27 dB,10.03 dB,and 16.73 dB at the receiving points from R1 to R3 under the heavy-loaded condition,while they increase by 0.67 dB at R1,and decrease by 14.93 dB at R2,and 22.20 dB at R3,for the light-loaded condition.The study also utilizes the pressure-time derivatives to visualize the numerical noise and to pinpoint the dynamics of the vortex cores,and the optimization of the grid design can significantly reduce the numerical noise.The computational accuracy of the P-FWH method can meet the noise requirements for the preliminary design of rim driven thrusters.
基金Supported by the National Natural Science Foundation of China under Grant No.51975138the High-Tech Ship Scientific Research Project from the Ministry of Industry and Information Technology under Grant No.CJ05N20the National Defense Basic Research Project under Grant No.JCKY2023604C006.
文摘Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.
基金Projects(51204040,U1202274)supported by the National Natural Science Foundation of ChinaProjects(2010AA03A405,2102AA062303)supported by the National High-tech Research and Development Program of China+1 种基金Project(2012BAE01B02)supported by the National Science and Technology Support Program of ChinaProject(N130702001)supported by the Fundamental Research Funds for the Central Universities,China
文摘Rare earth oxide was prepared via direct pyrolysis of rare earth chloride solution. Based on this technique, a new-type jet-flow pyrolysis reactor was designed, and then the fluid dynamics (pressure and velocity) inside the reactor was numerically simulated using a computational fluid dynamics method. The self-produced pressure (p) and the fuel inlet velocity (v) satisfied a quadratic function,p=0.06v2+0.23v?4.49. To fully utilize the combustion-generated heat in pyrolysis of rare earth chloride, an appropriate external pressure p=v2+3v?4.27 should be imposed at the feed inlet. The 1.25- and 1.5-fold increase of feed inlet diameter resulted in decline of adsorption dynamic pressure, but the intake of rare earth chloride increased by more than 30% and 60%, respectively. The fluid flow in the reactor was affected by the feeding rate; the fluid flow peaked near the throat of venturi and gradually smoothed down at the jet-flow reactor’s terminal along with the sharp decline of feeding rate.
基金Open Funds of State Key Laboratory of MillimeterWaves,China (No. K200401), Outstanding Teaching and ResearchAwards for Young Teachers of Nanjing Normal University (No.1320BL51)
文摘A numerical technique of the target-region locating （TRL） solver in conjunction with the wave-front method is presented for the application of the finite element method （FEM） for 3-D electromagnetic computation. First, the principle of TRL technique is described. Then, the availability of TRL solver for nonlinear application is particularly discussed demonstrating that this solver can be easily used while still remaining great efficiency. The implementation on how to apply this technique in FEM based on magnetic vector potential （MVP） is also introduced. Finally, a numerical example of 3-D magnetostatic modeling using the TRL solver and FEMLAB is given. It shows that a huge computer resource can be saved by employing the new solver.
文摘An efficient numerical model for wave refraction, diffraction and reflection is presented in this paper. In the model the modified time-dependent mild-slope equation is transformed into an evolution equation and an improved ADI method involving a relaxation factor is adopted to solve it. The method has the advantage of improving the numerical stability and convergence rate by properly determining the relaxation factor. The range of the relaxation factor making the differential scheme unconditionally stable is determined by stability analysis. Several verifications are performed to examine the accuracy of the present model. The numerical results coincide with the analytic solutions or experimental data very well and the computer time is reduced.
文摘The numerical analysis of heat transfer of laminar nanofluid flow over a fiat stretching sheet is presented. Two sets of boundary conditions （BCs） axe analyzed, i.e., a constant （Case 1） and a linear streamwise variation of nanopaxticle volume fraction and wall temperature （Case 2）. The governing equations and BCs axe reduced to a set of nonlinear ordinary differential equations （ODEs） and the corresponding BCs, respectively. The dependencies of solutions on Prandtl number Pr, Lewis number Le, Brownian motion number Nb, and thermophoresis number Nt are studied in detail. The results show that the reduced Nusselt number and the reduced Sherwood number increase for the BCs of Case 2 compared with Case 1. The increases of Nb, Nt, and Le numbers cause a decrease of the reduced Nusselt number, while the reduced Sherwood number increases with the increase of Nb and Le numbers. For low Prandtl numbers, an increase of Nt number can cause to decrease in the reduced Sherwood number, while it increases for high Prandtl numbers.
基金Supported by Russian Fund of Fund amental Investigations(Pr.990101064)and Russian Minister of Educatin
文摘In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarantee of computations with a given precision. The equations of programmed constraints and those of constraint perturbations are defined. The stability of the programmed manifold for numerical solutions of the kinematical and dynamical equations is obtained by corresponding construction of the constraint perturbation equations. The dynamical equations of system with programmed constraints are set up in the form of Lagrange’s equations in generalized coordinates. Certain inverse problems of rigid body dynamics are examined.
文摘The numerical solution for a type of quasilinear wave equation is studied.The three-level difference scheme for quasi-linear waver equation with strong dissipative term is constructed and the convergence is proved.The error of the difference solution is estimated.The theoretical results are controlled on a numerical example.
文摘This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.
文摘In this paper, authors discuss the numerical methods of general discontinuous boundary value problems for elliptic complex equations of first order, They first give the well posedness of general discontinuous boundary value problems, reduce the discontinuous boundary value problems to a variation problem, and then find the numerical solutions of above problem by the finite element method. Finally authors give some error-estimates of the foregoing numerical solutions.
基金supported by National Natural Science Foundation of China under Grant No.10735030Natural Science Foundation of Zhejiang Province of China under Grant No.Y604056Doctoral Science Foundation of Ningbo City under Grant No.2005A61030
文摘In this paper, the Adomian decomposition method is developed for the numerical solutions of a class of nonlinear evolution equations with nonlinear term of any order, utt＋auxx ＋ bu ＋ cu^p＋ du^2p-1=0, which contains some important famous equations. When setting the initial conditions in different forms, some new generalized numerical solutions： numerical hyperbolic solutions, numerical doubly periodic solutions are obtained. The numerical solutions are compared with exact solutions. The scheme is tested by choosing different values of p, positive and negative, integer and fraction, to illustrate the efficiency of the ADM method and the generalization of the solutions.
基金Supported by the National Natural Science Foundation of China(51008084)the Natural Science Foundation of Guangdong Province(9451009001002753)
文摘In this paper, T-stability of the Euler-Maruyama method is taken into account for linear stochastic delay differential equations with multiplicative noise and constant time lag in the Under a certain condition for coefficients, T-stability of the numerical scheme is researched. Moreover, some numerical examples will be presented to support the theoretical results.
基金National Natural Science Foundation of China(No.11371087)
文摘An initial value problem was considered for a coupled differential system with multi-term Caputo type fractional derivatives. By means of nonlinear alternative of Leray-Schauder and Banach contraction principle,the existence and uniqueness of solutions for the system were derived. Using a fractional predictorcorrector method, a numerical method was presented for the specified system. An example was given to illustrate the obtained results.
基金Project(KKSY201503006)supported by Scientific Research Foundation of Kunming University of Science and Technology,ChinaProject(2014FD009)supported by the Applied Basic Research Foundation(Youth Program)of ChinaProject(51090385)supported by the National Natural Science Foundation of China
文摘The exothermic efficiency of microwave heating an electrolyte/water solution is remarkably high due to the dielectric heating by orientation polarization of water and resistance heating by the Joule process occurred simultaneously compared with pure water.A three-dimensional finite element numerical model of multi-feed microwave heating industrial liquids continuously flowing in a meter-scale circular tube is presented.The temperature field inside the applicator tube in the cavity is solved by COMSOL Multiphysics and professional programming to describe the momentum,energy and Maxwell's equations.The evaluations of the electromagnetic field,the temperature distribution and the velocity field are simulated for the fluids dynamically heated by singleand multi-feed microwave system,respectively.Both the pilot experimental investigations and numerical results of microwave with single-feed heating for fluids with different effective permittivity and flow rates show that the presented numerical modeling makes it possible to analyze dynamic process of multi-feed microwave heating the industrial liquid.The study aids in enhancing the understanding and optimizing of dynamic process in the use of multi-feed microwave heating industrial continuous flow for a variety of material properties and technical parameters.
基金Supported by the National Natural Science Foundation of China.
文摘The control volume formulation with the QUICK finite difference scheme is used to solveincompressible liquid flow past a solid sphere in terms of stream function and vorticity.Several tech-nical points are addressed on improving the accuracy and efficiency of numerical simulation of similarproblems of fluid flow.In particular,the importance of suitable specification of the distortion func-tion to enforcing the far field boundarv conditions is emphasized.
基金Supported by the Natural Science Foundation of Anhui Province(1908085QA09)Higher Education Department of the Ministry of Education(201802358008)
文摘The boundary knot method(BKM) is a boundary-type meshfree method. Only non-singular general solutions are used during the whole solution procedures. The effective condition number(ECN), which depends on the right-hand side vector of a linear system, is considered as an alternative criterion to the traditional condition number. In this paper, the effective condition number is used to help determine the position and distribution of the collocation points as well as the quasi-optimal collocation point numbers. During the solution process, we propose an NMN-search algorithm. Numerical examples show that the ECN is reliable to measure the feasibility of the BKM.
基金This work was supported by the National Natural Science Foundation of China (No. 19604015) and Guangdong Provincial Natural Science Foundation (No. 980279 and 980346).
文摘Photoinduced birefringence is investigated in a new amorphous copolymer containing azobenzene groups. The levels of birefringence signal are found to depend an the polarization angle between the pump beam and the probe beam, and on the ellipticity of the pump beam. Both the growth and decay processes of the birefringence signal can be described by known biexponential equations. The rate constants and the amplitudes associated with the growth process of the photoinduced birefringence are observed to display a linear dependence with the pump beam intensity. A new dynamic model of the photoinduced birefringence is presented taking into account the contributions of both the bans and cis isomers of azobenzene groups and the local polymer segments. The numerical treatment of this model shows good agreement with the experimental data in the whole writing-erasing processes of the photoinduced birefringence conducted in our polymer samples.
文摘Practical resolution of consolidation problems that we often face requires an extensive and solid knowledge of the different parameters highlighted by the Terzaghi one-dimensional consolidation theory. This theory, with its assumptions, leads to a partial differential equation of second order in space and first order in time of pore water pressure. Analytical and numerical resolutions of this equation allow determining the water pressure variation before and after the application of a charge. Numerical modeling has enabled the simulation of the whole results obtained by the two methods of resolution (pressure, degree of consolidation, time factor, among others) to have a physical analysis and a lawful observation that lead to a suitable understanding of the phenomenon of Terzaghi one-dimensional consolidation.
基金Supported by Chinese National Programs for High Technology Research and Development(2008AA042207)
文摘A numerical solution of the weight function for a two-electrode electromagnetic flowmeter was proposed. The solution was obtained by using the finite element method based on the basic equation of a traditional two-electrode electromagnetic flowmeter. The two-dimensional distribution of the weight function of the electromagnetic flowmeter obtained was verified by the analytical solution. Three-dimensional distribution of the weight function was also presented in the paper. It can be employed to analyze the sensitivity and linearity of the electromagnetic flowmeter with non-uniform magnetic field, and even to assist the design of the excitation coil pair.
