The properties of the magnetic mold in magnetic mold casting directly determine the quality of the final cast parts.In this study,the magnetic mold properties in magnetic mold casting,were studied utilizing a coupled ...The properties of the magnetic mold in magnetic mold casting directly determine the quality of the final cast parts.In this study,the magnetic mold properties in magnetic mold casting,were studied utilizing a coupled electromagnetic-structural method through numerical simulation.This study investigated key factors including equivalent stress,the distribution of tensile and compressive stresses,and the area ratio of tensile stress.It compared molds made entirely of magnetic materials with those made partially of magnetic materials.Simulation results indicate that as current increases from 4 A to 8 A,both the initial magnetic mold and the material-replaced magnetic mold initially show an increasing trend in equivalent stress,tensile-compressive stress,and the area ratio of tensile stress,peaking at 6 A before declining.After material replacement,the area ratio of tensile stress at 6 A decreases to 19.84%,representing a reduction of 29.72%.Magnetic molds comprising a combination of magnetic and non-magnetic materials exhibit sufficient strength and a reduced area ratio of tensile stress compared to those made entirely from magnetic materials.This study provides valuable insights for optimizing magnetic mold casting processes and offers practical guidance for advancing the application of magnetic molds.展开更多
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.展开更多
In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement wa...In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement was simulated based on actual meteorological data of Nanjing. 24-hour rutting development under a transient temperature field was calculated in each month. The rutting depth accumulated under the static temperature field was also estimated and the relationship between constant temperature parameters was analyzed. Then the effective temperature for pavement rutting was determined based on the rutting equivalence principle. The results show that the monthly effective temperature is above 40 t in July and August, while in June and September it ranges from 30 to 40 Rutting development can be ignored when the monthly effective temperature is less than 30 t. The yearly effective temperature for rutting in Nanjing is around 38. 5 t. The long-term rutting prediction model based on the effective temperature can reflect the influences of meteorological factors and traffic time distribution.展开更多
Loess soils are widely distributed worldwide and typical in northwest China,and excessive agricultural irrigation has caused landslides in the area,specifically in the Heifangtai loess region in Lanzhou,Gansu,China.Ge...Loess soils are widely distributed worldwide and typical in northwest China,and excessive agricultural irrigation has caused landslides in the area,specifically in the Heifangtai loess region in Lanzhou,Gansu,China.Geophysical exploration is an essential method in landslide engineering geological surveys,and geological surveying,drilling,geophysical prospecting,monitoring,and other methods are used for performing engineering geological evaluation and obtaining comprehensive basic data for landslide protection design and construction.The theoretical feasibility of using geophysical methods in loess landslide detection is essential.On the basis of the shallow geological structure of the Heifangtai landslide region in Lanzhou,Gansu,China,a typical geoelectric model of the magnetotelluric method was established,and the loess landslide area was modeled through a two-dimensional finite element method,forward numerical simulation,and engineering geological analysis.The distribution characteristics of the magnetotelluric field were determined.This is a typical application of the geological process analysis method in geophysical exploration.This study provides the typical stratigraphic structure and electrical characteristics of different groundwater distributions in Heifangtai,Gansu,China,verifies the accuracy of forward modeling and calculation results,and provides a detailed theoretical basis for landslide detection through magnetotelluric methods.Through the numerical simulation of the forward modeling of the Heifangtai landslide region in Lanzhou,Gansu,China,this study can provide a detailed geophysical basis for landslide investigation,corroborate results of geological investigation and landslide design,and facilitate the sustainable development of agriculture in Heifangtai.展开更多
LS-DYNA program and the principle of ALE method were introduced, and the target features of the reinforced concrete penetration were analyzed by using the D material model and the ALE method. A numerical simulation ha...LS-DYNA program and the principle of ALE method were introduced, and the target features of the reinforced concrete penetration were analyzed by using the D material model and the ALE method. A numerical simulation has been done to show the penetration visually and veritably. The simulation results are analyzed carefully and explicitly prove their significance to the research of reinforced concrete penetration.展开更多
The modified Siemens method is the dominant process for the production of polysilicon,yet it is characterised by high energy consumption.Two models of laboratory-grade Siemens reduction furnace and 12 pairs of rods in...The modified Siemens method is the dominant process for the production of polysilicon,yet it is characterised by high energy consumption.Two models of laboratory-grade Siemens reduction furnace and 12 pairs of rods industrial-grade Siemens chemical vapor deposition(CVD)reduction furnace were established,and the effects of factors such as the diameter of silicon rods,the surface temperature of silicon rods,the air inlet velocity and temperature on the heat transfer process inside the reduction furnace were investigated by numerical simulation.The results show that the convective and radiant heat losses in the furnace increased with the diameter of the silicon rods.Furthermore,the radiant heat loss of the inner and outer rings of silicon rods was inconsistent for the industrial-grade reduction furnace.As the surface temperature of the silicon rods increases,the convective heat loss in the furnace increases,while the radiative heat loss remains relatively constant.When the inlet temperature and inlet velocity increase,the convective heat loss decreases,while the radiant heat loss remains relatively constant.Furthermore,the furnace wall surface emissivity increases,resulting in a significant increase in the amount of radiant heat loss in the furnace.In practice,this can be mitigated by polishing or adding coatings to reduce the furnace wall surface emissivity.展开更多
In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step ...In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.展开更多
Considering the limitation of computational capacity, a new finite element solution is used to simulate the welding deformation of the side sill of railroad car＇ s bogie frame based on the local-global method. Firstl...Considering the limitation of computational capacity, a new finite element solution is used to simulate the welding deformation of the side sill of railroad car＇ s bogie frame based on the local-global method. Firstly, a volumetric heat source defined by a double ellipsoid is adopted to simulate the thermal distributions of the arc welding process. And then, the local models extracted from the global model are computed with refined meshes. On these bases, the global distortions of the subject studied are ascertained by transferring the inner forces of computed local models to the global model. It indicates that the local-global method is feasible for simulating the large welded structures by comparing the computed results with the corresponding actual measured values. The work provides basis for optimizing the welding sequence and clamping conditions, and has theoretical values and engineering significance in the integral design, manufacturing technique selection of the bogie frame, as well as other kinds of large welded structures.展开更多
The mechanism of stress generation and propagation by detonation loading in five separate independent advance of ore breaking patterns is discussed in the paper. An elastic numerical model was developed using AN- SYS/...The mechanism of stress generation and propagation by detonation loading in five separate independent advance of ore breaking patterns is discussed in the paper. An elastic numerical model was developed using AN- SYS/LS-DYNA 3D Nonlinear Dynamic Finite Element Software. In this package ANSYS is the preprocessor and LS-DYNA is the postprocessor. Numerical models in the paper to actual were l：10 and the element mesh was dissected in scanning mode utilizing the symmetry characteristics of the numerical model. Five different advance rates were studied. Parameters, such as the time required to maximum stress, the action time of the available stress, the maximum velocity of the nodes, the stress penetration time, the magnitude of the stress peak and the time duration for high stress were numerically simulated. The 2.2 m advance appeared optimum from an analysis of the simulation results. The results from numerical simulation have been validated by tests with physical models.展开更多
In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models...In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with Mbeing a positive integer for regular grids are discussed and illustrated by constructing the second order （M= 1） and the fourth order （M = 2） recursion formulas.展开更多
Non-pillar continuous mining(NPCM) is regarded as a high-efficient, high-level and one-step mining technology, which can be divided into two substopes. Back fill stability status in substope I, which directly influenc...Non-pillar continuous mining(NPCM) is regarded as a high-efficient, high-level and one-step mining technology, which can be divided into two substopes. Back fill stability status in substope I, which directly influence the loss rate and dilution rate, etc, will determine whether the experimental research is successful or not. By employing energy method of limit analysis and finite element numerical simulation method, the critical backfill height was determined under the prerequisite condition of its stability, which put forward theoretical basis for reasonable and correct selection of backfill’s parameters. The result showed that the first backfill could not keep stable for NPCM, while the other was able to.展开更多
The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated includ...The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated including spherical, oblate and spherical. The bubble shape and drag coefficient were compared with experimental results. It is observed that the simulated results show good conformity to experimental results over a wide range of Reynolds number. In addition, the detailed flow field based on the reference coordinate system moving with the bubble is obtained, and the relationship among flow field, bubble shape and velocity is discussed.展开更多
The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method （LBM） was used to study the electrophysical property of rock...The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method （LBM） was used to study the electrophysical property of rock saturated with fluid because of its advantages over conventional numerical approaches in handling complex pore geometry and boundary conditions. The digital core model was constructed through the accumulation of matrix grains based on their radius distribution obtained by the measurements of core samples. The flow of electrical current through the core model saturated with oil and water was simulated on the mesoscopic scale to reveal the non-Archie relationship between resistivity index and water saturation （I-Sw）. The results from LBM simulation and laboratory measurements demonstrated that the I-Sw relation in the range of low water saturation was generally not a straight line in the log-log coordinates as described by the Archie equation. We thus developed a new equation based on numerical simulation and physical experiments. This new equation was used to fit the data from laboratory core measurements and previously published data. Determination of fluid saturation and reservoir evaluation could be significantly improved by using the new equation.展开更多
Numerical wave tanks are widely-acknowledged tools in studying waves and wave-structure interactions. They can generate waves under realistic scales and offers more information on the fluid field. However, most numeri...Numerical wave tanks are widely-acknowledged tools in studying waves and wave-structure interactions. They can generate waves under realistic scales and offers more information on the fluid field. However, most numerical wave tanks suffer from issues known as the numerical dissipation and numerical dispersion. The former causes wave energy to be slowly dissipated and the latter shifts wave frequencies during wave propagation. This paper proposes a simple method of depressing numerical dissipation effects on the basis of solving Euler equations using the finite difference method(FDM). The wave propagation solutions are solved analytically taking into account the influence of the damping terms. The main idea of the method is to append a source term to the momentum equation, whose strength is determined by how strong the numerical damping effect is. The method is verified by successfully depressing numerical effects during the simulation of regular linear waves, Stokes waves and irregular waves. By applying the method, wave energy is able to be close to its initial value after long distance of travel.展开更多
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.展开更多
Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid ...Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
Combined with the kinetic model of liquid film spreading, a new numerical method of solid-liquid-gas three-phase flow was developed for the moving of contact line, which was a hybrid method of computational fluid dyna...Combined with the kinetic model of liquid film spreading, a new numerical method of solid-liquid-gas three-phase flow was developed for the moving of contact line, which was a hybrid method of computational fluid dynamics and lattice Boltzmalm method （LBM）. By taking the effect of molecule force in droplet and the wall surface on liquid film into account, the changing law of contact angle with different surface tensions was analyzed on glass and aluminum foil surfaces. Compared with experimental results, the standard deviation by using LBM is less than 0.5°, which validates the feasibility of LBM simulation on the dynamic process of liquid film spreading. In addition, oscillations are discovered both at the initial and end phases. The phenomenon of retraction is also found and the maximum retraction angle is 7.58°. The obtained result shows that the retraction is proved to be correlative with precursor film by tracking the volume change of liquid film contour. Furthermore, non-dimensional coefficient 2 is introduced to measure the liquid film retraction capacity.展开更多
In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibi...In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.展开更多
Phase field method offers the prospect of being able to perform realistic numerical experiments on dendrite growth in a metallic system. In this paper, the equiaxed dendrite evolution during the solidification of a pu...Phase field method offers the prospect of being able to perform realistic numerical experiments on dendrite growth in a metallic system. In this paper, the equiaxed dendrite evolution during the solidification of a pure material was numerically simulated using the phase field model. The equiaxed dendrite growth in a two-dimensional square domain of undercooled melt (nickel) with four-fold anisotropy was simulated. The phase field model equations was solved using the explicit finite difference method on a uniform mesh. The formation of various equiaxed dendrite patterns was shown by a series of simulations, and the effect of anisotropy on equiaxed dendrite morphology was investigated.展开更多
In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform ...In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.展开更多
基金the National Natural Science Foundation of China(No.51875062,No.52205336)the China Postdoctoral Science Foundation(No.2021M700567).
文摘The properties of the magnetic mold in magnetic mold casting directly determine the quality of the final cast parts.In this study,the magnetic mold properties in magnetic mold casting,were studied utilizing a coupled electromagnetic-structural method through numerical simulation.This study investigated key factors including equivalent stress,the distribution of tensile and compressive stresses,and the area ratio of tensile stress.It compared molds made entirely of magnetic materials with those made partially of magnetic materials.Simulation results indicate that as current increases from 4 A to 8 A,both the initial magnetic mold and the material-replaced magnetic mold initially show an increasing trend in equivalent stress,tensile-compressive stress,and the area ratio of tensile stress,peaking at 6 A before declining.After material replacement,the area ratio of tensile stress at 6 A decreases to 19.84%,representing a reduction of 29.72%.Magnetic molds comprising a combination of magnetic and non-magnetic materials exhibit sufficient strength and a reduced area ratio of tensile stress compared to those made entirely from magnetic materials.This study provides valuable insights for optimizing magnetic mold casting processes and offers practical guidance for advancing the application of magnetic molds.
文摘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.
基金The National Natural Science Foundation of China(No.51378121)the Fok Ying Tung Education Foundation(No.141076)the Scientific Innovation Research of College Graduates in Jiangsu Province(No.KYLX_0164)
文摘In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement was simulated based on actual meteorological data of Nanjing. 24-hour rutting development under a transient temperature field was calculated in each month. The rutting depth accumulated under the static temperature field was also estimated and the relationship between constant temperature parameters was analyzed. Then the effective temperature for pavement rutting was determined based on the rutting equivalence principle. The results show that the monthly effective temperature is above 40 t in July and August, while in June and September it ranges from 30 to 40 Rutting development can be ignored when the monthly effective temperature is less than 30 t. The yearly effective temperature for rutting in Nanjing is around 38. 5 t. The long-term rutting prediction model based on the effective temperature can reflect the influences of meteorological factors and traffic time distribution.
文摘Loess soils are widely distributed worldwide and typical in northwest China,and excessive agricultural irrigation has caused landslides in the area,specifically in the Heifangtai loess region in Lanzhou,Gansu,China.Geophysical exploration is an essential method in landslide engineering geological surveys,and geological surveying,drilling,geophysical prospecting,monitoring,and other methods are used for performing engineering geological evaluation and obtaining comprehensive basic data for landslide protection design and construction.The theoretical feasibility of using geophysical methods in loess landslide detection is essential.On the basis of the shallow geological structure of the Heifangtai landslide region in Lanzhou,Gansu,China,a typical geoelectric model of the magnetotelluric method was established,and the loess landslide area was modeled through a two-dimensional finite element method,forward numerical simulation,and engineering geological analysis.The distribution characteristics of the magnetotelluric field were determined.This is a typical application of the geological process analysis method in geophysical exploration.This study provides the typical stratigraphic structure and electrical characteristics of different groundwater distributions in Heifangtai,Gansu,China,verifies the accuracy of forward modeling and calculation results,and provides a detailed theoretical basis for landslide detection through magnetotelluric methods.Through the numerical simulation of the forward modeling of the Heifangtai landslide region in Lanzhou,Gansu,China,this study can provide a detailed geophysical basis for landslide investigation,corroborate results of geological investigation and landslide design,and facilitate the sustainable development of agriculture in Heifangtai.
文摘LS-DYNA program and the principle of ALE method were introduced, and the target features of the reinforced concrete penetration were analyzed by using the D material model and the ALE method. A numerical simulation has been done to show the penetration visually and veritably. The simulation results are analyzed carefully and explicitly prove their significance to the research of reinforced concrete penetration.
基金funded by the Natural Science Foundation Projects in Sichuan Province(No.2022NSFSC0254).
文摘The modified Siemens method is the dominant process for the production of polysilicon,yet it is characterised by high energy consumption.Two models of laboratory-grade Siemens reduction furnace and 12 pairs of rods industrial-grade Siemens chemical vapor deposition(CVD)reduction furnace were established,and the effects of factors such as the diameter of silicon rods,the surface temperature of silicon rods,the air inlet velocity and temperature on the heat transfer process inside the reduction furnace were investigated by numerical simulation.The results show that the convective and radiant heat losses in the furnace increased with the diameter of the silicon rods.Furthermore,the radiant heat loss of the inner and outer rings of silicon rods was inconsistent for the industrial-grade reduction furnace.As the surface temperature of the silicon rods increases,the convective heat loss in the furnace increases,while the radiative heat loss remains relatively constant.When the inlet temperature and inlet velocity increase,the convective heat loss decreases,while the radiant heat loss remains relatively constant.Furthermore,the furnace wall surface emissivity increases,resulting in a significant increase in the amount of radiant heat loss in the furnace.In practice,this can be mitigated by polishing or adding coatings to reduce the furnace wall surface emissivity.
基金Project(41630642)supported by the Key Project of National Natural Science Foundation of ChinaProject(51974360)supported by the National Natural Science Foundation of ChinaProject(2018JJ3656)supported by the Natural Science Foundation of Hunan Province,China。
文摘In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.
文摘Considering the limitation of computational capacity, a new finite element solution is used to simulate the welding deformation of the side sill of railroad car＇ s bogie frame based on the local-global method. Firstly, a volumetric heat source defined by a double ellipsoid is adopted to simulate the thermal distributions of the arc welding process. And then, the local models extracted from the global model are computed with refined meshes. On these bases, the global distortions of the subject studied are ascertained by transferring the inner forces of computed local models to the global model. It indicates that the local-global method is feasible for simulating the large welded structures by comparing the computed results with the corresponding actual measured values. The work provides basis for optimizing the welding sequence and clamping conditions, and has theoretical values and engineering significance in the integral design, manufacturing technique selection of the bogie frame, as well as other kinds of large welded structures.
文摘The mechanism of stress generation and propagation by detonation loading in five separate independent advance of ore breaking patterns is discussed in the paper. An elastic numerical model was developed using AN- SYS/LS-DYNA 3D Nonlinear Dynamic Finite Element Software. In this package ANSYS is the preprocessor and LS-DYNA is the postprocessor. Numerical models in the paper to actual were l：10 and the element mesh was dissected in scanning mode utilizing the symmetry characteristics of the numerical model. Five different advance rates were studied. Parameters, such as the time required to maximum stress, the action time of the available stress, the maximum velocity of the nodes, the stress penetration time, the magnitude of the stress peak and the time duration for high stress were numerically simulated. The 2.2 m advance appeared optimum from an analysis of the simulation results. The results from numerical simulation have been validated by tests with physical models.
基金National Basic Research Program of China Under Grant No. 2007CB714200National Natural Science Foundation of China Under Grant No. 90715038
文摘In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with Mbeing a positive integer for regular grids are discussed and illustrated by constructing the second order （M= 1） and the fourth order （M = 2） recursion formulas.
文摘Non-pillar continuous mining(NPCM) is regarded as a high-efficient, high-level and one-step mining technology, which can be divided into two substopes. Back fill stability status in substope I, which directly influence the loss rate and dilution rate, etc, will determine whether the experimental research is successful or not. By employing energy method of limit analysis and finite element numerical simulation method, the critical backfill height was determined under the prerequisite condition of its stability, which put forward theoretical basis for reasonable and correct selection of backfill’s parameters. The result showed that the first backfill could not keep stable for NPCM, while the other was able to.
基金Project(21406141)supported by the National Natural Science Foundation of ChinaProject(20141078)supported by the Scientific Research Starting Foundation for Doctors of Liaoning Province,China+1 种基金Project(L2014060)supported by the Foundation of Department of Education of Liaoning Province,ChinaProject(157B21)supported by the Scientific Research Starting Foundation for Doctors of Shenyang Aerospace University,China
文摘The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated including spherical, oblate and spherical. The bubble shape and drag coefficient were compared with experimental results. It is observed that the simulated results show good conformity to experimental results over a wide range of Reynolds number. In addition, the detailed flow field based on the reference coordinate system moving with the bubble is obtained, and the relationship among flow field, bubble shape and velocity is discussed.
基金sponsored by the project No.50404001 from the National Natural Science Foundation of Chinathe National Key Fundamental Research & Development Project(Grant No.2007CB209601)+1 种基金the China National PetroleumCorporation Fundamental Research Program (Grant No.06A30102)the China Postdoctoral Science Foundation(Project No.2004035350)
文摘The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method （LBM） was used to study the electrophysical property of rock saturated with fluid because of its advantages over conventional numerical approaches in handling complex pore geometry and boundary conditions. The digital core model was constructed through the accumulation of matrix grains based on their radius distribution obtained by the measurements of core samples. The flow of electrical current through the core model saturated with oil and water was simulated on the mesoscopic scale to reveal the non-Archie relationship between resistivity index and water saturation （I-Sw）. The results from LBM simulation and laboratory measurements demonstrated that the I-Sw relation in the range of low water saturation was generally not a straight line in the log-log coordinates as described by the Archie equation. We thus developed a new equation based on numerical simulation and physical experiments. This new equation was used to fit the data from laboratory core measurements and previously published data. Determination of fluid saturation and reservoir evaluation could be significantly improved by using the new equation.
基金The National Natural Science Foundation of China under contract No.51609101 and 51909103the Natural Science Foundation of Fujian Province of China under contract Nos 2017J01701,2017J05085 and 2018J05090the Outstanding Young University Scientific Research Talents Cultivation Plan of Fujian Province of China
文摘Numerical wave tanks are widely-acknowledged tools in studying waves and wave-structure interactions. They can generate waves under realistic scales and offers more information on the fluid field. However, most numerical wave tanks suffer from issues known as the numerical dissipation and numerical dispersion. The former causes wave energy to be slowly dissipated and the latter shifts wave frequencies during wave propagation. This paper proposes a simple method of depressing numerical dissipation effects on the basis of solving Euler equations using the finite difference method(FDM). The wave propagation solutions are solved analytically taking into account the influence of the damping terms. The main idea of the method is to append a source term to the momentum equation, whose strength is determined by how strong the numerical damping effect is. The method is verified by successfully depressing numerical effects during the simulation of regular linear waves, Stokes waves and irregular waves. By applying the method, wave energy is able to be close to its initial value after long distance of travel.
基金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.
基金supported by the Natural Sciences and Engineering Research Council of Canada through Discovery Grant 341275 (G. Grasselli) and Engage EGP 461019-13
文摘Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
基金Project(U1261107)supported by the National Natural Science Foundation of China
文摘Combined with the kinetic model of liquid film spreading, a new numerical method of solid-liquid-gas three-phase flow was developed for the moving of contact line, which was a hybrid method of computational fluid dynamics and lattice Boltzmalm method （LBM）. By taking the effect of molecule force in droplet and the wall surface on liquid film into account, the changing law of contact angle with different surface tensions was analyzed on glass and aluminum foil surfaces. Compared with experimental results, the standard deviation by using LBM is less than 0.5°, which validates the feasibility of LBM simulation on the dynamic process of liquid film spreading. In addition, oscillations are discovered both at the initial and end phases. The phenomenon of retraction is also found and the maximum retraction angle is 7.58°. The obtained result shows that the retraction is proved to be correlative with precursor film by tracking the volume change of liquid film contour. Furthermore, non-dimensional coefficient 2 is introduced to measure the liquid film retraction capacity.
基金the National Natural Science Foundation of China(No.40774056)Program of Excellent Team in Harbin Institute of Technology
文摘In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.
文摘Phase field method offers the prospect of being able to perform realistic numerical experiments on dendrite growth in a metallic system. In this paper, the equiaxed dendrite evolution during the solidification of a pure material was numerically simulated using the phase field model. The equiaxed dendrite growth in a two-dimensional square domain of undercooled melt (nickel) with four-fold anisotropy was simulated. The phase field model equations was solved using the explicit finite difference method on a uniform mesh. The formation of various equiaxed dendrite patterns was shown by a series of simulations, and the effect of anisotropy on equiaxed dendrite morphology was investigated.
基金China Postdoctoral Science Foundation Under Grant No.20100480321National Basic Research Program of China Under Grant No. 2007CB714200
文摘In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.
