Aerodynamic performances of axial compressors are significantly affected by variation of Reynolds number in aero-engines.In the design and analysis of compressors,previous correction methods for cascades and stages ha...Aerodynamic performances of axial compressors are significantly affected by variation of Reynolds number in aero-engines.In the design and analysis of compressors,previous correction methods for cascades and stages have difficulties in predicting comprehensively Reynolds number effects on airfoils,matching and characteristics curves.This study proposes Re-correction models for loss,deviation angle and endwall blockage based on classical theories and cascade tests,and loss and deviation models show good agreement in test data of NACA65 and C4 cascades.Throughflow method considering Reynolds number effects is developed by integrating the correction models into a verified Streamline Curvature(SLC)tool.A three-stage axial compressor is investigated through SLC and CFD methods from design Reynolds number(Red=2106)to low Re=4104,and the numerical methods are validated with test data of characteristic curves and spanwise distributions at Red.With Re reduction,SLC method with correction models well predicts variation in overall performances compared with CFD calculations and Wassell's model.Streamwise and spanwise matching such as total pressure and loss distributions in SLC predictions are basically consistent with those in CFD results at near-stall points under design and low Reynolds numbers.SLC and CFD methods share similar detections of stall risks in the third stage(Stg3),and their analyses of diffusion processes deviate to some extent due to different predictions in separated endwall flow.The correction models can be adopted to consider Reynolds number effects in through-flow design and analysis of axial compressors.展开更多
A linear-correction least-squares（LCLS） estimation procedure is proposed for geolocation using frequency difference of arrival （FDOA） measurements only. We first analyze the measurements of FDOA, and further deriv...A linear-correction least-squares（LCLS） estimation procedure is proposed for geolocation using frequency difference of arrival （FDOA） measurements only. We first analyze the measurements of FDOA, and further derive the Cramer-Rao lower bound （CRLB） of geoloeation using FDOA measurements. For the localization model is a nonlinear least squares（LS） estimator with a nonlinear constrained, a linearizing method is used to convert the model to a linear least squares estimator with a nonlinear con- strained. The Gauss-Newton iteration method is developed to conquer the source localization problem. From the analysis of solving Lagrange multiplier, the algorithm is a generalization of linear-correction least squares estimation procedure under the condition of geolocation using FDOA measurements only. The algorithm is compared with common least squares estimation. Comparisons of their estimation accuracy and the CRLB are made, and the proposed method attains the CRLB. Simulation re- sults are included to corroborate the theoretical development.展开更多
This study discusses the challenges in logging the evaluation of low-resistivity oil reservoirs,especially the difficult problems involving their saturation calculation.A correction method for equivalent water conduct...This study discusses the challenges in logging the evaluation of low-resistivity oil reservoirs,especially the difficult problems involving their saturation calculation.A correction method for equivalent water conductivity is proposed,given the high conductivity caused by small amounts of water distributed in a network within the low-resistivity reservoir,which mimics the eff ects of high water saturation.This approach signifi cantly improves the accuracy of hydrocarbon saturation calculations in these low-resistivity reservoirs.The corrected hydrocarbon saturation values highly matched the porosity and are consistent with experimental results.This study also establishes a discrimination process to determine whether corrections are required,verifying the eff ectiveness and accuracy of the method through an application example.展开更多
In high-resolution cone-beam computed tomography （CBCT） using the flat-panel detector, imperfect or defect detector elements cause ring artifacts due to the none-uniformity of their X-ray response. They often distur...In high-resolution cone-beam computed tomography （CBCT） using the flat-panel detector, imperfect or defect detector elements cause ring artifacts due to the none-uniformity of their X-ray response. They often disturb the image quality. A dedicated fitting correction method for high-resolution micro-CT is presented. The method converts each elementary X-ray response curve to an average one, and eliminates response inconsistency among pixels. Other factors of the method are discussed, such as the correction factor variability by different sampling frames and nonlinear factors over the whole spectrum. Results show that the noise and artifacts are both reduced in reconstructed images展开更多
The efficacy of error correction and various kinds of correction approaches is one of the key issues in second language writing faced by both teachers and researchers. The current paper reviews the definition of error...The efficacy of error correction and various kinds of correction approaches is one of the key issues in second language writing faced by both teachers and researchers. The current paper reviews the definition of error correction and examines the different views on whether error correction in L2 writing should be corrected. In particular, the paper discusses and analyses the three common correction methods: direct correction, peer feedback and indirect correction. Teachers are encouraged to weigh and analyze the advantages and disadvantages of these methods according to the current literature, employ the most beneficial error correction method in L2 writing, and adapt its suitability to their teaching context.展开更多
Based on the hindcast results of summer rainfall anomalies over China for the period 1981-2000 by the Dynamical Climate Prediction System （IAP-DCP） developed by the Institute of Atmospheric Physics, a correction met...Based on the hindcast results of summer rainfall anomalies over China for the period 1981-2000 by the Dynamical Climate Prediction System （IAP-DCP） developed by the Institute of Atmospheric Physics, a correction method that can account for the dependence of model＇s systematic biases on SST anomalies is proposed. It is shown that this correction method can improve the hindcast skill of the IAP-DCP for summer rainfall anomalies over China, especially in western China and southeast China, which may imply its potential application to real-time seasonal prediction.展开更多
Surface low/distortion is one of the most challenging surface deflections that have a great effect on the exterior appearance of automobiles.Most studies on surface distortion/deflection have focused on evaluation and...Surface low/distortion is one of the most challenging surface deflections that have a great effect on the exterior appearance of automobiles.Most studies on surface distortion/deflection have focused on evaluation and visualization techniques,the research work on correction or prevention of surface low/distortion is limited,and there is no perfect surface low/distortion corrective method that can satisfy the needs of the engineering.A B-spline based geometry morphing algorithm is proposed and then a new program based on UG-NX platform is developed to modify the die face in the surface low/distortion areas.To verify this developed system,the experimental dies that can replicate the surface low/distortion phenomenon successfully is put to use.Five geometric variables are introduced to describe the basic geometry of typical depression features of automotive outer panels.The experimental dies are then designed to reflect various combinations of these geometric parameters.The stamping experiments are conducted on cold rolled grade 5（CR5） sheet steel and various static measurements,such as oil-stoning,laser scanner,etc,are performed to measure and record the surface low/distortions.Three approaches including good bearing,holds in blank and die face morphing that aim to correct low/distortions are tried out and surface low/distortions are observed in the specimen with reverse draw depth of 10 mm.The measurement results show that die morphing is a practical and effective method to correct the surface low/distortion.The correction method proposed can be used to minimize the occurrence of surface low/distortion in die manufacturing,which has certain reference significance to the correction of surface low/distortion.展开更多
In this paper, an analogue correction method of errors （ACE） based on a complicated atmospheric model is further developed and applied to numerical weather prediction （NWP）. The analysis shows that the ACE can eff...In this paper, an analogue correction method of errors （ACE） based on a complicated atmospheric model is further developed and applied to numerical weather prediction （NWP）. The analysis shows that the ACE can effectively reduce model errors by combining the statistical analogue method with the dynamical model together in order that the information of plenty of historical data is utilized in the current complicated NWP model, Furthermore, in the ACE, the differences of the similarities between different historical analogues and the current initial state are considered as the weights for estimating model errors. The results of daily, decad and monthly prediction experiments on a complicated T63 atmospheric model show that the performance of the ACE by correcting model errors based on the estimation of the errors of 4 historical analogue predictions is not only better than that of the scheme of only introducing the correction of the errors of every single analogue prediction, but is also better than that of the T63 model.展开更多
The blind-hole method is the most widely used approach to experimentally determine the distribution of residual stress. This paper aims to improve test accuracy of welding residual stress and conducts an experimental ...The blind-hole method is the most widely used approach to experimentally determine the distribution of residual stress. This paper aims to improve test accuracy of welding residual stress and conducts an experimental study on the strain release factors involved when using the blind-hole method for Q235 and Q345, two steels commonly used in building structures. The ranges of strain release factors A and B in the elastic stage, the effects of strain release factors on residual stress calculated values, and the plastic corrected strain release factors are analyzed considering of the effect of plastic deformation around the blind hole on measurement accuracy. Finally, a simplified calculation formula to determine strain release factors is proposed for use with the blind-hole method. Results show that in the elastic stage, strain release factor A for Q235 and Q345 ranges from-0.399 to-0.525 and strain release factor B from-0.791 to-0.960. Changing the strain release factors A and B shows that calculated residual tensile stress varies in relation to a decrease in both factor values. However, there is a increase in calculated residual compressive stress with a decrease in the strain release factor A value, but there is an decrease with a decrease in strain release factor B value. Calculated residual stress applied to elastic strain release factors is compared with that applied to amended plastic strain release factors for Q235 steel. The maximum deviation between calculated residual stress and test stress is reduced from 21.1 to 1.0%,and for Q345 steel from 26.5 to 1.2%. It is thus evident that the plastic correction formula proposed in this paper can be used in calculations when conducting a residual stress test.展开更多
A third-order correction was recently suggested to improve the accuracy of the half-power bandwidth method in estimating the damping of single DOF systems.This paper analyzes the accuracy of the half-power bandwidth m...A third-order correction was recently suggested to improve the accuracy of the half-power bandwidth method in estimating the damping of single DOF systems.This paper analyzes the accuracy of the half-power bandwidth method with the third-order correction in damping estimation for multi-DOF linear systems.Damping ratios in a two-DOF linear system are estimated using its displacement and acceleration frequency response curves,respectively.A wide range of important parameters that characterize the shape of these response curves are taken into account.Results show that the third-order correction may greatly improve the accuracy of the half-power bandwidth method in estimating damping in a two-DOF system.In spite of this,the half-power bandwidth method may significantly overestimate the damping ratios of two-DOF systems in some cases.展开更多
The main purpose of reverse engineering is to convert discrete data pointsinto piecewise smooth, continuous surface models. Before carrying out model reconstruction it issignificant to extract geometric features becau...The main purpose of reverse engineering is to convert discrete data pointsinto piecewise smooth, continuous surface models. Before carrying out model reconstruction it issignificant to extract geometric features because the quality of modeling greatly depends on therepresentation of features. Some fitting techniques of natural quadric surfaces with least-squaresmethod are described. And these techniques can be directly used to extract quadric surfaces featuresduring the process of segmentation for point cloud.展开更多
The Galerkin and least-squares methods are two classes of the most popular Krylov subspace methOds for solving large linear systems of equations. Unfortunately, both the methods may suffer from serious breakdowns of t...The Galerkin and least-squares methods are two classes of the most popular Krylov subspace methOds for solving large linear systems of equations. Unfortunately, both the methods may suffer from serious breakdowns of the same type: In a breakdown situation the Galerkin method is unable to calculate an approximate solution, while the least-squares method, although does not really break down, is unsucessful in reducing the norm of its residual. In this paper we first establish a unified theorem which gives a relationship between breakdowns in the two methods. We further illustrate theoretically and experimentally that if the coefficient matrix of a lienar system is of high defectiveness with the associated eigenvalues less than 1, then the restarted Galerkin and least-squares methods will be in great risks of complete breakdowns. It appears that our findings may help to understand phenomena observed practically and to derive treatments for breakdowns of this type.展开更多
Relaxation time spectra （RTS） derived from time domain induced polarization data （TDIP） are helpful to assess oil reservoir pore structures. However, due to the sensitivity to the signal-to-noise ratio （SNR）, th...Relaxation time spectra （RTS） derived from time domain induced polarization data （TDIP） are helpful to assess oil reservoir pore structures. However, due to the sensitivity to the signal-to-noise ratio （SNR）, the inversion accuracy of the traditional singular value decomposition （SVD） inversion method reduces with a decrease of SNR. In order to enhance the inversion accuracy and improve robustness of the inversion method to the SNR, an improved inversion method, based on damping factor and spectrum component residual correction, is proposed in this study. The numerical inversion results show that the oscillation of the RTS derived from the SVD method increased with a decrease of SNR, which makes it impossible to get accurate inversion components. However, the SNR has little influence on inversion components of the improved method, and the RTS has high inversion accuracy and robustness. Moreover, RTS derived from core sample data is basically in accord with the pore-size distribution curve, and the RTS derived from the actual induced polarization logging data is smooth and continuous, which indicates that the improved method is practicable.展开更多
A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercriti...A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.展开更多
In this work,synchronous cutting of concave and convex surfaces was achieved using the duplex helical method for the hypoid gear,and the problem of tooth surface error correction was studied.First,the mathematical mod...In this work,synchronous cutting of concave and convex surfaces was achieved using the duplex helical method for the hypoid gear,and the problem of tooth surface error correction was studied.First,the mathematical model of the hypoid gears machined by the duplex helical method was established.Second,the coordinates of discrete points on the tooth surface were obtained by measurement center,and the normal errors of the discrete points were calculated.Third,a tooth surface error correction model is established,and the tooth surface error was corrected using the Levenberg-Marquard algorithm with trust region strategy and least square method.Finally,grinding experiments were carried out on the machining parameters obtained by Levenberg-Marquard algorithm with trust region strategy,which had a better effect on tooth surface error correction than the least square method.After the tooth surface error is corrected,the maximum absolute error is reduced from 30.9μm before correction to 6.8μm,the root mean square of the concave error is reduced from 15.1 to 2.1μm,the root mean square of the convex error is reduced from 10.8 to 1.8μm,and the sum of squared errors of the concave and convex surfaces was reduced from 15471 to 358μm^(2).It is verified that the Levenberg-Marquard algorithm with trust region strategy has a good accuracy for the tooth surface error correction of hypoid gear machined by duplex helical method.展开更多
This is a study of the Durand-Kerner and Nourein methods for finding the roots of a given algebraic equation simultaneously. We consider the conditions under which the iterative methods fail. The numerical example is ...This is a study of the Durand-Kerner and Nourein methods for finding the roots of a given algebraic equation simultaneously. We consider the conditions under which the iterative methods fail. The numerical example is presented.展开更多
Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of s...Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.展开更多
The purpose of this article is to develop and analyze least-squares approximations for the incompressible magnetohydrodynamic equations. The major advantage of the least-squares finite element method is that it is not...The purpose of this article is to develop and analyze least-squares approximations for the incompressible magnetohydrodynamic equations. The major advantage of the least-squares finite element method is that it is not subjected to the so-called Ladyzhenskaya-Babuska-Brezzi （LBB） condition. The authors employ least-squares functionals which involve a discrete inner product which is related to the inner product in H^-1（Ω）.展开更多
With the development of computational power, there has been an increased focus on data-fitting related seismic inversion techniques for high fidelity seismic velocity model and image, such as full-waveform inversion a...With the development of computational power, there has been an increased focus on data-fitting related seismic inversion techniques for high fidelity seismic velocity model and image, such as full-waveform inversion and least squares migration. However, though more advanced than conventional methods, these data fitting methods can be very expensive in terms of computational cost. Recently, various techniques to optimize these data-fitting seismic inversion problems have been implemented to cater for the industrial need for much improved efficiency. In this study, we propose a general stochastic conjugate gradient method for these data-fitting related inverse problems. We first prescribe the basic theory of our method and then give synthetic examples. Our numerical experiments illustrate the potential of this method for large-size seismic inversion application.展开更多
For the second-order finite volume method,implicit schemes and reconstruction methods are two main algorithms which influence the robustness and efficiency of the numerical simulations of compressible turbulent flows....For the second-order finite volume method,implicit schemes and reconstruction methods are two main algorithms which influence the robustness and efficiency of the numerical simulations of compressible turbulent flows.In this paper,a compact least-squares reconstruction method is proposed to calculate the gradients for the distribution of flow field variables approximation.The compactness of the new reconstruction method is reflected in the gradient calculation process.The geometries of the face-neighboring elements are no longer utilized,and the weighted average values at the centroid of the interfaces are used to calculate the gradients instead of the values at the centroid of the face-neighboring elements.Meanwhile,an exact Jacobian solving strategy is developed for implicit temporal discretization.The accurate processing of Jacobian matrix can extensively improve the invertibility of the Jacobian matrix and avoid introducing extra numerical errors.In addition,a modified Venkatakrishnan limiter is applied to deal with the shock which may appear in transonic flows and the applicability of the mentioned methods is enhanced further.The combination of the proposed methods makes the numerical simulations of turbulent flow converge rapidly and steadily with an adaptive increasing CFL number.The numerical results of several benchmarks indicate that the proposed methods perform well in terms of robustness,efficiency and accuracy,and have good application potential in turbulent flow simulations of complex configurations.展开更多
基金supported by the National Science and Tech-nology Major Project of China(Nos.2017-II-0007-0021 and J2019-II-0017-0038)。
文摘Aerodynamic performances of axial compressors are significantly affected by variation of Reynolds number in aero-engines.In the design and analysis of compressors,previous correction methods for cascades and stages have difficulties in predicting comprehensively Reynolds number effects on airfoils,matching and characteristics curves.This study proposes Re-correction models for loss,deviation angle and endwall blockage based on classical theories and cascade tests,and loss and deviation models show good agreement in test data of NACA65 and C4 cascades.Throughflow method considering Reynolds number effects is developed by integrating the correction models into a verified Streamline Curvature(SLC)tool.A three-stage axial compressor is investigated through SLC and CFD methods from design Reynolds number(Red=2106)to low Re=4104,and the numerical methods are validated with test data of characteristic curves and spanwise distributions at Red.With Re reduction,SLC method with correction models well predicts variation in overall performances compared with CFD calculations and Wassell's model.Streamwise and spanwise matching such as total pressure and loss distributions in SLC predictions are basically consistent with those in CFD results at near-stall points under design and low Reynolds numbers.SLC and CFD methods share similar detections of stall risks in the third stage(Stg3),and their analyses of diffusion processes deviate to some extent due to different predictions in separated endwall flow.The correction models can be adopted to consider Reynolds number effects in through-flow design and analysis of axial compressors.
基金National High-tech Research and Development Program of China (2011AA7072043)National Defense Key Laboratory Foundation of China (9140C860304)Innovation Fund of Graduate School of NUDT (B120406)
文摘A linear-correction least-squares（LCLS） estimation procedure is proposed for geolocation using frequency difference of arrival （FDOA） measurements only. We first analyze the measurements of FDOA, and further derive the Cramer-Rao lower bound （CRLB） of geoloeation using FDOA measurements. For the localization model is a nonlinear least squares（LS） estimator with a nonlinear constrained, a linearizing method is used to convert the model to a linear least squares estimator with a nonlinear con- strained. The Gauss-Newton iteration method is developed to conquer the source localization problem. From the analysis of solving Lagrange multiplier, the algorithm is a generalization of linear-correction least squares estimation procedure under the condition of geolocation using FDOA measurements only. The algorithm is compared with common least squares estimation. Comparisons of their estimation accuracy and the CRLB are made, and the proposed method attains the CRLB. Simulation re- sults are included to corroborate the theoretical development.
文摘This study discusses the challenges in logging the evaluation of low-resistivity oil reservoirs,especially the difficult problems involving their saturation calculation.A correction method for equivalent water conductivity is proposed,given the high conductivity caused by small amounts of water distributed in a network within the low-resistivity reservoir,which mimics the eff ects of high water saturation.This approach signifi cantly improves the accuracy of hydrocarbon saturation calculations in these low-resistivity reservoirs.The corrected hydrocarbon saturation values highly matched the porosity and are consistent with experimental results.This study also establishes a discrimination process to determine whether corrections are required,verifying the eff ectiveness and accuracy of the method through an application example.
基金Supported by the National Basic Research Program of China ("973"Program)(2006CB601201)~~
文摘In high-resolution cone-beam computed tomography （CBCT） using the flat-panel detector, imperfect or defect detector elements cause ring artifacts due to the none-uniformity of their X-ray response. They often disturb the image quality. A dedicated fitting correction method for high-resolution micro-CT is presented. The method converts each elementary X-ray response curve to an average one, and eliminates response inconsistency among pixels. Other factors of the method are discussed, such as the correction factor variability by different sampling frames and nonlinear factors over the whole spectrum. Results show that the noise and artifacts are both reduced in reconstructed images
文摘The efficacy of error correction and various kinds of correction approaches is one of the key issues in second language writing faced by both teachers and researchers. The current paper reviews the definition of error correction and examines the different views on whether error correction in L2 writing should be corrected. In particular, the paper discusses and analyses the three common correction methods: direct correction, peer feedback and indirect correction. Teachers are encouraged to weigh and analyze the advantages and disadvantages of these methods according to the current literature, employ the most beneficial error correction method in L2 writing, and adapt its suitability to their teaching context.
文摘Based on the hindcast results of summer rainfall anomalies over China for the period 1981-2000 by the Dynamical Climate Prediction System （IAP-DCP） developed by the Institute of Atmospheric Physics, a correction method that can account for the dependence of model＇s systematic biases on SST anomalies is proposed. It is shown that this correction method can improve the hindcast skill of the IAP-DCP for summer rainfall anomalies over China, especially in western China and southeast China, which may imply its potential application to real-time seasonal prediction.
基金supported by Key Project of National Natural Science Foundation of China (Grant No. 10932003)National Hi-tech Research and Development Program of China (863 Program,Grant No.2009AA04Z101)National Basic Research Program of China (973Program,Grant No. 2010CB832700)
文摘Surface low/distortion is one of the most challenging surface deflections that have a great effect on the exterior appearance of automobiles.Most studies on surface distortion/deflection have focused on evaluation and visualization techniques,the research work on correction or prevention of surface low/distortion is limited,and there is no perfect surface low/distortion corrective method that can satisfy the needs of the engineering.A B-spline based geometry morphing algorithm is proposed and then a new program based on UG-NX platform is developed to modify the die face in the surface low/distortion areas.To verify this developed system,the experimental dies that can replicate the surface low/distortion phenomenon successfully is put to use.Five geometric variables are introduced to describe the basic geometry of typical depression features of automotive outer panels.The experimental dies are then designed to reflect various combinations of these geometric parameters.The stamping experiments are conducted on cold rolled grade 5（CR5） sheet steel and various static measurements,such as oil-stoning,laser scanner,etc,are performed to measure and record the surface low/distortions.Three approaches including good bearing,holds in blank and die face morphing that aim to correct low/distortions are tried out and surface low/distortions are observed in the specimen with reverse draw depth of 10 mm.The measurement results show that die morphing is a practical and effective method to correct the surface low/distortion.The correction method proposed can be used to minimize the occurrence of surface low/distortion in die manufacturing,which has certain reference significance to the correction of surface low/distortion.
基金Project supported by the National Natural Science Foundation of China （Grant Nos 40575036 and 40325015）.Acknowledgement The authors thank Drs Zhang Pei-Qun and Bao Ming very much for their valuable comments on the present paper.
文摘In this paper, an analogue correction method of errors （ACE） based on a complicated atmospheric model is further developed and applied to numerical weather prediction （NWP）. The analysis shows that the ACE can effectively reduce model errors by combining the statistical analogue method with the dynamical model together in order that the information of plenty of historical data is utilized in the current complicated NWP model, Furthermore, in the ACE, the differences of the similarities between different historical analogues and the current initial state are considered as the weights for estimating model errors. The results of daily, decad and monthly prediction experiments on a complicated T63 atmospheric model show that the performance of the ACE by correcting model errors based on the estimation of the errors of 4 historical analogue predictions is not only better than that of the scheme of only introducing the correction of the errors of every single analogue prediction, but is also better than that of the T63 model.
基金supported by the National Natural Science Foundation of China (no. 51478120)
文摘The blind-hole method is the most widely used approach to experimentally determine the distribution of residual stress. This paper aims to improve test accuracy of welding residual stress and conducts an experimental study on the strain release factors involved when using the blind-hole method for Q235 and Q345, two steels commonly used in building structures. The ranges of strain release factors A and B in the elastic stage, the effects of strain release factors on residual stress calculated values, and the plastic corrected strain release factors are analyzed considering of the effect of plastic deformation around the blind hole on measurement accuracy. Finally, a simplified calculation formula to determine strain release factors is proposed for use with the blind-hole method. Results show that in the elastic stage, strain release factor A for Q235 and Q345 ranges from-0.399 to-0.525 and strain release factor B from-0.791 to-0.960. Changing the strain release factors A and B shows that calculated residual tensile stress varies in relation to a decrease in both factor values. However, there is a increase in calculated residual compressive stress with a decrease in the strain release factor A value, but there is an decrease with a decrease in strain release factor B value. Calculated residual stress applied to elastic strain release factors is compared with that applied to amended plastic strain release factors for Q235 steel. The maximum deviation between calculated residual stress and test stress is reduced from 21.1 to 1.0%,and for Q345 steel from 26.5 to 1.2%. It is thus evident that the plastic correction formula proposed in this paper can be used in calculations when conducting a residual stress test.
基金National Natural Science Foundation under Grant No. 51179093National Basic Research Program of China under Grant No. 2011CB013602Program for New Century Excellent Talents in University under Grant No.NCET-10-0531
文摘A third-order correction was recently suggested to improve the accuracy of the half-power bandwidth method in estimating the damping of single DOF systems.This paper analyzes the accuracy of the half-power bandwidth method with the third-order correction in damping estimation for multi-DOF linear systems.Damping ratios in a two-DOF linear system are estimated using its displacement and acceleration frequency response curves,respectively.A wide range of important parameters that characterize the shape of these response curves are taken into account.Results show that the third-order correction may greatly improve the accuracy of the half-power bandwidth method in estimating damping in a two-DOF system.In spite of this,the half-power bandwidth method may significantly overestimate the damping ratios of two-DOF systems in some cases.
基金This project is supported by Research Foundation for Doctoral Program of Higher Education, China (No.98033532)
文摘The main purpose of reverse engineering is to convert discrete data pointsinto piecewise smooth, continuous surface models. Before carrying out model reconstruction it issignificant to extract geometric features because the quality of modeling greatly depends on therepresentation of features. Some fitting techniques of natural quadric surfaces with least-squaresmethod are described. And these techniques can be directly used to extract quadric surfaces featuresduring the process of segmentation for point cloud.
文摘The Galerkin and least-squares methods are two classes of the most popular Krylov subspace methOds for solving large linear systems of equations. Unfortunately, both the methods may suffer from serious breakdowns of the same type: In a breakdown situation the Galerkin method is unable to calculate an approximate solution, while the least-squares method, although does not really break down, is unsucessful in reducing the norm of its residual. In this paper we first establish a unified theorem which gives a relationship between breakdowns in the two methods. We further illustrate theoretically and experimentally that if the coefficient matrix of a lienar system is of high defectiveness with the associated eigenvalues less than 1, then the restarted Galerkin and least-squares methods will be in great risks of complete breakdowns. It appears that our findings may help to understand phenomena observed practically and to derive treatments for breakdowns of this type.
基金supported by a project from the Youth Science Foundation of the National Natural Science Foundation of China (11104089)
文摘Relaxation time spectra （RTS） derived from time domain induced polarization data （TDIP） are helpful to assess oil reservoir pore structures. However, due to the sensitivity to the signal-to-noise ratio （SNR）, the inversion accuracy of the traditional singular value decomposition （SVD） inversion method reduces with a decrease of SNR. In order to enhance the inversion accuracy and improve robustness of the inversion method to the SNR, an improved inversion method, based on damping factor and spectrum component residual correction, is proposed in this study. The numerical inversion results show that the oscillation of the RTS derived from the SVD method increased with a decrease of SNR, which makes it impossible to get accurate inversion components. However, the SNR has little influence on inversion components of the improved method, and the RTS has high inversion accuracy and robustness. Moreover, RTS derived from core sample data is basically in accord with the pore-size distribution curve, and the RTS derived from the actual induced polarization logging data is smooth and continuous, which indicates that the improved method is practicable.
基金the National Science Council ot Taiwan,China for funding this research(Project no.:NSC 94-2218-E-035-011)
文摘A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.
基金Projects(52075552,51575533,51805555,11662004)supported by the National Natural Science Foundation of China。
文摘In this work,synchronous cutting of concave and convex surfaces was achieved using the duplex helical method for the hypoid gear,and the problem of tooth surface error correction was studied.First,the mathematical model of the hypoid gears machined by the duplex helical method was established.Second,the coordinates of discrete points on the tooth surface were obtained by measurement center,and the normal errors of the discrete points were calculated.Third,a tooth surface error correction model is established,and the tooth surface error was corrected using the Levenberg-Marquard algorithm with trust region strategy and least square method.Finally,grinding experiments were carried out on the machining parameters obtained by Levenberg-Marquard algorithm with trust region strategy,which had a better effect on tooth surface error correction than the least square method.After the tooth surface error is corrected,the maximum absolute error is reduced from 30.9μm before correction to 6.8μm,the root mean square of the concave error is reduced from 15.1 to 2.1μm,the root mean square of the convex error is reduced from 10.8 to 1.8μm,and the sum of squared errors of the concave and convex surfaces was reduced from 15471 to 358μm^(2).It is verified that the Levenberg-Marquard algorithm with trust region strategy has a good accuracy for the tooth surface error correction of hypoid gear machined by duplex helical method.
文摘This is a study of the Durand-Kerner and Nourein methods for finding the roots of a given algebraic equation simultaneously. We consider the conditions under which the iterative methods fail. The numerical example is presented.
基金the National Science Council of Taiwan for funding this research (NSC 96-2221-E-019-061).
文摘Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.
基金supported by the National Basic Research Program of China (2005CB321701)NSF of mathematics research special fund of Hebei Province(08M005)
文摘The purpose of this article is to develop and analyze least-squares approximations for the incompressible magnetohydrodynamic equations. The major advantage of the least-squares finite element method is that it is not subjected to the so-called Ladyzhenskaya-Babuska-Brezzi （LBB） condition. The authors employ least-squares functionals which involve a discrete inner product which is related to the inner product in H^-1（Ω）.
基金partially supported by the National Natural Science Foundation of China (No.41230318)
文摘With the development of computational power, there has been an increased focus on data-fitting related seismic inversion techniques for high fidelity seismic velocity model and image, such as full-waveform inversion and least squares migration. However, though more advanced than conventional methods, these data fitting methods can be very expensive in terms of computational cost. Recently, various techniques to optimize these data-fitting seismic inversion problems have been implemented to cater for the industrial need for much improved efficiency. In this study, we propose a general stochastic conjugate gradient method for these data-fitting related inverse problems. We first prescribe the basic theory of our method and then give synthetic examples. Our numerical experiments illustrate the potential of this method for large-size seismic inversion application.
基金supported by the National Natural Science Foundation of China(Nos.11702329,12102247)the Shanghai Key Lab of Vehicle Aerodynamics and Vehicle Thermal Management Systems,China(No.VATLAB-2021-01)。
文摘For the second-order finite volume method,implicit schemes and reconstruction methods are two main algorithms which influence the robustness and efficiency of the numerical simulations of compressible turbulent flows.In this paper,a compact least-squares reconstruction method is proposed to calculate the gradients for the distribution of flow field variables approximation.The compactness of the new reconstruction method is reflected in the gradient calculation process.The geometries of the face-neighboring elements are no longer utilized,and the weighted average values at the centroid of the interfaces are used to calculate the gradients instead of the values at the centroid of the face-neighboring elements.Meanwhile,an exact Jacobian solving strategy is developed for implicit temporal discretization.The accurate processing of Jacobian matrix can extensively improve the invertibility of the Jacobian matrix and avoid introducing extra numerical errors.In addition,a modified Venkatakrishnan limiter is applied to deal with the shock which may appear in transonic flows and the applicability of the mentioned methods is enhanced further.The combination of the proposed methods makes the numerical simulations of turbulent flow converge rapidly and steadily with an adaptive increasing CFL number.The numerical results of several benchmarks indicate that the proposed methods perform well in terms of robustness,efficiency and accuracy,and have good application potential in turbulent flow simulations of complex configurations.
