In this paper,an improved computational fluid dynamic(CFD)model for gas-liquid flow in bubble column was developed using the one-equation Wary-Agarwal(WA)turbulence model coupled with the population balance model(PBM)...In this paper,an improved computational fluid dynamic(CFD)model for gas-liquid flow in bubble column was developed using the one-equation Wary-Agarwal(WA)turbulence model coupled with the population balance model(PBM).Through 18 orthogonal test cases,the optimal combination of interfacial force models,including drag force,lift force,turbulent dispersion force.The modified wall lubrication force model was proposed to improve the predictive ability for hydrodynamic behavior near the wall of the bubble column.The values simulated by optimized CFD model were in agreement with experimental data,and the errors were within±20%.In addition,the axial velocity,turbulent kinetic energy,bubble size distribution,and the dynamic characteristic of bubble plume were analyzed at different superficial gas velocities.This research work could provide a theoretical basis for the extension of the CFD-PBM coupled model to other multiphase reactors..展开更多
Surface waves have a considerable effect on vertical mixing in the upper ocean.In the past two decades,the vertical mixing induced through nonbreaking surface waves has been used in ocean and climate models to improve...Surface waves have a considerable effect on vertical mixing in the upper ocean.In the past two decades,the vertical mixing induced through nonbreaking surface waves has been used in ocean and climate models to improve the simulation of the upper ocean.Thus far,several nonbreaking wave-induced mixing parameterization schemes have been proposed;however,no quantitative comparison has been performed among them.In this paper,a one-dimensional ocean model was used to compare the performances of five schemes,including those of Qiao et al.(Q),Hu and Wang(HW),Huang and Qiao(HQ),Pleskachevsky et al.(P),and Ghantous and Babanin(GB).Similar to previous studies,all of these schemes can decrease the simulated sea surface temperature(SST),increase the subsurface temperature,and deepen the mixed layer,thereby alleviating the common thermal deviation problem of the ocean model for upper ocean simulation.Among these schemes,the HQ scheme exhibited the weakest wave-induced mixing effect,and the HW scheme exhibited the strongest effect;the other three schemes exhibited roughly the same effect.In particular,the Q and P schemes exhibited nearly the same effect.In the simulation based on observations from the Ocean Weather Station Papa,the HQ scheme exhibited the best performance,followed by the Q scheme.In the experiment with the HQ scheme,the root-mean-square deviation of the simulated SST from the observations was 0.43℃,and the mixed layer depth(MLD)was 2.0 m.As a contrast,the deviations of the SST and MLD reached 1.25℃ and 8.4 m,respectively,in the experiment without wave-induced mixing.展开更多
Machine-learned augmentations to turbulence models can be advantageous for flows within the training dataset but can often cause harm outside.This lack of generalizability arises because the constants(as well as the f...Machine-learned augmentations to turbulence models can be advantageous for flows within the training dataset but can often cause harm outside.This lack of generalizability arises because the constants(as well as the functions)in a Reynolds-averaged Navier–Stokes(RANS)model are coupled,and un-constrained re-calibration of these constants(and functions)can disrupt the calibrations of the baseline model,the preservation of which is critical to the model's generalizability.To safeguard the behaviors of the baseline model beyond the training dataset,machine learning must be constrained such that basic calibrations like the law of the wall are kept intact.This letter aims to identify such constraints in two-equation RANS models so that future machine learning work can be performed without violating these constraints.We demonstrate that the identified constraints are not limiting.Furthermore,they help preserve the generalizability of the baseline model.展开更多
A growing interest has been devoted to the contra-rotating propellers (CRPs) due to their high propulsive efficiency, torque balance, low fuel consumption, low cavitations, low noise performance and low hull vibrati...A growing interest has been devoted to the contra-rotating propellers (CRPs) due to their high propulsive efficiency, torque balance, low fuel consumption, low cavitations, low noise performance and low hull vibration. Compared with the single-screw system, it is more difficult for the open water performance prediction because forward and aft propellers interact with each other and generate a more complicated flow field around the CRPs system. The current work focuses on the open water performance prediction of contra-rotating propellers by RANS and sliding mesh method considering the effect of computational time step size and turbulence model. The validation study has been performed on two sets of contra-rotating propellers developed by David W Taylor Naval Ship R & D center. Compared with the experimental data, it shows that RANS with sliding mesh method and SST k-ω turbulence model has a good precision in the open water performance prediction of contra-rotating propellers, and small time step size can improve the level of accuracy for CRPs with the same blade number of forward and aft propellers, while a relatively large time step size is a better choice for CRPs with different blade numbers.展开更多
The mechanism of local scour around submarine pipelines is studied numerically based on a renormalized group (RNG) turbulence model. To validate the numerical model, the equilibrium profiles of local scour for two c...The mechanism of local scour around submarine pipelines is studied numerically based on a renormalized group (RNG) turbulence model. To validate the numerical model, the equilibrium profiles of local scour for two cases are simulated and compared with the experimental data. It shows that the RNG turbulence model can give an appropriate prediction for the configuration of equilibrium scour hole, and it is applicable to this situation. The local scour mechanism around submarine pipelines including the flow structure, shear stress distribution and pressure field is then analyzed and compared with experiments. For further comparison and validation, especially for the flow structure, a numerical calculation employing the large eddy simulation (LES) is also conducted. The numerical results of RNG demonstrate that the critical factor governing the equilibrium profile is the seabed shear stress distribution in the case of bed load sediment transport, and the two-equation RNG turbulence model coupled with the law of wall is capable of giving a satisfying estimation for the bed shear stress. Moreover, the piping phenomena due to the great difference of pressure between the upstream and downstream parts of pipelines and the vortex structure around submarine pipelines are also simulated successfully, which are believed to be the important factor that lead to the onset of local scour.展开更多
A steady three-dimensional fluid flow and solidification model was built based on CFD software by high-Reynolds-number and Lam-Bremhorst low-Reynolds-number k-ε model.During the simulation,the fixed-grid enthalpy-por...A steady three-dimensional fluid flow and solidification model was built based on CFD software by high-Reynolds-number and Lam-Bremhorst low-Reynolds-number k-ε model.During the simulation,the fixed-grid enthalpy-porosity technique was used to represent the solidification,and Darcy law was adopted to simulate the flow in mushy region.The prediction for steel flow and solidification was evaluated by the comparison of two turbulence models.It is found that both Lam-Bremhorst low-Reynolds-number and high-Reynolds-number k-ε models predict the same trend of the steel flow and temperature distribution.However,due to the effect of turbulent flow on heat transfer,the low-Reynolds-number turbulence model predicts longer penetration depth of molten steel in sub-mold region,less shell growth and higher shell surface temperature at the narrow face compared with standard k-ε model.展开更多
Cavitation typically occurs when the fluid pressure is lower than the vapor pressure at a local thermodynamic state, and the flow is frequently unsteady and turbulent. To assess the state-of-the-art of computational c...Cavitation typically occurs when the fluid pressure is lower than the vapor pressure at a local thermodynamic state, and the flow is frequently unsteady and turbulent. To assess the state-of-the-art of computational capabilities for unsteady cavitating flows, different cavitation and turbulence model combinations are conducted. The selected cavitation models include several widely-used models including one based on phenomenological argument and the other utilizing interface dynamics. The k-e turbulence model with additional implementation of the filter function and density correction function are considered to reduce the eddy viscosity according to the computed turbulence length scale and local fluid density respectively. We have also blended these alternative cavitation and lustrate that the eddy viscosity turbulence treatments, to ilnear the closure region can significantly influence the capture of detached cavity. From the experimental validations regarding the force analysis, frequency, and the cavity visualization, no single model combination performs best in all aspects. Furthermore, the implications of parameters contained in different cavitation models are investigated. The phase change process is more pronounced around the detached cavity, which is better illustrated by the interfacial dynamics model. Our study provides insight to aid further modeling development.展开更多
Data-driven turbulence modeling studies have reached such a stage that the basic framework is settled,but several essential issues remain that strongly affect the performance.Two problems are studied in the current re...Data-driven turbulence modeling studies have reached such a stage that the basic framework is settled,but several essential issues remain that strongly affect the performance.Two problems are studied in the current research:(1)the processing of the Reynolds stress tensor and(2)the coupling method between the machine learning model and flow solver.For the Reynolds stress processing issue,we perform the theoretical derivation to extend the relevant tensor arguments of Reynolds stress.Then,the tensor representation theorem is employed to give the complete irreducible invariants and integrity basis.An adaptive regularization term is employed to enhance the representation performance.For the coupling issue,an iterative coupling framework with consistent convergence is proposed and then applied to a canonical separated flow.The results have high consistency with the direct numerical simulation true values,which proves the validity of the current approach.展开更多
The emerging push of the differentiable programming paradigm in scientific computing is conducive to training deep learning turbulence models using indirect observations.This paper demonstrates the viability of this a...The emerging push of the differentiable programming paradigm in scientific computing is conducive to training deep learning turbulence models using indirect observations.This paper demonstrates the viability of this approach and presents an end-to-end differentiable framework for training deep neural networks to learn eddy viscosity models from indirect observations derived from the velocity and pressure fields.The framework consists of a Reynolds-averaged Navier–Stokes(RANS)solver and a neuralnetwork-represented turbulence model,each accompanied by its derivative computations.For computing the sensitivities of the indirect observations to the Reynolds stress field,we use the continuous adjoint equations for the RANS equations,while the gradient of the neural network is obtained via its built-in automatic differentiation capability.We demonstrate the ability of this approach to learn the true underlying turbulence closure when one exists by training models using synthetic velocity data from linear and nonlinear closures.We also train a linear eddy viscosity model using synthetic velocity measurements from direct numerical simulations of the Navier–Stokes equations for which no true underlying linear closure exists.The trained deep-neural-network turbulence model showed predictive capability on similar flows.展开更多
Five turbulence models of Reynolds average Navier-Stokes(RANS),including the standard k-ω model,the RNG k-e model taking into account the low Reynolds number effect,the realizable k-ω model,the SST k-ω model,and th...Five turbulence models of Reynolds average Navier-Stokes(RANS),including the standard k-ω model,the RNG k-e model taking into account the low Reynolds number effect,the realizable k-ω model,the SST k-ω model,and the Reynolds stress model(RSM),are employed in the numerical simulations of direct current(DC)arc plasma torches in the range of arc current from 80 A to 240 A and air gas flow rate from 10 m^3 h^-1 to 50 m^3 h^-1.The calculated voltage,electric field intensity,and the heat loss in the arc chamber are compared with the experiments.The results indicate that the arc voltage,the electric field,and the heat loss in the arc chamber calculated by using the standard k-ω model,the RNG k-ωmodel taking into account the low Reynolds number effect,and the realizable k-ω model are much larger than those in the experiments.The RSM predicts relatively close results to the experiments,but fails in the trend of heat loss varying with the gas flow rate.The calculated results of the SST k-ω model are in the best agreement with the experiments,which may be attributed to the reasonable predictions of the turbulence as well as its distribution.展开更多
Most current computations of trubulent flows with second-moment closure adopt the diffusion mo- dels which neglect the effect of pressure-velocity correlation.ln the present paper the importance of this correlation ef...Most current computations of trubulent flows with second-moment closure adopt the diffusion mo- dels which neglect the effect of pressure-velocity correlation.ln the present paper the importance of this correlation effect is elucidated the neglect of this effect accounts for some major defects in the wide application of the se- cond-moment closures.Through the relation between and ,established by Lumley,we propose here a new turbulence diffusion model which takes into consideration the pressure effect.Applications of this new model in the computation of shearless turbulence mixing layer and plane-and round-jet flows show that the spreading rate of these flows can be satisfactorily captured.展开更多
The main theoretical turbulence models for application to hydrodynamic lubrication problems were briefly reviewed, and the course of their development and their fundamentals were explained. Predictions by these models...The main theoretical turbulence models for application to hydrodynamic lubrication problems were briefly reviewed, and the course of their development and their fundamentals were explained. Predictions by these models on flow fields in turbulent Couette flows and shear-induced countercurrent flows were compared to existing measurements, and Zhang & Zhang's combined k-ε model was shown to have surpassingly satisfactory results. The method of application of this combined k-ε model to high speed journal bearings and annular seals was summarized, and the predicted results were shown to be satisfactory by comparisons with existing experiments of journal bearings and annular seals.展开更多
A parameter estimation algorithm is introduced and used to determine the parameters in the standard k-epsilon two equation turbulence model (SKE). It can be found from the estimation results that although the paramete...A parameter estimation algorithm is introduced and used to determine the parameters in the standard k-epsilon two equation turbulence model (SKE). It can be found from the estimation results that although the parameter estimation method is an effective method to determine model parameters, it is. difficult to obtain a set of parameters for SKE to suit all kinds of separated flow and a modification of the turbulence model structure should be considered. So, a new nonlinear k-e two-equation model (NNKE) is put forward in this paper and the corresponding parameter estimation technique is applied to determine the model parameters. By implementing the NNKE to solve some engineering turbulent flows, it is shown that NNKE is more accurate and versatile than SKE. Thus, the success of NNKE implies that the parameter estimation technique may have a bright prospect in engineering turbulence model research.展开更多
The research is motivated by the ongoing the electromagnetic continuous casting of molten metal. The revised k-ε model considering the effect of magnetic field application was derived. The specific model equations fo...The research is motivated by the ongoing the electromagnetic continuous casting of molten metal. The revised k-ε model considering the effect of magnetic field application was derived. The specific model equations for the electromagnetic braking were used to calculate the velocity distribution in the continuous casting mold of steel. The results show that the revised k-ε model considering the effect of magnetic field application tends to suppress the production of turbulence and difference between the conventional and revised k-e model is small.展开更多
A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carr...A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and εequations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model.展开更多
This study numerically and experimentally investigates the effects of wave loads on a monopile-type offshore wind turbine placed on a 1:25 slope at different water depths as well as the effect of choosing different tu...This study numerically and experimentally investigates the effects of wave loads on a monopile-type offshore wind turbine placed on a 1:25 slope at different water depths as well as the effect of choosing different turbulence models on the efficiency of the numerical model.The numerical model adopts a two-phase flow by solving Unsteady Reynolds-Averaged Navier−Stokes(URANS)equations using the Volume Of Fluid(VOF)method and three differentk-ωturbulence models.Typical environmental conditions from the East China Sea are studied.The wave run-up and the wave loads applied on the monopile are investigated and compared with relevant experimental data as well as with mathematical predictions based on relevant theories.The numerical model is well validated against the experimental data at model scale.The use of different turbulence models results in different predictions on the wave height but less differences on the wave period.The baseline k-ωturbulence model and Shear-Stress Transport(SST)k-ωturbulence model exhibit better performance on the prediction of hydrodynamic load,at a model-scale water depth of 0.42 m,while the laminar model provides better results for large water depths.The SST turbulence model performs better in predicting wave run-up for water depth 0.42 m,while the laminar model and standard k-ωmodel perform better at water depth 0.52 m and 0.62 m,respectively.展开更多
We investigate the role of extended intrinsic mean spin tensor introduced in this work for turbulence modelling in a non-inertial frame of reference. It is described by the Euclidean group of transformations and, in p...We investigate the role of extended intrinsic mean spin tensor introduced in this work for turbulence modelling in a non-inertial frame of reference. It is described by the Euclidean group of transformations and, in particular, its significance and importance in the approach of the algebraic Reynolds stress modelling, such as in a nonlinear K-ε model. To this end and for illustration of the effect of extended intrinsic spin tensor on turbulence modelling, we examine several recently developed nonlinear K-ε models and compare their performance in predicting the homogeneous turbulent shear flow in a rotating frame of reference with LES data. Our results and analysis indicate that, only if the deficiencies of these models and the like be well understood and properly corrected, may in the near future, more sophisticated nonlinear K-ε models be developed to better predict complex turbulent flows in a non-inertial frame of reference.展开更多
This paper describes a newly developednon-isotropic multiple-scale turbulence model (MS/ASM) for complex flow calculations. This model focuses on the direct modeling of Reynolds stresses and utilizes split-spectrum co...This paper describes a newly developednon-isotropic multiple-scale turbulence model (MS/ASM) for complex flow calculations. This model focuses on the direct modeling of Reynolds stresses and utilizes split-spectrum concepts to model multiple-scale effects in turbulence. Validation studies on free shear flows, rotating flows and recirculating flows show that the current model performs significantly better than the single-scale k-e model. The present model is relatively inexpensive in terms of CPU time which makes in suitable for broad engineering flow applications.展开更多
By analyzing the components of Reynolds stresses of implicit algebraic stress model(IASM) in this paper, that Reynolds stresses in buoyant turbulent flows were produced by both strain and buoyancy is considered. Conse...By analyzing the components of Reynolds stresses of implicit algebraic stress model(IASM) in this paper, that Reynolds stresses in buoyant turbulent flows were produced by both strain and buoyancy is considered. Consequently, a nonlinear anisotropy buoyant turbulence model was developed by applying linearity of equilibrium hypothesis to Reynolds stress transports. The model avoids numerical singularity and its reliability is verified by the comparisons between predictions and experimental data.展开更多
基金supported by the National Natural Science Foundation of China(22078009)National Key Research and Development Program of China(2021YFC3001102,2021YFC3001100)。
文摘In this paper,an improved computational fluid dynamic(CFD)model for gas-liquid flow in bubble column was developed using the one-equation Wary-Agarwal(WA)turbulence model coupled with the population balance model(PBM).Through 18 orthogonal test cases,the optimal combination of interfacial force models,including drag force,lift force,turbulent dispersion force.The modified wall lubrication force model was proposed to improve the predictive ability for hydrodynamic behavior near the wall of the bubble column.The values simulated by optimized CFD model were in agreement with experimental data,and the errors were within±20%.In addition,the axial velocity,turbulent kinetic energy,bubble size distribution,and the dynamic characteristic of bubble plume were analyzed at different superficial gas velocities.This research work could provide a theoretical basis for the extension of the CFD-PBM coupled model to other multiphase reactors..
基金supported by the Laoshan Laboratory(No.LSKJ202201600)the National Key Research and Development Program of China(No.2022YFC2808304).
文摘Surface waves have a considerable effect on vertical mixing in the upper ocean.In the past two decades,the vertical mixing induced through nonbreaking surface waves has been used in ocean and climate models to improve the simulation of the upper ocean.Thus far,several nonbreaking wave-induced mixing parameterization schemes have been proposed;however,no quantitative comparison has been performed among them.In this paper,a one-dimensional ocean model was used to compare the performances of five schemes,including those of Qiao et al.(Q),Hu and Wang(HW),Huang and Qiao(HQ),Pleskachevsky et al.(P),and Ghantous and Babanin(GB).Similar to previous studies,all of these schemes can decrease the simulated sea surface temperature(SST),increase the subsurface temperature,and deepen the mixed layer,thereby alleviating the common thermal deviation problem of the ocean model for upper ocean simulation.Among these schemes,the HQ scheme exhibited the weakest wave-induced mixing effect,and the HW scheme exhibited the strongest effect;the other three schemes exhibited roughly the same effect.In particular,the Q and P schemes exhibited nearly the same effect.In the simulation based on observations from the Ocean Weather Station Papa,the HQ scheme exhibited the best performance,followed by the Q scheme.In the experiment with the HQ scheme,the root-mean-square deviation of the simulated SST from the observations was 0.43℃,and the mixed layer depth(MLD)was 2.0 m.As a contrast,the deviations of the SST and MLD reached 1.25℃ and 8.4 m,respectively,in the experiment without wave-induced mixing.
基金supported by the Air Force Office of Scientific Research(Grant No.FA9550-23-1-0272)the National Natural Science Foundation of China(Grant Nos.11988102 and 91752202).
文摘Machine-learned augmentations to turbulence models can be advantageous for flows within the training dataset but can often cause harm outside.This lack of generalizability arises because the constants(as well as the functions)in a Reynolds-averaged Navier–Stokes(RANS)model are coupled,and un-constrained re-calibration of these constants(and functions)can disrupt the calibrations of the baseline model,the preservation of which is critical to the model's generalizability.To safeguard the behaviors of the baseline model beyond the training dataset,machine learning must be constrained such that basic calibrations like the law of the wall are kept intact.This letter aims to identify such constraints in two-equation RANS models so that future machine learning work can be performed without violating these constraints.We demonstrate that the identified constraints are not limiting.Furthermore,they help preserve the generalizability of the baseline model.
基金supported by the National Natural Science Foundation of China(Grant No.51079157)
文摘A growing interest has been devoted to the contra-rotating propellers (CRPs) due to their high propulsive efficiency, torque balance, low fuel consumption, low cavitations, low noise performance and low hull vibration. Compared with the single-screw system, it is more difficult for the open water performance prediction because forward and aft propellers interact with each other and generate a more complicated flow field around the CRPs system. The current work focuses on the open water performance prediction of contra-rotating propellers by RANS and sliding mesh method considering the effect of computational time step size and turbulence model. The validation study has been performed on two sets of contra-rotating propellers developed by David W Taylor Naval Ship R & D center. Compared with the experimental data, it shows that RANS with sliding mesh method and SST k-ω turbulence model has a good precision in the open water performance prediction of contra-rotating propellers, and small time step size can improve the level of accuracy for CRPs with the same blade number of forward and aft propellers, while a relatively large time step size is a better choice for CRPs with different blade numbers.
基金supported by the Program for Changjiang Scholars and Innovative Research Team in University of China under contract No,IRT0420the National Natural Science Foundation of China under contract No.50409015.
文摘The mechanism of local scour around submarine pipelines is studied numerically based on a renormalized group (RNG) turbulence model. To validate the numerical model, the equilibrium profiles of local scour for two cases are simulated and compared with the experimental data. It shows that the RNG turbulence model can give an appropriate prediction for the configuration of equilibrium scour hole, and it is applicable to this situation. The local scour mechanism around submarine pipelines including the flow structure, shear stress distribution and pressure field is then analyzed and compared with experiments. For further comparison and validation, especially for the flow structure, a numerical calculation employing the large eddy simulation (LES) is also conducted. The numerical results of RNG demonstrate that the critical factor governing the equilibrium profile is the seabed shear stress distribution in the case of bed load sediment transport, and the two-equation RNG turbulence model coupled with the law of wall is capable of giving a satisfying estimation for the bed shear stress. Moreover, the piping phenomena due to the great difference of pressure between the upstream and downstream parts of pipelines and the vortex structure around submarine pipelines are also simulated successfully, which are believed to be the important factor that lead to the onset of local scour.
文摘A steady three-dimensional fluid flow and solidification model was built based on CFD software by high-Reynolds-number and Lam-Bremhorst low-Reynolds-number k-ε model.During the simulation,the fixed-grid enthalpy-porosity technique was used to represent the solidification,and Darcy law was adopted to simulate the flow in mushy region.The prediction for steel flow and solidification was evaluated by the comparison of two turbulence models.It is found that both Lam-Bremhorst low-Reynolds-number and high-Reynolds-number k-ε models predict the same trend of the steel flow and temperature distribution.However,due to the effect of turbulent flow on heat transfer,the low-Reynolds-number turbulence model predicts longer penetration depth of molten steel in sub-mold region,less shell growth and higher shell surface temperature at the narrow face compared with standard k-ε model.
基金supported by the National Natural Science Foundation of China (10802026)
文摘Cavitation typically occurs when the fluid pressure is lower than the vapor pressure at a local thermodynamic state, and the flow is frequently unsteady and turbulent. To assess the state-of-the-art of computational capabilities for unsteady cavitating flows, different cavitation and turbulence model combinations are conducted. The selected cavitation models include several widely-used models including one based on phenomenological argument and the other utilizing interface dynamics. The k-e turbulence model with additional implementation of the filter function and density correction function are considered to reduce the eddy viscosity according to the computed turbulence length scale and local fluid density respectively. We have also blended these alternative cavitation and lustrate that the eddy viscosity turbulence treatments, to ilnear the closure region can significantly influence the capture of detached cavity. From the experimental validations regarding the force analysis, frequency, and the cavity visualization, no single model combination performs best in all aspects. Furthermore, the implications of parameters contained in different cavitation models are investigated. The phase change process is more pronounced around the detached cavity, which is better illustrated by the interfacial dynamics model. Our study provides insight to aid further modeling development.
基金This work was supported by the National Natural Science Foundation of China(91852108,11872230 and 92152301).
文摘Data-driven turbulence modeling studies have reached such a stage that the basic framework is settled,but several essential issues remain that strongly affect the performance.Two problems are studied in the current research:(1)the processing of the Reynolds stress tensor and(2)the coupling method between the machine learning model and flow solver.For the Reynolds stress processing issue,we perform the theoretical derivation to extend the relevant tensor arguments of Reynolds stress.Then,the tensor representation theorem is employed to give the complete irreducible invariants and integrity basis.An adaptive regularization term is employed to enhance the representation performance.For the coupling issue,an iterative coupling framework with consistent convergence is proposed and then applied to a canonical separated flow.The results have high consistency with the direct numerical simulation true values,which proves the validity of the current approach.
文摘The emerging push of the differentiable programming paradigm in scientific computing is conducive to training deep learning turbulence models using indirect observations.This paper demonstrates the viability of this approach and presents an end-to-end differentiable framework for training deep neural networks to learn eddy viscosity models from indirect observations derived from the velocity and pressure fields.The framework consists of a Reynolds-averaged Navier–Stokes(RANS)solver and a neuralnetwork-represented turbulence model,each accompanied by its derivative computations.For computing the sensitivities of the indirect observations to the Reynolds stress field,we use the continuous adjoint equations for the RANS equations,while the gradient of the neural network is obtained via its built-in automatic differentiation capability.We demonstrate the ability of this approach to learn the true underlying turbulence closure when one exists by training models using synthetic velocity data from linear and nonlinear closures.We also train a linear eddy viscosity model using synthetic velocity measurements from direct numerical simulations of the Navier–Stokes equations for which no true underlying linear closure exists.The trained deep-neural-network turbulence model showed predictive capability on similar flows.
基金National Natural Science Foundation of China(Nos.11675177,11875256)the Anhui Province Scientific and Technological Project(No.1604a0902145).
文摘Five turbulence models of Reynolds average Navier-Stokes(RANS),including the standard k-ω model,the RNG k-e model taking into account the low Reynolds number effect,the realizable k-ω model,the SST k-ω model,and the Reynolds stress model(RSM),are employed in the numerical simulations of direct current(DC)arc plasma torches in the range of arc current from 80 A to 240 A and air gas flow rate from 10 m^3 h^-1 to 50 m^3 h^-1.The calculated voltage,electric field intensity,and the heat loss in the arc chamber are compared with the experiments.The results indicate that the arc voltage,the electric field,and the heat loss in the arc chamber calculated by using the standard k-ω model,the RNG k-ωmodel taking into account the low Reynolds number effect,and the realizable k-ω model are much larger than those in the experiments.The RSM predicts relatively close results to the experiments,but fails in the trend of heat loss varying with the gas flow rate.The calculated results of the SST k-ω model are in the best agreement with the experiments,which may be attributed to the reasonable predictions of the turbulence as well as its distribution.
文摘Most current computations of trubulent flows with second-moment closure adopt the diffusion mo- dels which neglect the effect of pressure-velocity correlation.ln the present paper the importance of this correlation effect is elucidated the neglect of this effect accounts for some major defects in the wide application of the se- cond-moment closures.Through the relation between and ,established by Lumley,we propose here a new turbulence diffusion model which takes into consideration the pressure effect.Applications of this new model in the computation of shearless turbulence mixing layer and plane-and round-jet flows show that the spreading rate of these flows can be satisfactorily captured.
文摘The main theoretical turbulence models for application to hydrodynamic lubrication problems were briefly reviewed, and the course of their development and their fundamentals were explained. Predictions by these models on flow fields in turbulent Couette flows and shear-induced countercurrent flows were compared to existing measurements, and Zhang & Zhang's combined k-ε model was shown to have surpassingly satisfactory results. The method of application of this combined k-ε model to high speed journal bearings and annular seals was summarized, and the predicted results were shown to be satisfactory by comparisons with existing experiments of journal bearings and annular seals.
文摘A parameter estimation algorithm is introduced and used to determine the parameters in the standard k-epsilon two equation turbulence model (SKE). It can be found from the estimation results that although the parameter estimation method is an effective method to determine model parameters, it is. difficult to obtain a set of parameters for SKE to suit all kinds of separated flow and a modification of the turbulence model structure should be considered. So, a new nonlinear k-e two-equation model (NNKE) is put forward in this paper and the corresponding parameter estimation technique is applied to determine the model parameters. By implementing the NNKE to solve some engineering turbulent flows, it is shown that NNKE is more accurate and versatile than SKE. Thus, the success of NNKE implies that the parameter estimation technique may have a bright prospect in engineering turbulence model research.
文摘The research is motivated by the ongoing the electromagnetic continuous casting of molten metal. The revised k-ε model considering the effect of magnetic field application was derived. The specific model equations for the electromagnetic braking were used to calculate the velocity distribution in the continuous casting mold of steel. The results show that the revised k-ε model considering the effect of magnetic field application tends to suppress the production of turbulence and difference between the conventional and revised k-e model is small.
文摘A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and εequations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model.
基金the National Natural Science Foundation of China(Grant Nos.52071058 and 51939002)Liaoning Revitalization Talents Program(Grant No,XLYC1807208)the Special Funds for Promoting High Quality Development from Department of Natural Resources of Guangdong Province(Grant No.GDNRC[2020]015).
文摘This study numerically and experimentally investigates the effects of wave loads on a monopile-type offshore wind turbine placed on a 1:25 slope at different water depths as well as the effect of choosing different turbulence models on the efficiency of the numerical model.The numerical model adopts a two-phase flow by solving Unsteady Reynolds-Averaged Navier−Stokes(URANS)equations using the Volume Of Fluid(VOF)method and three differentk-ωturbulence models.Typical environmental conditions from the East China Sea are studied.The wave run-up and the wave loads applied on the monopile are investigated and compared with relevant experimental data as well as with mathematical predictions based on relevant theories.The numerical model is well validated against the experimental data at model scale.The use of different turbulence models results in different predictions on the wave height but less differences on the wave period.The baseline k-ωturbulence model and Shear-Stress Transport(SST)k-ωturbulence model exhibit better performance on the prediction of hydrodynamic load,at a model-scale water depth of 0.42 m,while the laminar model provides better results for large water depths.The SST turbulence model performs better in predicting wave run-up for water depth 0.42 m,while the laminar model and standard k-ωmodel perform better at water depth 0.52 m and 0.62 m,respectively.
文摘We investigate the role of extended intrinsic mean spin tensor introduced in this work for turbulence modelling in a non-inertial frame of reference. It is described by the Euclidean group of transformations and, in particular, its significance and importance in the approach of the algebraic Reynolds stress modelling, such as in a nonlinear K-ε model. To this end and for illustration of the effect of extended intrinsic spin tensor on turbulence modelling, we examine several recently developed nonlinear K-ε models and compare their performance in predicting the homogeneous turbulent shear flow in a rotating frame of reference with LES data. Our results and analysis indicate that, only if the deficiencies of these models and the like be well understood and properly corrected, may in the near future, more sophisticated nonlinear K-ε models be developed to better predict complex turbulent flows in a non-inertial frame of reference.
文摘This paper describes a newly developednon-isotropic multiple-scale turbulence model (MS/ASM) for complex flow calculations. This model focuses on the direct modeling of Reynolds stresses and utilizes split-spectrum concepts to model multiple-scale effects in turbulence. Validation studies on free shear flows, rotating flows and recirculating flows show that the current model performs significantly better than the single-scale k-e model. The present model is relatively inexpensive in terms of CPU time which makes in suitable for broad engineering flow applications.
文摘By analyzing the components of Reynolds stresses of implicit algebraic stress model(IASM) in this paper, that Reynolds stresses in buoyant turbulent flows were produced by both strain and buoyancy is considered. Consequently, a nonlinear anisotropy buoyant turbulence model was developed by applying linearity of equilibrium hypothesis to Reynolds stress transports. The model avoids numerical singularity and its reliability is verified by the comparisons between predictions and experimental data.