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.展开更多
This paper aims to investigate a dam break in a channel with a bend in the presence of several obstacles.To accurately determine the flood zones,it is necessary to take into account many factors such as terrain,reserv...This paper aims to investigate a dam break in a channel with a bend in the presence of several obstacles.To accurately determine the flood zones,it is necessary to take into account many factors such as terrain,reservoir volume.Numerical modeling was used to determine the flood zone.Numerical modeling based on the Navier-Stokes equations with a turbulent k-epsilon RNG model,the Volume of Fluid(VOF)method and the PISO algorithm were used to analyze the flow in a bend channel at an angle of 10 with the obstacles.To verify the numerical model,a test on dam break in the 450 channel was conducted.The simulation results were compared with experimental data and with the numerical data of existing data.Having been convinced of the correctness of the mathematical model,the authors carried out a numerical simulation of the main problem in three versions:without barriers,with one obstacle,with two obstacles.According to the obtained numerical results,it can be noted that irregular landforms held the flow,a decrease in water level and a slower time for water emergence could be seen.Thus,the water flow without an obstacle,with one obstacle and with two obstacles showed 4.2 s,4.4 s and 4.6 s of the time of water appearance,respectively.This time shift can give a certain advantage when conducting various events to evacuate people.展开更多
This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeabi...This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.展开更多
The resistivity difference between oil and gas layers and the water layers in low contrast tight sandstone reservoirs is subtle. Fluid identification and saturation calculation based on conventional logging methods ar...The resistivity difference between oil and gas layers and the water layers in low contrast tight sandstone reservoirs is subtle. Fluid identification and saturation calculation based on conventional logging methods are facing challenges in such reservoirs. In this paper, a new method is proposed for fluid identification and saturation calculation in low contrast tight sandstone reservoirs. First, a model for calculating apparent formation water resistivity is constructed, which takes into account the influence of shale on the resistivity calculation and avoids apparent formation water resistivity abnormal values.Based on the distribution of the apparent formation water resistivity obtained by the new model, the water spectrum is determined for fluid identification in low contrast tight sandstone reservoirs.Following this, according to the average, standard deviation, and endpoints of the water spectrum, a new four-parameter model for calculating reservoir oil and gas saturation is built. The methods proposed in this paper are applied to the low contrast tight sandstone reservoirs in the Q4 formation of the X53 block and X70 block in the south of Songliao Basin, China. The results show that the water spectrum method can effectively distinguish oil-water layers and water layers in the study area. The standard deviation of the water spectrum in the oil-water layer is generally greater than that in the water layer. The new four-parameter model yields more accurate oil and gas saturation. These findings verify the effectiveness of the proposed methods.展开更多
A stencil-like volume of fluid (VOF) method is proposed for tracking free interface. A stencil on a grid cell is worked out according to the normal direction of the interface, in which only three interface positions...A stencil-like volume of fluid (VOF) method is proposed for tracking free interface. A stencil on a grid cell is worked out according to the normal direction of the interface, in which only three interface positions are possible in 2D cases, and the interface can be reconstructed by only requiring the known local volume fraction information. On the other hand, the fluid-occupying-length is defined on each side of the stencil, through which a unified fluid-occupying volume model and a unified algorithm can be obtained to solve the interface advection equation. The method is suitable for the arbitrary geometry of the grid cell, and is extendible to 3D cases. Typical numerical examples show that the current method can give "sharp" results for tracking free interface.展开更多
Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise re...Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.展开更多
Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data imm...Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data immersion method(BDIM)for compressible flows,which were solved using a fully self-programmed Fortran code.Results show that the predicted thrust spectrum aligns closely with the experimental measurements.In addition,it captures the thrust dipole directivity of the noise around the rotating propeller due to random pressure pulsations on the blades,as well as the flow structures simultaneously.Furthermore,the differences in the statistical characteristics,flow structures,and low-frequency broadband thrust spectra due to different turbulence levels were investigated.This analysis indicates that the interaction between the upstream,which is characterized by a lower turbulence level and a higher turbulent length of scale,and the rotating propeller results in a lower amplitude in force spectra and a slight increase in the scale of tip vortices.展开更多
Characteristics of heat transfer and flow of Newtonian and non-Newtonian fluids through porous walls and in porous media are studied due to their wide range of applications including geothermal reservoirs,heat exchang...Characteristics of heat transfer and flow of Newtonian and non-Newtonian fluids through porous walls and in porous media are studied due to their wide range of applications including geothermal reservoirs,heat exchangers,marine propulsion,and aerodynamics.The current study investigates the characteristics of heat transport in a reactive third-grade fluid,moving through permeable parallel plates,with uniform suction/injection velocity.The two permeable,parallel plates are maintained at the same,constant temperature.After being transformed into its dimensionless equivalent,governing equations are solved by employing the Least Squares Method(LSM).The LSM results are further validated with numerical solutions for temperature and velocity.The impact of cross-flow Reynolds number,Peclet number,heat generation parameter,non-Newtonian parameter,and Brinkman number on entropy generation,velocity,temperature,and Bejan number are investigated.Theresults indicate that temperature distribution is significantly influenced by the third-grade fluid parameter.The maximum temperature drops from almost 0.12 to 0.10 as the third-grade fluid parameter increases from0.05 to 0.4.When the cross-flow Reynolds number is raised from 0.05 to 3,the maximum temperature drops from 0.12 to around 0.09.Temperature is strongly influenced by the heat generation parameter.A greater understanding of the thermal characteristics necessary for the design of a variety of systems,such as heat exchangers,marine propulsion,aerodynamic systems,etc.,may be gained from the findings of the current study.展开更多
This study rigorously examines the interplay between viscous dissipation,magnetic effects,and thermal radiation on the flow behavior of a non-Newtonian Carreau squeezed fluid passing by a sensor surface within a micro...This study rigorously examines the interplay between viscous dissipation,magnetic effects,and thermal radiation on the flow behavior of a non-Newtonian Carreau squeezed fluid passing by a sensor surface within a micro cantilever channel,aiming to deepen our understanding of heat transport processes in complex fluid dynamics scenarios.The primary objective is to elucidate how physical operational parameters influence both the velocity of fluid flow and its temperature distribution,utilizing a comprehensive numerical approach.Employing a combination of mathematical modeling techniques,including similarity transformation,this investigation transforms complex partial differential equations into more manageable ordinary ones,subsequently solving them using the homotopy perturbation method.By analyzing the obtained solutions and presenting them graphically,alongside detailed analysis,the study sheds light on the pivotal role of significant parameters in shaping fluid movement and energy distribution.Noteworthy observations reveal a substantial increase in fluid velocity with escalating magnetic parameters,while conversely,a contrasting trend emerges in the temperature distribution,highlighting the intricate relationship between magnetic effects,flow dynamics,and thermal behavior in non-Newtonian fluids.Further,the suction velocity enhance both the local skin friction and Nusselt numbers,whereas theWeissenberg number reduces them,opposite to the effect of the power-law index.展开更多
The Richtmyer–Meshkov(RM)instability plays an important role in various natural and engineering fields such as inertial confinement fusion.In this study,the effect of relaxation time on the RM instability under resho...The Richtmyer–Meshkov(RM)instability plays an important role in various natural and engineering fields such as inertial confinement fusion.In this study,the effect of relaxation time on the RM instability under reshock impact is investigated using a two-component discrete Boltzmann method.The hydrodynamic and thermodynamic characteristics of the fluid system are comprehensively analyzed from the perspectives of the density gradient,vorticity,kinetic energy,mixing degree,mixing width and non-equilibrium intensity.Simulation results indicate that for longer relaxation time,the diffusion and dissipation are enhanced,the physical gradients decrease,and the growth of the interface is suppressed.Furthermore,the non-equilibrium manifestations show complex patterns,driven by the competitive physical mechanisms of the diffusion,dissipation,shock wave,rarefaction wave,transverse wave and fluid instabilities.These findings provide valuable insight into the fundamental mechanism of compressible fluid flows.展开更多
We are pleased to introduce this special thematic section on Nanofluid Mechanics and Heat Transfer that is being included in Acta Mechanica Sinica(AMS).This thematic issue consists of 6 papers selected from papers tha...We are pleased to introduce this special thematic section on Nanofluid Mechanics and Heat Transfer that is being included in Acta Mechanica Sinica(AMS).This thematic issue consists of 6 papers selected from papers that were presented at the 18th International Symposium on Numer-ical Analysis of Fluid Flows,Heat and Mass Transfer-Nu-merical Fluids 2023,held in Heraklion,Crete Greece,11-17 September 2023,and invited through a general call.The symposium covers various subjects:from new numerical methods and fundamental research until engineering appli-cations,and it is a part of the International Conference of Numerical Analysis and Applied Mathematics(ICNAAM),held annually.展开更多
Physics-informed neural networks(PINNs)have shown considerable promise for performing numerical simulations in fluid mechanics.They provide mesh-free,end-to-end approaches by embedding physical laws into their loss fu...Physics-informed neural networks(PINNs)have shown considerable promise for performing numerical simulations in fluid mechanics.They provide mesh-free,end-to-end approaches by embedding physical laws into their loss functions.However,when addressing complex flow problems,PINNs still face some challenges such as activation saturation and vanishing gradients in deep network training,leading to slow convergence and insufficient prediction accuracy.We present physics-informed neural networks incorporating lattice Boltzmann method optimized by tanh robust weight initialization(T-PINN-LBM)to address these challenges.This approach fuses the mesoscopic lattice Boltzmann model with the automatic differentiation framework of PINNs.It also implements a tanh robust weight initialization method derived from fixed point analysis.This model effectively mitigates activation and gradient decay in deep networks,improving convergence speed and data efficiency in multiscale flow simulations.We validate the effectiveness of the model on the classical arithmetic example of lid-driven cavity flow.Compared to the traditional Xavier initialized PINN and PINN-LBM,T-PINNLBM reduces the mean absolute error(MAE)by one order of magnitude at the same network depth and maintains stable convergence in deeper networks.The results demonstrate that this model can accurately capture complex flow structures without prior data,providing a new feasible pathway for data-free driven fluid simulation.展开更多
Understanding the complex interaction between heat and mass transfer in non-Newtonian microflows is essential for the development and optimization of efficient microfluidic and thermal management systems.This study in...Understanding the complex interaction between heat and mass transfer in non-Newtonian microflows is essential for the development and optimization of efficient microfluidic and thermal management systems.This study investigates the magnetohydrodynamic(MHD)thermosolutal convection of a Casson fluid within an inclined,porous microchannel subjected to convective boundary conditions.The nonlinear,coupled equations governing momentum,energy,and species transport are solved numerically using the MATLAB bvp4c solver,ensuring high numerical accuracy and stability.To identify the dominant parameters influencing flow behavior and to optimize transport performance,a comprehensive hybrid optimization framework—combining a modified Taguchi design,Grey Relational Analysis(GRA),and Principal Component Analysis(PCA)—is proposed.This integrated strategy enables the simultaneous assessment of skin friction,Nusselt number,and Sherwood number,providing a rigorous multi-objective evaluation of system performance.Comparative validation with benchmark results from the literature confirms the accuracy and reliability of the present formulation and its numerical implementation.The results highlight the intricate coupling among flow slip,buoyancy effects,and convective transport mechanisms.Increased slip flow enhances axial velocity,while a higher solutal Biot number intensifies concentration gradients near the channel walls.Conversely,a lower thermal Biot number diminishes the temperature field,indicating weaker heat transfer across the boundaries.PCA results reveal that the first principal component(PC1)accounts for most of the system variance,demonstrating the dominant influence of coupled flow and transport parameters on overall system performance.展开更多
In this paper,the authors examine various slip effects on themagnetic field and thermal radiative impacts on the flow,mass and heat transfer of a Jeffrey nanofluid over a 2-dimensional inclined stretching sheet by a p...In this paper,the authors examine various slip effects on themagnetic field and thermal radiative impacts on the flow,mass and heat transfer of a Jeffrey nanofluid over a 2-dimensional inclined stretching sheet by a porous media.The offered work is modelled to be in the form of a combination of coupled highly nonlinear partial differential equations in dimensional contexts.Governing equations were obtained,dimensionless parameters were defined in terms of similarity parameters,and the solutionswere obtained by the Homotopy Analysis Method(HAM).The analysis is significant as the effects of viscosity are identified and the important parameters are to be determined that could eventually control a type of flowbehaviour,especially in promoting the flowand inhibiting flowof velocity,temperature,and concentrations.The findings show that such an increase in themagnetic parameter decreases the velocity profile by approximately 15%due to more Lorentz forces,and thermal radiation increases the temperature profile by up to 25%,therefore,enhancing the rate of heat transfer.The process of Brownian motion and thermophoresis increases the depth of the thermal boundary layer by 10–20 percent and reduces in concentration profiles by 12 percent when the Brownian motion parameter increases.A velocity slip parameter lowers the velocity field by about 18 percent,and a parameter of permeability lowers the momentum of flow by another 10 percent.The HAM solutions show very high accuracy levels,having an order of convergence at level 15 and errormargins are well below 0.01 percent compared to the earlier studies.All these findings can provide profound knowledge in improving heat transmission in non-Newtonian fluid systems and can be used in biomedical engineering,thermal insulation,and industrial processes such as polymer extrusion and cooling technology.Principles of heat and mass transfer give us the crucial foundation on which to study the behavior of heat andmaterial flows in other engineering and scientific disciplines.Such principles apply to various fields of study,including the following engineering fields:mechanical,chemical,aerospace,civil,and environmental.展开更多
VOF(Volume Of Fluid)方法能够通过在欧拉网格上使用离散的体积分数域表示光滑界面,在不可混合流体的数值模拟中得到广泛应用。针对多相流仿真中的液滴曲率计算问题,开发了一种计算界面曲率的算法。首先提出了一种新的数据生成方法,在...VOF(Volume Of Fluid)方法能够通过在欧拉网格上使用离散的体积分数域表示光滑界面,在不可混合流体的数值模拟中得到广泛应用。针对多相流仿真中的液滴曲率计算问题,开发了一种计算界面曲率的算法。首先提出了一种新的数据生成方法,在液滴界面上进行随机采样,增强网格内体积分数的信息量,并调整取值范围以覆盖正负曲率。然后改进了传统的深度神经网络(DNN)模型,使其在计算曲率时保持对称性。基于VOF方法与该模型,利用目标单元及邻近单元体积分数计算曲率。最后选取最优模型并应用于Basilisk软件中,以提高计算曲率的准确性和稳定性。测试结果表明,其计算曲率时准确稳定。在计算小半径液滴曲率时,误差减小了25%至50%,并能用于液滴融合仿真,证明了其应用价值。展开更多
基金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.
基金supported by the grant from the Ministry of science and Higher education of the Republic of Kazakhstan(AP23489948).
文摘This paper aims to investigate a dam break in a channel with a bend in the presence of several obstacles.To accurately determine the flood zones,it is necessary to take into account many factors such as terrain,reservoir volume.Numerical modeling was used to determine the flood zone.Numerical modeling based on the Navier-Stokes equations with a turbulent k-epsilon RNG model,the Volume of Fluid(VOF)method and the PISO algorithm were used to analyze the flow in a bend channel at an angle of 10 with the obstacles.To verify the numerical model,a test on dam break in the 450 channel was conducted.The simulation results were compared with experimental data and with the numerical data of existing data.Having been convinced of the correctness of the mathematical model,the authors carried out a numerical simulation of the main problem in three versions:without barriers,with one obstacle,with two obstacles.According to the obtained numerical results,it can be noted that irregular landforms held the flow,a decrease in water level and a slower time for water emergence could be seen.Thus,the water flow without an obstacle,with one obstacle and with two obstacles showed 4.2 s,4.4 s and 4.6 s of the time of water appearance,respectively.This time shift can give a certain advantage when conducting various events to evacuate people.
文摘This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.
基金funded by the National Natural Science Foundation of China (42174131)。
文摘The resistivity difference between oil and gas layers and the water layers in low contrast tight sandstone reservoirs is subtle. Fluid identification and saturation calculation based on conventional logging methods are facing challenges in such reservoirs. In this paper, a new method is proposed for fluid identification and saturation calculation in low contrast tight sandstone reservoirs. First, a model for calculating apparent formation water resistivity is constructed, which takes into account the influence of shale on the resistivity calculation and avoids apparent formation water resistivity abnormal values.Based on the distribution of the apparent formation water resistivity obtained by the new model, the water spectrum is determined for fluid identification in low contrast tight sandstone reservoirs.Following this, according to the average, standard deviation, and endpoints of the water spectrum, a new four-parameter model for calculating reservoir oil and gas saturation is built. The methods proposed in this paper are applied to the low contrast tight sandstone reservoirs in the Q4 formation of the X53 block and X70 block in the south of Songliao Basin, China. The results show that the water spectrum method can effectively distinguish oil-water layers and water layers in the study area. The standard deviation of the water spectrum in the oil-water layer is generally greater than that in the water layer. The new four-parameter model yields more accurate oil and gas saturation. These findings verify the effectiveness of the proposed methods.
基金Project supported by the National Natural Science Foundation of China (No.10672097)Shanghai Leading Academic Discipline Project (No.Y0103)
文摘A stencil-like volume of fluid (VOF) method is proposed for tracking free interface. A stencil on a grid cell is worked out according to the normal direction of the interface, in which only three interface positions are possible in 2D cases, and the interface can be reconstructed by only requiring the known local volume fraction information. On the other hand, the fluid-occupying-length is defined on each side of the stencil, through which a unified fluid-occupying volume model and a unified algorithm can be obtained to solve the interface advection equation. The method is suitable for the arbitrary geometry of the grid cell, and is extendible to 3D cases. Typical numerical examples show that the current method can give "sharp" results for tracking free interface.
基金supported by the National Natural Science Foundation of China(Grant Nos.12302435 and 12221002)。
文摘Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.
基金Supported by the National Key R&D Program of China(2022YFB3303500).
文摘Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data immersion method(BDIM)for compressible flows,which were solved using a fully self-programmed Fortran code.Results show that the predicted thrust spectrum aligns closely with the experimental measurements.In addition,it captures the thrust dipole directivity of the noise around the rotating propeller due to random pressure pulsations on the blades,as well as the flow structures simultaneously.Furthermore,the differences in the statistical characteristics,flow structures,and low-frequency broadband thrust spectra due to different turbulence levels were investigated.This analysis indicates that the interaction between the upstream,which is characterized by a lower turbulence level and a higher turbulent length of scale,and the rotating propeller results in a lower amplitude in force spectra and a slight increase in the scale of tip vortices.
文摘Characteristics of heat transfer and flow of Newtonian and non-Newtonian fluids through porous walls and in porous media are studied due to their wide range of applications including geothermal reservoirs,heat exchangers,marine propulsion,and aerodynamics.The current study investigates the characteristics of heat transport in a reactive third-grade fluid,moving through permeable parallel plates,with uniform suction/injection velocity.The two permeable,parallel plates are maintained at the same,constant temperature.After being transformed into its dimensionless equivalent,governing equations are solved by employing the Least Squares Method(LSM).The LSM results are further validated with numerical solutions for temperature and velocity.The impact of cross-flow Reynolds number,Peclet number,heat generation parameter,non-Newtonian parameter,and Brinkman number on entropy generation,velocity,temperature,and Bejan number are investigated.Theresults indicate that temperature distribution is significantly influenced by the third-grade fluid parameter.The maximum temperature drops from almost 0.12 to 0.10 as the third-grade fluid parameter increases from0.05 to 0.4.When the cross-flow Reynolds number is raised from 0.05 to 3,the maximum temperature drops from 0.12 to around 0.09.Temperature is strongly influenced by the heat generation parameter.A greater understanding of the thermal characteristics necessary for the design of a variety of systems,such as heat exchangers,marine propulsion,aerodynamic systems,etc.,may be gained from the findings of the current study.
文摘This study rigorously examines the interplay between viscous dissipation,magnetic effects,and thermal radiation on the flow behavior of a non-Newtonian Carreau squeezed fluid passing by a sensor surface within a micro cantilever channel,aiming to deepen our understanding of heat transport processes in complex fluid dynamics scenarios.The primary objective is to elucidate how physical operational parameters influence both the velocity of fluid flow and its temperature distribution,utilizing a comprehensive numerical approach.Employing a combination of mathematical modeling techniques,including similarity transformation,this investigation transforms complex partial differential equations into more manageable ordinary ones,subsequently solving them using the homotopy perturbation method.By analyzing the obtained solutions and presenting them graphically,alongside detailed analysis,the study sheds light on the pivotal role of significant parameters in shaping fluid movement and energy distribution.Noteworthy observations reveal a substantial increase in fluid velocity with escalating magnetic parameters,while conversely,a contrasting trend emerges in the temperature distribution,highlighting the intricate relationship between magnetic effects,flow dynamics,and thermal behavior in non-Newtonian fluids.Further,the suction velocity enhance both the local skin friction and Nusselt numbers,whereas theWeissenberg number reduces them,opposite to the effect of the power-law index.
基金supported by the National Natural Science Foundation of China(Grant No.U2242214)the Guangdong Basic and Applied Basic Research Foundation(Grant No.2024A1515010927)+3 种基金the Humanities and Social Science Foundation of the Ministry of Education in China(Grant No.24YJCZH163)the Fujian Provincial Units Special Funds for Education and Research(Grant No.2022639)Fundamental Research Funds for the Central Universities,Sun Yatsen University(Grant No.24qnpy044)partly supported by the Open Research Fund of Key Laboratory of Analytical Mathematics and Applications(Fujian Normal University),Ministry of Education,P.R.China(Grant No.JAM2405)。
文摘The Richtmyer–Meshkov(RM)instability plays an important role in various natural and engineering fields such as inertial confinement fusion.In this study,the effect of relaxation time on the RM instability under reshock impact is investigated using a two-component discrete Boltzmann method.The hydrodynamic and thermodynamic characteristics of the fluid system are comprehensively analyzed from the perspectives of the density gradient,vorticity,kinetic energy,mixing degree,mixing width and non-equilibrium intensity.Simulation results indicate that for longer relaxation time,the diffusion and dissipation are enhanced,the physical gradients decrease,and the growth of the interface is suppressed.Furthermore,the non-equilibrium manifestations show complex patterns,driven by the competitive physical mechanisms of the diffusion,dissipation,shock wave,rarefaction wave,transverse wave and fluid instabilities.These findings provide valuable insight into the fundamental mechanism of compressible fluid flows.
文摘We are pleased to introduce this special thematic section on Nanofluid Mechanics and Heat Transfer that is being included in Acta Mechanica Sinica(AMS).This thematic issue consists of 6 papers selected from papers that were presented at the 18th International Symposium on Numer-ical Analysis of Fluid Flows,Heat and Mass Transfer-Nu-merical Fluids 2023,held in Heraklion,Crete Greece,11-17 September 2023,and invited through a general call.The symposium covers various subjects:from new numerical methods and fundamental research until engineering appli-cations,and it is a part of the International Conference of Numerical Analysis and Applied Mathematics(ICNAAM),held annually.
基金supported by the National Natural Science Foundation of China(Grant Nos.12474453,12174085,and 12404530).
文摘Physics-informed neural networks(PINNs)have shown considerable promise for performing numerical simulations in fluid mechanics.They provide mesh-free,end-to-end approaches by embedding physical laws into their loss functions.However,when addressing complex flow problems,PINNs still face some challenges such as activation saturation and vanishing gradients in deep network training,leading to slow convergence and insufficient prediction accuracy.We present physics-informed neural networks incorporating lattice Boltzmann method optimized by tanh robust weight initialization(T-PINN-LBM)to address these challenges.This approach fuses the mesoscopic lattice Boltzmann model with the automatic differentiation framework of PINNs.It also implements a tanh robust weight initialization method derived from fixed point analysis.This model effectively mitigates activation and gradient decay in deep networks,improving convergence speed and data efficiency in multiscale flow simulations.We validate the effectiveness of the model on the classical arithmetic example of lid-driven cavity flow.Compared to the traditional Xavier initialized PINN and PINN-LBM,T-PINNLBM reduces the mean absolute error(MAE)by one order of magnitude at the same network depth and maintains stable convergence in deeper networks.The results demonstrate that this model can accurately capture complex flow structures without prior data,providing a new feasible pathway for data-free driven fluid simulation.
文摘Understanding the complex interaction between heat and mass transfer in non-Newtonian microflows is essential for the development and optimization of efficient microfluidic and thermal management systems.This study investigates the magnetohydrodynamic(MHD)thermosolutal convection of a Casson fluid within an inclined,porous microchannel subjected to convective boundary conditions.The nonlinear,coupled equations governing momentum,energy,and species transport are solved numerically using the MATLAB bvp4c solver,ensuring high numerical accuracy and stability.To identify the dominant parameters influencing flow behavior and to optimize transport performance,a comprehensive hybrid optimization framework—combining a modified Taguchi design,Grey Relational Analysis(GRA),and Principal Component Analysis(PCA)—is proposed.This integrated strategy enables the simultaneous assessment of skin friction,Nusselt number,and Sherwood number,providing a rigorous multi-objective evaluation of system performance.Comparative validation with benchmark results from the literature confirms the accuracy and reliability of the present formulation and its numerical implementation.The results highlight the intricate coupling among flow slip,buoyancy effects,and convective transport mechanisms.Increased slip flow enhances axial velocity,while a higher solutal Biot number intensifies concentration gradients near the channel walls.Conversely,a lower thermal Biot number diminishes the temperature field,indicating weaker heat transfer across the boundaries.PCA results reveal that the first principal component(PC1)accounts for most of the system variance,demonstrating the dominant influence of coupled flow and transport parameters on overall system performance.
文摘In this paper,the authors examine various slip effects on themagnetic field and thermal radiative impacts on the flow,mass and heat transfer of a Jeffrey nanofluid over a 2-dimensional inclined stretching sheet by a porous media.The offered work is modelled to be in the form of a combination of coupled highly nonlinear partial differential equations in dimensional contexts.Governing equations were obtained,dimensionless parameters were defined in terms of similarity parameters,and the solutionswere obtained by the Homotopy Analysis Method(HAM).The analysis is significant as the effects of viscosity are identified and the important parameters are to be determined that could eventually control a type of flowbehaviour,especially in promoting the flowand inhibiting flowof velocity,temperature,and concentrations.The findings show that such an increase in themagnetic parameter decreases the velocity profile by approximately 15%due to more Lorentz forces,and thermal radiation increases the temperature profile by up to 25%,therefore,enhancing the rate of heat transfer.The process of Brownian motion and thermophoresis increases the depth of the thermal boundary layer by 10–20 percent and reduces in concentration profiles by 12 percent when the Brownian motion parameter increases.A velocity slip parameter lowers the velocity field by about 18 percent,and a parameter of permeability lowers the momentum of flow by another 10 percent.The HAM solutions show very high accuracy levels,having an order of convergence at level 15 and errormargins are well below 0.01 percent compared to the earlier studies.All these findings can provide profound knowledge in improving heat transmission in non-Newtonian fluid systems and can be used in biomedical engineering,thermal insulation,and industrial processes such as polymer extrusion and cooling technology.Principles of heat and mass transfer give us the crucial foundation on which to study the behavior of heat andmaterial flows in other engineering and scientific disciplines.Such principles apply to various fields of study,including the following engineering fields:mechanical,chemical,aerospace,civil,and environmental.
文摘VOF(Volume Of Fluid)方法能够通过在欧拉网格上使用离散的体积分数域表示光滑界面,在不可混合流体的数值模拟中得到广泛应用。针对多相流仿真中的液滴曲率计算问题,开发了一种计算界面曲率的算法。首先提出了一种新的数据生成方法,在液滴界面上进行随机采样,增强网格内体积分数的信息量,并调整取值范围以覆盖正负曲率。然后改进了传统的深度神经网络(DNN)模型,使其在计算曲率时保持对称性。基于VOF方法与该模型,利用目标单元及邻近单元体积分数计算曲率。最后选取最优模型并应用于Basilisk软件中,以提高计算曲率的准确性和稳定性。测试结果表明,其计算曲率时准确稳定。在计算小半径液滴曲率时,误差减小了25%至50%,并能用于液滴融合仿真,证明了其应用价值。