Phonon polaritons(PhPs)exhibit directional in-plane propagation and ultralow losses in van der Waals(vdW)crystals,offering new possibilities for controlling the flow of light at the nanoscale.However,these PhPs,includ...Phonon polaritons(PhPs)exhibit directional in-plane propagation and ultralow losses in van der Waals(vdW)crystals,offering new possibilities for controlling the flow of light at the nanoscale.However,these PhPs,including their directional propagation,are inherently determined by the anisotropic crystal structure of the host materials.Although in-plane anisotropic PhPs can be manipulated by twisting engineering,such as twisting individual vdW slabs,dynamically adjusting their propagation presents a significant challenge.The limited application of the twisted bilayer structure in bare films further restricts its usage.In this study,we present a technique in which anisotropic PhPs supported by bare biaxial vdW slabs can be actively tuned by modifying their local dielectric environment.Excitingly,we predict that the iso-frequency contour of PhPs can be reoriented to enable propagation along forbidden directions when the crystal is placed on a substrate with a moderate negative permittivity.Besides,we systematically investigate the impact of polaritonic coupling on near-field radiative heat transfer(NFRHT)between heterostructures integrated with different substrates that have negative permittivity.Our main findings reveal that through the analysis of dispersion contour and photon transmission coefficient,the excitation and reorientation of the fundamental mode facilitate increased photon tunneling,thereby enhancing heat transfer between heterostructures.Conversely,the annihilation of the fundamental mode hinders heat transfer.Furthermore,we find the enhancement or suppression of radiative energy transport depends on the relative magnitude of the slab thickness and the vacuum gap width.Finally,the effect of negative permittivity substrates on NFRHT along the[001]crystalline direction ofα-MoO3 is considered.The spectral band where the excited fundamental mode resulting from the negative permittivity substrates is shifted to the first Reststrahlen Band(RB 1)ofα-MoO_(3) and is widened,resulting in more significant enhancement of heat flux from RB 1.We anticipate our results will motivate new direction for dynamical tunability of the PhPs in photonic devices.展开更多
The approximate compressible model is adopted to study the effects of strength and compressibility on the penetration by WHA long rod and copper jet into semi-infinite target in detail. For WHA rod penetrating PMMA at...The approximate compressible model is adopted to study the effects of strength and compressibility on the penetration by WHA long rod and copper jet into semi-infinite target in detail. For WHA rod penetrating PMMA at 2 km/s <V <5 km/s, the compressibility has a significant effect on the penetration efficiency. We clarify how compressibility affects the penetration efficiency by changing the stagnation pressures of the rod and target. For WHA rod penetrating 4340 Steel and 6061-T6 Al at 2 km/s < V < 10 km/s, the effect of strength is strong and the effect of compressibility is negligible at lower impact velocity, whilst the effect of strength is weak and the effect of compressibility becomes stronger at higher impact velocity. For the copper jet penetrating 4030 Steel, 6061-T6 Al and PMMA. the virtual origin model is adopted, and the compressibility and strength are implicitly considered by the linear relation between the penetration velocity and impact velocity. The effects of compressibility and target resistance on penetration efficiency are studied. The results show that the target resistance has a significant effect on the penetration efficiency. Howver PMMA is much more compressible than copper and the huge difference of compressibility has a significant effect on the penetration by hypervelocity copper jet into PMMA.展开更多
By solving the Boltzmann transport equation and considering the spin-dependent grain boundary scattering, the distribution of electrons in grains and the electrical transport properties in the applied magnetic field a...By solving the Boltzmann transport equation and considering the spin-dependent grain boundary scattering, the distribution of electrons in grains and the electrical transport properties in the applied magnetic field are studied. With regard to the dominant influence of grain boundary scattering which is taken as a boundary condition for the electrical transport, the grain size-dependent electrical conductivity is investigated. In addition, the reorientation of the relative magnetization between grains brings the change of the electron spin when the magnetonanocrystailine material is subjected to the magnetic field, resulting in the remarkable giant magnetoresistance effect.展开更多
We study the stability of vortices pinning and dynamics in a superconducting thin strip containing a square array of antidot triplets by using the nonlinear Ginzburg–Landau(GL)theory.Compared with the regular square ...We study the stability of vortices pinning and dynamics in a superconducting thin strip containing a square array of antidot triplets by using the nonlinear Ginzburg–Landau(GL)theory.Compared with the regular square array of circular holes,the vortices are no longer pinned inside the circular holes,but instead stabilized at the center of the antidot triplets depending on the geometry parameters.Moreover,the influences of the geometry parameters and the polarity of the applied current on the current–voltage(I–V)characteristics are also studied.The critical current for the sample turning into a normal state becomes smaller when the hole diameter D is smaller and the spacing B between the holes is larger.Due to the asymmetric pinning sites,our numerical simulations demonstrate that the positive and negative rectified voltages appear alternately in the resistive state of the sample under an ac current of square pulses.展开更多
Spectral and directional control of thermal emission based on excitation of confined electromagnetic resonant modes paves a viable way for the design and construction of microscale thermal emitters/absorbers. In this ...Spectral and directional control of thermal emission based on excitation of confined electromagnetic resonant modes paves a viable way for the design and construction of microscale thermal emitters/absorbers. In this paper, we present numerical simulation results of the thermal radiative properties of a silicon carbide（Si C） thermal emitter/absorber composed of periodic microstructures. We illustrate different electromagnetic resonant modes which can be excited with the structure,such as surface phonon polaritons, magnetic polaritons and photonic crystal modes, and the process of radiation spectrum optimization based on a non-linear optimization algorithm. We show that the spectral and directional control of thermal emission/absorption can be efficiently achieved by adjusting the geometrical parameters of the structure. Moreover, the optimized spectrum is insensitive to 3% dimension modification.展开更多
By using MTS815 rock mechanics test system,a series of acoustic emission(AE) location experiments were performed under unloading confining pressure,increasing the axial stress.The AE space-time evolution regularities ...By using MTS815 rock mechanics test system,a series of acoustic emission(AE) location experiments were performed under unloading confining pressure,increasing the axial stress.The AE space-time evolution regularities and energy releasing characteristics during deformation and failure process of coal of different loading rates are compared,the influence mechanism of loading rates on the microscopic crack evolution were studied,combining the AE characteristics and the macroscopic failure modes of the specimens,and the precursory characteristics of coal failure were also analyzed quantitatively.The results indicate that as the loading rate is higher,the AE activity and the main fracture will begin earlier.The destruction of coal body is mainly the function of shear strain at lower loading rate and tension strain at higher rate,and will transform from brittleness to ductility at critical velocities.When the deformation of the coal is mainly plasticity,the amplitude of the AE ringing counting rate increases largely and the AE energy curves appear an obvious ''step'',which can be defined as the first failure precursor point.Statics of AE information shows that the strongest AE activity begins when the axial stress level was 92-98%,which can be defined as the other failure precursor point.As the loading rate is smaller,the coal more easily reaches the latter precursor point after the first one,so attention should be aroused to prevent dynamic disaster in coal mining when the AE activity reaches the first precursor point.展开更多
A numerical method for simulating the motion and deformation of an axisymmetric bubble or drop rising or falling in another infinite and initially stationary fluid is developed based on the volume of fluid （VOF） met...A numerical method for simulating the motion and deformation of an axisymmetric bubble or drop rising or falling in another infinite and initially stationary fluid is developed based on the volume of fluid （VOF） method in the frame of two incompressible and immiscible viscous fluids under the action of gravity, taking into consideration of surface tension effects. A comparison of the numerical results by this method with those by other works indicates the validity of the method. In the frame of inviseid and incompressible fluids without taking into consideration of surface tension effects, the mechanisms of the generation of the liquid jet and the transition from spherical shape to toroidal shape during the bubble or drop deformation, the increase of the ring diameter of the toroidal bubble or drop and the decrease of its cross-section area during its motion, and the effects of the density ratio of the two fluids on the deformation of the bubble or drop are analysed both theoretically and numerically.展开更多
The rolling massage manipulation is a classic Chinese massage, which is expected to improve the circulation by pushing, pulling and kneading of the muscle. A model for the rolling massage manipulation is proposed and ...The rolling massage manipulation is a classic Chinese massage, which is expected to improve the circulation by pushing, pulling and kneading of the muscle. A model for the rolling massage manipulation is proposed and the lattice Boltzmann method is applied to study the blood flow in the blood vessels. The simulation results show that the blood flux is considerably modified by the rolling massage and the explicit value depends on the rolling frequency, the rolling depth, and the diameter of the vessel. The smaller the diameter of the blood vessel, the larger the enhancement of the blood flux by the rolling massage. The model, together with the simulation results, is expected to be helpful to understand the mechanism and further development of rolling massage techniques.展开更多
Liquid-filled compartment structure consists of a bulk steel plate with matrix blind holes which are filled with liquid and a steel front plate to seal up the liquid with rings and bolts.The liquid-filled compartment ...Liquid-filled compartment structure consists of a bulk steel plate with matrix blind holes which are filled with liquid and a steel front plate to seal up the liquid with rings and bolts.The liquid-filled compartment structure can resist the shaped charge warhead effectively.This paper presents experimental and theoretical investigations of the penetration ability of the residual shaped charge jet emerging from the liquid-filled compartment structure after the penetration process at different impact angles.On the basis of shock wave propagation theory,the influence of the liquid-filled compartment structure on jet stability is analysed.The interferences of the liquid backflow caused by a reflected shock wave and a back plate on jet stability under different impact angles are also examined.In addition,the range of the disturbed velocity segments of the jet at different impact angles and the penetration ability of the residual jet are obtained.A theoretical model is validated against the experimental penetration depths.展开更多
An experimental study on acceleration mechanism of flame propagation of propane-air mixture in ducts with obstacles was conducted. The acceleration mechanism of flame propagation is mainly due to the positive feedback...An experimental study on acceleration mechanism of flame propagation of propane-air mixture in ducts with obstacles was conducted. The acceleration mechanism of flame propagation is mainly due to the positive feedback of the turbulence region induced by obstacles for combustion process. It can be seen from the experimental results that the maximum explosion pressure can increase by 20%, the maximum rate of pressure rise can increase by 10 times and the flame propagation velocity can increase by 20 times when obstacles are present.展开更多
This paper studies the roughness effect combining with effects of rarefaction and compressibility by a lattice Boltzmann model for rarefied gas flows at high Knudsen numbers. By discussing the effect of the tangential...This paper studies the roughness effect combining with effects of rarefaction and compressibility by a lattice Boltzmann model for rarefied gas flows at high Knudsen numbers. By discussing the effect of the tangential momentum accommodation coefficient on the rough boundary condition, the lattice Boltzmann simulations of nitrogen and helium flows are performed in a two-dimensional microchannel with rough boundaries. The surface roughness effects in the microchannel on the velocity field, the mass flow rate and the friction coefficient are studied and analysed. Numerical results for the two gases in micro scale show different characteristics from macroscopic flows and demonstrate the feasibility of the lattice Boltzmann model in rarefied gas dynamics.展开更多
Molecular dynamics simulations are carried out to study atomic diffusion in the explosive welding process of NisoTis0-Cu （at.%）. By using a hybrid method which combines molecular dynamics simulation and classical di...Molecular dynamics simulations are carried out to study atomic diffusion in the explosive welding process of NisoTis0-Cu （at.%）. By using a hybrid method which combines molecular dynamics simulation and classical diffusion the- ory, the thickness of the diffusion layer and the atomic concentration distribution across the welding interface are obtained. The results indicate that the concentration distribution curves at different times have a geometric similarity. According to the geometric similarity, the atomic concentration distribution at any time in explosive welding can be calculated. NisoTis0- Cu explosive welding and scanning electron microscope experiments are done to verify the results. The simulation results and the experimental results are in good agreement.展开更多
The particle path tracking method is proposed and used in two-dimensional（2D） and three-dimensional（3D） numerical simulations of continuously rotating detonation engines（CRDEs）. This method is used to analyze th...The particle path tracking method is proposed and used in two-dimensional（2D） and three-dimensional（3D） numerical simulations of continuously rotating detonation engines（CRDEs）. This method is used to analyze the combustion and expansion processes of the fresh particles, and the thermodynamic cycle process of CRDE. In a 3D CRDE flow field, as the radius of the annulus increases, the no-injection area proportion increases, the non-detonation proportion decreases, and the detonation height decreases. The flow field parameters on the 3D mid annulus are different from in the 2D flow field under the same chamber size. The non-detonation proportion in the 3D flow field is less than in the 2D flow field. In the 2D and 3D CRDE, the paths of the flow particles have only a small fluctuation in the circumferential direction. The numerical thermodynamic cycle processes are qualitatively consistent with the three ideal cycle models, and they are right in between the ideal F–J cycle and ideal ZND cycle. The net mechanical work and thermal efficiency are slightly smaller in the 2D simulation than in the 3D simulation. In the 3D CRDE, as the radius of the annulus increases, the net mechanical work is almost constant, and the thermal efficiency increases. The numerical thermal efficiencies are larger than F–J cycle, and much smaller than ZND cycle.展开更多
Particulate composites are one of the widely used materials in producing numerous state-of-the-art components in biomedical,automobile,aerospace including defence technology.Variety of modelling techniques have been a...Particulate composites are one of the widely used materials in producing numerous state-of-the-art components in biomedical,automobile,aerospace including defence technology.Variety of modelling techniques have been adopted in the past to model mechanical behaviour of particulate composites.Due to their favourable properties,particle-based methods provide a convenient platform to model failure or fracture of these composites.Smooth particle hydrodynamics(SPH)is one of such methods which demonstrate excellent potential for modelling failure or fracture of particulate composites in a Lagrangian setting.One of the major challenges in using SPH method for modelling composite materials depends on accurate and efficient way to treat interface and boundary conditions.In this paper,a masterslave method based multi-freedom constraints is proposed to impose essential boundary conditions and interfacial displacement constraints in modelling mechanical behaviour of composite materials using SPH method.The proposed methodology enforces the above constraints more accurately and requires only smaller condition number for system stiffness matrix than the procedures based on typical penalty function approach.A minimum cut-off value-based error criteria is employed to improve the computational efficiency of the proposed methodology.In addition,the proposed method is further enhanced by adopting a modified numerical interpolation scheme along the boundary to increase the accuracy and computational efficiency.The numerical examples demonstrate that the proposed master-slave approach yields better accuracy in enforcing displacement constraints and requires approximately the same computational time as that of penalty method.展开更多
The relative scaling exponents and intermittency of three-dimensional compressible turbulent channel flow are investigated by using direct numerical simulation. One case is subsonic flow （Mα=0.8）, the other is supe...The relative scaling exponents and intermittency of three-dimensional compressible turbulent channel flow are investigated by using direct numerical simulation. One case is subsonic flow （Mα=0.8）, the other is supersonic （Mα=1.3）, and the Reynolds numbers based on the mean bulk velocity and channel half-width are 2826 and 3010, respectively. The analysis of the local slopes of sixth order velocity structure function to third order reveals that there is a well-defined scaling range for 10〈y^＋〈100. it is also noted that the intermittency of longitudinal velocity increments in this region is stronger than that of the transverse ones. Comparison with the incompressible case shows that the location of the most intensive intermittency moves toward the log-law region, which is related to the displacement of streamwise vortical structures in the near-wall region.展开更多
In this paper, a novel magnetoelectric（ME） composite structure is proposed, and the ME response in the structure is measured at the bias magnetic field up to 2000 Oe（1 Oe = 79.5775 A·m^（-1）） and the excitat...In this paper, a novel magnetoelectric（ME） composite structure is proposed, and the ME response in the structure is measured at the bias magnetic field up to 2000 Oe（1 Oe = 79.5775 A·m^（-1）） and the excitation frequency of alternating magnetic field ranging from 1 kHz to 200 kHz. The ME voltage of each PZT layer is detected. According to the measurement results, the phase differences are observed among three channels and the multi-peak phenomenon appears in each channel. Meanwhile, the results show that the ME structure can stay a relatively high ME response within a wide bandwidth.Besides, the hysteretic loops of three PZT layers are observed. When the frequency of alternating current（AC） magnetic field changes, the maximum value of ME coefficient appears in different layers due to the multiple vibration modes of the structure. Moreover, a finite element analysis is performed to evaluate the resonant frequency of the structure, and the theoretical calculating results accord well with the experimental results. The experiment results suggest that the proposed structure may be a good candidate for designing broadband magnetic field sensors.展开更多
The problem involving an edge-crack in a rectangular material under the anti-plane mechanical loading and in-plane electric loading is analyzed under the impermeable conditions. By using the series expansion, the gene...The problem involving an edge-crack in a rectangular material under the anti-plane mechanical loading and in-plane electric loading is analyzed under the impermeable conditions. By using the series expansion, the general solutions of electromechanical fields are obtained, which satisfied both governing equations and crack sufrace boundary conditions, and the unknown constants in which can be obtained by the boundary collocation method. Numerical results are given to show the effect of electromechanical interaction on energy release rate.展开更多
A new method for three-dimensional simulation of the interaction between the gas and the solid around is developed.The effects of the gas on the thermal-mechanical behaviors within the surrounded solid are performed b...A new method for three-dimensional simulation of the interaction between the gas and the solid around is developed.The effects of the gas on the thermal-mechanical behaviors within the surrounded solid are performed by replacing the internal gas with an equivalent solid in the modeling,which can make it convenient to simulate the thermal-mechanical coupling effects in the solid research objects with gases in them.The applied thermal expansion coefficient,Young's modulus and Poisson's ratio of the equivalent solid material are derived.A series of tests have been conducted;and the proposed equivalent solid method to simulate the gas effects is validated.展开更多
The multi-bifurcation effect of blood flow is investigated by lattice Boltzmann method at Re = 200 with six different bifurcation angles α which are 22.5°, 25°, 28°, 30°, 33°, 35°, respe...The multi-bifurcation effect of blood flow is investigated by lattice Boltzmann method at Re = 200 with six different bifurcation angles α which are 22.5°, 25°, 28°, 30°, 33°, 35°, respectively. The velocities and ratios of average velocity at various bifurcations are discussed. It is indicated that the maximum velocity at the section near the first divider increases and shifts towards the walls of branch with the increase of a. At the first bifurcation, the average horizontal velocities increase with the increase of α. The average horizontal velocities of outer branches at the secondary bifurcation decrease at 22.5°≤α≤ 30° and increase at 30° ≤α≤ 35°, whereas those of inner branches at the secondary bifurcation have the opposite variation, as the same as the above variations of the ratios of average horizontal velocities at various bifurcations. The ratios of average vertical velocities of branch at first bifurcation to that of outer branches at the secondary bifurcation increase at 22.5° ≤α≤ 30°and decrease at 30° ≤α≤ 35°, whereas the ratios of average vertical velocities of branch at first bifurcation to that of inner branches at the secondary bifurcation always decrease.展开更多
In the inviscid and incompressible fluid flow regime,surface tension effects on the behaviour of an initially spherical buoyancy-driven bubble rising in an infinite and initially stationary liquid are investigated num...In the inviscid and incompressible fluid flow regime,surface tension effects on the behaviour of an initially spherical buoyancy-driven bubble rising in an infinite and initially stationary liquid are investigated numerically by a volume of fluid （VOF） method. The ratio of the gas density to the liquid density is 0.001, which is close to the case of an air bubble rising in water. It is found by numerical experiment that there exist four critical Weber numbers We1,~We2,~We3 and We4, which distinguish five different kinds of bubble behaviours. It is also found that when 1≤We2, the bubble will finally reach a steady shape, and in this case after it rises acceleratedly for a moment, it will rise with an almost constant speed, and the lower the Weber number is, the higher the speed is. When We 〉We2, the bubble will not reach a steady shape, and in this case it will not rise with a constant speed. The mechanism of the above phenomena has been analysed theoretically and numerically.展开更多
基金supported by the National Natural Science Foundation of China(Nos.52106099 and 51576004)the Natural Science Foundation of Shandong Province(No.ZR2022YQ57)the Taishan Scholars Program.
文摘Phonon polaritons(PhPs)exhibit directional in-plane propagation and ultralow losses in van der Waals(vdW)crystals,offering new possibilities for controlling the flow of light at the nanoscale.However,these PhPs,including their directional propagation,are inherently determined by the anisotropic crystal structure of the host materials.Although in-plane anisotropic PhPs can be manipulated by twisting engineering,such as twisting individual vdW slabs,dynamically adjusting their propagation presents a significant challenge.The limited application of the twisted bilayer structure in bare films further restricts its usage.In this study,we present a technique in which anisotropic PhPs supported by bare biaxial vdW slabs can be actively tuned by modifying their local dielectric environment.Excitingly,we predict that the iso-frequency contour of PhPs can be reoriented to enable propagation along forbidden directions when the crystal is placed on a substrate with a moderate negative permittivity.Besides,we systematically investigate the impact of polaritonic coupling on near-field radiative heat transfer(NFRHT)between heterostructures integrated with different substrates that have negative permittivity.Our main findings reveal that through the analysis of dispersion contour and photon transmission coefficient,the excitation and reorientation of the fundamental mode facilitate increased photon tunneling,thereby enhancing heat transfer between heterostructures.Conversely,the annihilation of the fundamental mode hinders heat transfer.Furthermore,we find the enhancement or suppression of radiative energy transport depends on the relative magnitude of the slab thickness and the vacuum gap width.Finally,the effect of negative permittivity substrates on NFRHT along the[001]crystalline direction ofα-MoO3 is considered.The spectral band where the excited fundamental mode resulting from the negative permittivity substrates is shifted to the first Reststrahlen Band(RB 1)ofα-MoO_(3) and is widened,resulting in more significant enhancement of heat flux from RB 1.We anticipate our results will motivate new direction for dynamical tunability of the PhPs in photonic devices.
基金supported by the National Outstanding Young Scientist Foundation of China(11225213)the Key Subject “Computational solid mechanics” of China Academy of Engineering Physics
文摘The approximate compressible model is adopted to study the effects of strength and compressibility on the penetration by WHA long rod and copper jet into semi-infinite target in detail. For WHA rod penetrating PMMA at 2 km/s <V <5 km/s, the compressibility has a significant effect on the penetration efficiency. We clarify how compressibility affects the penetration efficiency by changing the stagnation pressures of the rod and target. For WHA rod penetrating 4340 Steel and 6061-T6 Al at 2 km/s < V < 10 km/s, the effect of strength is strong and the effect of compressibility is negligible at lower impact velocity, whilst the effect of strength is weak and the effect of compressibility becomes stronger at higher impact velocity. For the copper jet penetrating 4030 Steel, 6061-T6 Al and PMMA. the virtual origin model is adopted, and the compressibility and strength are implicitly considered by the linear relation between the penetration velocity and impact velocity. The effects of compressibility and target resistance on penetration efficiency are studied. The results show that the target resistance has a significant effect on the penetration efficiency. Howver PMMA is much more compressible than copper and the huge difference of compressibility has a significant effect on the penetration by hypervelocity copper jet into PMMA.
基金Supported by the National Natural Science Foundation of China under Grant No 90405005, the National Basic Research Program of China under Grant No 2007CB607506, the Specialized Research Fund for the Doctoral Programme of Higher Education of China under Grant No 20050730016, and the Fund of of Lanzhou University under Grant No WUT2005Z04.
文摘By solving the Boltzmann transport equation and considering the spin-dependent grain boundary scattering, the distribution of electrons in grains and the electrical transport properties in the applied magnetic field are studied. With regard to the dominant influence of grain boundary scattering which is taken as a boundary condition for the electrical transport, the grain size-dependent electrical conductivity is investigated. In addition, the reorientation of the relative magnetization between grains brings the change of the electron spin when the magnetonanocrystailine material is subjected to the magnetic field, resulting in the remarkable giant magnetoresistance effect.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11702034,11702218,and 11421062)Fundamental Research Funds for the Central Universities,China(Grant Nos.310812171011 and G2016KY0305)the National Key Project of Magneto-Constrained Fusion Energy Development Program,China(Grant No.2013GB110002)
文摘We study the stability of vortices pinning and dynamics in a superconducting thin strip containing a square array of antidot triplets by using the nonlinear Ginzburg–Landau(GL)theory.Compared with the regular square array of circular holes,the vortices are no longer pinned inside the circular holes,but instead stabilized at the center of the antidot triplets depending on the geometry parameters.Moreover,the influences of the geometry parameters and the polarity of the applied current on the current–voltage(I–V)characteristics are also studied.The critical current for the sample turning into a normal state becomes smaller when the hole diameter D is smaller and the spacing B between the holes is larger.Due to the asymmetric pinning sites,our numerical simulations demonstrate that the positive and negative rectified voltages appear alternately in the resistive state of the sample under an ac current of square pulses.
基金Project supported by the National Natural Science Foundation of China(Grant No.51076002)the National Basis Research Program of China(Grant No.2013CA328900)the Key Project of Complicated Electromagnetic Environment Laboratory of CAEP,China(Grant No.2015E0-01-1)
文摘Spectral and directional control of thermal emission based on excitation of confined electromagnetic resonant modes paves a viable way for the design and construction of microscale thermal emitters/absorbers. In this paper, we present numerical simulation results of the thermal radiative properties of a silicon carbide（Si C） thermal emitter/absorber composed of periodic microstructures. We illustrate different electromagnetic resonant modes which can be excited with the structure,such as surface phonon polaritons, magnetic polaritons and photonic crystal modes, and the process of radiation spectrum optimization based on a non-linear optimization algorithm. We show that the spectral and directional control of thermal emission/absorption can be efficiently achieved by adjusting the geometrical parameters of the structure. Moreover, the optimized spectrum is insensitive to 3% dimension modification.
文摘By using MTS815 rock mechanics test system,a series of acoustic emission(AE) location experiments were performed under unloading confining pressure,increasing the axial stress.The AE space-time evolution regularities and energy releasing characteristics during deformation and failure process of coal of different loading rates are compared,the influence mechanism of loading rates on the microscopic crack evolution were studied,combining the AE characteristics and the macroscopic failure modes of the specimens,and the precursory characteristics of coal failure were also analyzed quantitatively.The results indicate that as the loading rate is higher,the AE activity and the main fracture will begin earlier.The destruction of coal body is mainly the function of shear strain at lower loading rate and tension strain at higher rate,and will transform from brittleness to ductility at critical velocities.When the deformation of the coal is mainly plasticity,the amplitude of the AE ringing counting rate increases largely and the AE energy curves appear an obvious ''step'',which can be defined as the first failure precursor point.Statics of AE information shows that the strongest AE activity begins when the axial stress level was 92-98%,which can be defined as the other failure precursor point.As the loading rate is smaller,the coal more easily reaches the latter precursor point after the first one,so attention should be aroused to prevent dynamic disaster in coal mining when the AE activity reaches the first precursor point.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 10672043 and 10272032)
文摘A numerical method for simulating the motion and deformation of an axisymmetric bubble or drop rising or falling in another infinite and initially stationary fluid is developed based on the volume of fluid （VOF） method in the frame of two incompressible and immiscible viscous fluids under the action of gravity, taking into consideration of surface tension effects. A comparison of the numerical results by this method with those by other works indicates the validity of the method. In the frame of inviseid and incompressible fluids without taking into consideration of surface tension effects, the mechanisms of the generation of the liquid jet and the transition from spherical shape to toroidal shape during the bubble or drop deformation, the increase of the ring diameter of the toroidal bubble or drop and the decrease of its cross-section area during its motion, and the effects of the density ratio of the two fluids on the deformation of the bubble or drop are analysed both theoretically and numerically.
基金Supported by the 100 Persons Project of the Chinese Academy of Sciences, the National Natural Science Foundation of China under Grant No 10474109, the Foundation of Ministry of Personnel of China.
文摘The rolling massage manipulation is a classic Chinese massage, which is expected to improve the circulation by pushing, pulling and kneading of the muscle. A model for the rolling massage manipulation is proposed and the lattice Boltzmann method is applied to study the blood flow in the blood vessels. The simulation results show that the blood flux is considerably modified by the rolling massage and the explicit value depends on the rolling frequency, the rolling depth, and the diameter of the vessel. The smaller the diameter of the blood vessel, the larger the enhancement of the blood flux by the rolling massage. The model, together with the simulation results, is expected to be helpful to understand the mechanism and further development of rolling massage techniques.
基金This research was supported by the National Natural Science Foundation of China(Grant No.11472115,11872214)the China Scholarship Council(201706845026).
文摘Liquid-filled compartment structure consists of a bulk steel plate with matrix blind holes which are filled with liquid and a steel front plate to seal up the liquid with rings and bolts.The liquid-filled compartment structure can resist the shaped charge warhead effectively.This paper presents experimental and theoretical investigations of the penetration ability of the residual shaped charge jet emerging from the liquid-filled compartment structure after the penetration process at different impact angles.On the basis of shock wave propagation theory,the influence of the liquid-filled compartment structure on jet stability is analysed.The interferences of the liquid backflow caused by a reflected shock wave and a back plate on jet stability under different impact angles are also examined.In addition,the range of the disturbed velocity segments of the jet at different impact angles and the penetration ability of the residual jet are obtained.A theoretical model is validated against the experimental penetration depths.
文摘An experimental study on acceleration mechanism of flame propagation of propane-air mixture in ducts with obstacles was conducted. The acceleration mechanism of flame propagation is mainly due to the positive feedback of the turbulence region induced by obstacles for combustion process. It can be seen from the experimental results that the maximum explosion pressure can increase by 20%, the maximum rate of pressure rise can increase by 10 times and the flame propagation velocity can increase by 20 times when obstacles are present.
基金Project supported by the National Natural Science Foundation of China and NSAF (Grant No 10576010)the Creation Foundation of Fudan University (Grant No 2126003)
文摘This paper studies the roughness effect combining with effects of rarefaction and compressibility by a lattice Boltzmann model for rarefied gas flows at high Knudsen numbers. By discussing the effect of the tangential momentum accommodation coefficient on the rough boundary condition, the lattice Boltzmann simulations of nitrogen and helium flows are performed in a two-dimensional microchannel with rough boundaries. The surface roughness effects in the microchannel on the velocity field, the mass flow rate and the friction coefficient are studied and analysed. Numerical results for the two gases in micro scale show different characteristics from macroscopic flows and demonstrate the feasibility of the lattice Boltzmann model in rarefied gas dynamics.
基金supported by the National Natural Science Foundation of China(Grant Nos.10732010,10972010,11332002,and 11028206)
文摘Molecular dynamics simulations are carried out to study atomic diffusion in the explosive welding process of NisoTis0-Cu （at.%）. By using a hybrid method which combines molecular dynamics simulation and classical diffusion the- ory, the thickness of the diffusion layer and the atomic concentration distribution across the welding interface are obtained. The results indicate that the concentration distribution curves at different times have a geometric similarity. According to the geometric similarity, the atomic concentration distribution at any time in explosive welding can be calculated. NisoTis0- Cu explosive welding and scanning electron microscope experiments are done to verify the results. The simulation results and the experimental results are in good agreement.
文摘The particle path tracking method is proposed and used in two-dimensional（2D） and three-dimensional（3D） numerical simulations of continuously rotating detonation engines（CRDEs）. This method is used to analyze the combustion and expansion processes of the fresh particles, and the thermodynamic cycle process of CRDE. In a 3D CRDE flow field, as the radius of the annulus increases, the no-injection area proportion increases, the non-detonation proportion decreases, and the detonation height decreases. The flow field parameters on the 3D mid annulus are different from in the 2D flow field under the same chamber size. The non-detonation proportion in the 3D flow field is less than in the 2D flow field. In the 2D and 3D CRDE, the paths of the flow particles have only a small fluctuation in the circumferential direction. The numerical thermodynamic cycle processes are qualitatively consistent with the three ideal cycle models, and they are right in between the ideal F–J cycle and ideal ZND cycle. The net mechanical work and thermal efficiency are slightly smaller in the 2D simulation than in the 3D simulation. In the 3D CRDE, as the radius of the annulus increases, the net mechanical work is almost constant, and the thermal efficiency increases. The numerical thermal efficiencies are larger than F–J cycle, and much smaller than ZND cycle.
基金National Key R&D Program of China(No.2018YFC0809700,No.2017YFC0803300)National Natural Science Foundation of China(No.71673158,No.11702046).
文摘Particulate composites are one of the widely used materials in producing numerous state-of-the-art components in biomedical,automobile,aerospace including defence technology.Variety of modelling techniques have been adopted in the past to model mechanical behaviour of particulate composites.Due to their favourable properties,particle-based methods provide a convenient platform to model failure or fracture of these composites.Smooth particle hydrodynamics(SPH)is one of such methods which demonstrate excellent potential for modelling failure or fracture of particulate composites in a Lagrangian setting.One of the major challenges in using SPH method for modelling composite materials depends on accurate and efficient way to treat interface and boundary conditions.In this paper,a masterslave method based multi-freedom constraints is proposed to impose essential boundary conditions and interfacial displacement constraints in modelling mechanical behaviour of composite materials using SPH method.The proposed methodology enforces the above constraints more accurately and requires only smaller condition number for system stiffness matrix than the procedures based on typical penalty function approach.A minimum cut-off value-based error criteria is employed to improve the computational efficiency of the proposed methodology.In addition,the proposed method is further enhanced by adopting a modified numerical interpolation scheme along the boundary to increase the accuracy and computational efficiency.The numerical examples demonstrate that the proposed master-slave approach yields better accuracy in enforcing displacement constraints and requires approximately the same computational time as that of penalty method.
基金Supported by the National Natural Science Foundation of China under Grant Nos 10032020 and 10225210.
文摘The relative scaling exponents and intermittency of three-dimensional compressible turbulent channel flow are investigated by using direct numerical simulation. One case is subsonic flow （Mα=0.8）, the other is supersonic （Mα=1.3）, and the Reynolds numbers based on the mean bulk velocity and channel half-width are 2826 and 3010, respectively. The analysis of the local slopes of sixth order velocity structure function to third order reveals that there is a well-defined scaling range for 10〈y^＋〈100. it is also noted that the intermittency of longitudinal velocity increments in this region is stronger than that of the transverse ones. Comparison with the incompressible case shows that the location of the most intensive intermittency moves toward the log-law region, which is related to the displacement of streamwise vortical structures in the near-wall region.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11702120,11372120,11421062,and 11572143)the Fundamental Research Funds for the Central Universities,China(Grant No.lzujbky-2016-106)
文摘In this paper, a novel magnetoelectric（ME） composite structure is proposed, and the ME response in the structure is measured at the bias magnetic field up to 2000 Oe（1 Oe = 79.5775 A·m^（-1）） and the excitation frequency of alternating magnetic field ranging from 1 kHz to 200 kHz. The ME voltage of each PZT layer is detected. According to the measurement results, the phase differences are observed among three channels and the multi-peak phenomenon appears in each channel. Meanwhile, the results show that the ME structure can stay a relatively high ME response within a wide bandwidth.Besides, the hysteretic loops of three PZT layers are observed. When the frequency of alternating current（AC） magnetic field changes, the maximum value of ME coefficient appears in different layers due to the multiple vibration modes of the structure. Moreover, a finite element analysis is performed to evaluate the resonant frequency of the structure, and the theoretical calculating results accord well with the experimental results. The experiment results suggest that the proposed structure may be a good candidate for designing broadband magnetic field sensors.
文摘The problem involving an edge-crack in a rectangular material under the anti-plane mechanical loading and in-plane electric loading is analyzed under the impermeable conditions. By using the series expansion, the general solutions of electromechanical fields are obtained, which satisfied both governing equations and crack sufrace boundary conditions, and the unknown constants in which can be obtained by the boundary collocation method. Numerical results are given to show the effect of electromechanical interaction on energy release rate.
基金Supported by Natural Science Foundation of China(Nos.10772049 and 11072062)Research and Development Program of China(863 Program, No.2009AA04Z408)+1 种基金Natural Science Foundation of Shanghai(No.06ZR14009)Pujiang Scholar Program and the Wangdao Scholar Program(No.08076) of Fudan University
文摘A new method for three-dimensional simulation of the interaction between the gas and the solid around is developed.The effects of the gas on the thermal-mechanical behaviors within the surrounded solid are performed by replacing the internal gas with an equivalent solid in the modeling,which can make it convenient to simulate the thermal-mechanical coupling effects in the solid research objects with gases in them.The applied thermal expansion coefficient,Young's modulus and Poisson's ratio of the equivalent solid material are derived.A series of tests have been conducted;and the proposed equivalent solid method to simulate the gas effects is validated.
基金Supported by the National Natural Science Foundation of China under Grant No 10576010.
文摘The multi-bifurcation effect of blood flow is investigated by lattice Boltzmann method at Re = 200 with six different bifurcation angles α which are 22.5°, 25°, 28°, 30°, 33°, 35°, respectively. The velocities and ratios of average velocity at various bifurcations are discussed. It is indicated that the maximum velocity at the section near the first divider increases and shifts towards the walls of branch with the increase of a. At the first bifurcation, the average horizontal velocities increase with the increase of α. The average horizontal velocities of outer branches at the secondary bifurcation decrease at 22.5°≤α≤ 30° and increase at 30° ≤α≤ 35°, whereas those of inner branches at the secondary bifurcation have the opposite variation, as the same as the above variations of the ratios of average horizontal velocities at various bifurcations. The ratios of average vertical velocities of branch at first bifurcation to that of outer branches at the secondary bifurcation increase at 22.5° ≤α≤ 30°and decrease at 30° ≤α≤ 35°, whereas the ratios of average vertical velocities of branch at first bifurcation to that of inner branches at the secondary bifurcation always decrease.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10672043 and 10272032)
文摘In the inviscid and incompressible fluid flow regime,surface tension effects on the behaviour of an initially spherical buoyancy-driven bubble rising in an infinite and initially stationary liquid are investigated numerically by a volume of fluid （VOF） method. The ratio of the gas density to the liquid density is 0.001, which is close to the case of an air bubble rising in water. It is found by numerical experiment that there exist four critical Weber numbers We1,~We2,~We3 and We4, which distinguish five different kinds of bubble behaviours. It is also found that when 1≤We2, the bubble will finally reach a steady shape, and in this case after it rises acceleratedly for a moment, it will rise with an almost constant speed, and the lower the Weber number is, the higher the speed is. When We 〉We2, the bubble will not reach a steady shape, and in this case it will not rise with a constant speed. The mechanism of the above phenomena has been analysed theoretically and numerically.
