Search

The hybrid machining method called electrochemical mill-grinding is an advanced machining technology, which can achieve high-quality processing of various difficult-to-cut materials. However, achieving rapid transportation and removal of processed products within a small machining gap is a key challenge, which directly affects the final machined surface quality. In this paper, the product transport under different flushing modes was studied through numerical simulation of the flow field. By using the particle tracking method, the time required for the removal of machining products from the machining gap was dynamically simulated. The analysis results indicated that the removal speed of machining products could be significantly improved when the flushing pressure was 0.8 MPa and the electrolyte jet angle was 30°. In addition, machining experiments were conducted. The machining experiment results showed that the machined surface had a metallic luster with clear edge contours. The grinding marks on the machined surface indicated the grinding effect. The continuous machining of complex patterns demonstrated the reliability of this hybrid machining process. [ABSTRACT FROM AUTHOR]

The unified gas-kinetic wave–particle method (UGKWP) has been developed for the multiscale gas, plasma, and multiphase flow transport processes for the past years. In this work, we propose an implicit UGKWP (IUGKWP) method to remove the Courant–Friedrichs–Lewy time step constraint. Based on the local integral solution of the radiative transfer equation (RTE), the particle transport processes are categorized into the long-λ streaming process and the short-λ streaming process compared to a local physical characteristic time tp. In the construction of the IUGKWP method, the long-λ streaming process is tracked by the implicit Monte Carlo method; the short-λ streaming process is evolved by solving the implicit moment equations; and the photon distribution is closed by a local integral solution of RTE. In the IUGKWP method, the multiscale flux of radiation energy and the multiscale closure of photon distribution are constructed based on the local integral solution. The IUGKWP method preserves the second-order asymptotic expansion of RTE in the optically thick regime and adapts its computational complexity to the flow regime. The numerical dissipation is well controlled, and the teleportation error is significantly reduced in the optically thick regime. The computational complexity of the IUGKWP method decreases exponentially as the Knudsen number approaches zero, and the computational efficiency is remarkably improved in the optically thick regime. The IUGKWP is formulated on a generalized unstructured mesh, and multidimensional 2D and 3D algorithms are developed. Numerical tests are presented to validate the capability of IUGKWP in capturing the multiscale photon transport process. The algorithm and code will apply in the engineering applications of inertial confinement fusion. [ABSTRACT FROM AUTHOR]

To achieve rapid, efficient, and low-cost preparation of large-scale stable aluminum alloy superhydrophobic surfaces, a new preparation method is proposed. The outer surface of the array micro-protrusions was coated with a layer of armor, which was the molten spatter produced during picosecond laser processing. The molten sputters and micro-protrusions combined to form micro–nano composite multi-layer structures. Through these special array micro–nano composite multi-layer structures and chemical modification, the wear-resistant and superhydrophobic properties of aluminum alloy surfaces were realized. According to test results, the array micro–nano composite structures prepared by picosecond laser and chemical modification had a water drop contact angle of 154.6° and a water drop rolling angle of 2°, exhibiting excellent superhydrophobic and anti-adhesion properties. Its self-cleaning, corrosion resistance and friction and wear behavior were systematically analyzed. The analysis results showed that the rolling droplets on the prepared surface could easily take away contaminants. The corrosion voltage and corrosion current density of the prepared superhydrophobic surface are significantly lower than that of the raw surface. In addition, a water drop contact angle of the aluminum alloy sample maintained at 145.1° after five wear tests, indicating the prepared surface after wear testing still had hydrophobic performance. The innovative method proposed in this study provides a simple and effective method for preparing large-scale wear-resistant superhydrophobic surface of aluminum alloy. [ABSTRACT FROM AUTHOR]

This paper pioneers the direct numerical simulation (DNS) and physical analysis in supersonic three-temperature carbon dioxide (CO2) turbulent channel flow. CO2 is a linear and symmetric triatomic molecular, with the thermal non-equilibrium three-temperature effects arising from the interactions among translational, rotational, and vibrational modes at room temperature. Thus, the rotational and vibrational modes of CO2 are addressed. The thermal non-equilibrium effect of CO2 has been modeled in an extended three-temperature kinetic model, with the calibrated translational, rotational, and vibrational relaxation time. To solve the extended kinetic equation accurately and robustly, non-equilibrium high-accuracy gas-kinetic scheme is proposed within the well-established two-stage fourth-order framework. Compared with the one-temperature supersonic turbulent channel flow, supersonic three-temperature CO2 turbulence enlarges the ensemble heat transfer of the wall by approximate 20% and slightly decreases the ensemble frictional force. The ensemble density and temperature fields are greatly affected, and there is little change in Van Driest transformation of streamwise velocity. The thermal non-equilibrium three-temperature effects of CO2 also suppress the peak of normalized root mean square of density and temperature, normalized turbulent intensities and Reynolds stress. The vibrational modes of CO2 behave quite differently with rotational and translational modes. Compared with the vibrational temperature fields, the rotational temperature fields have the higher similarity with translational temperature fields, especially in temperature amplitude. Current thermal non-equilibrium models, high-accuracy DNS and physical analysis in supersonic CO2 turbulent flow can act as the benchmark for the long-term applicability of compressible CO2 turbulence. [ABSTRACT FROM AUTHOR]

A unified framework for particulate two-phase flow is presented with a wide range of solid particle concentration from dilute to dense limit. The two-phase flow is simulated by two coupled flow solvers, that is, the gas-kinetic scheme (GKS) for the gas phase and unified gas-kinetic wave–particle method (UGKWP) for the solid particle phase. The GKS is a second-order Navier–Stokes flow solver. The UGKWP is a multiscale method for all flow regimes. The wave and particle decomposition in UGKWP depends on the cell's Knudsen number (Kn). At a small Kn number, the highly concentrated solid particle phase will be modeled by the Eulerian hydrodynamic wave due to the intensive particle–particle collisions. At a large Kn number, the dilute solid particle will be followed by the Lagrangian particle to capture the non-equilibrium transport. In the transition regime, a smooth transition between the above limits is obtained according to the local Kn number. The distribution of solid particles in UGKWP is composed of analytical function and discrete particle, and both condensed and dilute phases can be automatically captured in the most efficient way. In the current scheme, the two-phase model improves the previous one in many aspects, such as drag force model, the frictional pressure formulation, and flux limiting model. The scheme is tested in many typical gas–particle two-phase problems, including the interaction of shock wave with solid particle layer, horizontal pneumatic conveying, bubble formation, and particle cluster phenomena in the fluidized bed. The results validate the GKS-UGKWP for the simulation of gas–particle flow. [ABSTRACT FROM AUTHOR]

In this study, the high-order gas kinetic scheme (GKS) is employed for explicit large eddy simulation (hereafter referred to simply as "LES") and implicit large eddy simulation (iLES) of turbulent channel flows. The main objective is to compare the performance of iLES and LES in the high-order finite volume framework, and study which is most suitable for turbulence simulation. The prediction abilities of different explicit LES models and iLES method on various statistics and flow structures are compared. Most results from both iLES and explicit LES are very close to those of direct numerical simulation. Moreover, iLES is generally superior to explicit LES in predicting several important flow properties, including the mean velocity profiles, Reynolds stress, and Q-criterion iso-surfaces. This superior performance of iLES may arise that the numerical dissipation of the high-order scheme is enough to replace the subgrid dissipation needed in large eddy simulation of turbulence. If the explicit LES model is further adopted, the dissipation will increase, so the results degrade. The overall satisfactory results show that the high-order GKS can provide appropriate numerical dissipation and is suitable for iLES of turbulence. [ABSTRACT FROM AUTHOR]

In this paper, a high-order gas-kinetic scheme is developed on three-dimensional unstructured meshes for compressible Euler and Navier–Stokes equations. To achieve the high-order spatial accuracy, the three-dimensional weighted essentially non-oscillatory (WENO) reconstruction is extended to the unstructured tetrahedral and hexahedral meshes. A simple strategy is adopted for the selection of candidate stencils, and the topologically independent linear weights are used for the spatial reconstruction. The efficiency and robustness of the classical WENO reconstruction are improved. In addition to the two-stage fourth-order temporal discretization and lower–upper symmetric Gauss–Seidel method, the explicit and implicit high-order gas-kinetic schemes are developed for unsteady and steady problems. Accuracy tests on hexahedral and tetrahedral grids validate the third-order of accuracy, and various three-dimensional incompressible and compressible numerical experiments are also presented. The results validate the accuracy and robustness of the proposed scheme for both inviscid and viscous flows. In the future, the current scheme will be extended to the hybrid unstructured meshes and Reynolds-averaged Navier–Stokes simulation with high Reynolds numbers. [ABSTRACT FROM AUTHOR]

A well-balanced compact high-order gas-kinetic scheme (GKS) on unstructured mesh is first developed for solving the shallow water equations with source terms. The distinguishable feature of the finite volume GKS is that based on the gas-kinetic formulation, a time-accurate gas distribution function can be constructed, from which both the fluxes and the flow variables can be explicitly evaluated at the cell interface. As a result, besides the update of cell-averaged conservative variables, the cell-averaged slopes of the flow variables can be updated as well. Equipped with both flow variables and their slopes, a fourth-order compact spatial reconstruction on unstructured mesh can be obtained as the initial condition at the beginning of each time step. For the shallow water flow, in order to preserve the well-balanced property, the advection and the source terms in the flux function have to be balanced properly. The current compact GKS achieves high-order accuracy, keeps the well-balanced property, and has super-robustness in the simulation of bore waves. The scheme is used in the shallow water flow studies, such as dam breaking and bore wave propagation. In addition, the pollution transport, morphodynamics, and bottom friction in the shallow water flow have been included in the scheme. In the end, the water discharge in the Pearl River estuary and the dam-break experiment with movable bed topography have been simulated. [ABSTRACT FROM AUTHOR]

Many kinetic relaxation models have been proposed for the study of rarefied flows. Based on the single relaxation time model, a discrete velocity method-based unified gas-kinetic scheme (UGKS) has been constructed. The UGKS models the gas dynamics on the discretized space directly on account of accumulating flow evolution from particle transport and collision within a time step. Under the UGKS framework, a unified gas-kinetic wave-particle (UGKWP) method has been further developed for non-equilibrium flow simulation, where the time evolution of the gas distribution function is composed of analytical wave and individual particles. In the highly rarefied regime, the flow evolution is mainly described by the particle transport and collision. Because of the use of single relaxation time for particle collision, there is a noticeable discrepancy between the UGKWP solution and the full Boltzmann or direct simulation Monte Carlo (DSMC) result, such as the temperature distribution inside a shock layer at high Mach numbers. In this Letter, a modification of the particle collision time according to the particle velocity will be implemented in the UGKWP. As a result, the new model greatly improves the performance of the UGKWP in the capturing of non-equilibrium flows. There is an excellent match between UGKWP and DSMC or Boltzmann solution in the highly rarefied regime. In the continuum flow regime, due to the absence of particles, the modification of the particle collision time will not take effect and the UGKWP will get back to the hydrodynamic Navier–Stokes flow solver with correct dissipative coefficients at small cell Knudsen numbers. [ABSTRACT FROM AUTHOR]

In this paper, the unified gas-kinetic wave-particle (UGKWP) method has been constructed on a three-dimensional unstructured mesh with parallel computing for multiscale flow simulation. Based on the direct modeling methodology, the unified gas-kinetic scheme (UGKS) models the flow dynamics directly on the numerical mesh size and time step scales, and it is able to capture the flow dynamics from the kinetic scale particle transport to the hydrodynamic wave propagation seamlessly according to the local cell Knudsen number. Instead of discretizing the particle velocity space in UGKS, the UGKWP method is composed of evolution of deterministic wave and stochastic particles. With dynamic wave-particle decomposition according to the cell Knudsen number, the UGKWP method is able to capture the continuum wave interaction and rarefied particle transport under a unified framework and achieves high efficiency in different flow regimes. The UGKWP flow solver is constructed in three-dimensional space and is validated by many test cases at different Mach and Knudsen numbers. The examples include a 3D shock tube problem, lid-driven cubic cavity flow, high-speed flow passing through a cubic object, and hypersonic flow around a space vehicle. The parallel performance has been tested on the Tianhe-2 supercomputer, and reasonable parallel performance has been observed up to 1000 cores. With the wave-particle formulation, the UGKWP method has great potential in solving three-dimensional multiscale transport problems with the co-existence of continuum and rarefied flow regimes, especially for the high-speed rarefied and continuum flow simulation around a space vehicle in near-space flight, where the local Knudsen number can vary significantly with five or six orders of magnitude differences. [ABSTRACT FROM AUTHOR]

In this paper, we present a unified gas-kinetic wave-particle (UGKWP) method on unstructured mesh for the multiscale simulation of continuum and rarefied flow. Inheriting from the multiscale transport in the unified gas-kinetic scheme (UGKS), the integral solution of the kinetic model equation is employed in the construction of the UGKWP method to model the flow physics on the scales of cell size and time step. A novel wave-particle adaptive formulation is introduced in the UGKWP method to describe the flow dynamics in each control volume. The local gas evolution is constructed through the dynamical interaction of the deterministic hydrodynamic wave and the stochastic kinetic particle. To model the gas dynamics on the scales of cell size and time step, the decomposition, interaction, and evolution of the hydrodynamic wave and the kinetic particle depend on the ratio of time step to local collision time. In the rarefied flow regime, the UGKWP method recovers the nonequilibrium flow physics by discrete particles and performs as a stochastic particle method. In the continuum flow regime, the UGKWP method captures the flow behavior solely by macroscopic variable evolution and becomes a gas-kinetic hydrodynamic flow solver, the same as the gas-kinetic scheme, for viscous and heat-conducting Navier–Stokes solutions. In the transition regime, both kinetic particle and hydrodynamic wave contribute adaptively in the UGKWP to capture the peculiar nonequilibrium flow physics in a most efficient way. In different flow regimes, the Sod shock tube, lid-driven cavity flow, laminar boundary layer, and high-speed flow around a circular cylinder are computed to validate the UGKWP method on unstructured mesh. The UGKWP method obtains the same UGKS solutions in all Knudsen regimes. However, with an automatic wave-particle decomposition, the UGKWP method becomes very efficient. For example, at Mach number 30 and Knudsen number 0.1, the UGKWP has several-order-of-magnitude reductions in computational cost and memory requirement in comparison with UGKS. Overall, the UGKWP can capture the gas dynamics in all flow regimes efficiently and accurately from the free molecular transport to the Navier-Stokes flow evolution. [ABSTRACT FROM AUTHOR]

For increasingly non-equilibrium flowfields, the Navier-Stokes equations lose accuracy partially due to the single temperature approximation. To overcome this barrier, a continuum multi-temperature model based on the Bhatnagar-Gross-Krook equation coupled with the Landau-Teller-Jeans relaxation model has been proposed for two-dimensional hypersonic non-equilibrium multi-temperature flow computation. In a recent study, a two-stage fourth-order gas-kinetic scheme (GKS) has been developed for equilibrium flows, which achieves a fourth-order accuracy in space and time as well as high efficiency and robustness. In this paper, targeting for accurate and efficient simulation of multi-temperature non-equilibrium flows, a high-order three-dimensional multi-temperature GKS is constructed under the two-stage fourth-order framework, with the fourth-order Simpson interpolation rule for the newly emerged source term. Simulations on decaying homogeneous isotropic turbulence, low-density nozzle flow, rarefied hypersonic flow over a flat plate, and type IV shock-shock interaction are used to validate the multi-temperature model through the comparison with experimental measurements. The unified gas kinetic scheme (UGKS) results and the Direct Simulation Monte Carlo (DSMC) solutions will be used as well in some cases for validation. Computational results not only confirm the high-order accuracy and quite robustness of this scheme but also show the significant improvement on computational efficiency compared with UGKS and DSMC, especially for the flow in the near continuum regime. [ABSTRACT FROM AUTHOR]

The dynamics of a 2D rotating Crookes radiometer is studied using a moving mesh unified gas kinetic scheme. The whole evolution process of a fan from an initial unsteady start-up to a final steady state rotational movement in a rarefied gas environment is simulated numerically. Through the numerical study, the unsteady force distribution along a vane which dynamically drives the fan movement is captured. And a quantitative connection between total torque and rotational speed of the fan in the Knudsen number regime of 10-3 < Kn < 102 is obtained. Based on the dimensional analysis, the total radiometric torque can be decomposed into a net radiometric driving torque and a rotational resistance. Based on the numerical data, the analytical functions of the torque and angular velocity of a rotating fan in terms of Knudsen number are quantitatively constructed. This relationship is used to explain the experimental observation of the Knudsen number shift for the appearance of the maximum torque and the maximum rotational speed in the transitional flow regime. [ABSTRACT FROM AUTHOR]

A gas-kinetic model with multiple translational temperature for the continuum and near continuum flow simulations is proposed. The main purpose for this work is to derive the generalized Navier-Stokes equations with multiple temperature. It is well recognized that for increasingly rarefied flowfields, the predictions from continuum formulation, such as the Navier-Stokes equations lose accuracy. These inaccuracies may be partially due to the single temperature assumption in the standard Navier-Stokes equations. Here, based on an extended Bhatnagar-Gross-Krook (BGK) model with multiple translational temperature, the numerical scheme for its corresponding Navier-Stokes equations is also constructed. In the current approach, the energy exchange between x, y, and z directions is modeled through the particle collision, and individual energy equation in different direction is obtained. The kinetic model, newly constructed is an enlarged system in comparison with Holway’s ellipsoid statistical BGK model (ES-BGK). The detailed difference is presented in this paper. In the newly derived “Navier-Stokes” equations from the current model, all viscous terms are replaced by the temperature relaxation terms. The relation between the stress and strain in the standard Navier-Stokes equations is recovered only in the limiting case when the flow is close to the equilibrium, such as small temperature differences in different directions. In order to validate the generalized Navier-Stokes equations, we apply them to the study of Couette and Poiseuille flows with a wide range of Knudsen numbers. In the continuum flow regime, the standard Navier-Stokes solutions are precisely recovered. In the near continuum flow regime, the simulation results are compared with the direct simulation Monte Carlo solutions. The anomalous phenomena in the pressure and temperature distributions from the standard Navier-Stokes equations in the Poiseuille flow case at Kn=0.1 are well resolved by the generalized Navier-Stokes equations. This paper clearly shows that many thermal nonequilibrium phenomena in the near continuum flow regime can be well captured by modifying some assumptions in the standard Navier-Stokes equations. [ABSTRACT FROM AUTHOR]

In the slip flow regime with Kn=0.1 for the force and pressure driven Poiseuille flow, Zheng et al. [Rarefied Gas Dynamics, Vol. 23 (Whistler, Canada, 2002)] found out that the Navier-Stokes equations with slip boundary condition could give qualitative different results in the cross-stream direction from the reliable direct simulation Monte Carlo (DSMC) solution. In this Brief Communication, we are going to show that the discrepancy between the Navier-Stokes and the DSMC results can be resolved based on the simulation results of higher-order equations, such as Burnett and super-Burnett ones. [ABSTRACT FROM AUTHOR]

In the continuum transition flow regime, we propose to truncate the Chapman-Enskog expansion of the Boltzmann equation to the Navier-Stokes order without going to the Burnett or super Burnett orders. However, the particle collision time has to be generalized to depend not only on the local macroscopic flow variables, but also their gradients in the rarefied gas regime. In other words, a generalized constitutive relation for the stress and heat flux is proposed. This new model is applied to the study of argon gas shock structure. There is good agreement between the predicted shock structure and experimental results for a wide range of Mach numbers. [ABSTRACT FROM AUTHOR]

The flow in a shock tube is extremely complex with dynamic multi-scale structures of sharp fronts, flow separation, and vortices due to the interaction of the shock wave, the contact surface, and the boundary layer over the side wall of the tube. Prediction and understanding of the complex fluid dynamics are of theoretical and practical importance. It is also an extremely challenging problem for numerical simulation, especially at relatively high Reynolds numbers. Daru and Tenaud [“Evaluation of TVD high resolution schemes for unsteady viscous shocked flows,” Comput. Fluids 30, 89–113 (2001)] proposed a two-dimensional model problem as a numerical test case for high-resolution schemes to simulate the flow field in a square closed shock tube. Though many researchers attempted this problem using a variety of computational methods, there is not yet an agreed-upon grid-converged solution of the problem at the Reynolds number of 1000. This paper presents a rigorous grid-convergence study and the resulting grid-converged solutions for this problem by using a newly developed, efficient, and high-order gas-kinetic scheme. Critical data extracted from the converged solutions are documented as benchmark data. The complex fluid dynamics of the flow at Re = 1000 are discussed and analyzed in detail. Major phenomena revealed by the numerical computations include the downward concentration of the fluid through the curved shock, the formation of the vortices, the mechanism of the shock wave bifurcation, the structure of the jet along the bottom wall, and the Kelvin-Helmholtz instability near the contact surface. Presentation and analysis of those flow processes provide important physical insight into the complex flow physics occurring in a shock tube. [ABSTRACT FROM AUTHOR]