Direct numerical simulation(DNS) of gas–solid flow at high resolution has been carried out by coupling the lattice Boltzmann method(LBM) for gas flow and the discrete element method(DEM) for solid particles. However,...Direct numerical simulation(DNS) of gas–solid flow at high resolution has been carried out by coupling the lattice Boltzmann method(LBM) for gas flow and the discrete element method(DEM) for solid particles. However,the body force periodic boundary condition(FPBC) commonly used to cut down the huge computational cost of such simulation has faced accuracy concerns. In this study, a novel two-region periodic boundary condition(TPBC) is presented to remedy this problem, with the flow driven in the region with body force and freely evolving in the other region. With simulation cases for simple circulating fluidized bed risers, the validity and advantages of TPBC are demonstrated with more reasonable heterogeneity of the particle distribution as compared to the corresponding case with FPBC.展开更多
The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifical...The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifically the capillary gravity waves generated by its motion at the surface. The study analyses the flow of an inviscid, barotropic, and compressible fluid around the stationary solid body. The dynamic behaviour of the fluid is analysed using a two-dimensional coupled Neumann-Kelvin model extended with capillarity and inertia terms. For computational purposes, it is necessary to truncate the unbounded spatial domain with artificial boundaries and then introduce appropriate absorbing boundary conditions. The propagation of short wavelength waves in a convective fluid medium with significant differences in properties between the interior and the surface of the fluid presents a number of difficulties in the design of these conditions. The results are illustrated numerically and commented upon.展开更多
An accurate description of the contact condition between the tool and the workpiece material is one of the most important issues for expounding the underlying multi-physics coupled mechanism during friction stir weldi...An accurate description of the contact condition between the tool and the workpiece material is one of the most important issues for expounding the underlying multi-physics coupled mechanism during friction stir welding(FSW)process.In the present study,a novel asymmetrical boundary condition around the tool-workpiece contact interface is proposed for the FSW of AA2024-T4 alloy.A three-dimensional computational fluid dynamics model is employed for the comparison of the coupled thermal and plastic material flow behavior between asymmetrical and symmetrical boundary conditions.Numerical results of heat generation,temperature distribution,tunnel defect formation and material flow streamline during the welding process are quantitatively analyzed.Besides,various experimental measuring methods are utilized to obtain information about temperature,thermal cycle,tool torque and horizontal cross-section around the exiting keyhole.It is revealed that the modeling results of heat flux density and temperature distribution around the pin,as well as material flow characteristics all change significantly for the two models with different boundary conditions.Particularly,the asymmetrical boundary condition is more capable of predicting temperature fluctuation,plastic material flow along the vertical direction,as well as tunnel defect formation during FSW.Therefore,the superiority of the model with asymmetrical boundary condition over the symmetrical one during the numerical simulation of FSW is elucidated.展开更多
Inverse method was used in single crystal superalloy DD6 processing simulation during solidification. Numerical modeling coupled with experiments has been used to estimate the interface heat transfer coefficient (IHT...Inverse method was used in single crystal superalloy DD6 processing simulation during solidification. Numerical modeling coupled with experiments has been used to estimate the interface heat transfer coefficient (IHTC) between the surface of slab casting and inner mold. Calculated temperature dependent values of IHTC were obtained from a numerical solution. The calculated temperatures agreed well with the measurement of cooling profile.展开更多
Levee or dam failure can cause a significant disaster in most cases. A good prediction of the flood process especially in a real complex terrain is necessary for working out emergency plans for levee or dam breaches. ...Levee or dam failure can cause a significant disaster in most cases. A good prediction of the flood process especially in a real complex terrain is necessary for working out emergency plans for levee or dam breaches. Numerical simulations of levee or dam breach flow were carried out often with constant flow parameters and in relatively simple channels rather than in natural rivers with complex boundaries. This article presents our dedicated studies on the 2-D numerical model of levee or dam breach hydraulics with finite difference schemes. The good performance of the model is demonstrated by comparisons with the theoretical solution of an idealized dam-break flow over a frictionless flat rectangular channel. The model is also validated through its stability and conservation properties. The model is applied to simulate the flood propagation under complex boundary conditions, and the unsteady flood process in a river and in the dry floodplain with a complex bed terrain simultaneously. Furthermore, with respect to engineering practice, the numerical solutions can give special guidance to the effects of parameters such as the flood depth at different sites and the inundated area at different time periods after the levee breach and the travel time of the flood waves, which may be very important for practicing engineers in an efficient flood management.展开更多
Direct numerical simulation of a spatially developing turbulent boundary layer over a compliant wall with anisotropic wall material properties is performed. The Reynolds number varies from 300 to approximately 860 alo...Direct numerical simulation of a spatially developing turbulent boundary layer over a compliant wall with anisotropic wall material properties is performed. The Reynolds number varies from 300 to approximately 860 along the streamwise direction, based on the external flow velocity and the momentum thickness. Eight typical cases are selected for numerical investigation under the guidance of the monoharmonic analysis. The instantaneous flow fields exhibit the traveling wavy motion of the compliant wall, and the frequency-wavenumber power spectrum of wall pressure fluctuation is computed to quantify the mutual influence of the wall compliance and the turbulent flow at different wave numbers. It is shown that the Reynolds shear stress and the pressure fluctuation are generally enhanced by the wall compliance with the parameters considered in the present study. A dynamical decomposition of the skin-friction coefficient is derived, and a new term (CW) appears due to the wall-induced Reynolds shear stress. The influence of the anisotropic compliant wall motion on the turbulent boundary layer through the wall-induced negative Reynolds shear stress is discussed. To elucidate the underlying mechanism, the budget analysis of the Reynolds stresses transportation is further carried out. The impact of the wall compliance on the turbulent flow is disclosed by examining the variations of the diffusion and velocity-pressure correlation terms. It is shown that increase of the Reynolds stresses inside the flow domain is caused by enhancement of the velocity-pressure correlation term, possibly through the long-range influence of the wall compliance on the pressure field, rather than diffusion of the wall-induced Reynolds shear stress into the fluid flow.展开更多
The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general condit...The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.展开更多
The cell model developed since 1950s is a useful tool forexploring the behavior of particle assemblages, but it demandsfurther careful development of the outer boundary conditions so thatinteraction in a particle swar...The cell model developed since 1950s is a useful tool forexploring the behavior of particle assemblages, but it demandsfurther careful development of the outer boundary conditions so thatinteraction in a particle swarm is better represented. In this paper,the cell model and its development were reviewed, and themodifications of outer cell boundary conditions were suggested. Atthe cell outer boundary, the restriction of uniform liquid flow wasremoved in our simulation conducted in the reference frame fixed withthe particle.展开更多
Based on water inrush accident of 1841 working face of Desheng Coal Mine in Wu'an, Hebei province, China, an evaluation model of hydrodynamic characteristics of the project is set up and simulated using Matlab. It...Based on water inrush accident of 1841 working face of Desheng Coal Mine in Wu'an, Hebei province, China, an evaluation model of hydrodynamic characteristics of the project is set up and simulated using Matlab. It is assumed that the pipe flow would transform into seepage flow when the aggregates are plugged into the water inrush channel and the seepage flow would disappear along with grouting process. The simulation results show that the flow velocity will increase with an increase in height of aggregates accumulation body during the aggregates filling process; the maximum seepage velocity occurs on the top of plugging zone; and the water flow decreases with increasing plugging height of water inrush channel. Finally, the field construction results show that the water inrush channel can be plugged effectively by the compacted body prepared with aggregate and cement slurry.展开更多
In this paper the explanation of the mechanism of high-frequency oscillation instability resulted from absorbing boundary conditions is further improved. And we analytically prove the proposition that for one dimensio...In this paper the explanation of the mechanism of high-frequency oscillation instability resulted from absorbing boundary conditions is further improved. And we analytically prove the proposition that for one dimensional discrete model of elastic wave motion, the module of reflection factor will be greater than 1 in high frequency band when artificial wave velocity is greater than 1.5 times the ratio of discrete space step to discrete time step. Based on the proof, the frequency band in which instability occurs is discussed in detail, showing such high-frequency waves are meaningless for the numerical simulation of wave motion.展开更多
We study response of a shear beam to seismic excitations at its base. The research is conducted using computer simulation of the wave propagation on a numerical model. The wave equation is solved using the method of f...We study response of a shear beam to seismic excitations at its base. The research is conducted using computer simulation of the wave propagation on a numerical model. The wave equation is solved using the method of finite differences (FD) where the spatial and temporal derivatives are approximated with FD. We used formulation of the wave equation via the particle velocities, strains, mid stresses. Integrating particle velocities in time, we obtained displacements at spatial points. The main goal in this research is to study phenomena occurring due to three different types of boundary conditions, Dirichlet, Neumann, and moving boundary when simple half-sine pulse propagates through 1D medium modeled as a shear beam.展开更多
For direct numerical simulation (DNS) of turbulent boundary layers, generation of an appropriate inflow condition needs to be considered. This paper proposes a method, with which the inflow condition for spatial-mod...For direct numerical simulation (DNS) of turbulent boundary layers, generation of an appropriate inflow condition needs to be considered. This paper proposes a method, with which the inflow condition for spatial-mode DNS of turbulent boundary layers on supersonic blunt cones with different Mach numbers, Reynolds numbers and wall temperature conditions can be generated. This is based only on a given instant flow field obtained by a temporal-mode DNS of a turbulent boundary layer on a flat plate. Effectiveness of the method is shown in three typical examples by comparing the results with those obtained by other methods.展开更多
This paper presents an efficient numerical method for solving the Euler equations on rectilinear grids. Wall boundary conditions on the surface of an airfoil are implemented by using their first order expansions on th...This paper presents an efficient numerical method for solving the Euler equations on rectilinear grids. Wall boundary conditions on the surface of an airfoil are implemented by using their first order expansions on the airfoil chord line, which is placed along a grid line. However, the method is not restricted to flows with small disturbances since there are no restrictions on the magnitude of the velocity or pressure perturbations. The mathematical formulation and the numerical implementation of the wall boundary conditions in a finite volume Euler code are described. Steady transonic flows are calculated about the NACA 0006, NACA 0012 and NACA 0015 airfoils, corresponding to thickness ratios of 6%, 12%, and 15%, respectively. The computed results, including surface pressure distributions, the lift coefficient, the wave drag coefficient, and the pitching moment coefficient, at angles of attack from 0° to 8° are compared with solutions at the same conditions by FLO52, a well established Euler code using body fitted curvilinear grids. Results demonstrate that the method yields acceptable accuracies even for the relatively thick NACA 0015 airfoil and at high angles of attack. This study establishes the potential of extending the method to computing unsteady fluid structure interaction problems, where the use of a stationary rectilinear grid offers substantial advantages in both computer time and human work since it would not require the generation of time dependent body fitted grids.展开更多
The inflow conditions for spatial direct numerical simulation(SDNS) of turbulent boundary layers should reflect the characteristics of upstream turbulence,which is a puzzle. In this paper a new method is suggested,in ...The inflow conditions for spatial direct numerical simulation(SDNS) of turbulent boundary layers should reflect the characteristics of upstream turbulence,which is a puzzle. In this paper a new method is suggested,in which the flow field obtained by using temporal direct numerical simulation(TDNS) for fully developed turbulent flow(only flow field for a single moment is sufficient) can be used as the inflow of SDNS with a proper transformation. The calculation results confirm that this method is feasible and effective. It is also found that,under a proper time-space transformation,all statistics of the fully developed turbulence obtained by both temporal mode and spatial mode DNS are in excellent agreement with each other,not only qualitatively,but also quantitatively. The normal-wise distributions of mean flow profile,turbulent Mach number and the root mean square(RMS) of the fluctuations of various variables,as well as the Reynolds stresses of the fully developed turbulence obtained by using SDNS,bear similarity in nature.展开更多
Numerical simulation based on computational fluid dynamics (CFD) is a useful approach for quantitatively investigating the underlying thermal-mechanical conditions during FSW, such as temperature field and material ...Numerical simulation based on computational fluid dynamics (CFD) is a useful approach for quantitatively investigating the underlying thermal-mechanical conditions during FSW, such as temperature field and material deformation field. One of the critical issues in CFD simulation of FSW is the use of the frictional boundary condition, which represents the friction between the welding tool and the workpiece in the numerical models. In this study, three-dimensional numerical simulation is conducted to analyze the heat transfer and plastic deformation behaviors during the FSW of AA2024. For comparison purposes, both the boundary velocity (BV) models and the boundary shear stress (BSS) models are employed in order to assess their performances in predicting the temperature and material deformation in FSW. It is interesting to note that different boundary conditions yield similar predictions on temperature, but quite different predictions on material deformation. The numerical predictions are compared with the experimental results. The predicted deformation zone geometry by the BSS model is consistent with the experimental results while there is large difference between the predictions by the BV models and the experimental measurements. The fact that the BSS model yields more reasonable predictions on the deformation zone geometry is attributed to its capacity to automatically adjust the contact state at the tool/workpiece interface. Based on the favorable predictions on both the temperature field and the material deformation field, the BSS model is suggested to have a better performance in numerical simulation of FSW than the BV model.展开更多
This paper describes the numerical simulation of unsteady flows due to incoming wakes and/or varying back pressure,The solution method is based upon the one-step finite-volume TVD Lax-Wendroff scheme.Dual time-step ap...This paper describes the numerical simulation of unsteady flows due to incoming wakes and/or varying back pressure,The solution method is based upon the one-step finite-volume TVD Lax-Wendroff scheme.Dual time-step approach and multigrid algorithm are adopted to improve the computational efficiency of the baseline scheme.Numerical results for the transonic unsteady flow in a channel bump and the unsteady flow in a flat plate cascade and the VKI cascade are presented.展开更多
In order to develop a wall function boundary condition for high-speed flows so as to reduce the grid-dependence of the simula- tion for the skin friction and heat flux, a research was performed to improve the compress...In order to develop a wall function boundary condition for high-speed flows so as to reduce the grid-dependence of the simula- tion for the skin friction and heat flux, a research was performed to improve the compressible wall function boundary condition proposed by Nichols. Values of parameters in the velocity law-of-the-wall were revised according to numerical experiments and the expression of temperature law-of-the-wall was modified based on theoretical analysis and numerical simulation. Be- sides, the formula of the heat conduction term in near-wall region was derived so that the coupling between the wall function boundary condition and CFD code was realized more accurately. Whereafter, the application study of the modified wall func- tion was carried out. The numerical case of supersonic turbulent boundary layer on a flat plate illustrated that the modified wall function produces reasonable results of skin friction and heat flux, and profiles of velocity, temperature and turbulent eddy viscosity for coarse grids with the initial wall spacing of y+〈400, and that the modifications to the original wall function can obviously improve the simulation precision. As for the application of separation flows, it was found from the numerical cases of supersonic cavity flow and hypersonic axisymmetric compression comer that the compressible velocity law-of-the-wall originally established based on the fully-developed attached turbulent boundary layer approximately holds in the near-wall re- gion inside the separation flows, which ensures that reliable skin friction and heat flux can be given by the wall function inside the separation flows, while for the region near separation and reattachment points, the wall function gives results with a rela- tively large error, because the velocity law-of-the-wall used in the wall function takes on obvious deviation from the real ve- locity profiles near the separation and reattachment points.展开更多
The three-dimensional Navier-Stokes characteristic boundary conditions(3D-NSCBC), although physically reasonable and popular in many applications, may encounter the instability problem in simulating complex flows, esp...The three-dimensional Navier-Stokes characteristic boundary conditions(3D-NSCBC), although physically reasonable and popular in many applications, may encounter the instability problem in simulating complex flows, especially for large Reynolds number reactive turbulence where locally the strong reversed flow appears at the outflow boundary surfaces. In the present work, a revised 3D-NSCBC strategy is proposed based on the kinematic relation in different moving coordinate systems. Following this strategy, a systematic formulation is presented for the outflow surface with local reversed flow and can be easily extended to the coupled edge and corner boundaries. Direct numerical simulation(DNS) tests of flow with different turbulence intensities are carried out. Compared with the conventional 3D-NSCBC, the newly proposed method exhibits satisfactory performance to confine numerical instability in the strong reversed flow region. The results confirm the robustness and effectiveness of this newly proposed algorithm.展开更多
As a basic problem in many engineering applications, transition from laminar to turbulence still remains a difficult problem in computational fluid dynamics (CFD). A numerical study of one transitional flow in two-d...As a basic problem in many engineering applications, transition from laminar to turbulence still remains a difficult problem in computational fluid dynamics (CFD). A numerical study of one transitional flow in two-dimensional is conducted by Reynolds averaged numerical simulation (RANS) in this paper. Turbulence model plays a significant role in the complex flows' simulation, and four advanced turbulence models are evaluated. Numerical solution of frictional resistance coefficient is compared with the measured one in the transitional zone, which indicates that Wilcox (2006) k-ω model with correction is the best candidate. Comparisons of numerical and analytical solutions for dimensionless velocity show that averaged streamwise dimensionless velocity profiles correct the shape rapidly in transitional region. Furthermore, turbulence quantities such as turbulence kinetic energy, eddy viscosity, and Reynolds stress are also studied, which are helpful to learn the transition's behavior.展开更多
The nonsplitting perfectly matched layer (NPML) absorbing boundary condition (ABC) was first provided by Wang and Tang (2003) for the finite-difference simulation of elastic wave propagation in solids. In this p...The nonsplitting perfectly matched layer (NPML) absorbing boundary condition (ABC) was first provided by Wang and Tang (2003) for the finite-difference simulation of elastic wave propagation in solids. In this paper, the method is developed to extend the NPML to simulating elastic wave propagation in poroelastic media. Biot's equations are discretized and approximated to a staggered-grid by applying a fourth-order accurate central difference in space and a second-order accurate central difference in time. A cylindrical twolayer seismic model and a borehole model are chosen to validate the effectiveness of the NPML. The results show that the numerical solutions agree well with the solutions of the discrete wavenumber (DW) method.展开更多
基金Supported by the National Natural Science Foundation of China(21821005,91834303)Science Challenge Project(TZ2016001)+1 种基金the Key Research Program of Frontier Science of the Chinese Academy of Sciences(QYZDJ-SSW-JSC029)the Strategic Priority Research Program of the CAS(XDA21030700).
文摘Direct numerical simulation(DNS) of gas–solid flow at high resolution has been carried out by coupling the lattice Boltzmann method(LBM) for gas flow and the discrete element method(DEM) for solid particles. However,the body force periodic boundary condition(FPBC) commonly used to cut down the huge computational cost of such simulation has faced accuracy concerns. In this study, a novel two-region periodic boundary condition(TPBC) is presented to remedy this problem, with the flow driven in the region with body force and freely evolving in the other region. With simulation cases for simple circulating fluidized bed risers, the validity and advantages of TPBC are demonstrated with more reasonable heterogeneity of the particle distribution as compared to the corresponding case with FPBC.
文摘The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifically the capillary gravity waves generated by its motion at the surface. The study analyses the flow of an inviscid, barotropic, and compressible fluid around the stationary solid body. The dynamic behaviour of the fluid is analysed using a two-dimensional coupled Neumann-Kelvin model extended with capillarity and inertia terms. For computational purposes, it is necessary to truncate the unbounded spatial domain with artificial boundaries and then introduce appropriate absorbing boundary conditions. The propagation of short wavelength waves in a convective fluid medium with significant differences in properties between the interior and the surface of the fluid presents a number of difficulties in the design of these conditions. The results are illustrated numerically and commented upon.
基金The authors are grateful for the financial support from the National Natural Science Foundation of China(Grant Nos.52005297 and 52035005).
文摘An accurate description of the contact condition between the tool and the workpiece material is one of the most important issues for expounding the underlying multi-physics coupled mechanism during friction stir welding(FSW)process.In the present study,a novel asymmetrical boundary condition around the tool-workpiece contact interface is proposed for the FSW of AA2024-T4 alloy.A three-dimensional computational fluid dynamics model is employed for the comparison of the coupled thermal and plastic material flow behavior between asymmetrical and symmetrical boundary conditions.Numerical results of heat generation,temperature distribution,tunnel defect formation and material flow streamline during the welding process are quantitatively analyzed.Besides,various experimental measuring methods are utilized to obtain information about temperature,thermal cycle,tool torque and horizontal cross-section around the exiting keyhole.It is revealed that the modeling results of heat flux density and temperature distribution around the pin,as well as material flow characteristics all change significantly for the two models with different boundary conditions.Particularly,the asymmetrical boundary condition is more capable of predicting temperature fluctuation,plastic material flow along the vertical direction,as well as tunnel defect formation during FSW.Therefore,the superiority of the model with asymmetrical boundary condition over the symmetrical one during the numerical simulation of FSW is elucidated.
基金supported by National Basic Research Program of China(No.2005CB724105)National Natural Science Foundation of China (No.10477010)National High Technical Research and Development Program of China(No.2007AA04Z141)
文摘Inverse method was used in single crystal superalloy DD6 processing simulation during solidification. Numerical modeling coupled with experiments has been used to estimate the interface heat transfer coefficient (IHTC) between the surface of slab casting and inner mold. Calculated temperature dependent values of IHTC were obtained from a numerical solution. The calculated temperatures agreed well with the measurement of cooling profile.
基金supported by the National Basic Research and Development Program of China (973 Program,Grant No.2007CB714100)
文摘Levee or dam failure can cause a significant disaster in most cases. A good prediction of the flood process especially in a real complex terrain is necessary for working out emergency plans for levee or dam breaches. Numerical simulations of levee or dam breach flow were carried out often with constant flow parameters and in relatively simple channels rather than in natural rivers with complex boundaries. This article presents our dedicated studies on the 2-D numerical model of levee or dam breach hydraulics with finite difference schemes. The good performance of the model is demonstrated by comparisons with the theoretical solution of an idealized dam-break flow over a frictionless flat rectangular channel. The model is also validated through its stability and conservation properties. The model is applied to simulate the flood propagation under complex boundary conditions, and the unsteady flood process in a river and in the dry floodplain with a complex bed terrain simultaneously. Furthermore, with respect to engineering practice, the numerical solutions can give special guidance to the effects of parameters such as the flood depth at different sites and the inundated area at different time periods after the levee breach and the travel time of the flood waves, which may be very important for practicing engineers in an efficient flood management.
基金the National Natural Science Foundation of China (Grants 11772172 and 11490551).
文摘Direct numerical simulation of a spatially developing turbulent boundary layer over a compliant wall with anisotropic wall material properties is performed. The Reynolds number varies from 300 to approximately 860 along the streamwise direction, based on the external flow velocity and the momentum thickness. Eight typical cases are selected for numerical investigation under the guidance of the monoharmonic analysis. The instantaneous flow fields exhibit the traveling wavy motion of the compliant wall, and the frequency-wavenumber power spectrum of wall pressure fluctuation is computed to quantify the mutual influence of the wall compliance and the turbulent flow at different wave numbers. It is shown that the Reynolds shear stress and the pressure fluctuation are generally enhanced by the wall compliance with the parameters considered in the present study. A dynamical decomposition of the skin-friction coefficient is derived, and a new term (CW) appears due to the wall-induced Reynolds shear stress. The influence of the anisotropic compliant wall motion on the turbulent boundary layer through the wall-induced negative Reynolds shear stress is discussed. To elucidate the underlying mechanism, the budget analysis of the Reynolds stresses transportation is further carried out. The impact of the wall compliance on the turbulent flow is disclosed by examining the variations of the diffusion and velocity-pressure correlation terms. It is shown that increase of the Reynolds stresses inside the flow domain is caused by enhancement of the velocity-pressure correlation term, possibly through the long-range influence of the wall compliance on the pressure field, rather than diffusion of the wall-induced Reynolds shear stress into the fluid flow.
文摘The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.
基金Supported by the National Natural Science Foundation of China (No. 29836130).
文摘The cell model developed since 1950s is a useful tool forexploring the behavior of particle assemblages, but it demandsfurther careful development of the outer boundary conditions so thatinteraction in a particle swarm is better represented. In this paper,the cell model and its development were reviewed, and themodifications of outer cell boundary conditions were suggested. Atthe cell outer boundary, the restriction of uniform liquid flow wasremoved in our simulation conducted in the reference frame fixed withthe particle.
基金Financial support for this work, provided by the National Natural Science Foundation of China (Nos. 41072031, 40172119)the Natural Science Foundation of Hebei Province of China(No. D2012402008)
文摘Based on water inrush accident of 1841 working face of Desheng Coal Mine in Wu'an, Hebei province, China, an evaluation model of hydrodynamic characteristics of the project is set up and simulated using Matlab. It is assumed that the pipe flow would transform into seepage flow when the aggregates are plugged into the water inrush channel and the seepage flow would disappear along with grouting process. The simulation results show that the flow velocity will increase with an increase in height of aggregates accumulation body during the aggregates filling process; the maximum seepage velocity occurs on the top of plugging zone; and the water flow decreases with increasing plugging height of water inrush channel. Finally, the field construction results show that the water inrush channel can be plugged effectively by the compacted body prepared with aggregate and cement slurry.
基金Basic Scientific Research-related Project from Institute of Engineering Mechanics (01180001 and 2007C01)
文摘In this paper the explanation of the mechanism of high-frequency oscillation instability resulted from absorbing boundary conditions is further improved. And we analytically prove the proposition that for one dimensional discrete model of elastic wave motion, the module of reflection factor will be greater than 1 in high frequency band when artificial wave velocity is greater than 1.5 times the ratio of discrete space step to discrete time step. Based on the proof, the frequency band in which instability occurs is discussed in detail, showing such high-frequency waves are meaningless for the numerical simulation of wave motion.
文摘We study response of a shear beam to seismic excitations at its base. The research is conducted using computer simulation of the wave propagation on a numerical model. The wave equation is solved using the method of finite differences (FD) where the spatial and temporal derivatives are approximated with FD. We used formulation of the wave equation via the particle velocities, strains, mid stresses. Integrating particle velocities in time, we obtained displacements at spatial points. The main goal in this research is to study phenomena occurring due to three different types of boundary conditions, Dirichlet, Neumann, and moving boundary when simple half-sine pulse propagates through 1D medium modeled as a shear beam.
基金the National Natural Science Foundation of China(Nos.10632050,90716007)the Special Foundation for the Authors of National Excellent Doctoral Dissertations(No.200328)
文摘For direct numerical simulation (DNS) of turbulent boundary layers, generation of an appropriate inflow condition needs to be considered. This paper proposes a method, with which the inflow condition for spatial-mode DNS of turbulent boundary layers on supersonic blunt cones with different Mach numbers, Reynolds numbers and wall temperature conditions can be generated. This is based only on a given instant flow field obtained by a temporal-mode DNS of a turbulent boundary layer on a flat plate. Effectiveness of the method is shown in three typical examples by comparing the results with those obtained by other methods.
文摘This paper presents an efficient numerical method for solving the Euler equations on rectilinear grids. Wall boundary conditions on the surface of an airfoil are implemented by using their first order expansions on the airfoil chord line, which is placed along a grid line. However, the method is not restricted to flows with small disturbances since there are no restrictions on the magnitude of the velocity or pressure perturbations. The mathematical formulation and the numerical implementation of the wall boundary conditions in a finite volume Euler code are described. Steady transonic flows are calculated about the NACA 0006, NACA 0012 and NACA 0015 airfoils, corresponding to thickness ratios of 6%, 12%, and 15%, respectively. The computed results, including surface pressure distributions, the lift coefficient, the wave drag coefficient, and the pitching moment coefficient, at angles of attack from 0° to 8° are compared with solutions at the same conditions by FLO52, a well established Euler code using body fitted curvilinear grids. Results demonstrate that the method yields acceptable accuracies even for the relatively thick NACA 0015 airfoil and at high angles of attack. This study establishes the potential of extending the method to computing unsteady fluid structure interaction problems, where the use of a stationary rectilinear grid offers substantial advantages in both computer time and human work since it would not require the generation of time dependent body fitted grids.
基金the National Natural Science Foundation of China (Grant No. 90205021)the China Postdoctoral Science Foundation (Grant No. 20060400707)the Foundation for the Author of National Excellent Doctoral Dissertation of China (Grant No. 200328)
文摘The inflow conditions for spatial direct numerical simulation(SDNS) of turbulent boundary layers should reflect the characteristics of upstream turbulence,which is a puzzle. In this paper a new method is suggested,in which the flow field obtained by using temporal direct numerical simulation(TDNS) for fully developed turbulent flow(only flow field for a single moment is sufficient) can be used as the inflow of SDNS with a proper transformation. The calculation results confirm that this method is feasible and effective. It is also found that,under a proper time-space transformation,all statistics of the fully developed turbulence obtained by both temporal mode and spatial mode DNS are in excellent agreement with each other,not only qualitatively,but also quantitatively. The normal-wise distributions of mean flow profile,turbulent Mach number and the root mean square(RMS) of the fluctuations of various variables,as well as the Reynolds stresses of the fully developed turbulence obtained by using SDNS,bear similarity in nature.
基金supported by the National Natural Science Foundation of China(Grant No.51375259 and Grant No.51705280)the Ministry of Science and Technology of China(Grant No.2012ZX04012-011)+1 种基金Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund(the second phase,Grant No.U1501501)the Tsinghua National Laboratory for Information Science and Technology
文摘Numerical simulation based on computational fluid dynamics (CFD) is a useful approach for quantitatively investigating the underlying thermal-mechanical conditions during FSW, such as temperature field and material deformation field. One of the critical issues in CFD simulation of FSW is the use of the frictional boundary condition, which represents the friction between the welding tool and the workpiece in the numerical models. In this study, three-dimensional numerical simulation is conducted to analyze the heat transfer and plastic deformation behaviors during the FSW of AA2024. For comparison purposes, both the boundary velocity (BV) models and the boundary shear stress (BSS) models are employed in order to assess their performances in predicting the temperature and material deformation in FSW. It is interesting to note that different boundary conditions yield similar predictions on temperature, but quite different predictions on material deformation. The numerical predictions are compared with the experimental results. The predicted deformation zone geometry by the BSS model is consistent with the experimental results while there is large difference between the predictions by the BV models and the experimental measurements. The fact that the BSS model yields more reasonable predictions on the deformation zone geometry is attributed to its capacity to automatically adjust the contact state at the tool/workpiece interface. Based on the favorable predictions on both the temperature field and the material deformation field, the BSS model is suggested to have a better performance in numerical simulation of FSW than the BV model.
文摘This paper describes the numerical simulation of unsteady flows due to incoming wakes and/or varying back pressure,The solution method is based upon the one-step finite-volume TVD Lax-Wendroff scheme.Dual time-step approach and multigrid algorithm are adopted to improve the computational efficiency of the baseline scheme.Numerical results for the transonic unsteady flow in a channel bump and the unsteady flow in a flat plate cascade and the VKI cascade are presented.
基金supported by the National Natural Science Foundation of China(Grant No.11202014)
文摘In order to develop a wall function boundary condition for high-speed flows so as to reduce the grid-dependence of the simula- tion for the skin friction and heat flux, a research was performed to improve the compressible wall function boundary condition proposed by Nichols. Values of parameters in the velocity law-of-the-wall were revised according to numerical experiments and the expression of temperature law-of-the-wall was modified based on theoretical analysis and numerical simulation. Be- sides, the formula of the heat conduction term in near-wall region was derived so that the coupling between the wall function boundary condition and CFD code was realized more accurately. Whereafter, the application study of the modified wall func- tion was carried out. The numerical case of supersonic turbulent boundary layer on a flat plate illustrated that the modified wall function produces reasonable results of skin friction and heat flux, and profiles of velocity, temperature and turbulent eddy viscosity for coarse grids with the initial wall spacing of y+〈400, and that the modifications to the original wall function can obviously improve the simulation precision. As for the application of separation flows, it was found from the numerical cases of supersonic cavity flow and hypersonic axisymmetric compression comer that the compressible velocity law-of-the-wall originally established based on the fully-developed attached turbulent boundary layer approximately holds in the near-wall re- gion inside the separation flows, which ensures that reliable skin friction and heat flux can be given by the wall function inside the separation flows, while for the region near separation and reattachment points, the wall function gives results with a rela- tively large error, because the velocity law-of-the-wall used in the wall function takes on obvious deviation from the real ve- locity profiles near the separation and reattachment points.
基金the funding support by National Science Foundation China(NSFC)under the grant No.91441116NSFC-CNRS joint research project(No.11611130099,NSFC China and PRC 2016-2018 LATUMAR"Turbulence lagrangienne:études numériques et applications environnementales marines",CNRS,France)
文摘The three-dimensional Navier-Stokes characteristic boundary conditions(3D-NSCBC), although physically reasonable and popular in many applications, may encounter the instability problem in simulating complex flows, especially for large Reynolds number reactive turbulence where locally the strong reversed flow appears at the outflow boundary surfaces. In the present work, a revised 3D-NSCBC strategy is proposed based on the kinematic relation in different moving coordinate systems. Following this strategy, a systematic formulation is presented for the outflow surface with local reversed flow and can be easily extended to the coupled edge and corner boundaries. Direct numerical simulation(DNS) tests of flow with different turbulence intensities are carried out. Compared with the conventional 3D-NSCBC, the newly proposed method exhibits satisfactory performance to confine numerical instability in the strong reversed flow region. The results confirm the robustness and effectiveness of this newly proposed algorithm.
基金Foundation item: Supported by the National Natural Science Foundation of China (Nos. 51309040, 51379025), and the Fundamental Research Funds for the Central Universities (Nos. 3132014224, 3132014318).
文摘As a basic problem in many engineering applications, transition from laminar to turbulence still remains a difficult problem in computational fluid dynamics (CFD). A numerical study of one transitional flow in two-dimensional is conducted by Reynolds averaged numerical simulation (RANS) in this paper. Turbulence model plays a significant role in the complex flows' simulation, and four advanced turbulence models are evaluated. Numerical solution of frictional resistance coefficient is compared with the measured one in the transitional zone, which indicates that Wilcox (2006) k-ω model with correction is the best candidate. Comparisons of numerical and analytical solutions for dimensionless velocity show that averaged streamwise dimensionless velocity profiles correct the shape rapidly in transitional region. Furthermore, turbulence quantities such as turbulence kinetic energy, eddy viscosity, and Reynolds stress are also studied, which are helpful to learn the transition's behavior.
基金This research was supported by Natural Science Foundation of China (No. 403740043).
文摘The nonsplitting perfectly matched layer (NPML) absorbing boundary condition (ABC) was first provided by Wang and Tang (2003) for the finite-difference simulation of elastic wave propagation in solids. In this paper, the method is developed to extend the NPML to simulating elastic wave propagation in poroelastic media. Biot's equations are discretized and approximated to a staggered-grid by applying a fourth-order accurate central difference in space and a second-order accurate central difference in time. A cylindrical twolayer seismic model and a borehole model are chosen to validate the effectiveness of the NPML. The results show that the numerical solutions agree well with the solutions of the discrete wavenumber (DW) method.