The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LW...The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LWD).Several numerical methods,including the method of moments(MoM),the electric field integral equation(EFIE)method,and the finite-element(FE)method have been developed for the simulation of EM telemetry systems.The computational process of these methods is complicated and time-consuming.To solve this problem,we introduce an axisymmetric semi-analytical FE method(SAFEM)in the cylindrical coordinate system with the virtual layering technique for rapid simulation of EM telemetry in a layered earth.The proposed method divides the computational domain into a series of homogeneous layers.For each layer,only its cross-section is discretized,and a high-precision integration method based on Riccati equations is employed for the calculation of longitudinally homogeneous sections.The block-tridiagonal structure of the global coefficient matrix enables the use of the block Thomas algorithm,facilitating the efficient simulation of EM telemetry problems in layered media.After the theoretical development,we validate the accuracy and efficiency of our algorithm through a series of numerical experiments and comparisons with the Multiphysics modeling software COMSOL.We also discussed the impact of system parameters on EM telemetry signal and demonstrated the applicability of our method by testing it on a field dataset acquired from Dezhou,Shandong Province,China.展开更多
This study conducted the experimental investigation of aerodynamic heating of Micro-scale Rotational Shearing Flow with Axial Limited-Length(MRSFALL).The temperature riseof the stator is captured by the high response ...This study conducted the experimental investigation of aerodynamic heating of Micro-scale Rotational Shearing Flow with Axial Limited-Length(MRSFALL).The temperature riseof the stator is captured by the high response thermocouples.The eccentricity ratio and clearanceheight are guaranteed by means of instantaneous trajectory and torsion monitoring of the rotator.The result shows that the maximum temperature rise takes place upstream of the minimum clear-ance height along circumferential direction.The distribution of temperature rise presents asymmet-ric curve along axial direction,and peak value occurs near the dimensionless axial position of-0.18.The effect of aerodynamic heating becomes notable as the rotational speed is larger than3×10^(4)r/min.The effect of end leakage and the viscous dissipation have great impact on temper-ature rise of MRSFALL.More specially,the peak value of temperature rise at dimensionless clear-ance height of 0.0080 is larger than the case at dimensionless clearance height of 0.0044.Furthermore,when the eccentricity ratio is too large,the viscous dissipation is induced,and theadditional temperature rise is achieved.The heat flux identification of shear flow has been realizedby Sequential Function Specification Method(SFSM)and its estimation of thermal load has been given.The heat flux induced by the aerodynamic heating in this study varies from 950 W/m^(2)to1330 W/m^(2).展开更多
Geomaterials with inferior hydraulic and strength characteristics often need improvement to enhance their engineering behaviors.Traditional ground improvement techniques require enormous mechanical effort or synthetic...Geomaterials with inferior hydraulic and strength characteristics often need improvement to enhance their engineering behaviors.Traditional ground improvement techniques require enormous mechanical effort or synthetic chemicals.Sustainable stabilization technique such as microbially induced calcite precipitation(MICP)utilizes bacterial metabolic processes to precipitate cementitious calcium carbonate.The reactive transport of biochemical species in the soil mass initiates the precipitation of biocement during the MICP process.The precipitated biocement alters the hydro-mechanical performance of the soil mass.Usually,the flow,deformation,and transport phenomena regulate the biocementation technique via coupled bio-chemo-hydro-mechanical(BCHM)processes.Among all,one crucial phenomenon controlling the precipitation mechanism is the encapsulation of biomass by calcium carbonate.Biomass encapsulation can potentially reduce the biochemical reaction rate and decelerate biocementation.Laboratory examination of the encapsulation process demands a thorough analysis of associated coupled effects.Despite this,a numerical model can assist in capturing the coupled processes influencing encapsulation during the MICP treatment.However,most numerical models did not consider biochemical reaction rate kinetics accounting for the influence of bacterial encapsulation.Given this,the current study developed a coupled BCHM model to evaluate the effect of encapsulation on the precipitated calcite content using a micro-scale semiempirical relationship.Firstly,the developed BCHM model was verified and validated using numerical and experimental observations of soil column tests.Later,the encapsulation phenomenon was investigated in the soil columns of variable maximum calcite crystal sizes.The results depict altered reaction rates due to the encapsulation phenomenon and an observable change in the precipitated calcite content for each maximum crystal size.Furthermore,the permeability and deformation of the soil mass were affected by the simultaneous precipitation of calcium carbonate.Overall,the present study comprehended the influence of the encapsulation of bacteria on cement morphology-induced permeability,biocement-induced stresses and displacements.展开更多
A finite-element model of the thermosetting epoxy asphalt mixture（EAM） microstructure is developed to simulate the indirect tension test（IDT）.Image techniques are used to capture the EAM microstructure which is di...A finite-element model of the thermosetting epoxy asphalt mixture（EAM） microstructure is developed to simulate the indirect tension test（IDT）.Image techniques are used to capture the EAM microstructure which is divided into two phases：aggregates and mastic.A viscoelastic constitutive relationship,which is obtained from the results of a creep test,is used to represent the mastic phase at intermittent temperatures.Model simulation results of the stiffness modulus in IDT compare favorably with experimental data.Different loading directions and velocities are employed in order to account for their influence on the modulus and the localized stress of the microstructure model.It is pointed out that the modulus is not consistent when the loading direction changes since the heterogeneous distribution of the mixture internal structure,and the loading velocity affects the localized stress as a result of the viscoelasticity of the mastic.The study results can provide a theoretical basis for the finite-element method,which can be extended to the numerical simulations of asphalt mixture micromechanical behavior.展开更多
The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-eleme...The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.展开更多
The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such...The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such as dipping interfaces and rough topography. We present an adaptive FE method for 2.5D forward modeling of induced polarization（IP）. In the presented method, an unstructured triangulation mesh that allows for local mesh refinement and flexible description of arbitrary model geometries is used. Furthermore, the mesh refinement process is guided by dual error estimate weighting to bias the refinement towards elements that affect the solution at the receiver locations. After the final mesh is generated, the Jacobian matrix is used to obtain the IP response on 2D structure models. We validate the adaptive FE algorithm using a vertical contact model. The validation shows that the elements near the receivers are highly refined and the average relative error of the potentials converges to 0.4 % and 1.2 % for the IP response. This suggests that the numerical solution of the adaptive FE algorithm converges to an accurate solution with the refined mesh. Finally, the accuracy and flexibility of the adaptive FE procedure are also validated using more complex models.展开更多
Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficien...Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.展开更多
A dry-gas seal system is a non-contact seal technology that is widely used in different industrial applications.Spiral-groove dry-gas seal utilizes fluid dynamic pressure effects to realize the seal and lubrication pr...A dry-gas seal system is a non-contact seal technology that is widely used in different industrial applications.Spiral-groove dry-gas seal utilizes fluid dynamic pressure effects to realize the seal and lubrication processes,while forming a high pressure gas film between two sealing faces due to the deceleration of the gas pumped in or out.There is little research into the effects and the influence on seal performance,if the grooves and the gas film are at the micro-scale.This paper investigates the micro-scale effects on spiral-groove dry-gas seal performance in a numerical solution of a corrected Reynolds equation.The Reynolds equation is discretized by means of the finite difference method with the second order scheme and solved by the successive-over-relaxation（SOR） iterative method.The Knudsen number of the flow in the sealing gas film is changed from 0.005 to 0.120 with a variation of film depth and sealing pressure.The numerical results show that the average pressure in the gas film and the sealed gas leakage increase due to micro-scale effects.The open force is enlarged,while the gas film stiffness is significantly decreased due to micro-scale effects.The friction torque and power consumption remain constant,even in low sealing pressure and spin speed conditions.In this paper,the seal performance at different rotor face spin speeds is also described.The proposed research clarifies the micro-scale effects in a spiral-groove dry-gas seal and their influence on seal performance,which is expected to be useful for the improvement of the design of dry-gas seal systems operating in the slip flow regime.展开更多
To understand the discharge characteristics under a gap of micrometers,the breakdown voltage and current–voltage curve are measured experimentally in a needle-to-plate electrode at a microscale gap of 3–50 μm in ai...To understand the discharge characteristics under a gap of micrometers,the breakdown voltage and current–voltage curve are measured experimentally in a needle-to-plate electrode at a microscale gap of 3–50 μm in air.The effect of the needle radius and the gas pressure on the discharge characteristics are tested.The results show that when the gap is larger than 10 μm,the relation between the breakdown voltage and the gap looks like the Paschen curve;while below 10 μm,the breakdown voltage is nearly constant in the range of the tested gap.However,at the same gap distance,the breakdown voltage is still affected by the pressure and shows a trend similar to Paschen＇s law.The current–voltage characteristic in all the gaps is similar and follows the trend of a typical Townsend-to-glow discharge.A simple model is used to explain the non-normality of breakdown in the micro-gaps.The Townsend mechanism is suggested to control the breakdown process in this configuration before the gap reduces much smaller in air.展开更多
A micro-scale finite element method(FEM) was proposed to precisely calculate the heat conduction between mortar and aggregate, and thus to accurately predict the non-uniformity of concrete pouring temperature. The con...A micro-scale finite element method(FEM) was proposed to precisely calculate the heat conduction between mortar and aggregate, and thus to accurately predict the non-uniformity of concrete pouring temperature. The concrete temperature field during vibration was also precisely calculated by accurate description of heat absorption characteristics of different parts of concrete when vibration. Based on the above method, the prediction model was used to predict the pouring temperature of a practical engineering. The comparison between actual results and simulated values shows that this method can be adopted to accurately predict the non-uniformity of concrete pouring temperature and the influence of mechanized vibration on concrete pouring temperature, and thus accurately predict pouring temperature. The control of casting temperature is crucial for preventing concrete fracture. The study provides a new method for predicting the pouring temperature of concrete structures, which has great practical value in engineering application.展开更多
Over the last decade, computational methods have been intensively applied to a variety of scientific researches and engineering designs. Although the computational fluid dynamics （CFD） method has played a dominant r...Over the last decade, computational methods have been intensively applied to a variety of scientific researches and engineering designs. Although the computational fluid dynamics （CFD） method has played a dominant role in studying and simulating transport phenomena involving fluid flow and heat and mass transfers, in recent years, other numerical methods for the simulations at meso- and micro-scales have also been actively applied to solve the physics of complex flow and fluid-interface interactions. This paper presents a review of recent advances in multi-scale computational simulation of biomimetics related fluid flow problems. The state-of-the-art numerical techniques, such as lattice Boltzmann method （LBM）, molecular dynamics （MD）, and conventional CFD, applied to different problems such as fish flow, electro-osmosis effect of earthworm motion, and self-cleaning hydrophobic surface, and the numerical approaches are introduced. The new challenging of modelling biomimetics problems in developing the physical conditions of self-clean hydrophobic surfaces is discussed.展开更多
The Luzon Island is a volcanic arc sandwiched by the eastward subducting South China Sea and the northwestward subducting Philippine Sea plate. Through experiments of plane-stress, elastic, and 2-dimensional finite-el...The Luzon Island is a volcanic arc sandwiched by the eastward subducting South China Sea and the northwestward subducting Philippine Sea plate. Through experiments of plane-stress, elastic, and 2-dimensional finite-element modeling, we evaluated the relationship between plate kinematics and present-day deformation of Luzon Island and adjacent sea areas. The concept of coupling rate was applied to define the boundary velocities along the subduction zones. The distribution of velocity fields calculated in our models was compared with the velocity field revealed by recent geodetic （GPS） observations. The best model was obtained that accounts for the observed velocity field within the limits of acceptable mechanical parameters and reasonable boundary conditions. Sensitivity of the selection of parameters and boundary conditions were evaluated. The model is sensitive to the direction of convergence between the South China Sea and the Philippine Sea plates, and to different coupling rates in the Manila trench, Philippine trench and eastern Luzon trough. We suggest that a change of ±15° of the direction of motion of the Philippine Sea plate can induce important changes in the distribution of the computed displacement trajectories, and the movement of the Philippine Sea plate toward azimuth 330° best explains the velocity pattern observed in Luzon Island. In addition, through sensitivity analysis we conclude that the coupling rate in the Manila trench is much smaller compared with the rates in the eastern Luzon trough and the Philippine trench. This indicates that a significant part of momentum of the Philippine Sea plate motion has been absorbed by the Manila trench; whereas, a part of the momentum has been transmitted into Luzon Island through the eastern Luzon trough and the Philippine trench.展开更多
In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can be...In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling （LWD）, in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.展开更多
Micro-scale Al-Zn-Mg/Fe composite powders (MAF) with high reactivity and good storage properties were prepared by reducing iron onto the surface of Al-Zn-Mg alloy powders. Experimental results show that MAF as advance...Micro-scale Al-Zn-Mg/Fe composite powders (MAF) with high reactivity and good storage properties were prepared by reducing iron onto the surface of Al-Zn-Mg alloy powders. Experimental results show that MAF as advanced zero-valent iron are highly effective for degradation of chlorinated organic compounds. The efficiency of degradation for carbon tetrachloride and perchloroethylene is higher than 99% within a period of 2 h. The efficiency of degradation for trichloroethylene by MAF after storing for one month is equivalent to that by freshly prepared nano-size zero-valent iron particles.展开更多
Micro-scale functionally graded material(FGM)pipes conveying fluid have many significant applications in engineering fields.In this work,the thermoelastic vibration of FGM fluid-conveying tubes in elastic medium is st...Micro-scale functionally graded material(FGM)pipes conveying fluid have many significant applications in engineering fields.In this work,the thermoelastic vibration of FGM fluid-conveying tubes in elastic medium is studied.Based on modified couple stress theory and Hamilton’s principle,the governing equation and boundary conditions are obtained.The differential quadrature method(DQM)is applied to investigating the thermoelastic vibration of the FGM pipes.The effect of temperature variation,scale effect of the microtubule,micro-fluid effect,material properties,elastic coefficient of elastic medium and outer radius on thermoelastic vibration of the FGM pipes conveying fluid are studied.The results show that in the condition of considering the scale effect and micro-fluid of the microtubule,the critical dimensionless velocity of the system is higher than that of the system which calculated using classical macroscopic model.The results also show that the variations of temperature,material properties,elastic coefficient and outer radius have significant influences on the first-order dimensionless natural frequency.展开更多
A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercriti...A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.展开更多
Fabrication of micro gratings using a femtosecond laser exposure system is experimentally investigated for the electron moire method. Micro holes and lines are firstly etched for parameter study. Grating profile is th...Fabrication of micro gratings using a femtosecond laser exposure system is experimentally investigated for the electron moire method. Micro holes and lines are firstly etched for parameter study. Grating profile is theoretically optimized to form high quality moire patterns. For a demonstration, a parallel grating is fabricated on a specimen of quartz glass. The minimum line width and the distance between two adjacent lines are both set to be 1 μm, and the frequency of grating is 500 lines/ram. The experimental results indicate that the quality of gratings is good and the relative error of the gratings pitch is about 1.5%. Based on molte method, scanning electron microscope （SEM） moire patterns are observed clearly, which manifests that gratings fabricated with the femtosecond laser exposure is suitable for micro scale deformation measurement.展开更多
A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been deve...A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.展开更多
Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of s...Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.展开更多
Stresses in a block around a dipping fracture simulating a damage zone of a fault are reconstructed by finite-element modeling. A fracture corresponding to a fault of different lengths, with its plane dipping at diffe...Stresses in a block around a dipping fracture simulating a damage zone of a fault are reconstructed by finite-element modeling. A fracture corresponding to a fault of different lengths, with its plane dipping at different angles, is assumed to follow a lithological interface and to experience either compression or shear. The stress associated with the destruction shows an asymmetrical pattern with different distances from the highest stress sites to the fault plane in the hanging and foot walls. As the dip angle decreases,the high-stress zone becomes wider in the hanging wall but its width changes negligibly in the foot wall.The length of the simulated fault and the deformation type affect only the magnitude of maximum stress,which remains asymmetrical relative to the fault plane. The Lh/Lfratio, where Lhand Lfare the widths of high-stress zones in the hanging and foot walls of the fault, respectively, is inversely proportional to the fault plane dip. The arithmetic mean of this ratio over different fault lengths in fractures subject to compression changes from 0.29 at a dip of 80°to 1.67 at 30°. In the case of shift displacement, ratios are increasing to 1.2 and 2.94, respectively.Usually they consider vertical fault planes and symmetry in a damage zone of faults. Following that assumption may cause errors in reconstructions of stress and fault patterns in areas of complex structural setting. According geological data, we know the structures are different and asymmetric in hanging and foot walls of fault. Thus, it is important to quantify zones of that asymmetry. The modeling results have to be taken into account in studies of natural faults, especially for practical applications in seismic risk mapping, engineering geology, hydrogeology, and tectonics.展开更多
基金supported by the Major Research Project on Scientific Instrument Development of the National Natural Science Foundation of China(42327901)National Natural Science Foundation of China(42030806,42074120,41904104,423B2405).
文摘The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LWD).Several numerical methods,including the method of moments(MoM),the electric field integral equation(EFIE)method,and the finite-element(FE)method have been developed for the simulation of EM telemetry systems.The computational process of these methods is complicated and time-consuming.To solve this problem,we introduce an axisymmetric semi-analytical FE method(SAFEM)in the cylindrical coordinate system with the virtual layering technique for rapid simulation of EM telemetry in a layered earth.The proposed method divides the computational domain into a series of homogeneous layers.For each layer,only its cross-section is discretized,and a high-precision integration method based on Riccati equations is employed for the calculation of longitudinally homogeneous sections.The block-tridiagonal structure of the global coefficient matrix enables the use of the block Thomas algorithm,facilitating the efficient simulation of EM telemetry problems in layered media.After the theoretical development,we validate the accuracy and efficiency of our algorithm through a series of numerical experiments and comparisons with the Multiphysics modeling software COMSOL.We also discussed the impact of system parameters on EM telemetry signal and demonstrated the applicability of our method by testing it on a field dataset acquired from Dezhou,Shandong Province,China.
基金supports from the National Natural Science Foundation of China(No.52206091)the Aeronautical Science Foundation of China(No.201928052008)the Natural Science Foundation of Jiangsu Province,China(No.BK20210303)。
文摘This study conducted the experimental investigation of aerodynamic heating of Micro-scale Rotational Shearing Flow with Axial Limited-Length(MRSFALL).The temperature riseof the stator is captured by the high response thermocouples.The eccentricity ratio and clearanceheight are guaranteed by means of instantaneous trajectory and torsion monitoring of the rotator.The result shows that the maximum temperature rise takes place upstream of the minimum clear-ance height along circumferential direction.The distribution of temperature rise presents asymmet-ric curve along axial direction,and peak value occurs near the dimensionless axial position of-0.18.The effect of aerodynamic heating becomes notable as the rotational speed is larger than3×10^(4)r/min.The effect of end leakage and the viscous dissipation have great impact on temper-ature rise of MRSFALL.More specially,the peak value of temperature rise at dimensionless clear-ance height of 0.0080 is larger than the case at dimensionless clearance height of 0.0044.Furthermore,when the eccentricity ratio is too large,the viscous dissipation is induced,and theadditional temperature rise is achieved.The heat flux identification of shear flow has been realizedby Sequential Function Specification Method(SFSM)and its estimation of thermal load has been given.The heat flux induced by the aerodynamic heating in this study varies from 950 W/m^(2)to1330 W/m^(2).
基金the funding support from the Ministry of Education,Government of India,under the Prime Minister Research Fellowship programme(Grant Nos.SB21221901CEPMRF008347 and SB22230217CEPMRF008347).
文摘Geomaterials with inferior hydraulic and strength characteristics often need improvement to enhance their engineering behaviors.Traditional ground improvement techniques require enormous mechanical effort or synthetic chemicals.Sustainable stabilization technique such as microbially induced calcite precipitation(MICP)utilizes bacterial metabolic processes to precipitate cementitious calcium carbonate.The reactive transport of biochemical species in the soil mass initiates the precipitation of biocement during the MICP process.The precipitated biocement alters the hydro-mechanical performance of the soil mass.Usually,the flow,deformation,and transport phenomena regulate the biocementation technique via coupled bio-chemo-hydro-mechanical(BCHM)processes.Among all,one crucial phenomenon controlling the precipitation mechanism is the encapsulation of biomass by calcium carbonate.Biomass encapsulation can potentially reduce the biochemical reaction rate and decelerate biocementation.Laboratory examination of the encapsulation process demands a thorough analysis of associated coupled effects.Despite this,a numerical model can assist in capturing the coupled processes influencing encapsulation during the MICP treatment.However,most numerical models did not consider biochemical reaction rate kinetics accounting for the influence of bacterial encapsulation.Given this,the current study developed a coupled BCHM model to evaluate the effect of encapsulation on the precipitated calcite content using a micro-scale semiempirical relationship.Firstly,the developed BCHM model was verified and validated using numerical and experimental observations of soil column tests.Later,the encapsulation phenomenon was investigated in the soil columns of variable maximum calcite crystal sizes.The results depict altered reaction rates due to the encapsulation phenomenon and an observable change in the precipitated calcite content for each maximum crystal size.Furthermore,the permeability and deformation of the soil mass were affected by the simultaneous precipitation of calcium carbonate.Overall,the present study comprehended the influence of the encapsulation of bacteria on cement morphology-induced permeability,biocement-induced stresses and displacements.
基金Program for New Century Excellent Talents in University(No. NCET-08-0118)Specialized Research Fund for the Doctoral Program of Higher Education (No. 20090092110049)
文摘A finite-element model of the thermosetting epoxy asphalt mixture（EAM） microstructure is developed to simulate the indirect tension test（IDT）.Image techniques are used to capture the EAM microstructure which is divided into two phases：aggregates and mastic.A viscoelastic constitutive relationship,which is obtained from the results of a creep test,is used to represent the mastic phase at intermittent temperatures.Model simulation results of the stiffness modulus in IDT compare favorably with experimental data.Different loading directions and velocities are employed in order to account for their influence on the modulus and the localized stress of the microstructure model.It is pointed out that the modulus is not consistent when the loading direction changes since the heterogeneous distribution of the mixture internal structure,and the loading velocity affects the localized stress as a result of the viscoelasticity of the mastic.The study results can provide a theoretical basis for the finite-element method,which can be extended to the numerical simulations of asphalt mixture micromechanical behavior.
基金sponsored by the National Natural Science Foundation of China Research(Grant No.41274138)the Science Foundation of China University of Petroleum(Beijing)(No.KYJJ2012-05-02)
文摘The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.
基金financially supported by the National Natural Science Foundation of China(No.41204055,41164003,and 41104074)Opening Project(No.SMIL-2014-06) of Hubei Subsurface Multi-scale Imaging Lab(SMIL),China University of Geosciences(Wuhan)
文摘The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such as dipping interfaces and rough topography. We present an adaptive FE method for 2.5D forward modeling of induced polarization（IP）. In the presented method, an unstructured triangulation mesh that allows for local mesh refinement and flexible description of arbitrary model geometries is used. Furthermore, the mesh refinement process is guided by dual error estimate weighting to bias the refinement towards elements that affect the solution at the receiver locations. After the final mesh is generated, the Jacobian matrix is used to obtain the IP response on 2D structure models. We validate the adaptive FE algorithm using a vertical contact model. The validation shows that the elements near the receivers are highly refined and the average relative error of the potentials converges to 0.4 % and 1.2 % for the IP response. This suggests that the numerical solution of the adaptive FE algorithm converges to an accurate solution with the refined mesh. Finally, the accuracy and flexibility of the adaptive FE procedure are also validated using more complex models.
基金National High Technology Research and Development Program(863 Program)(No.2006AA06Z105,2007AA06Z134)
文摘Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.
基金supported by Scientific Research Foundation for Returned Scholars of Ministry of Education of China
文摘A dry-gas seal system is a non-contact seal technology that is widely used in different industrial applications.Spiral-groove dry-gas seal utilizes fluid dynamic pressure effects to realize the seal and lubrication processes,while forming a high pressure gas film between two sealing faces due to the deceleration of the gas pumped in or out.There is little research into the effects and the influence on seal performance,if the grooves and the gas film are at the micro-scale.This paper investigates the micro-scale effects on spiral-groove dry-gas seal performance in a numerical solution of a corrected Reynolds equation.The Reynolds equation is discretized by means of the finite difference method with the second order scheme and solved by the successive-over-relaxation（SOR） iterative method.The Knudsen number of the flow in the sealing gas film is changed from 0.005 to 0.120 with a variation of film depth and sealing pressure.The numerical results show that the average pressure in the gas film and the sealed gas leakage increase due to micro-scale effects.The open force is enlarged,while the gas film stiffness is significantly decreased due to micro-scale effects.The friction torque and power consumption remain constant,even in low sealing pressure and spin speed conditions.In this paper,the seal performance at different rotor face spin speeds is also described.The proposed research clarifies the micro-scale effects in a spiral-groove dry-gas seal and their influence on seal performance,which is expected to be useful for the improvement of the design of dry-gas seal systems operating in the slip flow regime.
基金supported by National Natural Science Foundation of China(11475019)the Electrostatic Research Foundation of Liu Shanghe Academicians and Experts Workstation,Beijing Orient Institute of Measurement and Test(BOIMTLSHJD20161002)
文摘To understand the discharge characteristics under a gap of micrometers,the breakdown voltage and current–voltage curve are measured experimentally in a needle-to-plate electrode at a microscale gap of 3–50 μm in air.The effect of the needle radius and the gas pressure on the discharge characteristics are tested.The results show that when the gap is larger than 10 μm,the relation between the breakdown voltage and the gap looks like the Paschen curve;while below 10 μm,the breakdown voltage is nearly constant in the range of the tested gap.However,at the same gap distance,the breakdown voltage is still affected by the pressure and shows a trend similar to Paschen＇s law.The current–voltage characteristic in all the gaps is similar and follows the trend of a typical Townsend-to-glow discharge.A simple model is used to explain the non-normality of breakdown in the micro-gaps.The Townsend mechanism is suggested to control the breakdown process in this configuration before the gap reduces much smaller in air.
基金Supported by the National Key Research and Development Project of China(No.2018YFC0406703)the National Natural Science Foundation of China(Nos.51779277,51579252,51439005)+2 种基金the Special Scientific Research Project of the State Key Laboratory of Simulation and Regulation of Water Cycles in River Basins(No.2016ZY10)a Special Scientific Research Project of the China Institute of Water Resources and Hydropower Research(Nos.SS0145B392016,SS0145B612017)the Special Scientific Research Project of the China Institute of Water Resources and Hydropower Research(No.KY1799).
文摘A micro-scale finite element method(FEM) was proposed to precisely calculate the heat conduction between mortar and aggregate, and thus to accurately predict the non-uniformity of concrete pouring temperature. The concrete temperature field during vibration was also precisely calculated by accurate description of heat absorption characteristics of different parts of concrete when vibration. Based on the above method, the prediction model was used to predict the pouring temperature of a practical engineering. The comparison between actual results and simulated values shows that this method can be adopted to accurately predict the non-uniformity of concrete pouring temperature and the influence of mechanized vibration on concrete pouring temperature, and thus accurately predict pouring temperature. The control of casting temperature is crucial for preventing concrete fracture. The study provides a new method for predicting the pouring temperature of concrete structures, which has great practical value in engineering application.
文摘Over the last decade, computational methods have been intensively applied to a variety of scientific researches and engineering designs. Although the computational fluid dynamics （CFD） method has played a dominant role in studying and simulating transport phenomena involving fluid flow and heat and mass transfers, in recent years, other numerical methods for the simulations at meso- and micro-scales have also been actively applied to solve the physics of complex flow and fluid-interface interactions. This paper presents a review of recent advances in multi-scale computational simulation of biomimetics related fluid flow problems. The state-of-the-art numerical techniques, such as lattice Boltzmann method （LBM）, molecular dynamics （MD）, and conventional CFD, applied to different problems such as fish flow, electro-osmosis effect of earthworm motion, and self-cleaning hydrophobic surface, and the numerical approaches are introduced. The new challenging of modelling biomimetics problems in developing the physical conditions of self-clean hydrophobic surfaces is discussed.
基金supported by the National Scientific and Tech-nological Support Project of China (No. 2006BAB19B02)the Guangdong Natural Science Foundation (No. 07004206)the National Natural Science Foundation of China (No. 40476026)
文摘The Luzon Island is a volcanic arc sandwiched by the eastward subducting South China Sea and the northwestward subducting Philippine Sea plate. Through experiments of plane-stress, elastic, and 2-dimensional finite-element modeling, we evaluated the relationship between plate kinematics and present-day deformation of Luzon Island and adjacent sea areas. The concept of coupling rate was applied to define the boundary velocities along the subduction zones. The distribution of velocity fields calculated in our models was compared with the velocity field revealed by recent geodetic （GPS） observations. The best model was obtained that accounts for the observed velocity field within the limits of acceptable mechanical parameters and reasonable boundary conditions. Sensitivity of the selection of parameters and boundary conditions were evaluated. The model is sensitive to the direction of convergence between the South China Sea and the Philippine Sea plates, and to different coupling rates in the Manila trench, Philippine trench and eastern Luzon trough. We suggest that a change of ±15° of the direction of motion of the Philippine Sea plate can induce important changes in the distribution of the computed displacement trajectories, and the movement of the Philippine Sea plate toward azimuth 330° best explains the velocity pattern observed in Luzon Island. In addition, through sensitivity analysis we conclude that the coupling rate in the Manila trench is much smaller compared with the rates in the eastern Luzon trough and the Philippine trench. This indicates that a significant part of momentum of the Philippine Sea plate motion has been absorbed by the Manila trench; whereas, a part of the momentum has been transmitted into Luzon Island through the eastern Luzon trough and the Philippine trench.
基金supported by the National Natural Science Foundation of China(No.41204094)Science Foundation of China University of Petroleum,Beijing(No.2462015YQ0506)
文摘In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling （LWD）, in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.
文摘Micro-scale Al-Zn-Mg/Fe composite powders (MAF) with high reactivity and good storage properties were prepared by reducing iron onto the surface of Al-Zn-Mg alloy powders. Experimental results show that MAF as advanced zero-valent iron are highly effective for degradation of chlorinated organic compounds. The efficiency of degradation for carbon tetrachloride and perchloroethylene is higher than 99% within a period of 2 h. The efficiency of degradation for trichloroethylene by MAF after storing for one month is equivalent to that by freshly prepared nano-size zero-valent iron particles.
文摘Micro-scale functionally graded material(FGM)pipes conveying fluid have many significant applications in engineering fields.In this work,the thermoelastic vibration of FGM fluid-conveying tubes in elastic medium is studied.Based on modified couple stress theory and Hamilton’s principle,the governing equation and boundary conditions are obtained.The differential quadrature method(DQM)is applied to investigating the thermoelastic vibration of the FGM pipes.The effect of temperature variation,scale effect of the microtubule,micro-fluid effect,material properties,elastic coefficient of elastic medium and outer radius on thermoelastic vibration of the FGM pipes conveying fluid are studied.The results show that in the condition of considering the scale effect and micro-fluid of the microtubule,the critical dimensionless velocity of the system is higher than that of the system which calculated using classical macroscopic model.The results also show that the variations of temperature,material properties,elastic coefficient and outer radius have significant influences on the first-order dimensionless natural frequency.
基金the National Science Council ot Taiwan,China for funding this research(Project no.:NSC 94-2218-E-035-011)
文摘A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.
基金support from the National Natural Science Foundation of China (11372118 and 11302082)
文摘Fabrication of micro gratings using a femtosecond laser exposure system is experimentally investigated for the electron moire method. Micro holes and lines are firstly etched for parameter study. Grating profile is theoretically optimized to form high quality moire patterns. For a demonstration, a parallel grating is fabricated on a specimen of quartz glass. The minimum line width and the distance between two adjacent lines are both set to be 1 μm, and the frequency of grating is 500 lines/ram. The experimental results indicate that the quality of gratings is good and the relative error of the gratings pitch is about 1.5%. Based on molte method, scanning electron microscope （SEM） moire patterns are observed clearly, which manifests that gratings fabricated with the femtosecond laser exposure is suitable for micro scale deformation measurement.
基金Projects(41674080,41674079)supported by the National Natural Science Foundation of China
文摘A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.
基金the National Science Council of Taiwan for funding this research (NSC 96-2221-E-019-061).
文摘Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.
文摘Stresses in a block around a dipping fracture simulating a damage zone of a fault are reconstructed by finite-element modeling. A fracture corresponding to a fault of different lengths, with its plane dipping at different angles, is assumed to follow a lithological interface and to experience either compression or shear. The stress associated with the destruction shows an asymmetrical pattern with different distances from the highest stress sites to the fault plane in the hanging and foot walls. As the dip angle decreases,the high-stress zone becomes wider in the hanging wall but its width changes negligibly in the foot wall.The length of the simulated fault and the deformation type affect only the magnitude of maximum stress,which remains asymmetrical relative to the fault plane. The Lh/Lfratio, where Lhand Lfare the widths of high-stress zones in the hanging and foot walls of the fault, respectively, is inversely proportional to the fault plane dip. The arithmetic mean of this ratio over different fault lengths in fractures subject to compression changes from 0.29 at a dip of 80°to 1.67 at 30°. In the case of shift displacement, ratios are increasing to 1.2 and 2.94, respectively.Usually they consider vertical fault planes and symmetry in a damage zone of faults. Following that assumption may cause errors in reconstructions of stress and fault patterns in areas of complex structural setting. According geological data, we know the structures are different and asymmetric in hanging and foot walls of fault. Thus, it is important to quantify zones of that asymmetry. The modeling results have to be taken into account in studies of natural faults, especially for practical applications in seismic risk mapping, engineering geology, hydrogeology, and tectonics.
