We present a highly efficient lattice Boltzmann model for simulating compressible flows. This model is based on the combination of an appropriate finite difference scheme, a 16-discrete-velocity model [Kataoka and Tsu...We present a highly efficient lattice Boltzmann model for simulating compressible flows. This model is based on the combination of an appropriate finite difference scheme, a 16-discrete-velocity model [Kataoka and Tsutahara, Phys. Rev. E 69 (2004) 035701(R)] and reasonable dispersion and dissipation terms. The dispersion term effectively reduces the oscillation at the discontinuity and enhances numerical precision. The dissipation term makes the new model more easily meet with the yon Neumann stability condition. This model works for both high-speed and low-speed flows with arbitrary specific-heat-ratio. With the new model simulation results for the well-known benchmark problems get a high accuracy compared with the analytic or experimental ones. The used benchmark tests include (i) Shock tubes such as the Sod, Lax, Sjogreen, Colella explosion wave, and collision of two strong shocks, (ii) Regular and Mach shock reflections, and (iii) Shock wave reaction on cylindrical bubble problems. With a more realistic equation of state or free-energy functional, the new model has the potential tostudy the complex procedure of shock wave reaction on porous materials.展开更多
Topography around the Yellow River mouth has changed greatly in recent years, but studies on the current state of ma- rine dynamics off the Yellow River mouth are relatively scarce. This paper uses a two-dimension num...Topography around the Yellow River mouth has changed greatly in recent years, but studies on the current state of ma- rine dynamics off the Yellow River mouth are relatively scarce. This paper uses a two-dimension numerical model (MIKE 21) to reveal the tidal and wave dynamics in 2012, and conducts comparative analysis of the changes from 1996 to 2012. The results show that M2 amphidromic point moved southeastward by 11 kin. It further reveals that the tides around the Yellow River mouth are relatively stable due to the small variations in the tidal constituents. Over the study period, there is no noticeable change in the distribution of tidal types and tidal range, and the mean tidal range off the river mouth during the period studied is 0.5-1.1 m. However, the tidal currents changed greatly due to large change in topography. It is observed that the area with strong tidal currents shifted from the old river mouth (1976-1996) to the modem river mouth (1996-present). While the tidal current speeds decreased continually off the old river mouth, they increased off the modem river mouth. The Maximum Tidal Current Speed (MTCS) reached 1.4 m s-1, and the maximum current speed of 50-year return period reached 2.8 m s-1. Waves also changed greatly due to change in topography. The significant wave height (H1/3) of 50-year return period changed proportionately with the water depth, and the ratio of Hi/3 to depth being 0.4-0.6. H1/3 of the 50-year return period in erosion zone increased continually with increasing water depth, and the rate of change varied between 0.06 and 0.07myr-1. Based on the results of this study, we infer that in the future, the modem river mouth will protrude gradually northward, while the erosion zone, comprising the old river mouth and area between the modern river mouth and the old river mouth (Intermediate region) will continue to erode. As the modem river mouth protrudes towards the sea, there will be a gradual increase in the current speed and decrease in wave height. Conversely, the old river mouth will retreat, with gradual decrease in current speed and increase in wave height. As more coastal constructions spring up around the Yellow River mouth in the future, we recommend that variation in hydrodynamics over time should be taken into consideration when designing such coastal constructions.展开更多
Coal bump is a dynamic process,thus it is necessary to reveal the process validly.2D-ball code is an efficient approach developed by the authors based on the basicDEM rationale proposed by Cundall.Numerical simulation...Coal bump is a dynamic process,thus it is necessary to reveal the process validly.2D-ball code is an efficient approach developed by the authors based on the basicDEM rationale proposed by Cundall.Numerical simulations show that the coal bump experiencesa process of energy accumulation,sudden release of energy and energy decrease.The stiffness of coal particles has a great influence on the coal bump morphosisand particle velocity.Generally,the larger the stiffness of particles,the longer the shootingoff period and the larger the bump velocity.This is in agreement with the results of laboratoryexperiment and in-situ studies.However,the stiffness of particles has an influence onthe quantity value energy and no influence on the releasing energy pattern of the coalbump.展开更多
We represent a two-dimensional model of polydisperse granular mixtures with a power-law size distribution. The model consists of smooth hard disks in a rectangular box with inelastic collisions,driven by a homogeneous...We represent a two-dimensional model of polydisperse granular mixtures with a power-law size distribution. The model consists of smooth hard disks in a rectangular box with inelastic collisions,driven by a homogeneous heat bath at zero gravity.The width of particle size distribution is characterized by the only parameter,namely,the fractal dimension D.The energy dissipation of the mixture is increased as D increases or as e decreases.Furthermore,it is found that the steady-state properties of the mixture such as the collision rate,granular temperature,kinetic pressure and velocity distribution depend sensitively on size distribution parameter D.展开更多
The gradient model of two-dimensional defectless medium is formulated. A graphene sheet is examined as an example of such two-dimensional medium. The problem statement of a graphene sheet deforming in its plane and th...The gradient model of two-dimensional defectless medium is formulated. A graphene sheet is examined as an example of such two-dimensional medium. The problem statement of a graphene sheet deforming in its plane and the bending problem are examined. It is ascertained that the statement of the first problem is equivalent to the flat problem statement of Toupin gradient theory. The statement of the bending problem is equivalent to the plate bending theory of Timoshenko with certain reserves. The characteristic feature of both statements is the fact that the mechanical properties of the sheet of graphene are not defined by “volumetric” moduli but by adhesive ones which have different physical dimension that coincides with the dimension of the corresponding stiffness of classical and nonclassical plates.展开更多
The secondary structure of different Iβ cellulose was analyzed by a molecular dynamics sim- ulation with MARTINI coarse-grained force field, where each chain of the cellulose includes 40 D-glucoses units. Calculation...The secondary structure of different Iβ cellulose was analyzed by a molecular dynamics sim- ulation with MARTINI coarse-grained force field, where each chain of the cellulose includes 40 D-glucoses units. Calculation gives a satisfied description about the secondary structure of the cellulose. As the chain number increasing, the cellulose becomes the form of a helix, with the diameter of screw growing and spiral rising. Interestingly, the celluloses with chain number N of 4, 6, 24 and 36 do show right-hand twisting. On the contrast, the celluloses with N of 8, 12, 16 chains are left-hand twisting. These simulations indicate that the cellulose with chain number larger than 36 will break down to two parts. Besides, the result indicates that 36-chains cellulose model is the most stable among all models. Furthermore, the Lennard-Jones potential determines the secondary structure. In addition, an equation was set up to analyze the twisting structure.展开更多
The miniaturisation context leads to the rise of micro-machining processes. Micro-milling is one of the most flexible and fast of them. Although it is based on the same principles as macro-cutting, it is not a simple ...The miniaturisation context leads to the rise of micro-machining processes. Micro-milling is one of the most flexible and fast of them. Although it is based on the same principles as macro-cutting, it is not a simple scaling-down of it. This down-sizing involves new phenomena in the chip formation, such as the minimum chip thickness below which no chip is formed. This paper presents a review of the current state of the art in this field from an experimental and a numerical point of view. A 2D finite element model is then developed to study the influence of the depth of cut on the chip formation. After the model validation in macro-cutting, it highlights the phenomena reported in literature and allows to perform a minimum chip thickness estimation.展开更多
A numerical solution to a two-dimensional model of flow and transient heat transfer involving solidification in a pipe has been established.Where the temperature of pipe wall is below the freezing point of fluid,phase...A numerical solution to a two-dimensional model of flow and transient heat transfer involving solidification in a pipe has been established.Where the temperature of pipe wall is below the freezing point of fluid,phase change of flowing fluid and the influence of different boundary condition,such as pipe wall temperature,initial temperature and inlet velocity has been taken into account.Also it has been investigated to elicit proper non-dimensional numbers to show the solidification proceeding results.Additionally comparing the two acceptable inlet conditions,show distinctions between velocity inlet and pressure inlet boundary condition in such problems,which affect the whole freezing process.展开更多
In the present study the thermal model of skin and subdermal tissues (SST) of human breast have been developed. The human breast is assumed to be spherical in shape with upper hemisphere projecting out from the trun...In the present study the thermal model of skin and subdermal tissues (SST) of human breast have been developed. The human breast is assumed to be spherical in shape with upper hemisphere projecting out from the trunk of the body and lower hemisphere is considered to be a part of the body core. The upper hemisphere represents the breast and its SST region is divided into three layers namely epidermis, dermis and subdermal tissues. The inner part of the breast represents the core/shell of the breast. The outer surface of the breast is assumed to be exposed to the environment from where the heat loss takes place by conduction, convection, radiation and evaporation. The heat transfer from core to the surface takes place by thermal conduction and blood perfusion. Also metabolic activity takes place at different rates in different SST layers of the breast. Boundary conditions have been framed on the basis of physical conditions. A finite element model has been developed for a two-dimensional steady state case.展开更多
Temperature as an indicator of tissue response is widely used in clinical applications. In view of above a problem of temperature distribution in peripheral regions of extended spherical organs of a human body like, h...Temperature as an indicator of tissue response is widely used in clinical applications. In view of above a problem of temperature distribution in peripheral regions of extended spherical organs of a human body like, human breast involving uniformly perfused tumor is investigated in this paper. The human breast is assumed to be spherical in shape with upper hemisphere projecting out from the trunk of the body and lower hemisphere is considered to be a part of the body core. The outer surface of the breast is assumed to be exposed to the environment from where the heat loss takes place by conduction, convection, radiation and evaporation. The heat transfer from core to the surface takes place by thermal conduction and blood perfusion. Also metabolic activity takes place at different rates in different layers of the breast. An elliptical-shaped tumor is assumed to be present in the dermis region of human breast. A finite element model is developed for a two-dimensional steady state case incorporating the important parameters like blood flow, metabolic activity and thermal conductivity. The triangular ring elements are employed to discretize the region. Appropriate boundary conditions are framed using biophysical conditions. The numerical results are used to study the effect of tumor on temperature distribution in the region.展开更多
The Hooke's atoms with two or more than two electrons give rise to an interesting quantum mechanical model with valuable practical applications. In this work, we study the electronic properties near the origin of the...The Hooke's atoms with two or more than two electrons give rise to an interesting quantum mechanical model with valuable practical applications. In this work, we study the electronic properties near the origin of the harmonic potential. It is seen that the spherically averaged density, ~, exhibits an interesting character -- it has only even order terms in its small r expansion. The spherical average of the Hartree potential, ~, and the spherical average of the Kohn-Sham exchange-correlation potential, 9~c, are also shown to have the same property -- all odd order terms in their expansions vanish. F^rthermore, the analysis and results extend also to the case of two-dimensional models. While only models interacting via. the Coulomb potential are primarily considered in the article, the results also extend to models interacting via. other potentials (viz. Van der Waals potential).展开更多
In this paper, the lowtemperature properties of the spin1 twodimensionM frustrated Heisenberg antifer romagnet with the singleion anisotropy are investigated on a square lattice by using the spinwave theory. The influ...In this paper, the lowtemperature properties of the spin1 twodimensionM frustrated Heisenberg antifer romagnet with the singleion anisotropy are investigated on a square lattice by using the spinwave theory. The influence of the frustration and anisotropy is found in the thermodynamics of the model, such as the temperature dependence of the staggered magnetization and specific heat. For some selected values of the frustration and anisotropy parameters, the results for the specific heat are compared with those of existing theories and numerical estimates. Within a spinwave analysis, we have found the evidence for an intermediate magnetically disorder phase to separate the Nel and collinear phases.展开更多
基金Supported by the Science Foundations of LCP and CAEP under Grant Nos.2009A0102005 and 2009B0101012the National Basic Research Program (973 Program) under Grant No.2007CB815105the National Natural Science Foundation under Grant Nos.10775018,10702010,and 10775088
文摘We present a highly efficient lattice Boltzmann model for simulating compressible flows. This model is based on the combination of an appropriate finite difference scheme, a 16-discrete-velocity model [Kataoka and Tsutahara, Phys. Rev. E 69 (2004) 035701(R)] and reasonable dispersion and dissipation terms. The dispersion term effectively reduces the oscillation at the discontinuity and enhances numerical precision. The dissipation term makes the new model more easily meet with the yon Neumann stability condition. This model works for both high-speed and low-speed flows with arbitrary specific-heat-ratio. With the new model simulation results for the well-known benchmark problems get a high accuracy compared with the analytic or experimental ones. The used benchmark tests include (i) Shock tubes such as the Sod, Lax, Sjogreen, Colella explosion wave, and collision of two strong shocks, (ii) Regular and Mach shock reflections, and (iii) Shock wave reaction on cylindrical bubble problems. With a more realistic equation of state or free-energy functional, the new model has the potential tostudy the complex procedure of shock wave reaction on porous materials.
基金supported by the National Natural Science Foundation of China (Grant Nos. 41030856 and 41006024)the Foundation of Shandong Province (Grant No. BS2012HZ022)+1 种基金the Project of China Geological Survey (Grant No. GZH201100203)the Project of Taishan Scholar
文摘Topography around the Yellow River mouth has changed greatly in recent years, but studies on the current state of ma- rine dynamics off the Yellow River mouth are relatively scarce. This paper uses a two-dimension numerical model (MIKE 21) to reveal the tidal and wave dynamics in 2012, and conducts comparative analysis of the changes from 1996 to 2012. The results show that M2 amphidromic point moved southeastward by 11 kin. It further reveals that the tides around the Yellow River mouth are relatively stable due to the small variations in the tidal constituents. Over the study period, there is no noticeable change in the distribution of tidal types and tidal range, and the mean tidal range off the river mouth during the period studied is 0.5-1.1 m. However, the tidal currents changed greatly due to large change in topography. It is observed that the area with strong tidal currents shifted from the old river mouth (1976-1996) to the modem river mouth (1996-present). While the tidal current speeds decreased continually off the old river mouth, they increased off the modem river mouth. The Maximum Tidal Current Speed (MTCS) reached 1.4 m s-1, and the maximum current speed of 50-year return period reached 2.8 m s-1. Waves also changed greatly due to change in topography. The significant wave height (H1/3) of 50-year return period changed proportionately with the water depth, and the ratio of Hi/3 to depth being 0.4-0.6. H1/3 of the 50-year return period in erosion zone increased continually with increasing water depth, and the rate of change varied between 0.06 and 0.07myr-1. Based on the results of this study, we infer that in the future, the modem river mouth will protrude gradually northward, while the erosion zone, comprising the old river mouth and area between the modern river mouth and the old river mouth (Intermediate region) will continue to erode. As the modem river mouth protrudes towards the sea, there will be a gradual increase in the current speed and decrease in wave height. Conversely, the old river mouth will retreat, with gradual decrease in current speed and increase in wave height. As more coastal constructions spring up around the Yellow River mouth in the future, we recommend that variation in hydrodynamics over time should be taken into consideration when designing such coastal constructions.
基金Supported by the National Natural Science Fundation of China(50534080,50674063)Taishan Scholar Engineering Construction Fund of Shandong Province of China(J06N04)
文摘Coal bump is a dynamic process,thus it is necessary to reveal the process validly.2D-ball code is an efficient approach developed by the authors based on the basicDEM rationale proposed by Cundall.Numerical simulations show that the coal bump experiencesa process of energy accumulation,sudden release of energy and energy decrease.The stiffness of coal particles has a great influence on the coal bump morphosisand particle velocity.Generally,the larger the stiffness of particles,the longer the shootingoff period and the larger the bump velocity.This is in agreement with the results of laboratoryexperiment and in-situ studies.However,the stiffness of particles has an influence onthe quantity value energy and no influence on the releasing energy pattern of the coalbump.
基金Supported by the National Natural Science Foundation of China under Grant Nos.10675048 and 1068006
文摘We represent a two-dimensional model of polydisperse granular mixtures with a power-law size distribution. The model consists of smooth hard disks in a rectangular box with inelastic collisions,driven by a homogeneous heat bath at zero gravity.The width of particle size distribution is characterized by the only parameter,namely,the fractal dimension D.The energy dissipation of the mixture is increased as D increases or as e decreases.Furthermore,it is found that the steady-state properties of the mixture such as the collision rate,granular temperature,kinetic pressure and velocity distribution depend sensitively on size distribution parameter D.
文摘The gradient model of two-dimensional defectless medium is formulated. A graphene sheet is examined as an example of such two-dimensional medium. The problem statement of a graphene sheet deforming in its plane and the bending problem are examined. It is ascertained that the statement of the first problem is equivalent to the flat problem statement of Toupin gradient theory. The statement of the bending problem is equivalent to the plate bending theory of Timoshenko with certain reserves. The characteristic feature of both statements is the fact that the mechanical properties of the sheet of graphene are not defined by “volumetric” moduli but by adhesive ones which have different physical dimension that coincides with the dimension of the corresponding stiffness of classical and nonclassical plates.
文摘The secondary structure of different Iβ cellulose was analyzed by a molecular dynamics sim- ulation with MARTINI coarse-grained force field, where each chain of the cellulose includes 40 D-glucoses units. Calculation gives a satisfied description about the secondary structure of the cellulose. As the chain number increasing, the cellulose becomes the form of a helix, with the diameter of screw growing and spiral rising. Interestingly, the celluloses with chain number N of 4, 6, 24 and 36 do show right-hand twisting. On the contrast, the celluloses with N of 8, 12, 16 chains are left-hand twisting. These simulations indicate that the cellulose with chain number larger than 36 will break down to two parts. Besides, the result indicates that 36-chains cellulose model is the most stable among all models. Furthermore, the Lennard-Jones potential determines the secondary structure. In addition, an equation was set up to analyze the twisting structure.
文摘The miniaturisation context leads to the rise of micro-machining processes. Micro-milling is one of the most flexible and fast of them. Although it is based on the same principles as macro-cutting, it is not a simple scaling-down of it. This down-sizing involves new phenomena in the chip formation, such as the minimum chip thickness below which no chip is formed. This paper presents a review of the current state of the art in this field from an experimental and a numerical point of view. A 2D finite element model is then developed to study the influence of the depth of cut on the chip formation. After the model validation in macro-cutting, it highlights the phenomena reported in literature and allows to perform a minimum chip thickness estimation.
文摘A numerical solution to a two-dimensional model of flow and transient heat transfer involving solidification in a pipe has been established.Where the temperature of pipe wall is below the freezing point of fluid,phase change of flowing fluid and the influence of different boundary condition,such as pipe wall temperature,initial temperature and inlet velocity has been taken into account.Also it has been investigated to elicit proper non-dimensional numbers to show the solidification proceeding results.Additionally comparing the two acceptable inlet conditions,show distinctions between velocity inlet and pressure inlet boundary condition in such problems,which affect the whole freezing process.
文摘In the present study the thermal model of skin and subdermal tissues (SST) of human breast have been developed. The human breast is assumed to be spherical in shape with upper hemisphere projecting out from the trunk of the body and lower hemisphere is considered to be a part of the body core. The upper hemisphere represents the breast and its SST region is divided into three layers namely epidermis, dermis and subdermal tissues. The inner part of the breast represents the core/shell of the breast. The outer surface of the breast is assumed to be exposed to the environment from where the heat loss takes place by conduction, convection, radiation and evaporation. The heat transfer from core to the surface takes place by thermal conduction and blood perfusion. Also metabolic activity takes place at different rates in different SST layers of the breast. Boundary conditions have been framed on the basis of physical conditions. A finite element model has been developed for a two-dimensional steady state case.
文摘Temperature as an indicator of tissue response is widely used in clinical applications. In view of above a problem of temperature distribution in peripheral regions of extended spherical organs of a human body like, human breast involving uniformly perfused tumor is investigated in this paper. The human breast is assumed to be spherical in shape with upper hemisphere projecting out from the trunk of the body and lower hemisphere is considered to be a part of the body core. The outer surface of the breast is assumed to be exposed to the environment from where the heat loss takes place by conduction, convection, radiation and evaporation. The heat transfer from core to the surface takes place by thermal conduction and blood perfusion. Also metabolic activity takes place at different rates in different layers of the breast. An elliptical-shaped tumor is assumed to be present in the dermis region of human breast. A finite element model is developed for a two-dimensional steady state case incorporating the important parameters like blood flow, metabolic activity and thermal conductivity. The triangular ring elements are employed to discretize the region. Appropriate boundary conditions are framed using biophysical conditions. The numerical results are used to study the effect of tumor on temperature distribution in the region.
文摘The Hooke's atoms with two or more than two electrons give rise to an interesting quantum mechanical model with valuable practical applications. In this work, we study the electronic properties near the origin of the harmonic potential. It is seen that the spherically averaged density, ~, exhibits an interesting character -- it has only even order terms in its small r expansion. The spherical average of the Hartree potential, ~, and the spherical average of the Kohn-Sham exchange-correlation potential, 9~c, are also shown to have the same property -- all odd order terms in their expansions vanish. F^rthermore, the analysis and results extend also to the case of two-dimensional models. While only models interacting via. the Coulomb potential are primarily considered in the article, the results also extend to models interacting via. other potentials (viz. Van der Waals potential).
文摘In this paper, the lowtemperature properties of the spin1 twodimensionM frustrated Heisenberg antifer romagnet with the singleion anisotropy are investigated on a square lattice by using the spinwave theory. The influence of the frustration and anisotropy is found in the thermodynamics of the model, such as the temperature dependence of the staggered magnetization and specific heat. For some selected values of the frustration and anisotropy parameters, the results for the specific heat are compared with those of existing theories and numerical estimates. Within a spinwave analysis, we have found the evidence for an intermediate magnetically disorder phase to separate the Nel and collinear phases.