Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise re...Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.展开更多
The difference of sintering crunodes of metal powders and fibers is discussed. The mathematical model of the surface diffusion described by the difference in mean curvature is defined as a Hamilton-Jacobi-type equatio...The difference of sintering crunodes of metal powders and fibers is discussed. The mathematical model of the surface diffusion described by the difference in mean curvature is defined as a Hamilton-Jacobi-type equation, and the model is numerically solved by the level set method. The three-dimensional numerical simulations of two metal powders and fibers(the fiber angle is 0° or 90°) are implemented by this mathematical model, respectively. The numerical simulation results accord with the experimental ones. The sintering neck growth trends of metal powders and metal fibers are similar. The sintering neck radius of metal fibers is larger than that of metal powders. The difference of the neck radius is caused by the difference of geometric structure which makes an important influence on the curvature affecting the migration rate of atoms.展开更多
Based on the analysis of the principles of frequency-phase scanning 3-D （three-dimensional） radar and the scattering mechanism of 3-D radar, the target and clutter IF （intermediate frequency） signals model of freq...Based on the analysis of the principles of frequency-phase scanning 3-D （three-dimensional） radar and the scattering mechanism of 3-D radar, the target and clutter IF （intermediate frequency） signals model of frequencyphase scanning 3-D radar is presented. The IF signals model of different channels of 3-D radar is presented in one simple formula in which complex waveform effects are considered. The simulation results obtained during tests are also provided in the end.展开更多
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.展开更多
Two-dimensional images of the granular ore media with different grain sizes were obtained from the X-ray computed tomography.Combined with the digital image processing and finite element techniques,the original graysc...Two-dimensional images of the granular ore media with different grain sizes were obtained from the X-ray computed tomography.Combined with the digital image processing and finite element techniques,the original grayscale images were transformed into the finite element models directly.By using these models,the simulations of pore scale fluid flow among particles were conducted with the COMSOL Multiphysics,and the distribution characteristics of fluid flow velocity and pressure were analyzed.The simulation results show that there exist obvious preferential flow and leaching blind zone in each granular medium.The flow velocity at pore throat is larger than that of pore body and the largest velocity reaches 0.22 m/s.The velocity decreases gradually from the center of pore throat and body to the surface of particles.The flow paths of granular media with larger grain size distribute equally,while the fluid flow velocities in most of areas of granular media with smaller grain size are lower,and some of them approach to zero,so the permeability is very low.There exist some pore clusters with different pressures,which is the basic reason for the uneven flow velocity distribution.展开更多
According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in...According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.展开更多
The magnetohydrodynamic(MHD) three-dimensional flow of Jeffrey fluid in the presence of Newtonian heating is investigated. Flow is caused by a bidirectional stretching surface. Series solutions are constructed for the...The magnetohydrodynamic(MHD) three-dimensional flow of Jeffrey fluid in the presence of Newtonian heating is investigated. Flow is caused by a bidirectional stretching surface. Series solutions are constructed for the velocity and temperature fields. Convergence of series solutions is ensured graphically and numerically. The variations of key parameters on the physical quantities are shown and discussed in detail. Constructed series solutions are compared with the existing solutions in the limiting case and an excellent agreement is noticed. Nusselt numbers are computed with and without magnetic fields. It is observed that the Nusselt number decreases in the presence of magnetic field.展开更多
Heat and mass transfer effects in three-dimensional flow of Maxwell fluid over a stretching surface were addressed.Analysis was performed in the presence of internal heat generation/absorption. Concentration and therm...Heat and mass transfer effects in three-dimensional flow of Maxwell fluid over a stretching surface were addressed.Analysis was performed in the presence of internal heat generation/absorption. Concentration and thermal buoyancy effects were accounted. Convective boundary conditions for heat and mass transfer analysis were explored. Series solutions of the resulting problem were developed. Effects of mixed convection, internal heat generation/absorption parameter and Biot numbers on the dimensionless velocity, temperature and concentration distributions were illustrated graphically. Numerical values of local Nusselt and Sherwood numbers were obtained and analyzed for all the physical parameters. It is found that both thermal and concentration boundary layer thicknesses are decreasing functions of stretching ratio. Variations of mixed convection parameter and concentration buoyancy parameter on the velocity profiles and associated boundary layer thicknesses are enhanced. Velocity profiles and temperature increase in the case of internal heat generation while they reduce for heat absorption. Heat transfer Biot number increases the thermal boundary layer thickness and temperature. Also concentration and its associated boundary layer are enhanced with an increase in mass transfer Biot number. The local Nusselt and Sherwood numbers have quite similar behaviors for increasing values of mixed convection parameter, concentration buoyancy parameter and Deborah number.展开更多
Fatigue and tensile behaviors of homogenized WE 54 magnesium alloy before and after immersion in simulated body fluid(SBF)were investigated.According to the tensile test,the alloy without immersion in SBF solution has...Fatigue and tensile behaviors of homogenized WE 54 magnesium alloy before and after immersion in simulated body fluid(SBF)were investigated.According to the tensile test,the alloy without immersion in SBF solution has the highest tensile strength of 278 MPa,which decreased to 190 MPa after 336 h of immersion..The fatigue life of the homogenized WE 54 magnesium alloy before immersion in the SBF solution under a constant stress of 15 MPa is 3598 cycles.However,the fatigue life of the alloy decreased to 453 cycles after 336 h of immersion in the SBF solution under the same stress.Examination of the fracture surface of the samples by SEM reveals that the origin of the fatigue crack before immersion is micro-pores and defects.While corrosion pits and cracks are the main reasons for forming the initial fatigue crack after immersion.Moreover,the results obtained from practical work were evaluated and compared to theoretical calculations.The area of the hysteresis loops of the samples after the fatigue test,determined using Triangles and Monte Carlo methods,decreased from 4954.5 MPa and 4842.9 MPa before immersion to 192.0 MPa and 175.8 MPa after 336 h of immersion,respectively.展开更多
Probabilistic analysis is a rational approach for engineering design because it provides more insight than traditional deterministic analysis. Probabilistic evaluation on seismic stability of three dimensional (3D) sl...Probabilistic analysis is a rational approach for engineering design because it provides more insight than traditional deterministic analysis. Probabilistic evaluation on seismic stability of three dimensional (3D) slopes is studied in this paper. The slope safety factor is computed by combining the kinematic approach of limit analysis using a three-dimensional rotational failure mechanism with the pseudo-dynamic approach. The variability of input parameters, including six pseudo-dynamic parameters and two soil shear strength parameters, are taken into account by means of Monte-Carlo Simulations (MCS) method. The influences of pseudo-dynamic input variables on the computed failure probabilities are investigated and discussed. It is shown that the obtained failure probabilities increase with the pseudo-dynamic input variables and the pseudo-dynamic approach gives more conservative failure probability estimates compared with the pseudo-static approach.展开更多
A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-indu...A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.展开更多
Proper room and pillar sizes are both critical factors for safe mining and high ore recovery rate in shrinkage stoping mining of underground metal mines.The rock masses of Tangdan copper mine of China are fractured,wh...Proper room and pillar sizes are both critical factors for safe mining and high ore recovery rate in shrinkage stoping mining of underground metal mines.The rock masses of Tangdan copper mine of China are fractured,which needs much reinforcement and support prior to mining.Cement-sodium silicate grout technology was selected,then its related parameters such as grout pressure,diffusion radius and time were calculated and proposed.In order to test the effect of the pressured grout in the fractured No.4 ore block,field experiments were conducted.To optimize stoping configuration,three-dimensional numerical simulation with ANSYS and FLAC 3 D softwares was proposed.The results show that the drilling porosity and mechanical properties of the rock masses are increased obviously.After grout,ore recovery rate is increased by 10.2%employing the newly designed stoping configuration compared with the previous.Last,analyzed from the surface movements,roof subsidence and the maximum principal stress of the pillars,the mining safety is probable of being ensured.展开更多
Thermal comfort and indoor air quality as well as the energy efficiency have been recognized as essential parts of sustainable building assessment. This work aims to analyze the energy conservation of the heat recover...Thermal comfort and indoor air quality as well as the energy efficiency have been recognized as essential parts of sustainable building assessment. This work aims to analyze the energy conservation of the heat recovery ventilator and to investigate the effect of the air supply arrangement. Three types of mixing ventilation are chosen for the analysis of coupling ANSYS/FLUENT （a computational fluid dynamics （CFD） program） with TRNSYS （a building energy simulation （BES） software）. The adoption of mutual complementary boundary conditions for CFD and BES provides more accurate and complete information of indoor air distribution and thermal performance in buildings. A typical office-space situated in a middle storey is chosen for the analysis. The office-space is equipped with air-conditioners on the ceiling. A heat recovery ventilation system directly supplies flesh air to the office space. Its thermal performance and indoor air distribution predicted by the coupled method are compared under three types of ventilation system. When the supply and return openings for ventilation are arranged on the ceiling, there is no critical difference between the predictions of the coupled method and BES on the energy consumption of HVAC because PID control is adopted for the supply air temperature of the occupied zone. On the other hand, approximately 21% discrepancy for the heat recovery estimation in the maximum between the simulated results of coupled method and BES-only can be obviously found in the floor air supply ventilation case. The discrepancy emphasizes the necessity of coupling CFD with BES when vertical air temperature gradient exists. Our future target is to estimate the optimum design of heat recovery ventilation system to control CO2 concentration by adjusting flow rate of flesh air.展开更多
A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂F...A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂Fish＂ language provided by PFC3D and was employed to simulate the three-point bending beam test at two temperature levels： -10 ℃ and 15℃. The AC beam was modeled with the consideration of the microstructural features of asphalt mixtures. Uniaxial complex modulus test and indirect tensile strength test were conducted to obtain material input parameters for numerical modeling. The 3D predictions were validated using laboratory experimental measurements of AC beams prepared by the same mixture design. Effects of mastic stiffness, cohesive and adhesive strength on AC fracture behavior were investigated using the DEM model. The results show that the 3D DEM fracture model can accurately predict the fracture patterns of asphalt concrete. The ratio of stress at interfaces to the stress in mastics increases as the mastic stiffness decreases; however, the increase in the cohesive strength or adhesive strength shows no significant influence on the tensile strength.展开更多
A multi-domain nonlinear dynamic model of a proportional solenoid valve was presented.The electro-magnetic,mechanical and fluid subsystems of the valve were investigated,including their interactions.Governing equation...A multi-domain nonlinear dynamic model of a proportional solenoid valve was presented.The electro-magnetic,mechanical and fluid subsystems of the valve were investigated,including their interactions.Governing equations of the valve were derived in the form of nonlinear state equations.By comparing the simulated and measured data,the simulation model is validated with a deviation less than 15%,which can be used for the structural design and control algorithm optimization of proportional solenoid valves.展开更多
The frequent occurrence of safety accidents during the calendering process is caused by the flammable and explosive properties of composite modified double-base(CMDB)propellant.Optimization of process parameters with ...The frequent occurrence of safety accidents during the calendering process is caused by the flammable and explosive properties of composite modified double-base(CMDB)propellant.Optimization of process parameters with the aid of fluid simulation technology could effectively ensure the safety of the calendering process.To improve the accuracy of the simulation results,material parameters and model structure were corrected based on actual conditions,and adaptive grid technology was applied in the local mesh refinement.In addition,the rheological behavior,motion trajectories and heat transfer mechanisms of CMDB propellant slurry were studied with different gaps,rotational rates and temperatures of two rollers.The results indicated that the refined mesh could significantly improve the contour clarity of boundaries and simulate the characteristics of CMDB propellant slurry reflux movement caused by the convergent flow near the outlet.Compared with the gap,the increased rotational rate of roller could promote the reflux movement and intensify the shear flow of slurry inside the flow region by viscous shear dragging.Meanwhile,under the synergistic effect of contact heat transfer as well as convective heat exchange,heat accumulated near the outlet and diffused along the reflux movement,which led to the countercurrent heat dissipation behavior of CMDB propellant slurry.The plasticizing mechanism of slurry and the safety of calendering under different conditions were explored,which provided theoretical guidance and reference data for the optimization of calendering process conditions.Based on the simulation results,the safety of the CMDB propellant calendering process could be significantly improved with a few tests conducted during a short research and development cycle.展开更多
The finite element method was used to solve fluid dynamic interaction problems between the crust and mantle of the Earth. To consider different mechanical behaviours, the lithosphere consisting of the crust and upper ...The finite element method was used to solve fluid dynamic interaction problems between the crust and mantle of the Earth. To consider different mechanical behaviours, the lithosphere consisting of the crust and upper mantle was simulated as fluid-saturated porous rocks, while the upper aesthenospheric part of the mantle was simulated as viscous fluids. Since the whole lithosphere was computationally simulated, the dynamic interaction between the crust and the upper mantle was appropriately considered. In particular, the mixing of mantle fluids and crustal fluids was simulated in the corresponding computational model. The related computational simulation results from an example problem demonstrate that the mantle fluids can flow into the crust and mix with the crustal fluids due to the resulting convective flows in the crust-mantle system. Likewise, the crustal fluids can also flow into the upper mantle and mix with the mantle fluids. This kind of fluids mixing and exchange is very important to the better understanding of the governing processes that control the ore body formation and mineralization in the upper crust of the Earth.展开更多
Porous and dense TiNi alloys were successfully fabricated by powder metallurgy(P/M) method, and to further improve their surface biocompatibility, surface modification techniques including grind using silicon-carbide(...Porous and dense TiNi alloys were successfully fabricated by powder metallurgy(P/M) method, and to further improve their surface biocompatibility, surface modification techniques including grind using silicon-carbide(SiC) paper, acid etching and alkali treatment were employed to produce either irregularly rough surface or micro-porous surface roughness. X-ray diffractometry(XRD), scanning electron microscopy(SEM) and energy dispersive X-ray spectroscopy(EDX) attached to SEM were used to characterize surface structure and the Ca-P coatings. Effects of the above surface treatments on the surface morphology, apatite forming ability were systematically investigated. Results indicate that all the above surface treatments increase the apatite forming ability of TiNi alloys in varying degrees when soaked in simulated body fluid(SBF). More apatite coatings formed on TiNi samples sintered at 1050℃ and 1100℃ due to their high porosity and pure TiNi phase that is beneficial to heterogeneous nucleation. Furthermore, more uniform apatite was fabricated on the sample sintered from the mixture of Ni and Ti powders.展开更多
In order to improve the benefits of base bleed in base flow field,the base flow with hot base bleed for two jet models is studied.Twodimensional axisymmetric Navier-Stokes equations are computed by using a finite volu...In order to improve the benefits of base bleed in base flow field,the base flow with hot base bleed for two jet models is studied.Twodimensional axisymmetric Navier-Stokes equations are computed by using a finite volume scheme.The base flow of a cylinder afterbody with base bleed is simulated.The simulation results are validated with the experimental data,and the experimental results are well reproduced.On this basis,the base flow fields with base bleed for a circular jet model and an annulus jet model are investigated by selecting the injection temperature from 830 K to 2200 K.The results show that the base pressure of the annular jet model is higher than that of the circular jet model with the changes of the injection parameter and the injection temperature.For the circular jet model,the hot gases are concentrated in the vicinity of the base.For the annular jet model,the bleed gases flow into the shear layer directly so that the hot gases are concentrated in the shear layer.The latter temperature distribution is better for the increase of base pressure.展开更多
The Ailaoshan aquamarine-bearing pegmatites are associated with Proterozoic metamorphic rocks in the southern portion of the Ailaoshan fault-folded complex.The gem-bearing pegmatite mineralization zones of the region ...The Ailaoshan aquamarine-bearing pegmatites are associated with Proterozoic metamorphic rocks in the southern portion of the Ailaoshan fault-folded complex.The gem-bearing pegmatite mineralization zones of the region occur in areas generally consistent with the regional tectonic trend.The pegmatites are found in metamorphic rocks,migmatites and in the inner/outer contact zones of gneissoid granites. The Rb-Sr isochron drawn for the pegmatites is 26～31 Ma,(i.e.in Himalayan).The homogenization temperatures of melt and liquid inclusions in minerals vary from 185 to 920℃,which are comparable to the inclusions observed in banded migmatites and ptygmatic quartz veins in the surrounding metamorphic rocks. The mineralization fluids of the pegmatite were rich in HCO_3 and CO_2,and their compositional assemblages are comparable to metamorphic fluids.Results of H,O,C,Si etc.isotopic analyses and REE,and Be analyses indicates that the sources of mineralization components that formed the pegmatites are closely associated with metamorphic fluids and the enclosing metamorphic rocks. A pegmatite structure simulation experiment was conducted at high temperature and pressure(840℃and 1,500×105Pa.),with various metamorphic rock samples in a water-rich and volatile-rich environment.When the liquidus was reached,the temperature was gradually decreased at the rate of 5～10℃/day over a time period of three months.SEM energy-dispersive spectrum analyses were performed on the experimental products.A series of pegmatoid textures were observed including zonal texture,megacryst texture,drusy cavities,crystal druses,and vesicular texture along with more than ten types of minerals including plagioclase,microcline,quartz and biotite.Different metamorphic rock melts generated different mineral assemblages.Experiment results revealed that the partial melting of metamorphic rocks could form melts similar to pegmatite magmas. Based upon the geological characteristics,geochemistry,and pegmatite texture simulation experimental results,it is concluded that the mineralization components of Ailaoshan aquamarine-bearing pegmatites came from metamorphic rocks.The petrogenetic model for the origin of pegmatites is related to ultrametamorphism and metamorphic anatexis.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.12302435 and 12221002)。
文摘Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.
基金Projects(51174236,51134003)supported by the National Natural Science Foundation of ChinaProject(2011CB606306)supported by the National Basic Research Program of ChinaProject(PMM-SKL-4-2012)supported by the Opening Project of State Key Laboratory of Porous Metal Materials(Northwest Institute for Nonferrous Metal Research),China
文摘The difference of sintering crunodes of metal powders and fibers is discussed. The mathematical model of the surface diffusion described by the difference in mean curvature is defined as a Hamilton-Jacobi-type equation, and the model is numerically solved by the level set method. The three-dimensional numerical simulations of two metal powders and fibers(the fiber angle is 0° or 90°) are implemented by this mathematical model, respectively. The numerical simulation results accord with the experimental ones. The sintering neck growth trends of metal powders and metal fibers are similar. The sintering neck radius of metal fibers is larger than that of metal powders. The difference of the neck radius is caused by the difference of geometric structure which makes an important influence on the curvature affecting the migration rate of atoms.
文摘Based on the analysis of the principles of frequency-phase scanning 3-D （three-dimensional） radar and the scattering mechanism of 3-D radar, the target and clutter IF （intermediate frequency） signals model of frequencyphase scanning 3-D radar is presented. The IF signals model of different channels of 3-D radar is presented in one simple formula in which complex waveform effects are considered. The simulation results obtained during tests are also provided in the end.
基金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.
基金Projects(50934002,51074013) supported by the National Natural Science Foundation of China
文摘Two-dimensional images of the granular ore media with different grain sizes were obtained from the X-ray computed tomography.Combined with the digital image processing and finite element techniques,the original grayscale images were transformed into the finite element models directly.By using these models,the simulations of pore scale fluid flow among particles were conducted with the COMSOL Multiphysics,and the distribution characteristics of fluid flow velocity and pressure were analyzed.The simulation results show that there exist obvious preferential flow and leaching blind zone in each granular medium.The flow velocity at pore throat is larger than that of pore body and the largest velocity reaches 0.22 m/s.The velocity decreases gradually from the center of pore throat and body to the surface of particles.The flow paths of granular media with larger grain size distribute equally,while the fluid flow velocities in most of areas of granular media with smaller grain size are lower,and some of them approach to zero,so the permeability is very low.There exist some pore clusters with different pressures,which is the basic reason for the uneven flow velocity distribution.
基金Project(51074027)supported by the National Natural Science Foundation of China
文摘According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.
文摘The magnetohydrodynamic(MHD) three-dimensional flow of Jeffrey fluid in the presence of Newtonian heating is investigated. Flow is caused by a bidirectional stretching surface. Series solutions are constructed for the velocity and temperature fields. Convergence of series solutions is ensured graphically and numerically. The variations of key parameters on the physical quantities are shown and discussed in detail. Constructed series solutions are compared with the existing solutions in the limiting case and an excellent agreement is noticed. Nusselt numbers are computed with and without magnetic fields. It is observed that the Nusselt number decreases in the presence of magnetic field.
文摘Heat and mass transfer effects in three-dimensional flow of Maxwell fluid over a stretching surface were addressed.Analysis was performed in the presence of internal heat generation/absorption. Concentration and thermal buoyancy effects were accounted. Convective boundary conditions for heat and mass transfer analysis were explored. Series solutions of the resulting problem were developed. Effects of mixed convection, internal heat generation/absorption parameter and Biot numbers on the dimensionless velocity, temperature and concentration distributions were illustrated graphically. Numerical values of local Nusselt and Sherwood numbers were obtained and analyzed for all the physical parameters. It is found that both thermal and concentration boundary layer thicknesses are decreasing functions of stretching ratio. Variations of mixed convection parameter and concentration buoyancy parameter on the velocity profiles and associated boundary layer thicknesses are enhanced. Velocity profiles and temperature increase in the case of internal heat generation while they reduce for heat absorption. Heat transfer Biot number increases the thermal boundary layer thickness and temperature. Also concentration and its associated boundary layer are enhanced with an increase in mass transfer Biot number. The local Nusselt and Sherwood numbers have quite similar behaviors for increasing values of mixed convection parameter, concentration buoyancy parameter and Deborah number.
文摘Fatigue and tensile behaviors of homogenized WE 54 magnesium alloy before and after immersion in simulated body fluid(SBF)were investigated.According to the tensile test,the alloy without immersion in SBF solution has the highest tensile strength of 278 MPa,which decreased to 190 MPa after 336 h of immersion..The fatigue life of the homogenized WE 54 magnesium alloy before immersion in the SBF solution under a constant stress of 15 MPa is 3598 cycles.However,the fatigue life of the alloy decreased to 453 cycles after 336 h of immersion in the SBF solution under the same stress.Examination of the fracture surface of the samples by SEM reveals that the origin of the fatigue crack before immersion is micro-pores and defects.While corrosion pits and cracks are the main reasons for forming the initial fatigue crack after immersion.Moreover,the results obtained from practical work were evaluated and compared to theoretical calculations.The area of the hysteresis loops of the samples after the fatigue test,determined using Triangles and Monte Carlo methods,decreased from 4954.5 MPa and 4842.9 MPa before immersion to 192.0 MPa and 175.8 MPa after 336 h of immersion,respectively.
文摘Probabilistic analysis is a rational approach for engineering design because it provides more insight than traditional deterministic analysis. Probabilistic evaluation on seismic stability of three dimensional (3D) slopes is studied in this paper. The slope safety factor is computed by combining the kinematic approach of limit analysis using a three-dimensional rotational failure mechanism with the pseudo-dynamic approach. The variability of input parameters, including six pseudo-dynamic parameters and two soil shear strength parameters, are taken into account by means of Monte-Carlo Simulations (MCS) method. The influences of pseudo-dynamic input variables on the computed failure probabilities are investigated and discussed. It is shown that the obtained failure probabilities increase with the pseudo-dynamic input variables and the pseudo-dynamic approach gives more conservative failure probability estimates compared with the pseudo-static approach.
文摘A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.
基金Projects(51374034,51674012)supported by the National Natural Science Foundation of ChinaProject(2013BAB02B05)supported by the China National Science and Technology Support Program during the 12th Five-Year Plan Period
文摘Proper room and pillar sizes are both critical factors for safe mining and high ore recovery rate in shrinkage stoping mining of underground metal mines.The rock masses of Tangdan copper mine of China are fractured,which needs much reinforcement and support prior to mining.Cement-sodium silicate grout technology was selected,then its related parameters such as grout pressure,diffusion radius and time were calculated and proposed.In order to test the effect of the pressured grout in the fractured No.4 ore block,field experiments were conducted.To optimize stoping configuration,three-dimensional numerical simulation with ANSYS and FLAC 3 D softwares was proposed.The results show that the drilling porosity and mechanical properties of the rock masses are increased obviously.After grout,ore recovery rate is increased by 10.2%employing the newly designed stoping configuration compared with the previous.Last,analyzed from the surface movements,roof subsidence and the maximum principal stress of the pillars,the mining safety is probable of being ensured.
基金Project supported by Grant-in-Aid for Scientific Research (JSPS KAKENHI for Young Scientists (S), 21676005)
文摘Thermal comfort and indoor air quality as well as the energy efficiency have been recognized as essential parts of sustainable building assessment. This work aims to analyze the energy conservation of the heat recovery ventilator and to investigate the effect of the air supply arrangement. Three types of mixing ventilation are chosen for the analysis of coupling ANSYS/FLUENT （a computational fluid dynamics （CFD） program） with TRNSYS （a building energy simulation （BES） software）. The adoption of mutual complementary boundary conditions for CFD and BES provides more accurate and complete information of indoor air distribution and thermal performance in buildings. A typical office-space situated in a middle storey is chosen for the analysis. The office-space is equipped with air-conditioners on the ceiling. A heat recovery ventilation system directly supplies flesh air to the office space. Its thermal performance and indoor air distribution predicted by the coupled method are compared under three types of ventilation system. When the supply and return openings for ventilation are arranged on the ceiling, there is no critical difference between the predictions of the coupled method and BES on the energy consumption of HVAC because PID control is adopted for the supply air temperature of the occupied zone. On the other hand, approximately 21% discrepancy for the heat recovery estimation in the maximum between the simulated results of coupled method and BES-only can be obviously found in the floor air supply ventilation case. The discrepancy emphasizes the necessity of coupling CFD with BES when vertical air temperature gradient exists. Our future target is to estimate the optimum design of heat recovery ventilation system to control CO2 concentration by adjusting flow rate of flesh air.
基金Project(51208178)supported by the National Natural Science Foundation of ChinaProject(2012M520991)supported by China Postdoctoral Science Foundation
文摘A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂Fish＂ language provided by PFC3D and was employed to simulate the three-point bending beam test at two temperature levels： -10 ℃ and 15℃. The AC beam was modeled with the consideration of the microstructural features of asphalt mixtures. Uniaxial complex modulus test and indirect tensile strength test were conducted to obtain material input parameters for numerical modeling. The 3D predictions were validated using laboratory experimental measurements of AC beams prepared by the same mixture design. Effects of mastic stiffness, cohesive and adhesive strength on AC fracture behavior were investigated using the DEM model. The results show that the 3D DEM fracture model can accurately predict the fracture patterns of asphalt concrete. The ratio of stress at interfaces to the stress in mastics increases as the mastic stiffness decreases; however, the increase in the cohesive strength or adhesive strength shows no significant influence on the tensile strength.
基金Project(2008ZHZX1A0502) supported by the Independence Innovation Achievements Transformation Crucial Special Program of Shandong Province,China
文摘A multi-domain nonlinear dynamic model of a proportional solenoid valve was presented.The electro-magnetic,mechanical and fluid subsystems of the valve were investigated,including their interactions.Governing equations of the valve were derived in the form of nonlinear state equations.By comparing the simulated and measured data,the simulation model is validated with a deviation less than 15%,which can be used for the structural design and control algorithm optimization of proportional solenoid valves.
文摘The frequent occurrence of safety accidents during the calendering process is caused by the flammable and explosive properties of composite modified double-base(CMDB)propellant.Optimization of process parameters with the aid of fluid simulation technology could effectively ensure the safety of the calendering process.To improve the accuracy of the simulation results,material parameters and model structure were corrected based on actual conditions,and adaptive grid technology was applied in the local mesh refinement.In addition,the rheological behavior,motion trajectories and heat transfer mechanisms of CMDB propellant slurry were studied with different gaps,rotational rates and temperatures of two rollers.The results indicated that the refined mesh could significantly improve the contour clarity of boundaries and simulate the characteristics of CMDB propellant slurry reflux movement caused by the convergent flow near the outlet.Compared with the gap,the increased rotational rate of roller could promote the reflux movement and intensify the shear flow of slurry inside the flow region by viscous shear dragging.Meanwhile,under the synergistic effect of contact heat transfer as well as convective heat exchange,heat accumulated near the outlet and diffused along the reflux movement,which led to the countercurrent heat dissipation behavior of CMDB propellant slurry.The plasticizing mechanism of slurry and the safety of calendering under different conditions were explored,which provided theoretical guidance and reference data for the optimization of calendering process conditions.Based on the simulation results,the safety of the CMDB propellant calendering process could be significantly improved with a few tests conducted during a short research and development cycle.
基金Project(10872219) supported by the National Natural Science Foundation of China
文摘The finite element method was used to solve fluid dynamic interaction problems between the crust and mantle of the Earth. To consider different mechanical behaviours, the lithosphere consisting of the crust and upper mantle was simulated as fluid-saturated porous rocks, while the upper aesthenospheric part of the mantle was simulated as viscous fluids. Since the whole lithosphere was computationally simulated, the dynamic interaction between the crust and the upper mantle was appropriately considered. In particular, the mixing of mantle fluids and crustal fluids was simulated in the corresponding computational model. The related computational simulation results from an example problem demonstrate that the mantle fluids can flow into the crust and mix with the crustal fluids due to the resulting convective flows in the crust-mantle system. Likewise, the crustal fluids can also flow into the upper mantle and mix with the mantle fluids. This kind of fluids mixing and exchange is very important to the better understanding of the governing processes that control the ore body formation and mineralization in the upper crust of the Earth.
基金Project(51274247) supported by the National Natural Science Foundation of ChinaProject(2014zzts177) support by the Fundamental Research Funds for the Central Universities,China
文摘Porous and dense TiNi alloys were successfully fabricated by powder metallurgy(P/M) method, and to further improve their surface biocompatibility, surface modification techniques including grind using silicon-carbide(SiC) paper, acid etching and alkali treatment were employed to produce either irregularly rough surface or micro-porous surface roughness. X-ray diffractometry(XRD), scanning electron microscopy(SEM) and energy dispersive X-ray spectroscopy(EDX) attached to SEM were used to characterize surface structure and the Ca-P coatings. Effects of the above surface treatments on the surface morphology, apatite forming ability were systematically investigated. Results indicate that all the above surface treatments increase the apatite forming ability of TiNi alloys in varying degrees when soaked in simulated body fluid(SBF). More apatite coatings formed on TiNi samples sintered at 1050℃ and 1100℃ due to their high porosity and pure TiNi phase that is beneficial to heterogeneous nucleation. Furthermore, more uniform apatite was fabricated on the sample sintered from the mixture of Ni and Ti powders.
基金supported by National Nature Science Foundation of China(Grant No.51176076)
文摘In order to improve the benefits of base bleed in base flow field,the base flow with hot base bleed for two jet models is studied.Twodimensional axisymmetric Navier-Stokes equations are computed by using a finite volume scheme.The base flow of a cylinder afterbody with base bleed is simulated.The simulation results are validated with the experimental data,and the experimental results are well reproduced.On this basis,the base flow fields with base bleed for a circular jet model and an annulus jet model are investigated by selecting the injection temperature from 830 K to 2200 K.The results show that the base pressure of the annular jet model is higher than that of the circular jet model with the changes of the injection parameter and the injection temperature.For the circular jet model,the hot gases are concentrated in the vicinity of the base.For the annular jet model,the bleed gases flow into the shear layer directly so that the hot gases are concentrated in the shear layer.The latter temperature distribution is better for the increase of base pressure.
基金the National Natural Science Foundation of Ghina （Funds Nos.49673189 , 49873018 ） the Doctoral Scientific Program Special Fund of the Education Department （Fund No.98055807） The State Key Laboratory for Mineral Deposits. ...Research, Nanjing University.
文摘The Ailaoshan aquamarine-bearing pegmatites are associated with Proterozoic metamorphic rocks in the southern portion of the Ailaoshan fault-folded complex.The gem-bearing pegmatite mineralization zones of the region occur in areas generally consistent with the regional tectonic trend.The pegmatites are found in metamorphic rocks,migmatites and in the inner/outer contact zones of gneissoid granites. The Rb-Sr isochron drawn for the pegmatites is 26～31 Ma,(i.e.in Himalayan).The homogenization temperatures of melt and liquid inclusions in minerals vary from 185 to 920℃,which are comparable to the inclusions observed in banded migmatites and ptygmatic quartz veins in the surrounding metamorphic rocks. The mineralization fluids of the pegmatite were rich in HCO_3 and CO_2,and their compositional assemblages are comparable to metamorphic fluids.Results of H,O,C,Si etc.isotopic analyses and REE,and Be analyses indicates that the sources of mineralization components that formed the pegmatites are closely associated with metamorphic fluids and the enclosing metamorphic rocks. A pegmatite structure simulation experiment was conducted at high temperature and pressure(840℃and 1,500×105Pa.),with various metamorphic rock samples in a water-rich and volatile-rich environment.When the liquidus was reached,the temperature was gradually decreased at the rate of 5～10℃/day over a time period of three months.SEM energy-dispersive spectrum analyses were performed on the experimental products.A series of pegmatoid textures were observed including zonal texture,megacryst texture,drusy cavities,crystal druses,and vesicular texture along with more than ten types of minerals including plagioclase,microcline,quartz and biotite.Different metamorphic rock melts generated different mineral assemblages.Experiment results revealed that the partial melting of metamorphic rocks could form melts similar to pegmatite magmas. Based upon the geological characteristics,geochemistry,and pegmatite texture simulation experimental results,it is concluded that the mineralization components of Ailaoshan aquamarine-bearing pegmatites came from metamorphic rocks.The petrogenetic model for the origin of pegmatites is related to ultrametamorphism and metamorphic anatexis.
