Based on non-Darcian flow law described by exponent and threshold gradient within a double-layered soil, the classic theory of one-dimensional consolidation of double-layered soil was modified to consider the change o...Based on non-Darcian flow law described by exponent and threshold gradient within a double-layered soil, the classic theory of one-dimensional consolidation of double-layered soil was modified to consider the change of vertical total stress with depth and time together. Because of the complexity of governing equations, the numerical solutions were obtained in detail by finite difference method. Then, the numerical solutions were compared with the analytical solutions in condition that non-Darcian flow law was degenerated to Dary's law, and the comparison results show that numerical solutions are reliable. Finally, consolidation behavior of double-layered soil with different parameters was analyzed, and the results show that the consolidation rate of double-layered soil decreases with increasing the value of exponent and threshold of non-Darcian flow, and the exponent and threshold gradient of the first soil layer greatly influence the consolidation rate of double-layered soil. The larger the ratio of the equivalent water head of external load to the total thickness of double-layered soil, the larger the rate of the consolidation, and the similitude relationship in classical consolidation theory of double-layered soil is not satisfied. The other consolidation behavior of double-layered soil with non-Darcian flow is the same as that with Darcy's law.展开更多
The influence of initial groove angle on strain rate inside and outside groove of Ti6Al4V alloy was investigated.Based on the evolution of strain rate inside and outside groove,the effect of strain rate difference on ...The influence of initial groove angle on strain rate inside and outside groove of Ti6Al4V alloy was investigated.Based on the evolution of strain rate inside and outside groove,the effect of strain rate difference on the evolution of normal stress and effective stress inside and outside groove was also analyzed.The results show that when linear loading path changes from uniaxial tension to equi-biaxial tension,the initial groove angle plays a weaker role in the evolution of strain rate in the M-K model.Due to the constraint of force equilibrium between inside and outside groove,the strain rate difference makes the normal stress inside groove firstly decrease and then increase during calculation,which makes the prediction algorithm of forming limit convergent at elevated temperature.The decrease of normal stress inside groove is mainly caused by high temperature softening effect and the rotation of groove,while the increase of normal stress inside groove is mainly due to strain rate hardening effect.展开更多
The anisotropic properties of subsurface media cause waveform distortions in seismic wave propagation,resulting in a negative infl uence on seismic imaging.In addition,wavefields simulated by the conventional coupled ...The anisotropic properties of subsurface media cause waveform distortions in seismic wave propagation,resulting in a negative infl uence on seismic imaging.In addition,wavefields simulated by the conventional coupled pseudo-acoustic equation are not only aff ected by SV-wave artifacts but are also limited by anisotropic parameters.We propose a least-squares reverse time migration(LSRTM)method based on the pure q P-wave equation in vertically transverse isotropic media.A fi nite diff erence and fast Fourier transform method,which can improve the effi ciency of the numerical simulation compared to a pseudo-spectral method,is used to solve the pure q P-wave equation.We derive the corresponding demigration operator,migration operator,and gradient updating formula to implement the LSRTM.Numerical tests on the Hess model and field data confirm that the proposed method has a good correction eff ect for the travel time deviation caused by underground anisotropic media.Further,it signifi cantly suppresses the migration noise,balances the imaging amplitude,and improves the imaging resolution.展开更多
Supercavitating flow around a slender symmetric wedge moving at variable velocity in static fluid has been studied. Singular integral equation for the flow has been founded through distributing the sources and sinks o...Supercavitating flow around a slender symmetric wedge moving at variable velocity in static fluid has been studied. Singular integral equation for the flow has been founded through distributing the sources and sinks on the symmetrical axis. The supereavity length at each moment is determined by solving the singular integral equation with finite difference method. The supercavity shape at each moment is obtained by solving the partial differential equation with variable coefficient. For the case that the wedge takes the impulse and uniformly variable motion, numerical results of time history of the supercavity length and shape are presented. The calculated results indicate that the shape and the length of the supercavity vary in a similar way to the case that the wedge takes variable motion, and there is a time lag in unsteady supercavitating flow induced by the variation of wedge velocity.展开更多
Gas-particle two-phase flow is a very important consideration in designing various machines. Although a great deal of theoretical, experimental, and numerical research has been carried out, particle motion in a supers...Gas-particle two-phase flow is a very important consideration in designing various machines. Although a great deal of theoretical, experimental, and numerical research has been carried out, particle motion in a supersonic flow has not been sufficiently clarified. Hence, in order to clarify the interactions between flow and particles, the authors consider the characteristics of particle motion, especially at high temperatures. In the present study, the flow of a gas with a diluted particle load is to be simulated in a conventional converging-diverging supersonic nozzle. The turbulent gas flow in the nozzle is computed with the finite difference and RANS (raynolds averaged navier-stokes simulation) methods. The particle motion is simulated in a Lagrangian manner. In addition, taking into account the light particle loading, a weak coupling method is used. Through this investigation, it is shown that the particle velocity increases monotonically from the nozzle throat to the outlet. And it is shown that particles can be accelerated to higher velocities in helium than in nitrogen, and smaller particles tend to attain higher speed and lower static temperature.展开更多
The Effects of pressure stress work and viscous dissipation in mixed convection flow along a vertical fiat plate have been investigated. The results are obtained by transforming the governing system of boundary layer ...The Effects of pressure stress work and viscous dissipation in mixed convection flow along a vertical fiat plate have been investigated. The results are obtained by transforming the governing system of boundary layer equations into a system of non-dimensional equations and by applying implicit finite difference method together with Newton's linearization approximation. Numerical results for different values of pressure stress work parameter, viscous dissipation parameter and Prandtl number have been obtained. The velocity profiles, temperature distributions, skin friction co-efficient and the rate of heat transfer have been presented graphically for the effects of the aforementioned parameters.展开更多
Using the hybrid finite difference method, we solve the Fokker-Planck equation to study the effect of seed electron injection on acceleration of radiation belt electrons driven by chorus waves. Numerical results show ...Using the hybrid finite difference method, we solve the Fokker-Planck equation to study the effect of seed electron injection on acceleration of radiation belt electrons driven by chorus waves. Numerical results show that in the absence of injection chorus waves can accelerate electrons at large pitch angles (ae〉60°), producing enhancements in the phase space density (PSD) of (1-2 MeV ) electrons by a factor of 100-1000 within 1-2 days. In the presence of injection, chorus waves yield increase in PSD of electrons by accelerating the injected seed electrons. Meanwhile, the PSD evolution increases as the pitch angle in- creases but decreases as electron energy increases. Moreover, the PSD evolution can extend to higher energies with a time scale of 1-2 days for 1-2 MeV energies. When the injection increases by a factor of 10 higher than the initial value and re- mains for about two days, maximum values of PSD for 1 or 2 MeV increase to 6 or 3 times respectively higher than those without injection in two days. The current results suggest that the injected seed electrons play an important role in the evolu- tion of the radiation belt electrons.展开更多
In this study a mathematical model for two-dimensional pulsatile blood flow through overlapping constricted tapered vessels is presented. In order to establish resemblance to the in vivo conditions, an improved shape ...In this study a mathematical model for two-dimensional pulsatile blood flow through overlapping constricted tapered vessels is presented. In order to establish resemblance to the in vivo conditions, an improved shape of the time-variant overlapping stenosis in the elastic tapered artery subject to pulsatile pressure gradient is considered. Because it contains a suspension of all erythrocytes, the flowing blood is represented by micropolar fluid. By applying a suitable coordinate transformation, tapered cosine-shaped artery turned into non-tapered rectangular and a rigid artery. The governing nonlinear partial differential equations under the imposed realistic boundary conditions are solved using the finite difference method. The effects of vessel tapering on flow characteristics consid- ering their dependencies with time are investigated. The results show that by increasing the taper angle the axial velocity and volumetric flow rate increase and the microrota- tional velocity and resistive impedance reduce. It has been shown that the results are in agreement with similar data from the literature.展开更多
基金Projects(50878191,51109092)supported by the National Natural Science Foundation of China
文摘Based on non-Darcian flow law described by exponent and threshold gradient within a double-layered soil, the classic theory of one-dimensional consolidation of double-layered soil was modified to consider the change of vertical total stress with depth and time together. Because of the complexity of governing equations, the numerical solutions were obtained in detail by finite difference method. Then, the numerical solutions were compared with the analytical solutions in condition that non-Darcian flow law was degenerated to Dary's law, and the comparison results show that numerical solutions are reliable. Finally, consolidation behavior of double-layered soil with different parameters was analyzed, and the results show that the consolidation rate of double-layered soil decreases with increasing the value of exponent and threshold of non-Darcian flow, and the exponent and threshold gradient of the first soil layer greatly influence the consolidation rate of double-layered soil. The larger the ratio of the equivalent water head of external load to the total thickness of double-layered soil, the larger the rate of the consolidation, and the similitude relationship in classical consolidation theory of double-layered soil is not satisfied. The other consolidation behavior of double-layered soil with non-Darcian flow is the same as that with Darcy's law.
基金Project(51775023)supported by the National Natural Science Foundation of ChinaProject(YWF-18-BJ-J-75)supported by the Fundamental Research Funds for the Central Universities,China
文摘The influence of initial groove angle on strain rate inside and outside groove of Ti6Al4V alloy was investigated.Based on the evolution of strain rate inside and outside groove,the effect of strain rate difference on the evolution of normal stress and effective stress inside and outside groove was also analyzed.The results show that when linear loading path changes from uniaxial tension to equi-biaxial tension,the initial groove angle plays a weaker role in the evolution of strain rate in the M-K model.Due to the constraint of force equilibrium between inside and outside groove,the strain rate difference makes the normal stress inside groove firstly decrease and then increase during calculation,which makes the prediction algorithm of forming limit convergent at elevated temperature.The decrease of normal stress inside groove is mainly caused by high temperature softening effect and the rotation of groove,while the increase of normal stress inside groove is mainly due to strain rate hardening effect.
基金financially supported by the National Key R&D Program of China (No. 2019YFC0605503)the Major Scientific and Technological Projects of CNPC (No. ZD2019-183-003)the National Natural Science Foundation of China (No. 41922028,41874149)。
文摘The anisotropic properties of subsurface media cause waveform distortions in seismic wave propagation,resulting in a negative infl uence on seismic imaging.In addition,wavefields simulated by the conventional coupled pseudo-acoustic equation are not only aff ected by SV-wave artifacts but are also limited by anisotropic parameters.We propose a least-squares reverse time migration(LSRTM)method based on the pure q P-wave equation in vertically transverse isotropic media.A fi nite diff erence and fast Fourier transform method,which can improve the effi ciency of the numerical simulation compared to a pseudo-spectral method,is used to solve the pure q P-wave equation.We derive the corresponding demigration operator,migration operator,and gradient updating formula to implement the LSRTM.Numerical tests on the Hess model and field data confirm that the proposed method has a good correction eff ect for the travel time deviation caused by underground anisotropic media.Further,it signifi cantly suppresses the migration noise,balances the imaging amplitude,and improves the imaging resolution.
基金Sponsored by the National Natural Science Foundation of China(Grant No.10832007)
文摘Supercavitating flow around a slender symmetric wedge moving at variable velocity in static fluid has been studied. Singular integral equation for the flow has been founded through distributing the sources and sinks on the symmetrical axis. The supereavity length at each moment is determined by solving the singular integral equation with finite difference method. The supercavity shape at each moment is obtained by solving the partial differential equation with variable coefficient. For the case that the wedge takes the impulse and uniformly variable motion, numerical results of time history of the supercavity length and shape are presented. The calculated results indicate that the shape and the length of the supercavity vary in a similar way to the case that the wedge takes variable motion, and there is a time lag in unsteady supercavitating flow induced by the variation of wedge velocity.
文摘Gas-particle two-phase flow is a very important consideration in designing various machines. Although a great deal of theoretical, experimental, and numerical research has been carried out, particle motion in a supersonic flow has not been sufficiently clarified. Hence, in order to clarify the interactions between flow and particles, the authors consider the characteristics of particle motion, especially at high temperatures. In the present study, the flow of a gas with a diluted particle load is to be simulated in a conventional converging-diverging supersonic nozzle. The turbulent gas flow in the nozzle is computed with the finite difference and RANS (raynolds averaged navier-stokes simulation) methods. The particle motion is simulated in a Lagrangian manner. In addition, taking into account the light particle loading, a weak coupling method is used. Through this investigation, it is shown that the particle velocity increases monotonically from the nozzle throat to the outlet. And it is shown that particles can be accelerated to higher velocities in helium than in nitrogen, and smaller particles tend to attain higher speed and lower static temperature.
文摘The Effects of pressure stress work and viscous dissipation in mixed convection flow along a vertical fiat plate have been investigated. The results are obtained by transforming the governing system of boundary layer equations into a system of non-dimensional equations and by applying implicit finite difference method together with Newton's linearization approximation. Numerical results for different values of pressure stress work parameter, viscous dissipation parameter and Prandtl number have been obtained. The velocity profiles, temperature distributions, skin friction co-efficient and the rate of heat transfer have been presented graphically for the effects of the aforementioned parameters.
基金supported by the Initiative Project of Chinese Academy of Sciences (Grant No. YYYJ-1110)
文摘Using the hybrid finite difference method, we solve the Fokker-Planck equation to study the effect of seed electron injection on acceleration of radiation belt electrons driven by chorus waves. Numerical results show that in the absence of injection chorus waves can accelerate electrons at large pitch angles (ae〉60°), producing enhancements in the phase space density (PSD) of (1-2 MeV ) electrons by a factor of 100-1000 within 1-2 days. In the presence of injection, chorus waves yield increase in PSD of electrons by accelerating the injected seed electrons. Meanwhile, the PSD evolution increases as the pitch angle in- creases but decreases as electron energy increases. Moreover, the PSD evolution can extend to higher energies with a time scale of 1-2 days for 1-2 MeV energies. When the injection increases by a factor of 10 higher than the initial value and re- mains for about two days, maximum values of PSD for 1 or 2 MeV increase to 6 or 3 times respectively higher than those without injection in two days. The current results suggest that the injected seed electrons play an important role in the evolu- tion of the radiation belt electrons.
文摘In this study a mathematical model for two-dimensional pulsatile blood flow through overlapping constricted tapered vessels is presented. In order to establish resemblance to the in vivo conditions, an improved shape of the time-variant overlapping stenosis in the elastic tapered artery subject to pulsatile pressure gradient is considered. Because it contains a suspension of all erythrocytes, the flowing blood is represented by micropolar fluid. By applying a suitable coordinate transformation, tapered cosine-shaped artery turned into non-tapered rectangular and a rigid artery. The governing nonlinear partial differential equations under the imposed realistic boundary conditions are solved using the finite difference method. The effects of vessel tapering on flow characteristics consid- ering their dependencies with time are investigated. The results show that by increasing the taper angle the axial velocity and volumetric flow rate increase and the microrota- tional velocity and resistive impedance reduce. It has been shown that the results are in agreement with similar data from the literature.