In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering comp...In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.展开更多
Wire arc additive manufacture(WAAM) is a new technique to fabricate large-scale complex aluminum alloy components.However, the performance of the parts is critically influenced by residual stresses and deformation. A ...Wire arc additive manufacture(WAAM) is a new technique to fabricate large-scale complex aluminum alloy components.However, the performance of the parts is critically influenced by residual stresses and deformation. A sequentially thermal-mechanical coupled model of residual stress and deformation for aluminum alloy WAAM parts was established based on commercial FE software ABAQUS. The temperature field was calculated by the moving heat source(MHS) method. The temperature function was obtained according to the distribution of the peak temperature. Furthermore, the MHS method and segmented temperature function(STF) method were used to calculate the residual stress and deformation. The results show that the STF method satisfies both the efficiency and accuracy requirements. 1-segment, 3-segment, and 5-segment STF methods can shorten the time for mechanical analysis by 91%, 79%, 63%, respectively.The error of the residual stress and deformation are all less than 20%. STF method provides an effective way to predict the residual stress and deformation of large-scale WAAM parts.展开更多
With the applications of Nd-Fe-B material extending in recent years, the materials of neodymium metal and other rare earth metal alloy confront the increased demand and the high quality request at the same time.These ...With the applications of Nd-Fe-B material extending in recent years, the materials of neodymium metal and other rare earth metal alloy confront the increased demand and the high quality request at the same time.These factors stimulated greatly to perfect the producing craft of RE metals and improve the equipments.The rare earth electrolysis cell is developing towards large-scale way.Notwithstanding the present electrolysis cell of Nd metal, include 6 kA and 10 kA cell, exists some insurmountable problems during operation and these problems lead to lower electric efficiency and higher operating costs.So it is significant to study the physical fields of rare earth electrolysis cell.In this paper,a numerical flow mode is established using vortex- flowing function method and the fluid flow field of 3000A Nd electrolysis cell is computed using MATLAB.The results of the study will be important reference in theory for improving and enlarging rare earth fluoride system cell.展开更多
A numerical model for determining fishing net configuration and load is developed, based on the lumped masses method. The model is used to analyze the influences on the net of gravity sea cage by different sinker weig...A numerical model for determining fishing net configuration and load is developed, based on the lumped masses method. The model is used to analyze the influences on the net of gravity sea cage by different sinker weights in uniform flow. The net of gravity sea cage is simulated under different current velocities and sinker weights. In order to verify the validity of the numerical results, model test results made by Lader and Enerhaug are cited and compared with the numerical ones. The results of numerical simulation agree well with the experimental ones and the agreement is within an error range of 13% under different velocities and sinker weights. The numerical results indicate that increasing the sinker weight can effectively reduce the net deformation, but the total drag force on the net will increase accordingly. Results will give references for better knowledge of the hydrodynamic behavior of gravity sea cage.展开更多
In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through poro...In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through porous media. A gas dispersion equation in a 3-D field is calculated by use of numerical method on a weighted upstream multi-element balance. Based on data of an example with a U type ventilation mode,surface charts of air pressure distribution and gas concentration are drawn by Graphtool software. Finally,a comparison between actually measured results in the model test and the numerical simulation results is made to proves the numerical implementation feasible.展开更多
A numerical solution of the weight function for a two-electrode electromagnetic flowmeter was proposed. The solution was obtained by using the finite element method based on the basic equation of a traditional two-ele...A numerical solution of the weight function for a two-electrode electromagnetic flowmeter was proposed. The solution was obtained by using the finite element method based on the basic equation of a traditional two-electrode electromagnetic flowmeter. The two-dimensional distribution of the weight function of the electromagnetic flowmeter obtained was verified by the analytical solution. Three-dimensional distribution of the weight function was also presented in the paper. It can be employed to analyze the sensitivity and linearity of the electromagnetic flowmeter with non-uniform magnetic field, and even to assist the design of the excitation coil pair.展开更多
Due to the difficulties associated with preprocessing activities and poor grid convergence when simulating shale reservoirs in the context of traditional grid methods,in this study an innovative two-phase oil-water se...Due to the difficulties associated with preprocessing activities and poor grid convergence when simulating shale reservoirs in the context of traditional grid methods,in this study an innovative two-phase oil-water seepage model is elaborated.The modes is based on the radial basis meshless approach and is used to determine the pressure and water saturation in a sample reservoir.Two-dimensional examples demonstrate that,when compared to the finite difference method,the radial basis function method produces less errors and is more accurate in predicting daily oil production.The radial basis function and finite difference methods provide errors of 5.78 percent and 7.5 percent,respectively,when estimating the daily oil production data for a sample well.A sensitivity analysis of the key parameters that affect the radial basis function’s computation outcomes is also presented.展开更多
In this paper,a simple direct space-time semi-analytical meshless scheme is proposed for the numerical approximation of the coupled Burgers'equations.During the whole solution procedure,two different schemes are c...In this paper,a simple direct space-time semi-analytical meshless scheme is proposed for the numerical approximation of the coupled Burgers'equations.During the whole solution procedure,two different schemes are considered in terms of radial and non-radial basis functions.The time-dependent variable in the first radial scheme is directly considered as the normal space variables to formulate an"isotropic"space-time radial basis function.The second non-radial scheme considered relationship between time-dependent and spacedependent variables.Under such circumstance,we can get a one-step space-time meshless scheme.The numerical findings demonstrate that the proposed meshless schemes are precise,user-friendly,and effective in solving the coupled Burgers'equations.展开更多
This study investigates the dynamics of pneumococcal pneumonia using a novel fractal-fractional Susceptible-Carrier-Infected-Recovered model formulated with the Atangana-Baleanu in Caputo(ABC)sense.Unlike traditional ...This study investigates the dynamics of pneumococcal pneumonia using a novel fractal-fractional Susceptible-Carrier-Infected-Recovered model formulated with the Atangana-Baleanu in Caputo(ABC)sense.Unlike traditional epidemiological models that rely on classical or Caputo fractional derivatives,the proposed model incorporates nonlocal memory effects,hereditary properties,and complex transmission dynamics through fractalfractional calculus.The Atangana-Baleanu operator,with its non-singular Mittag-Leffler kernel,ensures a more realistic representation of disease progression compared to classical integer-order models and singular kernel-based fractional models.The study establishes the existence and uniqueness of the proposed system and conducts a comprehensive stability analysis,including local and global stability.Furthermore,numerical simulations illustrate the effectiveness of the ABC operator in capturing long-memory effects and nonlocal interactions in disease transmission.The results provide valuable insights into public health interventions,particularly in optimizing vaccination strategies,treatment approaches,and mitigation measures.By extending epidemiological modeling through fractal-fractional derivatives,this study offers an advanced framework for analyzing infectious disease dynamics with enhanced accuracy and predictive capabilities.展开更多
Sensor management plays an important role in data fusion system, and this paper presents an algorithm of sensor management that can be used in target detection, identification and tracking. First, the basic concept, r...Sensor management plays an important role in data fusion system, and this paper presents an algorithm of sensor management that can be used in target detection, identification and tracking. First, the basic concept, rule, range and function of sensor management are introduced. Then, the quantifying problems of target priority and sensor (or combination)-target pairing in multisensor management are discussed and the efficiency and waste functions are established based on the functions of target priority and sensor-target pairing. On this basis, a distribution algorithm of multi-sensor resources is given, which is optimized by the principle of maximum synthesis efficiency in the multisensor system and constrained by sensor maximum tracking power and what target must be scanned. In addition, the waste measure of sensor resources is introduced to improve the algorithm. Finally, a tactical task that includes three sensors and ten targets is set, and the simulation results show that the algorithm is feasible and effective.展开更多
Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved met...Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded interface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.展开更多
(2 + 1) dimensional Boussinesq and Kadomtsev-Petviashvili equation are investigated by employing Jacobi elliptic function expansion method in this paper. As a result, some new forms traveling wave solutions of the equ...(2 + 1) dimensional Boussinesq and Kadomtsev-Petviashvili equation are investigated by employing Jacobi elliptic function expansion method in this paper. As a result, some new forms traveling wave solutions of the equation are reported. Numerical simulation results are shown. These new solutions may be important for the explanation of some practical physical problems. The results of this paper show that Jacobi elliptic function method can be a useful tool in obtaining evolution solutions of nonlinear system.展开更多
Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved met...Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded in- terface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.展开更多
The element-free method is a new numerical technique presented in recent years.It uses the moving least square(MLS) approximation as its shape function,and it is determined by the basic function and weight function.Th...The element-free method is a new numerical technique presented in recent years.It uses the moving least square(MLS) approximation as its shape function,and it is determined by the basic function and weight function.The weight function is the mainly determining factor,so it greatly affects the accuracy of the computational results.The element-free Galerkin method(EFGM) was applied for the solution to plastic large deformation.The simulation of metal rheological forming was successfully done by programming and its results were visualized by using the plotting and data analyses software Tecplot.Then plastic strain under different stages during rheological forming and the three principal stresses at the last deformation were obtained.The example shows the feasibility of EFGM used for metal rheological forming and provides a new method for numerical simulation of rheological forming of complex parts.展开更多
Two computational cases that have analytic solutions are employed for studying the adaptive grid tech-nique based on the variational principle. The results show that for the computational case of traveling shock waves...Two computational cases that have analytic solutions are employed for studying the adaptive grid tech-nique based on the variational principle. The results show that for the computational case of traveling shock waves the weight function, with the 2nd-order derivation terms taken into consideration, can more effectively reduce the error than one with gradient terms. For the case of cyclonic frontogenesis, weight func-tions only related to the gradient are unable to enhance the computational accuracy while ones with the wind field and frontogenesis function taken into consideration can more reasonably arrange the grid. Com-pared with analytic solutions, the adaptive grid technique suggested in this paper can improve computational accuracy and it displays the prominent advantage of saving memory.展开更多
文摘In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.
基金supported by the National Key Technologies R&D Program (Grant No. 2018YFB1106000)Innovation Funds of China Academy of Launch Vehicle Technology (CALT) for Universities (Grant No.CALT201709)Tsinghua Grants for Autonomous Research。
文摘Wire arc additive manufacture(WAAM) is a new technique to fabricate large-scale complex aluminum alloy components.However, the performance of the parts is critically influenced by residual stresses and deformation. A sequentially thermal-mechanical coupled model of residual stress and deformation for aluminum alloy WAAM parts was established based on commercial FE software ABAQUS. The temperature field was calculated by the moving heat source(MHS) method. The temperature function was obtained according to the distribution of the peak temperature. Furthermore, the MHS method and segmented temperature function(STF) method were used to calculate the residual stress and deformation. The results show that the STF method satisfies both the efficiency and accuracy requirements. 1-segment, 3-segment, and 5-segment STF methods can shorten the time for mechanical analysis by 91%, 79%, 63%, respectively.The error of the residual stress and deformation are all less than 20%. STF method provides an effective way to predict the residual stress and deformation of large-scale WAAM parts.
文摘With the applications of Nd-Fe-B material extending in recent years, the materials of neodymium metal and other rare earth metal alloy confront the increased demand and the high quality request at the same time.These factors stimulated greatly to perfect the producing craft of RE metals and improve the equipments.The rare earth electrolysis cell is developing towards large-scale way.Notwithstanding the present electrolysis cell of Nd metal, include 6 kA and 10 kA cell, exists some insurmountable problems during operation and these problems lead to lower electric efficiency and higher operating costs.So it is significant to study the physical fields of rare earth electrolysis cell.In this paper,a numerical flow mode is established using vortex- flowing function method and the fluid flow field of 3000A Nd electrolysis cell is computed using MATLAB.The results of the study will be important reference in theory for improving and enlarging rare earth fluoride system cell.
基金This study was financially supported by the National Natural Science Foundation of China under contract No. 50279002the National "863" High Technology Project of China under contract Nos 2001AA623010 and 2003AA623010the Program for Changjiang Scholars and Innovative Research Team in University of China under contract No. IRT0420.
文摘A numerical model for determining fishing net configuration and load is developed, based on the lumped masses method. The model is used to analyze the influences on the net of gravity sea cage by different sinker weights in uniform flow. The net of gravity sea cage is simulated under different current velocities and sinker weights. In order to verify the validity of the numerical results, model test results made by Lader and Enerhaug are cited and compared with the numerical ones. The results of numerical simulation agree well with the experimental ones and the agreement is within an error range of 13% under different velocities and sinker weights. The numerical results indicate that increasing the sinker weight can effectively reduce the net deformation, but the total drag force on the net will increase accordingly. Results will give references for better knowledge of the hydrodynamic behavior of gravity sea cage.
基金Projects 50534090 and 50674090 supported by the National Natural Science Foundation of China2005CB221503 by the National Key Basic ResearchDevelopment Program (973 Program)
文摘In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through porous media. A gas dispersion equation in a 3-D field is calculated by use of numerical method on a weighted upstream multi-element balance. Based on data of an example with a U type ventilation mode,surface charts of air pressure distribution and gas concentration are drawn by Graphtool software. Finally,a comparison between actually measured results in the model test and the numerical simulation results is made to proves the numerical implementation feasible.
基金Supported by Chinese National Programs for High Technology Research and Development(2008AA042207)
文摘A numerical solution of the weight function for a two-electrode electromagnetic flowmeter was proposed. The solution was obtained by using the finite element method based on the basic equation of a traditional two-electrode electromagnetic flowmeter. The two-dimensional distribution of the weight function of the electromagnetic flowmeter obtained was verified by the analytical solution. Three-dimensional distribution of the weight function was also presented in the paper. It can be employed to analyze the sensitivity and linearity of the electromagnetic flowmeter with non-uniform magnetic field, and even to assist the design of the excitation coil pair.
基金supported by The China Postdoctoral Science Foundation(2021M702304)Natural Science Foundation of Shandong Province(ZR2021QE260).
文摘Due to the difficulties associated with preprocessing activities and poor grid convergence when simulating shale reservoirs in the context of traditional grid methods,in this study an innovative two-phase oil-water seepage model is elaborated.The modes is based on the radial basis meshless approach and is used to determine the pressure and water saturation in a sample reservoir.Two-dimensional examples demonstrate that,when compared to the finite difference method,the radial basis function method produces less errors and is more accurate in predicting daily oil production.The radial basis function and finite difference methods provide errors of 5.78 percent and 7.5 percent,respectively,when estimating the daily oil production data for a sample well.A sensitivity analysis of the key parameters that affect the radial basis function’s computation outcomes is also presented.
基金the Science and Technology Research Project of Henan Province (242102231052)the Key Scientific Research Plan of Colleges and Universities in Henan Province (23B140006)the Natural Science Foundation of Jiangxi Province (20224BAB201018)。
文摘In this paper,a simple direct space-time semi-analytical meshless scheme is proposed for the numerical approximation of the coupled Burgers'equations.During the whole solution procedure,two different schemes are considered in terms of radial and non-radial basis functions.The time-dependent variable in the first radial scheme is directly considered as the normal space variables to formulate an"isotropic"space-time radial basis function.The second non-radial scheme considered relationship between time-dependent and spacedependent variables.Under such circumstance,we can get a one-step space-time meshless scheme.The numerical findings demonstrate that the proposed meshless schemes are precise,user-friendly,and effective in solving the coupled Burgers'equations.
基金funded by the Research,Development,and Innovation Authority(RDIA)-Kingdom of Saudi Arabia-with grant number 12803-baha-2023-BU-R-3-1-EI.
文摘This study investigates the dynamics of pneumococcal pneumonia using a novel fractal-fractional Susceptible-Carrier-Infected-Recovered model formulated with the Atangana-Baleanu in Caputo(ABC)sense.Unlike traditional epidemiological models that rely on classical or Caputo fractional derivatives,the proposed model incorporates nonlocal memory effects,hereditary properties,and complex transmission dynamics through fractalfractional calculus.The Atangana-Baleanu operator,with its non-singular Mittag-Leffler kernel,ensures a more realistic representation of disease progression compared to classical integer-order models and singular kernel-based fractional models.The study establishes the existence and uniqueness of the proposed system and conducts a comprehensive stability analysis,including local and global stability.Furthermore,numerical simulations illustrate the effectiveness of the ABC operator in capturing long-memory effects and nonlocal interactions in disease transmission.The results provide valuable insights into public health interventions,particularly in optimizing vaccination strategies,treatment approaches,and mitigation measures.By extending epidemiological modeling through fractal-fractional derivatives,this study offers an advanced framework for analyzing infectious disease dynamics with enhanced accuracy and predictive capabilities.
基金National Natural Science Foundation of China! ( No.6 9772 0 3 1)Natural Science Foundation of Henan Province
文摘Sensor management plays an important role in data fusion system, and this paper presents an algorithm of sensor management that can be used in target detection, identification and tracking. First, the basic concept, rule, range and function of sensor management are introduced. Then, the quantifying problems of target priority and sensor (or combination)-target pairing in multisensor management are discussed and the efficiency and waste functions are established based on the functions of target priority and sensor-target pairing. On this basis, a distribution algorithm of multi-sensor resources is given, which is optimized by the principle of maximum synthesis efficiency in the multisensor system and constrained by sensor maximum tracking power and what target must be scanned. In addition, the waste measure of sensor resources is introduced to improve the algorithm. Finally, a tactical task that includes three sensors and ten targets is set, and the simulation results show that the algorithm is feasible and effective.
基金Supported by the National Natural Science Foundation of China (No.20490206) and the Special Funds for Major State Basic Research Program of China (973 Program, 2004CB217604).
文摘Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded interface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.
文摘(2 + 1) dimensional Boussinesq and Kadomtsev-Petviashvili equation are investigated by employing Jacobi elliptic function expansion method in this paper. As a result, some new forms traveling wave solutions of the equation are reported. Numerical simulation results are shown. These new solutions may be important for the explanation of some practical physical problems. The results of this paper show that Jacobi elliptic function method can be a useful tool in obtaining evolution solutions of nonlinear system.
基金the National Natural Science Foundation of China (No.20490206) the Special Funds for Major State BasicResearch Program of China (973 Program, 2004CB217604).
文摘Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded in- terface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.
基金Key project(02103) supported by National Education Department of ChinaKey project(02A008) supported by the Education Department of Hunan Province,China+3 种基金Project(2005090) supported by Central South University of Forestry and TechnologyProject(03JJY3007) supported by the Natural Science Foundation of Hunan Province,ChinaProject supported by the Rewarding Project for Excellent PhD Thesis of Hunan Province,ChinaProject(07031B) supported by Scientific Research Fund of Central South University of Forestry and Technology
文摘The element-free method is a new numerical technique presented in recent years.It uses the moving least square(MLS) approximation as its shape function,and it is determined by the basic function and weight function.The weight function is the mainly determining factor,so it greatly affects the accuracy of the computational results.The element-free Galerkin method(EFGM) was applied for the solution to plastic large deformation.The simulation of metal rheological forming was successfully done by programming and its results were visualized by using the plotting and data analyses software Tecplot.Then plastic strain under different stages during rheological forming and the three principal stresses at the last deformation were obtained.The example shows the feasibility of EFGM used for metal rheological forming and provides a new method for numerical simulation of rheological forming of complex parts.
基金Acknowledgments. The research report has been supported jointly by the National Natural Science Foundation of China under Grant Nos.40075024 and 49945009, and by the National Key Basic Research and Development Proj-ect under Grant No. G1998040911.
文摘Two computational cases that have analytic solutions are employed for studying the adaptive grid tech-nique based on the variational principle. The results show that for the computational case of traveling shock waves the weight function, with the 2nd-order derivation terms taken into consideration, can more effectively reduce the error than one with gradient terms. For the case of cyclonic frontogenesis, weight func-tions only related to the gradient are unable to enhance the computational accuracy while ones with the wind field and frontogenesis function taken into consideration can more reasonably arrange the grid. Com-pared with analytic solutions, the adaptive grid technique suggested in this paper can improve computational accuracy and it displays the prominent advantage of saving memory.