To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Un...To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Undersea Tunnel.To evaluate the discharging performance of short screw conveyor in different cases,the full-scale transient slurry-rock two-phase model for a short screw conveyor actively discharging rocks was established using computational fluid dynamics-discrete element method(CFD-DEM)coupling approach.In the fluid domain of coupling model,the sliding mesh technology was utilized to describe the rotations of the atmospheric composite cutterhead and the short screw conveyor.In the particle domain of coupling model,the dynamic particle factories were established to produce rock particles with the rotation of the cutterhead.And the accuracy and reliability of the CFD-DEM simulation results were validated via the field test and model test.Furthermore,a comprehensive parameter analysis was conducted to examine the effects of TBM operating parameters,the geometric design of screw conveyor and the size of rocks on the discharging performance of short screw conveyor.Accordingly,a reasonable rotational speed of screw conveyor was suggested and applied to Jiaozhou Bay Second Undersea Tunnel project.The findings in this paper could provide valuable references for addressing the excavation chamber clogging during ultra-large-diameter slurry TBM tunneling in hard rock for similar future.展开更多
Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the productio...Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the production efficiency.In this paper,these interactions were revealed by computational fluid dynamics–discrete element method(CFD-DEM)method.Based on simulation results,both liquid phase flow hydrodynamics and particle motion behavior have been analyzed,which gave the general information of the multiphase flow behavior inside the stirred tank reactor as to graphene production.By calculating the threshold at the beginning of graphite exfoliation process,the shear force from the slip velocity was determined as the active force.These results can support the optimization of the graphene production process.展开更多
To solve the problem of temperature rise caused by the high power density of high-speed permanent magnet synchronous traction motors,the temperature rise of various components in the motor is analyzed by coupling the ...To solve the problem of temperature rise caused by the high power density of high-speed permanent magnet synchronous traction motors,the temperature rise of various components in the motor is analyzed by coupling the equivalent thermal circuit method and computational fluid dynamics.Also,a cooling strategy is proposed to solve the problem of temperature rise,which is expected to prolong the service life of these devices.First,the theoretical bases of the approaches used to study heat transfer and fluid mechanics are discussed,then the fluid flow for the considered motor is analyzed,and the equivalent thermal circuit method is introduced for the calculation of the temperature rise.Finally,the stator,rotor loss,motor temperature rise,and the proposed cooling method are also explored through experiments.According to the results,the stator temperature at 50,000 r/min and 60,000 r/min at no-load operation is 68℃ and 76℃,respectively.By monitoring the temperature of the air outlets inside and outside the motor at different speeds,it is also found that the motor reaches a stable temperature rise after 65 min of operation.Coupling of the thermal circuit method and computational fluid dynamics is a strategy that can provide the average temperature rise of each component and can also comprehensively calculate the temperature of each local point.We conclude that a hybrid cooling strategy based on axial air cooling of the inner air duct of the motor and water cooling of the stator can meet the design requirements for the ventilation and cooling of this type of motors.展开更多
The present paper reviews the recent developments of a high⁃order⁃spectral method(HOS)and the combination with computational fluid dynamics(CFD)method for wave⁃structure interactions.As the numerical simulations of wa...The present paper reviews the recent developments of a high⁃order⁃spectral method(HOS)and the combination with computational fluid dynamics(CFD)method for wave⁃structure interactions.As the numerical simulations of wave⁃structure interaction require efficiency and accuracy,as well as the ability in calculating in open sea states,the HOS method has its strength in both generating extreme waves in open seas and fast convergence in simulations,while computational fluid dynamics(CFD)method has its advantages in simulating violent wave⁃structure interactions.This paper provides the new thoughts for fast and accurate simulations,as well as the future work on innovations in fine fluid field of numerical simulations.展开更多
To find out and improve the flow characteristics inside the intake system of cylinder head,the application of computational fluid dynamics(CFD)in the evaluation and optimization of the reconstructed intake system base...To find out and improve the flow characteristics inside the intake system of cylinder head,the application of computational fluid dynamics(CFD)in the evaluation and optimization of the reconstructed intake system based on slicing reverse method was proposed.The flow characteristics were found out through CFD,and the velocity vector field,pressure field and turbulent kinetic energy field for different valve lifts were discussed,which were in good agreement with experimental data,and the quality of reconstruction was evaluated.In order to improve its flow characteristic,an optimization plan was proposed.The results show that the flow characteristics after optimization are obviously improved.The results can provide a reference for the design and optimization of the intake system of cylinder head.展开更多
To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters of nasal structure, three-dimensional, anatomically accurate representations of 30 adult nasal cavity models...To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters of nasal structure, three-dimensional, anatomically accurate representations of 30 adult nasal cavity models were recons- tructed based on processed tomography images collected from normal people. The airflow fields in nasal cavities were simulated by fluid dynamics with finite element software ANSYS. The results showed that the difference of human nasal cavity structure led to different airflow distribution in the nasal cavities and variation of the main airstream passing through the common nasal meatus. The nasal resistance in the regions of nasal valve and nasal vestibule accounted for more than half of the overall resistance. The characteristic model of nasal cavity was extracted on the basis of characteristic points and dimensions deduced from the original models. It showed that either the geometric structure or the airflow field of the two kinds of models was similar. The characteristic dimensions were the characteristic parameters of nasal cavity that could properly represent the original model in model studies on nasal cavity.展开更多
A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using A USMPW scheme...A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using A USMPW scheme in conjunction with the third-order MUSCL scheme with Van Leer limiter. The present method was applied to solve the multidimensional compressible Navier-Stokes equations in curvilinear coordinates. Characteristic boundary conditions based on the eigensystem of the preconditioned equations were employed. In order to examine the performance of present method, driven-cavity flow at various Reynolds numbers and viscous flow through a convergent-divergent nozzle at supersonic were selected to rest this method. The computed results were compared with the experimental data or the other numerical results available in literature and good agreements between them are obtained. The results show that the present method is accurate, self-adaptive and stable for a wide range of flow conditions from low speed to supersonic flows.展开更多
The roll motions are influenced by significant viscous effects such as the flow separation.The 3D simulations of free decay roll motions for the ship model DTMB 5512 are carried out by Reynold averaged NavierStokes(RA...The roll motions are influenced by significant viscous effects such as the flow separation.The 3D simulations of free decay roll motions for the ship model DTMB 5512 are carried out by Reynold averaged NavierStokes(RANS) method based on the dynamic mesh technique.A new moving mesh technique is adopted and discussed in details for the present simulations.The purpose of the research is to obtain accurate numerical prediction for roll motions with their respective numerical/modeling errors and uncertainties.Errors and uncertainties are estimated by performing the modern verification and validation(V&V) procedures.Simulation results for the free-floating surface combatant are used to calculate the linear,nonlinear damping coefficients and resonant frequencies including a wide range of forward speed.The present work can provide a useful reference to calculate roll damping by computational fluid dynamics(CFD) method and simulate a general ship motions in waves.展开更多
A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-unifor...A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-uniform stress field and to analyse its phase drift phenomenon. In the model, the drift-inhibition angle and the expansion-inhibition angle are also deduced and used as evaluating indexes to describe the drifting trend of different ingredients among the mixed fluids. For solving above two indexes of the model, a new calculation method is developed and used to compute the phase distributions of multiphase fluid at peak stress and gradient area stress, respectively. As an example, the flow process of grease in a pipe is analysed by simulation method and used to verify the validity of the model.展开更多
While finite volume methodologies (FVM) have predominated in fluid flow computations, many flow problems, including groundwater models, would benefit from the use of boundary methods, such as the Complex Variable Boun...While finite volume methodologies (FVM) have predominated in fluid flow computations, many flow problems, including groundwater models, would benefit from the use of boundary methods, such as the Complex Variable Boundary Element Method (CVBEM). However, to date, there has been no reporting of a comparison of computational results between the FVM and the CVBEM in the assessment of flow field characteristics. In this work, the CVBEM is used to develop a flow field vector outcome of ideal fluid flow in a 90-degree bend which is then compared to the computational results from a finite volume model of the same situation. The focus of the modelling comparison in the current work is flow field trajectory vectors of the fluid flow, with respect to vector magnitude and direction. Such a comparison is necessary to validate the development of flow field vectors from the CVBEM and is of interest to many engineering flow problems, specifically groundwater modelling. Comparison of the CVBEM and FVM flow field trajectory vectors for the target problem of ideal flow in a 90-degree bend shows good agreement between the considered methodologies.展开更多
This paper presents an adapted stabilisation method for the equal-order mixed scheme of finite elements on convex polygonal meshes to analyse the high velocity and pressure gradient of incompressible fluid flows that ...This paper presents an adapted stabilisation method for the equal-order mixed scheme of finite elements on convex polygonal meshes to analyse the high velocity and pressure gradient of incompressible fluid flows that are governed by Stokes equations system.This technique is constructed by a local pressure projection which is extremely simple,yet effective,to eliminate the poor or even non-convergence as well as the instability of equal-order mixed polygonal technique.In this research,some numerical examples of incompressible Stokes fluid flow that is coded and programmed by MATLAB will be presented to examine the effectiveness of the proposed stabilised method.展开更多
In this paper,the author recounts his forty-year plus struggle to find a sound basis for understanding the computational fluid dynamics of compressible flow.
This paper presents and proves the mixed compatible finite element variationalprinciples in dynamics of viscous barotropic fluids. When the principles are proved, itis found that the compatibility conditions of stress...This paper presents and proves the mixed compatible finite element variationalprinciples in dynamics of viscous barotropic fluids. When the principles are proved, itis found that the compatibility conditions of stress can be naturally satisfied. The gene-rallzed variational principles with mixed hybrid incompatible finite elements are alsopresented and proved, and they can reduce the computation of incompatible elements indynamics of viscous barotropic flows.展开更多
This work presents a new application for the Hierarchical Function Expansion Method for the solution of the Navier-Stokes equations for compressible fluids in two dimensions and in high velocity. This method is based ...This work presents a new application for the Hierarchical Function Expansion Method for the solution of the Navier-Stokes equations for compressible fluids in two dimensions and in high velocity. This method is based on the finite elements method using the Petrov-Galerkin formulation, know as SUPG (Streamline Upwind Petrov-Galerkin), applied with the expansion of the variables into hierarchical functions. To test and validate the numerical method proposed as well as the computational program developed simulations are performed for some cases whose theoretical solutions are known. These cases are the following: continuity test, stability and convergence test, temperature step problem, and several oblique shocks. The objective of the last cases is basically to verify the capture of the shock wave by the method developed. The results obtained in the simulations with the proposed method were good both qualitatively and quantitatively when compared with the theoretical solutions. This allows concluding that the objectives of this work are reached.展开更多
Three recent breakthroughs due to AI in arts and science serve as motivation:An award winning digital image,protein folding,fast matrix multiplication.Many recent developments in artificial neural networks,particularl...Three recent breakthroughs due to AI in arts and science serve as motivation:An award winning digital image,protein folding,fast matrix multiplication.Many recent developments in artificial neural networks,particularly deep learning(DL),applied and relevant to computational mechanics(solid,fluids,finite-element technology)are reviewed in detail.Both hybrid and pure machine learning(ML)methods are discussed.Hybrid methods combine traditional PDE discretizations with ML methods either(1)to help model complex nonlinear constitutive relations,(2)to nonlinearly reduce the model order for efficient simulation(turbulence),or(3)to accelerate the simulation by predicting certain components in the traditional integration methods.Here,methods(1)and(2)relied on Long-Short-Term Memory(LSTM)architecture,with method(3)relying on convolutional neural networks.Pure ML methods to solve(nonlinear)PDEs are represented by Physics-Informed Neural network(PINN)methods,which could be combined with attention mechanism to address discontinuous solutions.Both LSTM and attention architectures,together with modern and generalized classic optimizers to include stochasticity for DL networks,are extensively reviewed.Kernel machines,including Gaussian processes,are provided to sufficient depth for more advanced works such as shallow networks with infinite width.Not only addressing experts,readers are assumed familiar with computational mechanics,but not with DL,whose concepts and applications are built up from the basics,aiming at bringing first-time learners quickly to the forefront of research.History and limitations of AI are recounted and discussed,with particular attention at pointing out misstatements or misconceptions of the classics,even in well-known references.Positioning and pointing control of a large-deformable beam is given as an example.展开更多
A simulation was carried out on an unsteady flow of a viscous, incompressible and electrically conducting fluid past an infinite vertical porous plate. A generic computer program using the Galerkin finite element meth...A simulation was carried out on an unsteady flow of a viscous, incompressible and electrically conducting fluid past an infinite vertical porous plate. A generic computer program using the Galerkin finite element method is employed to solve the coupled non-linear differential equations for velocity and temperature fields. The diffusion equation, the energy equation, the momentum equations and other relevant parameters are transformed into interpretable postfix codes. Numerical calculations are carried out on the flow fields both in the presence of cooling and heating of the plate by free convection currents. The effects of the dimensionless parameters, namely, the Prandtl number, the Eckert number, the modified Grashof number, the Schmidt number and the time on the temperature and velocity distributions are discussed.展开更多
An advanced Actuator Surface Method(ASM)coupled with Computational Fluid Dynamics(CFD)is developed and applied to the complex unsteady aerodynamic simulation of helicopter.By introducing an improved three-dimensional ...An advanced Actuator Surface Method(ASM)coupled with Computational Fluid Dynamics(CFD)is developed and applied to the complex unsteady aerodynamic simulation of helicopter.By introducing an improved three-dimensional anisotropic Gaussian kernel,this method effectively addresses the severe aerodynamic load fluctuations commonly associated with traditional Virtual Blade Method(VBM)due to turbulent flow around blade elements.To manage the issues of regional shape and grid cell quantity variations caused by virtual blade sweeping,a universal hybrid grid generation strategy is established without body-fitted and disk interpolation grids,which enhances the computational stability at both blade elements and blade edges.Aerodynamic numerical simulations of helicopter are performed using this method,focusing on rotor/fuselage interaction dominated by rotor wake motion and fuselage blockage effects,Blade-Vortex Interaction(BVI)induced by tip vortices,and maneuvering flights involving collective pitch ramp increases.The results indicate that the advanced ASM demonstrates reliability and robustness in the simulation of complex unsteady flow fields around helicopter.Under similar computational accuracy,the advanced ASM improves computational efficiency by nearly 40 times compared to the oversetgrid-based full Blade-Resolved(B-R)method,and by 6 times compared to the VBM.It shows significant advantages when applied to complex full-aircraft interaction and maneuvering flight conditions that require substantial computational resources.展开更多
A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underex...A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.展开更多
The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S）...The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S） equation. This paper proposes a novel immersed boundary-lattice Boltzmann method （IB-LBM） based on the feedback law. The method uses the immersed boundary concept in the LBM framework to capture the coupling between a body with complex geometry and a uniform fluid, Then, the flows around a stationary circular cylinder and two circular cylinders in a side by side arrangement are simulated by using the method. Results are agreed well with the benchmark data, so, the capability of the method for complex geometry is demonstrated. Different from the conventional IB-LBM, which uses the Hook＇s law or the direct forcing method to compute the interae- tion force, the method uses the feedback law--the feedback of velocity field and displacement information to calculate the force, thus ensuring the method has advantages of easy implementation and full parallelism.展开更多
基金supported by the Fundamental Research Funds for the Central Universities(Grant No.2023YJS053)the National Natural Science Foundation of China(Grant No.52278386).
文摘To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Undersea Tunnel.To evaluate the discharging performance of short screw conveyor in different cases,the full-scale transient slurry-rock two-phase model for a short screw conveyor actively discharging rocks was established using computational fluid dynamics-discrete element method(CFD-DEM)coupling approach.In the fluid domain of coupling model,the sliding mesh technology was utilized to describe the rotations of the atmospheric composite cutterhead and the short screw conveyor.In the particle domain of coupling model,the dynamic particle factories were established to produce rock particles with the rotation of the cutterhead.And the accuracy and reliability of the CFD-DEM simulation results were validated via the field test and model test.Furthermore,a comprehensive parameter analysis was conducted to examine the effects of TBM operating parameters,the geometric design of screw conveyor and the size of rocks on the discharging performance of short screw conveyor.Accordingly,a reasonable rotational speed of screw conveyor was suggested and applied to Jiaozhou Bay Second Undersea Tunnel project.The findings in this paper could provide valuable references for addressing the excavation chamber clogging during ultra-large-diameter slurry TBM tunneling in hard rock for similar future.
基金National Natural Science Foundation of China(U2004176,22008055)Technology Research Project of Henan Province(232102240034)are gratefully acknowledged.
文摘Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the production efficiency.In this paper,these interactions were revealed by computational fluid dynamics–discrete element method(CFD-DEM)method.Based on simulation results,both liquid phase flow hydrodynamics and particle motion behavior have been analyzed,which gave the general information of the multiphase flow behavior inside the stirred tank reactor as to graphene production.By calculating the threshold at the beginning of graphite exfoliation process,the shear force from the slip velocity was determined as the active force.These results can support the optimization of the graphene production process.
文摘To solve the problem of temperature rise caused by the high power density of high-speed permanent magnet synchronous traction motors,the temperature rise of various components in the motor is analyzed by coupling the equivalent thermal circuit method and computational fluid dynamics.Also,a cooling strategy is proposed to solve the problem of temperature rise,which is expected to prolong the service life of these devices.First,the theoretical bases of the approaches used to study heat transfer and fluid mechanics are discussed,then the fluid flow for the considered motor is analyzed,and the equivalent thermal circuit method is introduced for the calculation of the temperature rise.Finally,the stator,rotor loss,motor temperature rise,and the proposed cooling method are also explored through experiments.According to the results,the stator temperature at 50,000 r/min and 60,000 r/min at no-load operation is 68℃ and 76℃,respectively.By monitoring the temperature of the air outlets inside and outside the motor at different speeds,it is also found that the motor reaches a stable temperature rise after 65 min of operation.Coupling of the thermal circuit method and computational fluid dynamics is a strategy that can provide the average temperature rise of each component and can also comprehensively calculate the temperature of each local point.We conclude that a hybrid cooling strategy based on axial air cooling of the inner air duct of the motor and water cooling of the stator can meet the design requirements for the ventilation and cooling of this type of motors.
基金National Natural Science Foundation of China(Grant No.51879159)the National Key Research and Development Program of China(Grant Nos.2019YFB1704200 and 2019YFC0312400)+2 种基金the Chang Jiang Scholars Program(Grant No.T2014099)the Shanghai Excellent Academic Leaders Program(Grant No.17XD1402300)the Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China(Grant No.2016-23/09).
文摘The present paper reviews the recent developments of a high⁃order⁃spectral method(HOS)and the combination with computational fluid dynamics(CFD)method for wave⁃structure interactions.As the numerical simulations of wave⁃structure interaction require efficiency and accuracy,as well as the ability in calculating in open sea states,the HOS method has its strength in both generating extreme waves in open seas and fast convergence in simulations,while computational fluid dynamics(CFD)method has its advantages in simulating violent wave⁃structure interactions.This paper provides the new thoughts for fast and accurate simulations,as well as the future work on innovations in fine fluid field of numerical simulations.
基金“Strategic Cooperation of Science and Technology between Nanchong City and Southwest Petroleum University 2018” Special Fund Project,China(Nos.18SXHZ0030,18SXHZ0054)
文摘To find out and improve the flow characteristics inside the intake system of cylinder head,the application of computational fluid dynamics(CFD)in the evaluation and optimization of the reconstructed intake system based on slicing reverse method was proposed.The flow characteristics were found out through CFD,and the velocity vector field,pressure field and turbulent kinetic energy field for different valve lifts were discussed,which were in good agreement with experimental data,and the quality of reconstruction was evaluated.In order to improve its flow characteristic,an optimization plan was proposed.The results show that the flow characteristics after optimization are obviously improved.The results can provide a reference for the design and optimization of the intake system of cylinder head.
基金the National Natural Science Foundation of China (1047202510672036)the Natural Science Foundation of Liaoning Province,China (20032109)
文摘To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters of nasal structure, three-dimensional, anatomically accurate representations of 30 adult nasal cavity models were recons- tructed based on processed tomography images collected from normal people. The airflow fields in nasal cavities were simulated by fluid dynamics with finite element software ANSYS. The results showed that the difference of human nasal cavity structure led to different airflow distribution in the nasal cavities and variation of the main airstream passing through the common nasal meatus. The nasal resistance in the regions of nasal valve and nasal vestibule accounted for more than half of the overall resistance. The characteristic model of nasal cavity was extracted on the basis of characteristic points and dimensions deduced from the original models. It showed that either the geometric structure or the airflow field of the two kinds of models was similar. The characteristic dimensions were the characteristic parameters of nasal cavity that could properly represent the original model in model studies on nasal cavity.
文摘A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using A USMPW scheme in conjunction with the third-order MUSCL scheme with Van Leer limiter. The present method was applied to solve the multidimensional compressible Navier-Stokes equations in curvilinear coordinates. Characteristic boundary conditions based on the eigensystem of the preconditioned equations were employed. In order to examine the performance of present method, driven-cavity flow at various Reynolds numbers and viscous flow through a convergent-divergent nozzle at supersonic were selected to rest this method. The computed results were compared with the experimental data or the other numerical results available in literature and good agreements between them are obtained. The results show that the present method is accurate, self-adaptive and stable for a wide range of flow conditions from low speed to supersonic flows.
基金the National Natural Science Foundation of China(No.51579147)
文摘The roll motions are influenced by significant viscous effects such as the flow separation.The 3D simulations of free decay roll motions for the ship model DTMB 5512 are carried out by Reynold averaged NavierStokes(RANS) method based on the dynamic mesh technique.A new moving mesh technique is adopted and discussed in details for the present simulations.The purpose of the research is to obtain accurate numerical prediction for roll motions with their respective numerical/modeling errors and uncertainties.Errors and uncertainties are estimated by performing the modern verification and validation(V&V) procedures.Simulation results for the free-floating surface combatant are used to calculate the linear,nonlinear damping coefficients and resonant frequencies including a wide range of forward speed.The present work can provide a useful reference to calculate roll damping by computational fluid dynamics(CFD) method and simulate a general ship motions in waves.
基金Project supported by the National Natural Science Foundation of China (Grant No. 51075311)
文摘A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-uniform stress field and to analyse its phase drift phenomenon. In the model, the drift-inhibition angle and the expansion-inhibition angle are also deduced and used as evaluating indexes to describe the drifting trend of different ingredients among the mixed fluids. For solving above two indexes of the model, a new calculation method is developed and used to compute the phase distributions of multiphase fluid at peak stress and gradient area stress, respectively. As an example, the flow process of grease in a pipe is analysed by simulation method and used to verify the validity of the model.
文摘While finite volume methodologies (FVM) have predominated in fluid flow computations, many flow problems, including groundwater models, would benefit from the use of boundary methods, such as the Complex Variable Boundary Element Method (CVBEM). However, to date, there has been no reporting of a comparison of computational results between the FVM and the CVBEM in the assessment of flow field characteristics. In this work, the CVBEM is used to develop a flow field vector outcome of ideal fluid flow in a 90-degree bend which is then compared to the computational results from a finite volume model of the same situation. The focus of the modelling comparison in the current work is flow field trajectory vectors of the fluid flow, with respect to vector magnitude and direction. Such a comparison is necessary to validate the development of flow field vectors from the CVBEM and is of interest to many engineering flow problems, specifically groundwater modelling. Comparison of the CVBEM and FVM flow field trajectory vectors for the target problem of ideal flow in a 90-degree bend shows good agreement between the considered methodologies.
基金The authors would like to present our gratitude to the Flemish Government financially supporting for the VLIR-OUS TEAM Project,VN2017TEA454A103‘An innovative solution to protect Vietnamese coastal riverbanks from floods and erosion’.
文摘This paper presents an adapted stabilisation method for the equal-order mixed scheme of finite elements on convex polygonal meshes to analyse the high velocity and pressure gradient of incompressible fluid flows that are governed by Stokes equations system.This technique is constructed by a local pressure projection which is extremely simple,yet effective,to eliminate the poor or even non-convergence as well as the instability of equal-order mixed polygonal technique.In this research,some numerical examples of incompressible Stokes fluid flow that is coded and programmed by MATLAB will be presented to examine the effectiveness of the proposed stabilised method.
文摘In this paper,the author recounts his forty-year plus struggle to find a sound basis for understanding the computational fluid dynamics of compressible flow.
文摘This paper presents and proves the mixed compatible finite element variationalprinciples in dynamics of viscous barotropic fluids. When the principles are proved, itis found that the compatibility conditions of stress can be naturally satisfied. The gene-rallzed variational principles with mixed hybrid incompatible finite elements are alsopresented and proved, and they can reduce the computation of incompatible elements indynamics of viscous barotropic flows.
文摘This work presents a new application for the Hierarchical Function Expansion Method for the solution of the Navier-Stokes equations for compressible fluids in two dimensions and in high velocity. This method is based on the finite elements method using the Petrov-Galerkin formulation, know as SUPG (Streamline Upwind Petrov-Galerkin), applied with the expansion of the variables into hierarchical functions. To test and validate the numerical method proposed as well as the computational program developed simulations are performed for some cases whose theoretical solutions are known. These cases are the following: continuity test, stability and convergence test, temperature step problem, and several oblique shocks. The objective of the last cases is basically to verify the capture of the shock wave by the method developed. The results obtained in the simulations with the proposed method were good both qualitatively and quantitatively when compared with the theoretical solutions. This allows concluding that the objectives of this work are reached.
文摘Three recent breakthroughs due to AI in arts and science serve as motivation:An award winning digital image,protein folding,fast matrix multiplication.Many recent developments in artificial neural networks,particularly deep learning(DL),applied and relevant to computational mechanics(solid,fluids,finite-element technology)are reviewed in detail.Both hybrid and pure machine learning(ML)methods are discussed.Hybrid methods combine traditional PDE discretizations with ML methods either(1)to help model complex nonlinear constitutive relations,(2)to nonlinearly reduce the model order for efficient simulation(turbulence),or(3)to accelerate the simulation by predicting certain components in the traditional integration methods.Here,methods(1)and(2)relied on Long-Short-Term Memory(LSTM)architecture,with method(3)relying on convolutional neural networks.Pure ML methods to solve(nonlinear)PDEs are represented by Physics-Informed Neural network(PINN)methods,which could be combined with attention mechanism to address discontinuous solutions.Both LSTM and attention architectures,together with modern and generalized classic optimizers to include stochasticity for DL networks,are extensively reviewed.Kernel machines,including Gaussian processes,are provided to sufficient depth for more advanced works such as shallow networks with infinite width.Not only addressing experts,readers are assumed familiar with computational mechanics,but not with DL,whose concepts and applications are built up from the basics,aiming at bringing first-time learners quickly to the forefront of research.History and limitations of AI are recounted and discussed,with particular attention at pointing out misstatements or misconceptions of the classics,even in well-known references.Positioning and pointing control of a large-deformable beam is given as an example.
文摘A simulation was carried out on an unsteady flow of a viscous, incompressible and electrically conducting fluid past an infinite vertical porous plate. A generic computer program using the Galerkin finite element method is employed to solve the coupled non-linear differential equations for velocity and temperature fields. The diffusion equation, the energy equation, the momentum equations and other relevant parameters are transformed into interpretable postfix codes. Numerical calculations are carried out on the flow fields both in the presence of cooling and heating of the plate by free convection currents. The effects of the dimensionless parameters, namely, the Prandtl number, the Eckert number, the modified Grashof number, the Schmidt number and the time on the temperature and velocity distributions are discussed.
基金co-supported by the Foundation of the State Key Laboratory of Aerodynamics(No.RAL202203)the National Key Laboratory Foundation(No.6142202202)the China Postdoctoral Science Foundation(No.2024M754133)。
文摘An advanced Actuator Surface Method(ASM)coupled with Computational Fluid Dynamics(CFD)is developed and applied to the complex unsteady aerodynamic simulation of helicopter.By introducing an improved three-dimensional anisotropic Gaussian kernel,this method effectively addresses the severe aerodynamic load fluctuations commonly associated with traditional Virtual Blade Method(VBM)due to turbulent flow around blade elements.To manage the issues of regional shape and grid cell quantity variations caused by virtual blade sweeping,a universal hybrid grid generation strategy is established without body-fitted and disk interpolation grids,which enhances the computational stability at both blade elements and blade edges.Aerodynamic numerical simulations of helicopter are performed using this method,focusing on rotor/fuselage interaction dominated by rotor wake motion and fuselage blockage effects,Blade-Vortex Interaction(BVI)induced by tip vortices,and maneuvering flights involving collective pitch ramp increases.The results indicate that the advanced ASM demonstrates reliability and robustness in the simulation of complex unsteady flow fields around helicopter.Under similar computational accuracy,the advanced ASM improves computational efficiency by nearly 40 times compared to the oversetgrid-based full Blade-Resolved(B-R)method,and by 6 times compared to the VBM.It shows significant advantages when applied to complex full-aircraft interaction and maneuvering flight conditions that require substantial computational resources.
文摘A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.
基金Supported by the Aeronautical Science Foundation of China(20111453012)the National Defense Pre-Research Foundation of China(9140A13040111HK0329)~~
文摘The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S） equation. This paper proposes a novel immersed boundary-lattice Boltzmann method （IB-LBM） based on the feedback law. The method uses the immersed boundary concept in the LBM framework to capture the coupling between a body with complex geometry and a uniform fluid, Then, the flows around a stationary circular cylinder and two circular cylinders in a side by side arrangement are simulated by using the method. Results are agreed well with the benchmark data, so, the capability of the method for complex geometry is demonstrated. Different from the conventional IB-LBM, which uses the Hook＇s law or the direct forcing method to compute the interae- tion force, the method uses the feedback law--the feedback of velocity field and displacement information to calculate the force, thus ensuring the method has advantages of easy implementation and full parallelism.
