期刊文献+
共找到1,652篇文章
< 1 2 83 >
每页显示 20 50 100
A High-Order Conservative Semi-Lagrangian Solver for 3D Free Surface Flows with Sediment Transport on Voronoi Meshes
1
作者 Matteo Bergami Walter Boscheri Giacomo Dimarco 《Communications on Applied Mathematics and Computation》 2023年第2期596-637,共42页
In this paper,we present a conservative semi-Lagrangian scheme designed for the numeri-cal solution of 3D hydrostatic free surface flows involving sediment transport on unstruc-tured Voronoi meshes.A high-order recons... In this paper,we present a conservative semi-Lagrangian scheme designed for the numeri-cal solution of 3D hydrostatic free surface flows involving sediment transport on unstruc-tured Voronoi meshes.A high-order reconstruction procedure is employed for obtaining a piecewise polynomial representation of the velocity field and sediment concentration within each control volume.This is subsequently exploited for the numerical integration of the Lagrangian trajectories needed for the discretization of the nonlinear convective and viscous terms.The presented method is fully conservative by construction,since the transported quantity or the vector field is integrated for each cell over the deformed vol-ume obtained at the foot of the characteristics that arises from all the vertexes defining the computational element.The semi-Lagrangian approach allows the numerical scheme to be unconditionally stable for what concerns the advection part of the governing equations.Furthermore,a semi-implicit discretization permits to relax the time step restriction due to the acoustic impedance,hence yielding a stability condition which depends only on the explicit discretization of the viscous terms.A decoupled approach is then employed for the hydrostatic fluid solver and the transport of suspended sediment,which is assumed to be passive.The accuracy and the robustness of the resulting conservative semi-Lagrangian scheme are assessed through a suite of test cases and compared against the analytical solu-tion whenever is known.The new numerical scheme can reach up to fourth order of accu-racy on general orthogonal meshes composed by Voronoi polygons. 展开更多
关键词 Conservative semi-Lagrangian free surface flows Sediment transport High-order reconstruction Hydrostatic model
下载PDF
Validation of a CIP-based tank for numerical simulation of free surface flows 被引量:2
2
作者 Xi-Zeng Zhao 《Acta Mechanica Sinica》 SCIE EI CAS CSCD 2011年第6期877-890,共14页
A constrained interpolation profile CIP-based numerical tank is developed to simulate violent free surface flows.The numerical simulation is performed by the CIP-based Cartesian grid method,which is described in the p... A constrained interpolation profile CIP-based numerical tank is developed to simulate violent free surface flows.The numerical simulation is performed by the CIP-based Cartesian grid method,which is described in the present paper.The tangent of hyperbola for interface capturing(THINC) scheme is applied for capturing complex free surfaces.The new model is capable of simulating a flow with violently varied free surface.A series of computations are conducted to assess the developed algorithm and its versatility.These tests include the collapse of water column with and without an obstacle,sloshing in a fixed tank,the generation of regular waves in a tank,the generation of extreme waves in a tank.Excellent agreements are obtained when numerical results are compared with available analytical,experimental,and other numerical results. 展开更多
关键词 CIP method. THINC scheme free surface flow Wave breaking Wave tank
下载PDF
A Higher-Efficient Non-Hydrostatic Finite Volume Model for Strong Three-Dimensional Free Surface Flows and Sediment Transport 被引量:1
3
作者 LIU Xin MA Dian-guang ZHANG Qing-he 《China Ocean Engineering》 SCIE EI CSCD 2017年第6期736-746,共11页
In order to accurately simulate strong three-dimensional (3-D) free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equati... In order to accurately simulate strong three-dimensional (3-D) free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equations and convection-diffusion equation of sediment concentration with the mixing triangle and quadrilateral grids. The governing equations are discretized with the unstructured finite volume method in order to provide conservation properties of mass and momentum, and flexibility with practical application. It is shown that it is first-order accurate on nonuniform plane two-dimensional (2-D) grids and second-order accurate on uniform plane grids. A third-order approximation of the vertical velocity at the top-layer is applied. In such a way, free surface zero stress boundary condition is satisfied maturely, and very few vertical layers are needed to give an accurate solution even for complex discontinuous flow and short wave simulation. The model is applied to four examples to simulate strong 3-D free surface flows and sediment transport where non-hydrostatic pressures have a considerable effect on the velocity field. The newly developed model is verified against analytical solutions with an excellent agreement. 展开更多
关键词 higher-efficient NON-HYDROSTATIC strong 3-D free surface flows sediment transport 3-D numerical model
下载PDF
Effect of Viscosity on Free-Surface Waves in Oseen Flows 被引量:2
4
作者 卢东强 《Journal of Shanghai University(English Edition)》 CAS 2004年第1期8-12,共5页
Based on the complex dispersion relation for the two-dimensional free-surface waves generated by a moving body in the steady Oseen flows, the effect of viscosity on wavelength and wave amplitude was investigated by me... Based on the complex dispersion relation for the two-dimensional free-surface waves generated by a moving body in the steady Oseen flows, the effect of viscosity on wavelength and wave amplitude was investigated by means of an asymptotic method and a numerical analysis. A comparison between the asymptotic and numerical analysis for the viscous decay factor demonstrates the validity of the perturbation expansions for the wave profile. The numerical result shows that the wavelength of viscous wave is slightly elongated in comparison with that of inviscid wave. 展开更多
关键词 free-surface wave Oseen flow VISCOSITY complex dispersion relation.
下载PDF
Simulation of Free Surface Flows Based on MapReduce in Cloud Computing
5
作者 刘旭 班晓娟 王墨涵 《China Communications》 SCIE CSCD 2011年第6期28-35,共8页
Simulating turbulent liquids with breaking waves and splashes is among the most desired features in fluid animation. Lagrangian methods such as Smoothed Particle Hydrodynamics method (SPH) are a promising way to captu... Simulating turbulent liquids with breaking waves and splashes is among the most desired features in fluid animation. Lagrangian methods such as Smoothed Particle Hydrodynamics method (SPH) are a promising way to capture such properties. However, the Particle-based liquid surface simulation has not been applied very well since its consumption is way too large. This paper derives the governing equations in SPH approaches and parallelizes the dynamics-based surface simulation with the MapReduce program models which apply the SPH approach in Cloud Computing. Compared to the serial methods, this approach obtained a 3.11 times speedup on the experimental platform. 展开更多
关键词 MAPREDUCE SPH cloud computing free surface flows parallelize
下载PDF
Modeling of multiphase flow in low permeability porous media:Effect of wettability and pore structure properties
6
作者 Xiangjie Qin Yuxuan Xia +3 位作者 Juncheng Qiao Jiaheng Chen Jianhui Zeng Jianchao Cai 《Journal of Rock Mechanics and Geotechnical Engineering》 SCIE CSCD 2024年第4期1127-1139,共13页
Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the ef... Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the effects of complex pore structures and wettability.To address this issue,based on the digital rock of low permeability sandstone,a direct numerical simulation is performed considering the interphase drag and boundary slip to clarify the microscopic water-oil displacement process.In addition,a dual-porosity pore network model(PNM)is constructed to obtain the water-oil relative permeability of the sample.The displacement efficiency as a recovery process is assessed under different wetting and pore structure properties.Results show that microscopic displacement mechanisms explain the corresponding macroscopic relative permeability.The injected water breaks through the outlet earlier with a large mass flow,while thick oil films exist in rough hydrophobic surfaces and poorly connected pores.The variation of water-oil relative permeability is significant,and residual oil saturation is high in the oil-wet system.The flooding is extensive,and the residual oil is trapped in complex pore networks for hydrophilic pore surfaces;thus,water relative permeability is lower in the water-wet system.While the displacement efficiency is the worst in mixed-wetting systems for poor water connectivity.Microporosity negatively correlates with invading oil volume fraction due to strong capillary resistance,and a large microporosity corresponds to low residual oil saturation.This work provides insights into the water-oil flow from different modeling perspectives and helps to optimize the development plan for enhanced recovery. 展开更多
关键词 Low permeability porous media Water-oil flow WETTABILITY Pore structures Dual porosity pore network model(PNM) free surface model
下载PDF
A Higher-Efficient Three-Dimensional Numerical Model for Small Amplitude Free Surface Flows
7
作者 吕彪 《China Ocean Engineering》 SCIE EI CSCD 2014年第5期617-628,共12页
A higher-efficient three-dimensional non-hydrostatic model is developed to simulate small amplitude free surface flows based on a staggered unstructured grid. In this model, a fractional step algorithm is adopted to s... A higher-efficient three-dimensional non-hydrostatic model is developed to simulate small amplitude free surface flows based on a staggered unstructured grid. In this model, a fractional step algorithm is adopted to solve the Navier-Stokes equations in two major steps. A top-layer pressure method is proposed to minimize the number of vertical layers and subsequently the computational cost. Three classical examples of small amplitude free surface flows are used to demonstrate the capability and efficiency of the model. The satisfactory results demonstrated the capability and efficiency of modelling a range of small amplitude free surface flows with only a small number of vertical layers. 展开更多
关键词 higher-efficient NON-HYDROSTATIC small amplitude free surface flows 3D numerical model
下载PDF
A New Efficient Finite Volume Modeling of Small Amplitude Free Surface Flows with Unstructured Grid
8
作者 吕彪 《China Ocean Engineering》 SCIE EI CSCD 2013年第4期509-522,共14页
A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing t... A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing time. The advection and horizontal diffusion terms in the momentum equation are discretized by an integral interpolation method on the orthogonal unstructured staggered mesh and, while it has the attractive property of being conservative. The pressure-correction algorithm is employed for the non-hydrostatic pressure in order to achieve second-order temporal accuracy. A conservative scalar transport algorithm is also applied to discretize k - c equations in this model. The eddy viscosity is calculated from the k-c turbulent model. The resulting model is mass and momentum conservative. The model is verified by two examples to simulate unsteady small amplitude free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field, and then applied to simulate the tidal flow in the Bohai Sea. 展开更多
关键词 orthogonal unstructured grid NON-HYDROSTATIC small amplitude free surface flows 3-D numerical model k - ~ turbulent model
下载PDF
Effect of Internal Flow on Vortex-Induced Vibration of Submarine Free Spanning Pipelines 被引量:10
9
作者 娄敏 丁坚 +1 位作者 郭海燕 董晓林 《China Ocean Engineering》 SCIE EI 2005年第1期147-154,共8页
At present, most researches on the vortex-induced vibration of submarine free spanning pipelines ignore the effect of internal flowing fluid; furthermore, there are no research reports considering the coupling effect ... At present, most researches on the vortex-induced vibration of submarine free spanning pipelines ignore the effect of internal flowing fluid; furthermore, there are no research reports considering the coupling effect of internal and external fluid with the free span. In this paper, combining Iwan's wake oscillator model with the differential equation derived for the dynamic response of submarine free spanning pipelines with inclusion of internal flow, the pipe-fluid coupling equations are developed to investigate the effect of internal flow on the vortex-induced vibration of the free spans. The finite element approximation is implemented to derive the matrix equations of equilibrium. The Newmark method combined with simple iteration is used to solve the system of equations. The results indicate that the internal fluid flow may cause the shift of resonance band to the lower frequency and a slight decrease in the peak value; the effect will be more pronounced with the increase of the span length and can be weakened in the presence of the axial tension. 展开更多
关键词 vortex-induced vibration internal flow free span submarine pipelines
下载PDF
Vertical 2D Modeling of Free Surface Flow with Hydrodynamic Pressure Using SIMPLE Arithmetic in σ Coordinates 被引量:3
10
作者 吴修广 沈永明 郑永红 《海洋工程:英文版》 2004年第1期79-92,共14页
A numerical model for shallow water flow has been developed based on the unsteady Reynolds-averaged Navier-Stokes equations with the hydrodynamic pressure instead of hydrostatic pressure assumption. The equations are ... A numerical model for shallow water flow has been developed based on the unsteady Reynolds-averaged Navier-Stokes equations with the hydrodynamic pressure instead of hydrostatic pressure assumption. The equations are transformed into the σ-coordinate system and the eddy viscosity is calculated with the standard k-ε turbulence model. The control volume method is used to discrete the equations, and the boundary conditions at the bed for shallow water models only include vertical diffusion terms expressed with wall functions. And the semi-implicit method for pressure linked equation arithmetic is adopted to solve the equations. The model is applied to the 2D vertical plane flow of a current over two steep-sided trenches for which experiment data are available for comparison and good agreement is obtained. And the model is used to predicting the flow in a channel with a steep-sided submerged breakwater at the bottom, and the streamline is drawn. 展开更多
关键词 free surface flow hydrodynamic pressure kinematic boundary condition Σ-COORDINATE SIMPLE arithmetic
下载PDF
Model for seawater fouling and effects of temperature,flow velocity and surface free energy on seawater fouling 被引量:3
11
作者 Dazhang Yang Jianhua Liu +1 位作者 Xiaoxue E Linlin Jiang 《Chinese Journal of Chemical Engineering》 SCIE EI CAS CSCD 2016年第5期658-664,共7页
A kinetic model was proposed to predict the seawater fouling process in the seawater heat exchangers.The new model adopted an expression combining depositional and removal behaviors for seawater fouling based on the K... A kinetic model was proposed to predict the seawater fouling process in the seawater heat exchangers.The new model adopted an expression combining depositional and removal behaviors for seawater fouling based on the Kern–Seaton model.The present model parameters include the integrated kinetic rate of deposition(k d)and the integrated kinetic rate of removal(k r),which have clear physical signi ficance.A seawater-fouling monitoring device was established to validate the model.The experimental data were well fitted to the model,and the parameters were obtained in different conditions.SEM and EDX analyses were performed after the experiments,and the results show that the main components of seawater fouling are magnesium hydroxide and aluminum hydroxide.The effects of surface temperature,flow velocity and surface free energy were assessed by the model and the experimental data.The results indicate that the seawater fouling becomes aggravated as the surface temperature increased in a certain range,and the seawater fouling resistance reduced as the flow velocity of seawater increased.Furthermore,the effect of the surface free energy of metals was analyzed,showing that the lower surface free energy mitigates the seawater fouling accumulation. 展开更多
关键词 FOULING Seawater MODEL Surface temperature flow velocity Surface free energy
下载PDF
Computational modeling of free-surface slurry flow problems using particle simulation method 被引量:2
12
作者 赵崇斌 彭省临 +2 位作者 刘亮明 B.E.Hobbs A.Ord 《Journal of Central South University》 SCIE EI CAS 2013年第6期1653-1660,共8页
The particle simulation method is used to solve free-surface slurry flow problems that may be encountered in several scientific and engineering fields.The main idea behind the use of the particle simulation method is ... The particle simulation method is used to solve free-surface slurry flow problems that may be encountered in several scientific and engineering fields.The main idea behind the use of the particle simulation method is to treat granular or other materials as an assembly of many particles.Compared with the continuum-mechanics-based numerical methods such as the finite element and finite volume methods,the movement of each particle is accurately described in the particle simulation method so that the free surface of a slurry flow problem can be automatically obtained.The major advantage of using the particle simulation method is that only a simple numerical algorithm is needed to solve the governing equation of a particle simulation system.For the purpose of illustrating how to use the particle simulation method to solve free-surface flow problems,three examples involving slurry flow on three different types of river beds have been considered.The related particle simulation results obtained from these three examples have demonstrated that:1) The particle simulation method is a promising and useful method for solving free-surface flow problems encountered in both the scientific and engineering fields;2) The shape and irregular roughness of a river bed can have a significant effect on the free surface morphologies of slurry flow when it passes through the river bed. 展开更多
关键词 particle simulation free surface slurry flow numerical method
下载PDF
MHD Effect of Liquid Metal Film Flows as Plasma-Facing Components 被引量:1
13
作者 张秀杰 许增裕 潘传杰 《Plasma Science and Technology》 SCIE EI CAS CSCD 2008年第6期685-689,共5页
Stability of liquid metal film flow under gradient magnetic field is investigated. Three dimensional numerical simulations on magnetohydrodynamics (MHD) effect of free surface film flow were carried out, with emphas... Stability of liquid metal film flow under gradient magnetic field is investigated. Three dimensional numerical simulations on magnetohydrodynamics (MHD) effect of free surface film flow were carried out, with emphasis on the film thickness variation and its surface stability. Three different MHD phenomena of film flow were observed in the experiment, namely, retardant, rivulet and flat film flow. From our experiment and numerical simulation it can be concluded that flat film flow is a good choice for plasma-facing components (PFCs) 展开更多
关键词 liquid metal free surface film flow MHD effect numerical simulation
下载PDF
A Study of Two Dimensional Unsteady MHD Free Convection Flow over a Vertical Plate 被引量:1
14
作者 Absana Tarammim Mohammad Sharif Ullah Mohammed Jahir Uddin 《Open Journal of Fluid Dynamics》 2020年第4期342-355,共14页
In this paper, unsteady free convection heat transfer flow over a vertical plate in the presence of a magnetic field is discussed in detail. The dimensionless partial differential equations of continuity, momentum alo... In this paper, unsteady free convection heat transfer flow over a vertical plate in the presence of a magnetic field is discussed in detail. The dimensionless partial differential equations of continuity, momentum along energy are analyzed with suitable transformations. For numerical calculation, an implicit finite difference method is applied to solve a set of nonlinear dimensionless partial differential equations. Dimensionless velocity and temperature profile are also investigated due to the effects of assumed parameters in the concerned problem. An explicit finite difference technique is used to compute velocity and temperature profiles. The stability conditions are also examined. 展开更多
关键词 free Convection Heat Transfer flow Magnetic Field Explicit Finite Difference Method
下载PDF
Effect of Variable Viscosity, Dufour, Soret and Thermal Conductivity on Free Convective Heat and Mass Transfer of Non-Darcian Flow past Porous Flat Surface 被引量:1
15
作者 Isaac L. Animasaun Anselm O. Oyem 《American Journal of Computational Mathematics》 2014年第4期357-365,共9页
The motion of incompressible fluid of a variable fluid viscosity and variable thermal conductivity with thermal radiation, Dufour, Soret with heat and mass transfer over a linearly moving porous vertical semi-infinite... The motion of incompressible fluid of a variable fluid viscosity and variable thermal conductivity with thermal radiation, Dufour, Soret with heat and mass transfer over a linearly moving porous vertical semi-infinite plate with suction is investigated. The governing equations are transformed into a system of coupled nonlinear ordinary differential equations using similarity transformations with dimensionless variables and solved numerically using shooting method with Runge-Kutta fourth-order method and Newton-Raphson’s interpolation scheme implemented in MATLAB. The result showed that with increase in Dufour and Soret parameter, fluid velocity increases and temperature increases with increase in variation of Dufour while, temperature decreases with increase in Soret. The effects of variable fluid viscosity, variable thermal conductivity, thermal radiation, Soret, Dufour, Prandtl and Schmidt parameters on the dimensionless velocity, temperature and concentration profiles are shown graphically. 展开更多
关键词 VARIABLE Viscosity VARIABLE Thermal Conductivity Non-Darcian flow free Convection Dufour SORET POROUS Flat Surface
下载PDF
Phase transition in evolution of traffic flow with scale-free property 被引量:2
16
作者 沈波 高自友 《Chinese Physics B》 SCIE EI CAS CSCD 2008年第9期3284-3288,共5页
This paper investigates the behaviour of traffic flow in traffic systems with a new model based on the NaSch model and cluster approximation of mean-field theory. The proposed model aims at constructing a mapping rela... This paper investigates the behaviour of traffic flow in traffic systems with a new model based on the NaSch model and cluster approximation of mean-field theory. The proposed model aims at constructing a mapping relationship between the microcosmic behaviour and the macroscopic property of traffic flow. Results demonstrate that scale-free phenomenon of the evolution network becomes obvious when the density value of traffic flow reaches at the critical point of phase transition from free flow to traffic congestion, and jamming is limited in this scale-free structure. 展开更多
关键词 phase transition traffic flow scale-free property cellular automata
下载PDF
Kinetic-spectrophotometric determination of free cyanide by stopped-flow reversed flow injection analysis
17
《Journal of Environmental Sciences》 SCIE EI CAS CSCD 1993年第3期367-372,共6页
A stopped-flow reversed flow injection method for the determination of free cyanide is proposed. Pyridine-barbituric acid mixture is injected in the flow system as reagent to form the colour species with cyanide. The ... A stopped-flow reversed flow injection method for the determination of free cyanide is proposed. Pyridine-barbituric acid mixture is injected in the flow system as reagent to form the colour species with cyanide. The flow is stopped when the reagent zone comes in the flow cell, where absorbance-time data are collected at 580nm wavelength. The linear range of the determination is 0.1 -10μg/ml CN-. The sampling rate is 60h-1 and the relative standard deviation is 1.6% (n=16) at 5.0 μg/ml CN-1 level. With satisfactory results, the proposed method was applied to the determination of free cyanide in wastewater without sample pretreatment. 展开更多
关键词 kinetic analysis STOPPED-flow reversed flow injection system free cyanide wastewater.
下载PDF
Transient free convection flow of a viscoelastic fluid over vertical surface
18
作者 H. M. DUWAIRI R. A. DAMSEH +1 位作者 A. J. CHAMKHA M. S. ABDEL-JABER 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 2010年第5期557-564,共8页
In this paper, the viscoelsatic boundary layer flow and the heat transfer near a vertical isothermal impermeable surface and in a quiescent fluid are examined. The governing equations are formulated and solved numeric... In this paper, the viscoelsatic boundary layer flow and the heat transfer near a vertical isothermal impermeable surface and in a quiescent fluid are examined. The governing equations are formulated and solved numerically using MackCormak's technique. The results show excellent agreement with previously published results by a comparision. Representative results for the velocity and temperature profiles, boundary layer thicknesses, Nusselt numbers, and local skin friction coefficients are shown graphically for different values of viscoelsatic parameters. In general, it is found that the velocities increase inside the hydrodynamic boundary layers and the temperatures decrease inside the thermal boundary layers for the viscoelsatic fluid as compared with the Newtonian fluid due to favorable tensile stresses. Consequently, the coefficients of friction and heat transfer enhance for higher viscoelsatic parameters. 展开更多
关键词 viscoelastic flows TRANSIENT free convection heat transfer
下载PDF
Analysis on MHD Stability of Free Surface Jet flow in a Gradient Magnetic Fields
19
作者 许增裕 康伟山 潘传杰 《Plasma Science and Technology》 SCIE EI CAS CSCD 2004年第3期2301-2306,共6页
The simplified modeling for analysis on MHD stability of free surface jet flow in a gradient magnetic fields is based on the theoretical and experimental results on channel liquid metal MHD flow, especially, the resul... The simplified modeling for analysis on MHD stability of free surface jet flow in a gradient magnetic fields is based on the theoretical and experimental results on channel liquid metal MHD flow, especially, the results of MHD flow velocity distribution in cross-section of channels (rectangular duct and circular pipe), and the expected results from the modeling are well agreed with the recent experimental data obtained. It is the first modeling which can efficiently explain the experimental results of liquid-metal free surface jet flow. 展开更多
关键词 free surface jet flow MHD stability MODELING
下载PDF
Three dimensional free convection couette flow with transpiration cooling
20
作者 JAIN N. C. GUPTA P. 《Journal of Zhejiang University-Science A(Applied Physics & Engineering)》 SCIE EI CAS CSCD 2006年第3期340-346,共7页
Free convection flow between two vertical parallel plates with transverse sinusoidal injection of the fluid at the stationary plate and its corresponding removal by constant suction through the plate in uniform motion... Free convection flow between two vertical parallel plates with transverse sinusoidal injection of the fluid at the stationary plate and its corresponding removal by constant suction through the plate in uniform motion has been analyzed. Due to this type of injection velocity, the flow becomes three-dimensional. Analytical expressions for the velocity, temperature, skin friction and rate of heat transfer were obtained. The important characteristics of the problem, namely the skin friction and the rate of heat transfer are discussed in detail with the help of graphs. 展开更多
关键词 Couette flow free convection Three-dimensional flow Transpiration cooling Transverse sinusoidal injection
下载PDF
上一页 1 2 83 下一页 到第
使用帮助 返回顶部