Background: Studies have shown a strong correlation between the growth of E2 in serum and estrone-3-glucuronide (E1-3G) in urine during ovarian stimulation. Thus, we developed theoretical models for using urinary E1-3...Background: Studies have shown a strong correlation between the growth of E2 in serum and estrone-3-glucuronide (E1-3G) in urine during ovarian stimulation. Thus, we developed theoretical models for using urinary E1-3G in ovarian stimulation and focused on their experimental verification and analysis. Methods: A prospective, observational pilot study was conducted involving 54 patients who underwent 54 cycles of ovarian stimulation. The goal was to establish the growth rate of urinary E1-3G during the course of stimulation and to determine the daily upper and lower limits of growth rates at which stimulation is appropriate and safe. Controlled ovarian stimulation was performed using two different stimulation protocols—an antagonist protocol in 25 cases and a progestin-primed ovarian stimulation protocol (PPOS) in 29 cases, with fixed doses of gonadotropins. From the second day of stimulation, patients self-measured their daily urine E1-3G levels at home using a portable analyzer. In parallel, a standard ultrasound follow-up protocol accompanied by a determination of E2, LH, and P levels was applied to optimally control stimulation. Results: The average daily growth rates in both groups were about 50%. The daily increase in E1-3G for the antagonist protocol ranged from 14% to 79%, while they were 28% to 79% for the PPOS protocol. Conclusion: This is the first study to analyze the dynamics of E1-3G in two different protocols and to estimate the limits of its increase during the entire course of the stimulation. The results confirm our theoretical model for the viability of using urinary E1-3G for monitoring ovarian stimulation.展开更多
A hyperbolic Lindstedt-Poincare method is presented to determine the homoclinic solutions of a kind of nonlinear oscillators, in which critical value of the homoclinic bifurcation parameter can be determined. The gene...A hyperbolic Lindstedt-Poincare method is presented to determine the homoclinic solutions of a kind of nonlinear oscillators, in which critical value of the homoclinic bifurcation parameter can be determined. The generalized Lienard oscillator is studied in detail, and the present method's predictions are compared with those of Runge-Kutta method to illustrate its accuracy.展开更多
With the development of large-scale spectral surveys, fiber positioning technology has been developing rapidly. Because of the performance advantages of a four-quadrant(4Q) detector, a fiber positioning and real-tim...With the development of large-scale spectral surveys, fiber positioning technology has been developing rapidly. Because of the performance advantages of a four-quadrant(4Q) detector, a fiber positioning and real-time monitoring system based on the 4Q detector is proposed. The detection accuracy of this system is directly determined by the precision of the center of the spot. A Gaussian fitting algorithm based on the 4Q detector is studied and applied in the fiber positioning process to improve the calculated accuracy of the spot center. The relationship between the center position of the incident spot and the detector output signal is deduced. An experimental platform is built to complete the simulated experiment. Then we use the Gaussian fitting method to process experimental data, compare the fitting value with the theoretical one and calculate the corresponding error.展开更多
The three-way catalysts (TWCs) promoters Ce0.6Zr0.4- x TbxO2-y were prepared by sol-gel method. BET surface areas analysis indicated that an increase of the dopant Tb content from x = 0.05 to x = 0.15 favors an incr...The three-way catalysts (TWCs) promoters Ce0.6Zr0.4- x TbxO2-y were prepared by sol-gel method. BET surface areas analysis indicated that an increase of the dopant Tb content from x = 0.05 to x = 0.15 favors an increase of surface area from 66.8 to 80.4 m^2· g^-1 compared with the undoped sample Ce0 .6oZr0.40O2 65.1 m^2·g^- 1 after calcination at 650℃. Transmission electron microscopy (TEM) observation indicated that the doped samples have a higher thermal stability. The XRD and Raman spectra confirmed that the Ce0.6Zr0.4-xTbxO2-y cubic solid solution is formed. XPS analysis revealed that Ce and Tb mainly existed in the form of Ce^4+ and Tb^3 + , and Zr existed in the form of Zr^4+ on the surface of the samples. The doped samples were homogenous in composition ; the introduction of Tb into the CeO2-ZrO2 promoters resuited in the formation of a solid solution, and the concentration of surface lattice oxygen was increased.展开更多
In general, triangular and quadrilateral elements are commonly applied in two-dimensional finite element methods. If they are used to compute polycrystalline materials, the cost of computation can be quite significant...In general, triangular and quadrilateral elements are commonly applied in two-dimensional finite element methods. If they are used to compute polycrystalline materials, the cost of computation can be quite significant. Polygonal elements can do well in simulation of the materials behavior and provide greater flexibility for the meshing of complex geometries. Hence, the study on the polygonal element is a very useful and necessary part in the finite element method. In this paper, an n-sided polygonal element based on quadratic spline interpolant, denoted by PS2 element, is presented using the triangular area coordinates and the B-net method. The PS2 element is conforming and can exactly model the quadratic field. It is valid for both convex and non-convex polygonal element, and insensitive to mesh distortions. In addition, no mapping or coordinate transformation is required and thus no Jacobian matrix and its inverse are evaluated. Some appropriate examples are employed to evaluate the performance of the proposed element.展开更多
Multi-train modeling and simulation plays a vital role in railway electrification during operation and planning phase. Study of peak power demand and energy consumed by each traction substation needs to be deter- mine...Multi-train modeling and simulation plays a vital role in railway electrification during operation and planning phase. Study of peak power demand and energy consumed by each traction substation needs to be deter- mined to verify that electrical energy flowing in its railway power feeding system is appropriate or not. Gauss-Seidel, conventional Newton-Raphson, and current injection methods are well-known and widely accepted as a tool for electrical power network solver in DC railway power supply study. In this paper, a simplified Newton-Raphson method has been proposed. The proposed method employs a set of current-balance equations at each electrical node instead of the conventional power-balance equation used in the conventional Newton-Raphson method. This concept can remarkably reduce execution time and computing complexity for multi-train simulation. To evaluate its use, Sukhumvit line of Bangkok transit system (BTS) of Thai- land with 21.6-km line length and 22 passenger stopping stations is set as a test system. The multi-train simulation integrated with the proposed power network solver is developed to simulate 1-h operation service of selected 5-min headway. From the obtained results, the proposed method is more efficient with approximately 18 % faster than the conventional Newton-Raphson method and just over 6 % faster than the current injection method.展开更多
BACKGROUND: Our group previously reported that right-sided vagus nerve stimulation(RVNS) significantly improved outcomes after cardiopulmonary resuscitation(CPR) in a rat model of cardiac arrest(CA). However, whether ...BACKGROUND: Our group previously reported that right-sided vagus nerve stimulation(RVNS) significantly improved outcomes after cardiopulmonary resuscitation(CPR) in a rat model of cardiac arrest(CA). However, whether left-sided vagus nerve stimulation(LVNS) could achieve the same effect as RVNS in CPR outcomes remains unknown.METHODS: A rat model of CA was established using modified percutaneous epicardial electrical stimulation to induce ventricular fibrillation(VF). Rats were treated with LVNS or RVNS for 30 minutes before the induction of VF. All animals were observed closely within 72 hours after return of spontaneous circulation(ROSC), and their health and behavior were evaluated every 24 hours.RESULTS: Compared with those in the RVNS group, the hemodynamic measurements in the LVNS group decreased more notably. Vagus nerve stimulation(VNS) decreased the serum levels of tumor necrosis factor-alpha(TNF-α) and the arrhythmia score, and attenuated inflammatory infiltration in myocardial tissue after ROSC, regardless of the side of stimulation, compared with findings in the CPR group. Both LVNS and RVNS ameliorated myocardial function and increased the expression of α-7 nicotinic acetylcholine receptor in the myocardium after ROSC. Moreover, a clear improvement in 72-hour survival was shown with VNS pre-treatment, with no significant difference in efficacy when comparing the laterality of stimulation. CONCLUSIONS: LVNS may have similar effects as RVNS on improving outcomes after CPR.展开更多
This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional (2D) incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fra...This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional (2D) incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling. The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by join- ing the centroid of cells sharing the common vertex. For the temporal integration of the momentum equations, an im- plicit second-order scheme is utilized to enhance the com- putational stability and eliminate the time step limit due to the diffusion term. The momentum equations are discretized by the vertex-centered finite volume method (FVM) and the pressure Poisson equation is solved by the Galerkin finite el- ement method (FEM). The momentum interpolation is used to damp out the spurious pressure wiggles. The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both veloc- ity and pressure. The classic test cases, the lid-driven cavity flow, the skew cavity flow and the backward-facing step flow, show that numerical results are in good agreement with the published benchmark solutions.展开更多
Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes w...Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes with- out connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approxima- tion. Then local discrete equations can be simplified by con- densation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by as- sembling all local discrete equations and are solved by using the standard implicit Newmark's time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is imple- mented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.展开更多
A multi-model integration method is proposed to develop a multi-source and heterogeneous model for short-term solar flare prediction. Different prediction models are constructed on the basis of extracted predictors fr...A multi-model integration method is proposed to develop a multi-source and heterogeneous model for short-term solar flare prediction. Different prediction models are constructed on the basis of extracted predictors from a pool of observation databases. The outputs of the base models are normal- ized first because these established models extract predictors from many data resources using different prediction methods. Then weighted integration of the base models is used to develop a multi-model integrated model (MIM). The weight set that single models assign is optimized by a genetic algorithm. Seven base models and data from Solar and Heliospheric Observatory/Michelson Doppler Imager lon- gitudinal magnetograms are used to construct the MIM, and then its performance is evaluated by cross validation. Experimental results showed that the MIM outperforms any individual model in nearly every data group, and the richer the diversity of the base models, the better the performance of the MIM. Thus, integrating more diversified models, such as an expert system, a statistical model and a physical model, will greatly improve the performance of the MIM.展开更多
In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We the...In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.展开更多
This paper extends Le van's work to the case of nonlinear problem and the complicated configuration. The wrinkling stress distribution and the pressure effects are also included in our analysis. Pseudo-beam method is...This paper extends Le van's work to the case of nonlinear problem and the complicated configuration. The wrinkling stress distribution and the pressure effects are also included in our analysis. Pseudo-beam method is presented based on the inflatable beam theory to model the inflatable structures as a set of inflatable beam elements with a prestressed state. In this method, the discretized nonlinear equations are given based upon the virtual work principle with a 3-node Timoshenko's beam model. Finite element simulation is performed by using a 3-node BEAM189 element incorporating ANSYS nonlinear program. The pressure effect is equivalent included in our method by modifying beam element cross-section parameters related to pressure. A benchmark example, the bending case of an inflatable cantilever beam, is performed to verify the accuracy of our proposed method. The comparisons reveal that the numerical results obtained with our method are close to open published analytical and membrane finite element results. The method is then used to evaluate the whole buckling and the loadcarrying characteristics of an inflatable support frame subjected to a compression force. The wrinkling stress and region characteristics are also shown in the end. This method gives better convergence characteristics, and requires much less computation time. It is very effective to deal with the whole load-carrying ability analytical problems for large scale inflatable structures with complex configuration.展开更多
Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions...Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions have some properties of simplicity and conformality. Two 8-node quadrilateral elements have been developed using the trian- gular area coordinates and the B-net method, which can ex- actly model the quadratic field for both convex and concave quadrangles. Some appropriate examples are employed to evaluate the performance of the proposed elements. The nu- merical results show that the two spline elements can obtain solutions which are highly accurate and insensitive to mesh distortions.展开更多
A precursor of BaCe0.5Zr0.4Y0. 1O3-α electrolytes was synthesized by the sol-gel method and sintered at temperature which were 150 - 250 ℃ lower than by solid state reaction. The AC impedance spectrums of electrolyt...A precursor of BaCe0.5Zr0.4Y0. 1O3-α electrolytes was synthesized by the sol-gel method and sintered at temperature which were 150 - 250 ℃ lower than by solid state reaction. The AC impedance spectrums of electrolytes were measured by AUTOLAB PGSTA30 electrochemical measuring device at different temperatures. The conductivities of the electrolytes are 1.62×10^-4 - 6.43×10^-3, 2.52×10^-5 - 3.73×10^-3S·cm^-1 in the temperature range of 350-800℃. The activity energies are 0.54 and 0.84 eV. At the same time BaCe0.9Y0.1O3-α was prepared by direct solid state reaction. The conductivity of BaCe0.9Y0.1O3-α is 1 × 10^-4- 4×10^-3 S·cm^-1 and the activation energy is 0.50 eV at the same condition. The results show that conduction of electrolyte prepared by sol-gel method is higher than the one by solid state reaction. As far as BaCe0.9Y0.1O3-α concerned, its conductivity of the Zr-substituted specimens is decreased.展开更多
As concepts closely related to microwave absorption properties,impedance matching and phase matching were rarely combined with material parameters to regulate properties and explore related mechanisms.In this work,red...As concepts closely related to microwave absorption properties,impedance matching and phase matching were rarely combined with material parameters to regulate properties and explore related mechanisms.In this work,reduction–diffusion method was innovatively applied to synthesize rare earth alloy Y_(2)Fe_(17).In order to regulate the electromagnetic parameters of absorbers,the Y_(2)Fe_(17)N_(3-δ)particles were coated with silica(Y_(2)Fe_(17)N_(3-δ)@SiO_(2))and absorbers with different volume fractions were prepared.The relationship between impedance matching,matching thickness,and the strongest reflection loss peak(RLmin)was presented obviously.Compared to the microwave absorption properties of Y_(2)Fe_(17)N_(3-δ)/PU absorber,Y_(2)Fe_(17)N_(3-δ)@SiO_(2)/PU absorbers are more conducive to the realization of microwave absorption material standards which are thin thickness,light weight,strong absorbing intensity,and broad bandwidth.Based on microwave frequency bands,the microwave absorption properties of the absorbers were analyzed and the related parameters were listed.As an important parameter related to perfect matching,reflection factor(√ε_(r)/μ_(r))was discussed combined with microwave amplitude attenuation.According to the origin and mathematical model of bandwidth,the formula of EAB(RL<-10 dB)was derived and simplified.The calculated bandwidths agreed well with experimental results.展开更多
In this article, we propose a generalized exp(-Φ(ξ))-expansion method and successfully implement it to find exact traveling wave solutions to the fifth order standard Sawada-Kotera (SK) equation. The exact traveling...In this article, we propose a generalized exp(-Φ(ξ))-expansion method and successfully implement it to find exact traveling wave solutions to the fifth order standard Sawada-Kotera (SK) equation. The exact traveling wave solutions are established in the form of trigonometric, hyperbolic, exponential and rational functions with some free parameters. It is shown that this method is standard, effective and easily applicable mathematical tool for solving nonlinear partial differential equations arises in the field of mathematical physics and engineering.展开更多
Applying calculation method in alloy design should be an important tendency due to its characters of inexpensive cost, high efficiency and prediction. DOS calculations of AuSn, AsSn and SbSn Sn- based alloys have ...Applying calculation method in alloy design should be an important tendency due to its characters of inexpensive cost, high efficiency and prediction. DOS calculations of AuSn, AsSn and SbSn Sn- based alloys have been investigated by employing DV - Xa method, in which different cluster models were adopted to calculate electron structure.It is proved that some regulations must be taken into ac- count in order to carry out alloy design calculation successfully,which are described in this paper in detail.展开更多
文摘Background: Studies have shown a strong correlation between the growth of E2 in serum and estrone-3-glucuronide (E1-3G) in urine during ovarian stimulation. Thus, we developed theoretical models for using urinary E1-3G in ovarian stimulation and focused on their experimental verification and analysis. Methods: A prospective, observational pilot study was conducted involving 54 patients who underwent 54 cycles of ovarian stimulation. The goal was to establish the growth rate of urinary E1-3G during the course of stimulation and to determine the daily upper and lower limits of growth rates at which stimulation is appropriate and safe. Controlled ovarian stimulation was performed using two different stimulation protocols—an antagonist protocol in 25 cases and a progestin-primed ovarian stimulation protocol (PPOS) in 29 cases, with fixed doses of gonadotropins. From the second day of stimulation, patients self-measured their daily urine E1-3G levels at home using a portable analyzer. In parallel, a standard ultrasound follow-up protocol accompanied by a determination of E2, LH, and P levels was applied to optimally control stimulation. Results: The average daily growth rates in both groups were about 50%. The daily increase in E1-3G for the antagonist protocol ranged from 14% to 79%, while they were 28% to 79% for the PPOS protocol. Conclusion: This is the first study to analyze the dynamics of E1-3G in two different protocols and to estimate the limits of its increase during the entire course of the stimulation. The results confirm our theoretical model for the viability of using urinary E1-3G for monitoring ovarian stimulation.
基金supported by the National Natural Science Foundation of China (10672193)Sun Yat-sen University (Fu Lan Scholarship)the University of Hong Kong (CRGC grant).
文摘A hyperbolic Lindstedt-Poincare method is presented to determine the homoclinic solutions of a kind of nonlinear oscillators, in which critical value of the homoclinic bifurcation parameter can be determined. The generalized Lienard oscillator is studied in detail, and the present method's predictions are compared with those of Runge-Kutta method to illustrate its accuracy.
基金support by the Fundamental Research Funds for the Central Universities of China (2013/B15020271)the National Natural Science Foundation of China (1014/515029111)the National Undergraduate Training Program for Innovation and Entrepreneurship (201610294069)
文摘With the development of large-scale spectral surveys, fiber positioning technology has been developing rapidly. Because of the performance advantages of a four-quadrant(4Q) detector, a fiber positioning and real-time monitoring system based on the 4Q detector is proposed. The detection accuracy of this system is directly determined by the precision of the center of the spot. A Gaussian fitting algorithm based on the 4Q detector is studied and applied in the fiber positioning process to improve the calculated accuracy of the spot center. The relationship between the center position of the incident spot and the detector output signal is deduced. An experimental platform is built to complete the simulated experiment. Then we use the Gaussian fitting method to process experimental data, compare the fitting value with the theoretical one and calculate the corresponding error.
基金Project supported by the Natural Science Basic Research Plan in Shaanxi Province of China (2004B13)
文摘The three-way catalysts (TWCs) promoters Ce0.6Zr0.4- x TbxO2-y were prepared by sol-gel method. BET surface areas analysis indicated that an increase of the dopant Tb content from x = 0.05 to x = 0.15 favors an increase of surface area from 66.8 to 80.4 m^2· g^-1 compared with the undoped sample Ce0 .6oZr0.40O2 65.1 m^2·g^- 1 after calcination at 650℃. Transmission electron microscopy (TEM) observation indicated that the doped samples have a higher thermal stability. The XRD and Raman spectra confirmed that the Ce0.6Zr0.4-xTbxO2-y cubic solid solution is formed. XPS analysis revealed that Ce and Tb mainly existed in the form of Ce^4+ and Tb^3 + , and Zr existed in the form of Zr^4+ on the surface of the samples. The doped samples were homogenous in composition ; the introduction of Tb into the CeO2-ZrO2 promoters resuited in the formation of a solid solution, and the concentration of surface lattice oxygen was increased.
基金supported by the National Natural Science Foundation of China (60533060, 10672032, 10726067)Science Foundation of Dalian University of Technology (SFDUT07001)
文摘In general, triangular and quadrilateral elements are commonly applied in two-dimensional finite element methods. If they are used to compute polycrystalline materials, the cost of computation can be quite significant. Polygonal elements can do well in simulation of the materials behavior and provide greater flexibility for the meshing of complex geometries. Hence, the study on the polygonal element is a very useful and necessary part in the finite element method. In this paper, an n-sided polygonal element based on quadratic spline interpolant, denoted by PS2 element, is presented using the triangular area coordinates and the B-net method. The PS2 element is conforming and can exactly model the quadratic field. It is valid for both convex and non-convex polygonal element, and insensitive to mesh distortions. In addition, no mapping or coordinate transformation is required and thus no Jacobian matrix and its inverse are evaluated. Some appropriate examples are employed to evaluate the performance of the proposed element.
文摘Multi-train modeling and simulation plays a vital role in railway electrification during operation and planning phase. Study of peak power demand and energy consumed by each traction substation needs to be deter- mined to verify that electrical energy flowing in its railway power feeding system is appropriate or not. Gauss-Seidel, conventional Newton-Raphson, and current injection methods are well-known and widely accepted as a tool for electrical power network solver in DC railway power supply study. In this paper, a simplified Newton-Raphson method has been proposed. The proposed method employs a set of current-balance equations at each electrical node instead of the conventional power-balance equation used in the conventional Newton-Raphson method. This concept can remarkably reduce execution time and computing complexity for multi-train simulation. To evaluate its use, Sukhumvit line of Bangkok transit system (BTS) of Thai- land with 21.6-km line length and 22 passenger stopping stations is set as a test system. The multi-train simulation integrated with the proposed power network solver is developed to simulate 1-h operation service of selected 5-min headway. From the obtained results, the proposed method is more efficient with approximately 18 % faster than the conventional Newton-Raphson method and just over 6 % faster than the current injection method.
基金supported by research grants from the National Natural Science Foundation of China (81571866, 82072137)。
文摘BACKGROUND: Our group previously reported that right-sided vagus nerve stimulation(RVNS) significantly improved outcomes after cardiopulmonary resuscitation(CPR) in a rat model of cardiac arrest(CA). However, whether left-sided vagus nerve stimulation(LVNS) could achieve the same effect as RVNS in CPR outcomes remains unknown.METHODS: A rat model of CA was established using modified percutaneous epicardial electrical stimulation to induce ventricular fibrillation(VF). Rats were treated with LVNS or RVNS for 30 minutes before the induction of VF. All animals were observed closely within 72 hours after return of spontaneous circulation(ROSC), and their health and behavior were evaluated every 24 hours.RESULTS: Compared with those in the RVNS group, the hemodynamic measurements in the LVNS group decreased more notably. Vagus nerve stimulation(VNS) decreased the serum levels of tumor necrosis factor-alpha(TNF-α) and the arrhythmia score, and attenuated inflammatory infiltration in myocardial tissue after ROSC, regardless of the side of stimulation, compared with findings in the CPR group. Both LVNS and RVNS ameliorated myocardial function and increased the expression of α-7 nicotinic acetylcholine receptor in the myocardium after ROSC. Moreover, a clear improvement in 72-hour survival was shown with VNS pre-treatment, with no significant difference in efficacy when comparing the laterality of stimulation. CONCLUSIONS: LVNS may have similar effects as RVNS on improving outcomes after CPR.
基金supported by the Natural Science Foundation of China (11061021)the Program of Higher-level talents of Inner Mongolia University (SPH-IMU,Z200901004)the Scientific Research Projection of Higher Schools of Inner Mongolia(NJ10016,NJ10006)
文摘This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional (2D) incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling. The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by join- ing the centroid of cells sharing the common vertex. For the temporal integration of the momentum equations, an im- plicit second-order scheme is utilized to enhance the com- putational stability and eliminate the time step limit due to the diffusion term. The momentum equations are discretized by the vertex-centered finite volume method (FVM) and the pressure Poisson equation is solved by the Galerkin finite el- ement method (FEM). The momentum interpolation is used to damp out the spurious pressure wiggles. The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both veloc- ity and pressure. The classic test cases, the lid-driven cavity flow, the skew cavity flow and the backward-facing step flow, show that numerical results are in good agreement with the published benchmark solutions.
基金supported by the National Natural Science Founda-tion of China(11272118)Open Found of State Key Laboratory of Explosion Science and Technology(KFJJ12-5M)
文摘Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes with- out connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approxima- tion. Then local discrete equations can be simplified by con- densation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by as- sembling all local discrete equations and are solved by using the standard implicit Newmark's time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is imple- mented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.
基金supported by the National Natural Science Foundation of China(Grant No.11078010)SOHO is a project of international cooperation between the European Space Agency(ESA) and NASA
文摘A multi-model integration method is proposed to develop a multi-source and heterogeneous model for short-term solar flare prediction. Different prediction models are constructed on the basis of extracted predictors from a pool of observation databases. The outputs of the base models are normal- ized first because these established models extract predictors from many data resources using different prediction methods. Then weighted integration of the base models is used to develop a multi-model integrated model (MIM). The weight set that single models assign is optimized by a genetic algorithm. Seven base models and data from Solar and Heliospheric Observatory/Michelson Doppler Imager lon- gitudinal magnetograms are used to construct the MIM, and then its performance is evaluated by cross validation. Experimental results showed that the MIM outperforms any individual model in nearly every data group, and the richer the diversity of the base models, the better the performance of the MIM. Thus, integrating more diversified models, such as an expert system, a statistical model and a physical model, will greatly improve the performance of the MIM.
基金The project supported by the China NKBRSF(2001CB409604)
文摘In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.
基金supported by the Specialized Fund for the Doctoral Program of Higher Education of China (200802131046)China Postdoctoral Science Foundation Funded Major Project (200801290)+1 种基金Development Program of Outstanding Young Teachers in Harbin Institute of Technology (HITQNJS.2008.004)Specialized Fund for Innovation Talents of Science and Technology in Harbin (2008RFQXG057).
文摘This paper extends Le van's work to the case of nonlinear problem and the complicated configuration. The wrinkling stress distribution and the pressure effects are also included in our analysis. Pseudo-beam method is presented based on the inflatable beam theory to model the inflatable structures as a set of inflatable beam elements with a prestressed state. In this method, the discretized nonlinear equations are given based upon the virtual work principle with a 3-node Timoshenko's beam model. Finite element simulation is performed by using a 3-node BEAM189 element incorporating ANSYS nonlinear program. The pressure effect is equivalent included in our method by modifying beam element cross-section parameters related to pressure. A benchmark example, the bending case of an inflatable cantilever beam, is performed to verify the accuracy of our proposed method. The comparisons reveal that the numerical results obtained with our method are close to open published analytical and membrane finite element results. The method is then used to evaluate the whole buckling and the loadcarrying characteristics of an inflatable support frame subjected to a compression force. The wrinkling stress and region characteristics are also shown in the end. This method gives better convergence characteristics, and requires much less computation time. It is very effective to deal with the whole load-carrying ability analytical problems for large scale inflatable structures with complex configuration.
基金supported by the National Natural Science Foundation of China(11001037,11102037 and 11290143)the Fundamental Research Funds for the Central Universities
文摘Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions have some properties of simplicity and conformality. Two 8-node quadrilateral elements have been developed using the trian- gular area coordinates and the B-net method, which can ex- actly model the quadratic field for both convex and concave quadrangles. Some appropriate examples are employed to evaluate the performance of the proposed elements. The nu- merical results show that the two spline elements can obtain solutions which are highly accurate and insensitive to mesh distortions.
文摘A precursor of BaCe0.5Zr0.4Y0. 1O3-α electrolytes was synthesized by the sol-gel method and sintered at temperature which were 150 - 250 ℃ lower than by solid state reaction. The AC impedance spectrums of electrolytes were measured by AUTOLAB PGSTA30 electrochemical measuring device at different temperatures. The conductivities of the electrolytes are 1.62×10^-4 - 6.43×10^-3, 2.52×10^-5 - 3.73×10^-3S·cm^-1 in the temperature range of 350-800℃. The activity energies are 0.54 and 0.84 eV. At the same time BaCe0.9Y0.1O3-α was prepared by direct solid state reaction. The conductivity of BaCe0.9Y0.1O3-α is 1 × 10^-4- 4×10^-3 S·cm^-1 and the activation energy is 0.50 eV at the same condition. The results show that conduction of electrolyte prepared by sol-gel method is higher than the one by solid state reaction. As far as BaCe0.9Y0.1O3-α concerned, its conductivity of the Zr-substituted specimens is decreased.
基金Project supported by the National Key Research and Development Program of China(Grant No.2021YFB3501300)the National Natural Science Foundation of China(Grant No.51731001)the Fund from the State Key Laboratory of Baiyunobo Rare Earth Resource Researches and Comprehensive Utilization’s Key Research and Development Projects。
文摘As concepts closely related to microwave absorption properties,impedance matching and phase matching were rarely combined with material parameters to regulate properties and explore related mechanisms.In this work,reduction–diffusion method was innovatively applied to synthesize rare earth alloy Y_(2)Fe_(17).In order to regulate the electromagnetic parameters of absorbers,the Y_(2)Fe_(17)N_(3-δ)particles were coated with silica(Y_(2)Fe_(17)N_(3-δ)@SiO_(2))and absorbers with different volume fractions were prepared.The relationship between impedance matching,matching thickness,and the strongest reflection loss peak(RLmin)was presented obviously.Compared to the microwave absorption properties of Y_(2)Fe_(17)N_(3-δ)/PU absorber,Y_(2)Fe_(17)N_(3-δ)@SiO_(2)/PU absorbers are more conducive to the realization of microwave absorption material standards which are thin thickness,light weight,strong absorbing intensity,and broad bandwidth.Based on microwave frequency bands,the microwave absorption properties of the absorbers were analyzed and the related parameters were listed.As an important parameter related to perfect matching,reflection factor(√ε_(r)/μ_(r))was discussed combined with microwave amplitude attenuation.According to the origin and mathematical model of bandwidth,the formula of EAB(RL<-10 dB)was derived and simplified.The calculated bandwidths agreed well with experimental results.
文摘In this article, we propose a generalized exp(-Φ(ξ))-expansion method and successfully implement it to find exact traveling wave solutions to the fifth order standard Sawada-Kotera (SK) equation. The exact traveling wave solutions are established in the form of trigonometric, hyperbolic, exponential and rational functions with some free parameters. It is shown that this method is standard, effective and easily applicable mathematical tool for solving nonlinear partial differential equations arises in the field of mathematical physics and engineering.
文摘Applying calculation method in alloy design should be an important tendency due to its characters of inexpensive cost, high efficiency and prediction. DOS calculations of AuSn, AsSn and SbSn Sn- based alloys have been investigated by employing DV - Xa method, in which different cluster models were adopted to calculate electron structure.It is proved that some regulations must be taken into ac- count in order to carry out alloy design calculation successfully,which are described in this paper in detail.