This paper concerns the implementation of the orthogonal polynomials using the Galerkin method for solving Volterra integro-differential and Fredholm integro-differential equations. The constructed orthogonal polynomi...This paper concerns the implementation of the orthogonal polynomials using the Galerkin method for solving Volterra integro-differential and Fredholm integro-differential equations. The constructed orthogonal polynomials are used as basis functions in the assumed solution employed. Numerical examples for some selected problems are provided and the results obtained show that the Galerkin method with orthogonal polynomials as basis functions performed creditably well in terms of absolute errors obtained.展开更多
Sodium silicate modification can improve the overall performance of wood.The modification process has a great influence on the properties of modified wood.In this study,a new method was introduced to analyze the wood ...Sodium silicate modification can improve the overall performance of wood.The modification process has a great influence on the properties of modified wood.In this study,a new method was introduced to analyze the wood modification process,and the properties of modified wood were studied.Poplar wood was modified with sodium silicate by vacuum-pressure impregnation.After screening using single-factor experiments,an orthogonal experiment was carried out with solution concentration,impregnation time,impregnation pressure,and the cycle times as experimental factors.The modified poplar with the best properties was selected by fuzzy mathematics and characterized by SEM,FT-IR,XRD and TG.The results showed that some lignin and hemicellulose were removed from the wood due to the alkaline action of sodium silicate,and the orderly crystal area of poplar became disorderly,resulting in the reduction of crystallinity of the modified poplar wood.FT-IR analysis showed that sodium silicate was hydrolyzed to form polysilicic acid in wood,and structural analysis revealed the formation of Si-O-Si and Si-O-C,indicating that sodium silicate reacted with fibers on the wood cell wall.TG-DTG curves showed that the final residual mass of modified poplar wood increased from 25%to 67%,and the temperature of the maximum loss rate decreased from 343℃ to 276℃.The heat release and smoke release of modified poplar wood decreased obviously.This kind of material with high strength and fire resistance can be used in the outdoor building and indoor furniture.展开更多
Finite element method (FEM) is an efficient numerical tool for the solution of partial differential equations (PDEs). It is one of the most general methods when compared to other numerical techniques. PDEs posed in a ...Finite element method (FEM) is an efficient numerical tool for the solution of partial differential equations (PDEs). It is one of the most general methods when compared to other numerical techniques. PDEs posed in a variational form over a given space, say a Hilbert space, are better numerically handled with the FEM. The FEM algorithm is used in various applications which includes fluid flow, heat transfer, acoustics, structural mechanics and dynamics, electric and magnetic field, etc. Thus, in this paper, the Finite Element Orthogonal Collocation Approach (FEOCA) is established for the approximate solution of Time Fractional Telegraph Equation (TFTE) with Mamadu-Njoseh polynomials as grid points corresponding to new basis functions constructed in the finite element space. The FEOCA is an elegant mixture of the Finite Element Method (FEM) and the Orthogonal Collocation Method (OCM). Two numerical examples are experimented on to verify the accuracy and rate of convergence of the method as compared with the theoretical results, and other methods in literature.展开更多
The proper orthogonal decomposition (POD) method for the instationary Navier-Stokes equations is considered. Several numerical approaches to evaluating the POD eigenfunctions are presented. The POD eigenfunctions are ...The proper orthogonal decomposition (POD) method for the instationary Navier-Stokes equations is considered. Several numerical approaches to evaluating the POD eigenfunctions are presented. The POD eigenfunctions are applied as a basis for a Galerkin projection of the instationary Navier-Stokes equations. And a low-dimensional ordinary differential models for fluid flows governed by the instationary Navier-Stokes equations are constructed. The numerical examples show that the method is feasible and efficient for optimal control of fluids.展开更多
This work is aim at providing a numerical technique for the Volterra integral equations using Galerkin method. For this purpose, an effective matrix formulation is proposed to solve linear Volterra integral equations ...This work is aim at providing a numerical technique for the Volterra integral equations using Galerkin method. For this purpose, an effective matrix formulation is proposed to solve linear Volterra integral equations of the first and second kind respectively using orthogonal polynomials as trial functions which are constructed in the interval [-1,1] with respect to the weight function w(x)=1+x<sup>2</sup>. The efficiency of the proposed method is tested on several numerical examples and compared with the analytic solutions available in the literature.展开更多
In this paper, the fully discrete orthogonal collocation method for Sobolev equations is considered, and the equivalence for discrete Garlerkin method is proved. Optimal order error estimate is obtained.
In this Paper, we have proposed a new weighted residual method known as orthogonal collocation-based on mixed interpolation (OCMI). Mixed interpolation uses the classical polynomial approximation with two correction t...In this Paper, we have proposed a new weighted residual method known as orthogonal collocation-based on mixed interpolation (OCMI). Mixed interpolation uses the classical polynomial approximation with two correction terms given in the form of sine and cosine function. By these correction terms, we can control the error in the solution. We have applied this approach to a non-linear boundary value problem (BVP) in ODE which governs the electrohydrodynamic flow in a cylindrical conduit. The solution profiles shown in the figures are in good agreement with the work of Paullet (1999) and Ghasemi et al. (2014). Our solution is monotonic decreasing and satisfies , where, α governs the strength of non-linearity and for large values of α solutions are . The residual errors are given in Table 1 and Table 2 which are significantly small. Comparison of residual errors between our proposed method, Least square method and Homotopy analysis method is also given and shown via the Table 3 where as the profiles of the residual error are depicted in Figures 4-8. Table and graphs show that efficiency of the proposed method. The error bound and its L2-norm with relevant theorems for mixed interpolation are also given.展开更多
The inverse heat conduction method is one of methods to identify the casting simulation parameters. A new inverse method was presented according to the Tikhonov regularization theory. One appropriate regularized funct...The inverse heat conduction method is one of methods to identify the casting simulation parameters. A new inverse method was presented according to the Tikhonov regularization theory. One appropriate regularized functional was established, and the functional was solved by the sensitivity coefficient and Newtonaphson iteration method. Moreover, the orthogonal experimental design was used to estimate the appropriate initial value and variation domain of each variable to decrease the number of iteration and improve the identification accuracy and efficiency. It illustrated a detailed case of AlSiTMg sand mold casting and the temperature measurement experiment was done. The physical properties of sand mold and the interracial heat transfer coefficient were identified at the meantime. The results indicated that the new regularization method was efficient in overcoming the ill-posedness of the inverse heat conduction problem and improving the stability and accuracy of the solutions.展开更多
The influence and signification of casting parameters on the solidification process of steel ingot were discussed based on the finite element method (FEM) results by orthogonal experiment method. The range analysis, a...The influence and signification of casting parameters on the solidification process of steel ingot were discussed based on the finite element method (FEM) results by orthogonal experiment method. The range analysis, analysis of variance (ANOVA) and optimization project were used to investigate the FEM results. In order to decrease the ingot riser head and improve the utilization ratio of ingot, the casting parameters involved casting temperature, pouring velocity and interface heat transfer were optimized to decrease shrinkage pore and microporosity. The results show that the heat transfer coefficient between melt and heated board is a more sensitive factor. It is favor to decrease the shrinkage pore and microporosity under the conditions of low temperature, high pouring velocity and high heat transfer between melt and mold. If heat transfer in the ingot body is quicker than that in the riser, the position of shrinkage pore and microporosity will be closer to riser top. The results of optimization project show that few of shrinkage pore and microporosity reach into ingot body with the rational parameters, so the riser size can be reduced.展开更多
Due to the current trend towards lightweight design in automotive industry,hollow stepped gear shafts for automobile and its radial forging process are widely investigated.Utilizing coupled finite element thermo-mecha...Due to the current trend towards lightweight design in automotive industry,hollow stepped gear shafts for automobile and its radial forging process are widely investigated.Utilizing coupled finite element thermo-mechanical model,radial forging process of a hollow stepped gear shaft for automobile was simulated.The optimal combination of three process parameters including initial temperature,rotation rate and radial reduction was also selected using orthogonal design method.To examine the strain inhomogeneity of the forging workpiece,the strain inhomogeneity factor was introduced.The results reveal that the maximum effective strain and the minimum effective strain appeared in the outermost and innermost zones of different cross sections for the hollow stepped gear shaft,respectively.Optimal forging parameters are determined as a combination of initial temperature of 780°C,rotation rate of 21°/stroke and radial reduction of 3 mm.展开更多
For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eig...For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eigenfunctions) of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors (SIFs), so any order SIF can be extracted based on a certain known solution of displacement (an analytic result or a numerical result). Many numerical examples based on the finite element method of lines (FEMOL) show that the present method is very powerful and efficient.展开更多
Accurate real-time simulations of nuclear reactor circuit systems are particularly important for system safety analysis and design.To effectively improve computational efficiency without reducing accuracy,this study e...Accurate real-time simulations of nuclear reactor circuit systems are particularly important for system safety analysis and design.To effectively improve computational efficiency without reducing accuracy,this study establishes a thermal-hydraulics reduced-order model(ROM)for nuclear reactor circuit systems.The full-order circuit system calculation model is first established and verified and then used to calculate the thermal-hydraulic properties of the circuit system under different states as snapshots.The proper orthogonal decomposition method is used to extract the basis functions from snapshots,and the ROM is constructed using the least-squares method,effectively reducing the difficulty in constructing the ROM.A comparison between the full-order simulation and ROM prediction results of the AP1000 circuit system shows that the proposed ROM can improve computational efficiency by 1500 times while achieving a maximum relative error of 0.223%.This research develops a new direction and perspective for the digital twin modeling of nuclear reactor system circuits.展开更多
The effect of temperature and metal ion components on the viscosity of HPAM solution was studied by means of orthogonal method.Five factors and 4 levels were considered for the orthogonal design.The five factors inclu...The effect of temperature and metal ion components on the viscosity of HPAM solution was studied by means of orthogonal method.Five factors and 4 levels were considered for the orthogonal design.The five factors included temperature,contents of Na+,Mg2+,Ca2+ and S2-,and the four levels mainly considered the level of produced water components in main oilfields in China.The experiment results show that temperature is the most important factor to control HPAM solution viscosity prepared by produced water.The effect of ions was in order of Na+>S2->Mg2+>Ca2+.The mechanism of each factor influencing viscosity was also discussed.Based on the actual condition of Nanyang oilfield,the desulfurization was used to improve the solution viscosity,and satisfactory result was obtained.展开更多
It is necessary to reduce the currents of poloidal field(PF) coils as small as possible, during the static equilibrium design procedure of Experimental Advanced Superconductive Tokamak(EAST). The quasi-snowflake(QSF) ...It is necessary to reduce the currents of poloidal field(PF) coils as small as possible, during the static equilibrium design procedure of Experimental Advanced Superconductive Tokamak(EAST). The quasi-snowflake(QSF) divertor configuration is studied in this paper. Starting from a standard QSF plasma equilibrium, a new QSF equilibrium with 300 kA total plasma current is designed. In order to reduce the currents of PF6 and PF14, the influence of plasma shape on PF coil current distribution is analyzed. A fixed boundary equilibrium solver based on a non-rigid plasma model is used to calculate the flux distribution and PF coil current distribution. Then the plasma shape parameters are studied by the orthogonal method. According to the result, the plasma shape is redefined, and the calculated equilibrium shows that the currents of PF6 and PF14 are reduced by 3.592 kA and 2.773 kA, respectively.展开更多
文摘This paper concerns the implementation of the orthogonal polynomials using the Galerkin method for solving Volterra integro-differential and Fredholm integro-differential equations. The constructed orthogonal polynomials are used as basis functions in the assumed solution employed. Numerical examples for some selected problems are provided and the results obtained show that the Galerkin method with orthogonal polynomials as basis functions performed creditably well in terms of absolute errors obtained.
基金This work was financially supported by National Natural Science Foundation of China(32201485)Natural Science Foundation of Hunan Province,China(2022JJ40863)+1 种基金Scientific Research Project of Hunan Provincial Education Department,China(21B0238)The Science and Technology Innovation Program of Hunan Province(2021RC4062).
文摘Sodium silicate modification can improve the overall performance of wood.The modification process has a great influence on the properties of modified wood.In this study,a new method was introduced to analyze the wood modification process,and the properties of modified wood were studied.Poplar wood was modified with sodium silicate by vacuum-pressure impregnation.After screening using single-factor experiments,an orthogonal experiment was carried out with solution concentration,impregnation time,impregnation pressure,and the cycle times as experimental factors.The modified poplar with the best properties was selected by fuzzy mathematics and characterized by SEM,FT-IR,XRD and TG.The results showed that some lignin and hemicellulose were removed from the wood due to the alkaline action of sodium silicate,and the orderly crystal area of poplar became disorderly,resulting in the reduction of crystallinity of the modified poplar wood.FT-IR analysis showed that sodium silicate was hydrolyzed to form polysilicic acid in wood,and structural analysis revealed the formation of Si-O-Si and Si-O-C,indicating that sodium silicate reacted with fibers on the wood cell wall.TG-DTG curves showed that the final residual mass of modified poplar wood increased from 25%to 67%,and the temperature of the maximum loss rate decreased from 343℃ to 276℃.The heat release and smoke release of modified poplar wood decreased obviously.This kind of material with high strength and fire resistance can be used in the outdoor building and indoor furniture.
文摘Finite element method (FEM) is an efficient numerical tool for the solution of partial differential equations (PDEs). It is one of the most general methods when compared to other numerical techniques. PDEs posed in a variational form over a given space, say a Hilbert space, are better numerically handled with the FEM. The FEM algorithm is used in various applications which includes fluid flow, heat transfer, acoustics, structural mechanics and dynamics, electric and magnetic field, etc. Thus, in this paper, the Finite Element Orthogonal Collocation Approach (FEOCA) is established for the approximate solution of Time Fractional Telegraph Equation (TFTE) with Mamadu-Njoseh polynomials as grid points corresponding to new basis functions constructed in the finite element space. The FEOCA is an elegant mixture of the Finite Element Method (FEM) and the Orthogonal Collocation Method (OCM). Two numerical examples are experimented on to verify the accuracy and rate of convergence of the method as compared with the theoretical results, and other methods in literature.
基金National Natural Science Foundation of China (No.10671153)
文摘The proper orthogonal decomposition (POD) method for the instationary Navier-Stokes equations is considered. Several numerical approaches to evaluating the POD eigenfunctions are presented. The POD eigenfunctions are applied as a basis for a Galerkin projection of the instationary Navier-Stokes equations. And a low-dimensional ordinary differential models for fluid flows governed by the instationary Navier-Stokes equations are constructed. The numerical examples show that the method is feasible and efficient for optimal control of fluids.
文摘This work is aim at providing a numerical technique for the Volterra integral equations using Galerkin method. For this purpose, an effective matrix formulation is proposed to solve linear Volterra integral equations of the first and second kind respectively using orthogonal polynomials as trial functions which are constructed in the interval [-1,1] with respect to the weight function w(x)=1+x<sup>2</sup>. The efficiency of the proposed method is tested on several numerical examples and compared with the analytic solutions available in the literature.
文摘In this paper, the fully discrete orthogonal collocation method for Sobolev equations is considered, and the equivalence for discrete Garlerkin method is proved. Optimal order error estimate is obtained.
文摘In this Paper, we have proposed a new weighted residual method known as orthogonal collocation-based on mixed interpolation (OCMI). Mixed interpolation uses the classical polynomial approximation with two correction terms given in the form of sine and cosine function. By these correction terms, we can control the error in the solution. We have applied this approach to a non-linear boundary value problem (BVP) in ODE which governs the electrohydrodynamic flow in a cylindrical conduit. The solution profiles shown in the figures are in good agreement with the work of Paullet (1999) and Ghasemi et al. (2014). Our solution is monotonic decreasing and satisfies , where, α governs the strength of non-linearity and for large values of α solutions are . The residual errors are given in Table 1 and Table 2 which are significantly small. Comparison of residual errors between our proposed method, Least square method and Homotopy analysis method is also given and shown via the Table 3 where as the profiles of the residual error are depicted in Figures 4-8. Table and graphs show that efficiency of the proposed method. The error bound and its L2-norm with relevant theorems for mixed interpolation are also given.
文摘The inverse heat conduction method is one of methods to identify the casting simulation parameters. A new inverse method was presented according to the Tikhonov regularization theory. One appropriate regularized functional was established, and the functional was solved by the sensitivity coefficient and Newtonaphson iteration method. Moreover, the orthogonal experimental design was used to estimate the appropriate initial value and variation domain of each variable to decrease the number of iteration and improve the identification accuracy and efficiency. It illustrated a detailed case of AlSiTMg sand mold casting and the temperature measurement experiment was done. The physical properties of sand mold and the interracial heat transfer coefficient were identified at the meantime. The results indicated that the new regularization method was efficient in overcoming the ill-posedness of the inverse heat conduction problem and improving the stability and accuracy of the solutions.
基金Projects(50435010 50705080 50675187) supported by the National Natural Science Foundation of China
文摘The influence and signification of casting parameters on the solidification process of steel ingot were discussed based on the finite element method (FEM) results by orthogonal experiment method. The range analysis, analysis of variance (ANOVA) and optimization project were used to investigate the FEM results. In order to decrease the ingot riser head and improve the utilization ratio of ingot, the casting parameters involved casting temperature, pouring velocity and interface heat transfer were optimized to decrease shrinkage pore and microporosity. The results show that the heat transfer coefficient between melt and heated board is a more sensitive factor. It is favor to decrease the shrinkage pore and microporosity under the conditions of low temperature, high pouring velocity and high heat transfer between melt and mold. If heat transfer in the ingot body is quicker than that in the riser, the position of shrinkage pore and microporosity will be closer to riser top. The results of optimization project show that few of shrinkage pore and microporosity reach into ingot body with the rational parameters, so the riser size can be reduced.
基金Projects(51774054,51974050)supported by the National Natural Science Foundation of China。
文摘Due to the current trend towards lightweight design in automotive industry,hollow stepped gear shafts for automobile and its radial forging process are widely investigated.Utilizing coupled finite element thermo-mechanical model,radial forging process of a hollow stepped gear shaft for automobile was simulated.The optimal combination of three process parameters including initial temperature,rotation rate and radial reduction was also selected using orthogonal design method.To examine the strain inhomogeneity of the forging workpiece,the strain inhomogeneity factor was introduced.The results reveal that the maximum effective strain and the minimum effective strain appeared in the outermost and innermost zones of different cross sections for the hollow stepped gear shaft,respectively.Optimal forging parameters are determined as a combination of initial temperature of 780°C,rotation rate of 21°/stroke and radial reduction of 3 mm.
基金Project supported by the National Natural Science Foundation of China (Nos. 59525813 and 19872066).
文摘For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eigenfunctions) of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors (SIFs), so any order SIF can be extracted based on a certain known solution of displacement (an analytic result or a numerical result). Many numerical examples based on the finite element method of lines (FEMOL) show that the present method is very powerful and efficient.
基金supported by the National Natural Science Foundation of China(No.12205389)Guangdong Basic and Applied Basic Research Foundation(No.2022A1515011735)Science and Technology on Reactor System Design Technology Laboratory(No.KFKT-05-FWHT-WU-2023014).
文摘Accurate real-time simulations of nuclear reactor circuit systems are particularly important for system safety analysis and design.To effectively improve computational efficiency without reducing accuracy,this study establishes a thermal-hydraulics reduced-order model(ROM)for nuclear reactor circuit systems.The full-order circuit system calculation model is first established and verified and then used to calculate the thermal-hydraulic properties of the circuit system under different states as snapshots.The proper orthogonal decomposition method is used to extract the basis functions from snapshots,and the ROM is constructed using the least-squares method,effectively reducing the difficulty in constructing the ROM.A comparison between the full-order simulation and ROM prediction results of the AP1000 circuit system shows that the proposed ROM can improve computational efficiency by 1500 times while achieving a maximum relative error of 0.223%.This research develops a new direction and perspective for the digital twin modeling of nuclear reactor system circuits.
基金Project(2007AA06Z214) supported by the National High-Tech Research and Development Program of ChinaProject(ts20070704) supported by Taishan Scholars Construction EngineeringProject(200873181) supported by National Natural Science Foundation of China
文摘The effect of temperature and metal ion components on the viscosity of HPAM solution was studied by means of orthogonal method.Five factors and 4 levels were considered for the orthogonal design.The five factors included temperature,contents of Na+,Mg2+,Ca2+ and S2-,and the four levels mainly considered the level of produced water components in main oilfields in China.The experiment results show that temperature is the most important factor to control HPAM solution viscosity prepared by produced water.The effect of ions was in order of Na+>S2->Mg2+>Ca2+.The mechanism of each factor influencing viscosity was also discussed.Based on the actual condition of Nanyang oilfield,the desulfurization was used to improve the solution viscosity,and satisfactory result was obtained.
基金Project supported by the National Natural Science Foundation of China(Grant No.51677051)the Institute of Plasma Physics,Chinese Academy of Sciences
文摘It is necessary to reduce the currents of poloidal field(PF) coils as small as possible, during the static equilibrium design procedure of Experimental Advanced Superconductive Tokamak(EAST). The quasi-snowflake(QSF) divertor configuration is studied in this paper. Starting from a standard QSF plasma equilibrium, a new QSF equilibrium with 300 kA total plasma current is designed. In order to reduce the currents of PF6 and PF14, the influence of plasma shape on PF coil current distribution is analyzed. A fixed boundary equilibrium solver based on a non-rigid plasma model is used to calculate the flux distribution and PF coil current distribution. Then the plasma shape parameters are studied by the orthogonal method. According to the result, the plasma shape is redefined, and the calculated equilibrium shows that the currents of PF6 and PF14 are reduced by 3.592 kA and 2.773 kA, respectively.