To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced v...To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced velocities and gap spacing to observe flow-induced vibrations(FIVs).The immersed boundary–lattice Boltzmann flux solver(IB–LBFS)was applied as a numerical solution method,allowing for straightforward application on a simple Cartesian mesh.The accuracy and rationality of this method have been verified through comparisons with previous numerical results,including studies on flow past three stationary circular cylinders arranged in a similar pattern and vortex-induced vibrations of a single cylinder across different reduced velocities.When examining the FIVs of three cylinders,numerical simulations were carried out across a range of reduced velocities(3.0≤Ur≤13.0)and gap spacing(L=3D,4D,and 5D).The observed vibration response included several regimes:the desynchronization regime,the initial branch,and the lower branch.Notably,the transverse amplitude peaked,and a double vortex street formed in the wake when the reduced velocity reached the lower branch.This arrangement of three cylinders proved advantageous for energy capture as the upstream cylinder’s vibration response mirrored that of an isolated cylinder,while the response of each downstream cylinder was significantly enhanced.Compared to a single cylinder,the vibration and flow characteristics of this system are markedly more complex.The maximum transverse amplitudes of the downstream cylinders are nearly identical and exceed those observed in a single-cylinder set-up.Depending on the gap spacing,the flow pattern varied:it was in-phase for L=3D,antiphase for L=4D,and exhibited vortex shedding for L=5D.The wake configuration mainly featured double vortex streets for L=3D and evolved into two pairs of double vortex streets for L=5D.Consequently,it well illustrates the coupling mechanism that dynamics characteristics and wake vortex change with gap spacing and reduced velocities.展开更多
The diffuse-interface immersed boundary method(IBM)possesses excellent capabilities for simulating flows around complex geometries and moving boundaries.In this method,the flow field is solved on a fixed Cartesian mes...The diffuse-interface immersed boundary method(IBM)possesses excellent capabilities for simulating flows around complex geometries and moving boundaries.In this method,the flow field is solved on a fixed Cartesian mesh,while the solid boundary is discretized into a series of Lagrangian points immersed in the flow field.The boundary condition is implemented by introducing a force term into the momentum equation,and the interaction between the immersed boundary and the fluid domain is achieved via an interpolation process.Over the past decades,the diffuse-interface IBM has gained popularity and spawned many variants,effectively handling a wide range of flow problems from isothermal to thermal flows,from laminar to turbulent flows,and from complex geometries to fluidstructure interaction scenarios.This paper first outlines the basic principles of the diffuse-interface IBM,then highlights recent advancements achieved by the authors’research group,and finally shows the method’s excellent numerical performance and wide applicability through several case studies involving complex moving boundary problems.展开更多
Jet ventilation is widely used in the ventilation design of highway and railway tunnels as an important air supply method during tunnel operation and disaster periods.This ventilation method has also been applied for ...Jet ventilation is widely used in the ventilation design of highway and railway tunnels as an important air supply method during tunnel operation and disaster periods.This ventilation method has also been applied for fire control in immersed tunnels.We conduct numerical simulations using computational fluid dynamics(CFD)to study positive ventilation in the upstream and reverse ventilation in the downstream(P-R)for an extra-wide immersed tunnel.The effects of fire source location and jet fan air velocity response strategy on the ceiling temperature decay,carbon monoxide(CO)distribution,and smoke exhaust efficiency were investigated for varying fire source locations.The results show that flames will be tilted to the side of the jet fan with a smaller air velocity.Additionally,the jet fan air velocity should be adjusted based on the relative distance between the fire source and the smoke vent.Among the studied scenarios,the most effective outcome was achieved when the air velocity was adjusted to 25 m/s on the side near the smoke vent.Also in this scenario,the phenomenon of smoke deposition was effectively mitigated and the average smoke exhaust efficiency reached 87%.Moreover,we found that the temperature decay of the tunnel follows an exponential decay law.The temperature decay rate is significantly higher on the side closest to the smoke vent compared to the farther side.This research provides a theoretical basis for smoke control strategies for fires that occur in immersed tunnels.展开更多
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.展开更多
Safe operation of underground reservoirs in coal mines is crucial for the coordinated exploitation of coal and water resources in western China.Mine water infiltration significantly influences the stability of the coa...Safe operation of underground reservoirs in coal mines is crucial for the coordinated exploitation of coal and water resources in western China.Mine water infiltration significantly influences the stability of the coal pillar.Therefore,laboratory tests were systematically carried out on coal from the Daliuta Coal Mine in Northwest China.Samples were taken in the vertical and parallel bedding directions and soaked for 0 d,2 d,4 d,or 16 d.In this study,atomic absorption spectroscopy(AAS),X-ray diffraction(XRD),and scanning electron microscopy(SEM)were used to analyze the variations in the water absorption characteristics and corresponding internal structure of the coal.Uniaxial compression tests and synchronous acoustic emission(AE)monitoring revealed the sample failure process and mechanical properties of the samples.Finally,the time-and frequency-domain characteristics of the AE signal were comprehensively analyzed using fractal dimension,fast Fourier transform,and cluster analysis.The strength and elastic modulus demonstrate significant anisotropy with different bedding planes and reveal the existence of the optimum water content.Specifically,the sample strength increases after 2 d of immersion,with increments of 23.3% and 0.6% for the vertical and parallel bedding samples,and decreases after 16 d of immersion,with decreases of 29% and 45% for the vertical and parallel samples,respectively.Additionally,shear cracks account for over 60% during the damage development of the samples.The proportion of tensile cracks is higher for samples with longer immersion times and parallel bedding planes.This research provides a theoretical basis for stability evaluation and protection of coal pillars in underground reservoirs using the AE technique.展开更多
In industrial applications,plate-like structures such as steel strips in continuous hot-dip galvanizing and papers under fan action are ubiquitous.The vibration issues that arise when these structures are in axial mot...In industrial applications,plate-like structures such as steel strips in continuous hot-dip galvanizing and papers under fan action are ubiquitous.The vibration issues that arise when these structures are in axial motion,and are influenced by fluids and thermal fields,have attracted significant attention from the academic community.This study focuses on the nonlinear dynamic behavior of axially transporting immersed viscoelastic plates with particular emphasis on internal resonance and speed-dependent tension.The governing equation and the related boundary conditions for the axially transporting viscoelastic immersed plate are derived with Hamilton's principle,prioritizing the impact of time-varying tension induced by speed perturbations.Based on the second-order Galerkin truncation,the governing equation is discretized into a system of second-order ordinary differential equations.The multi-scale method is used to analyze the stable steady-state response of the immersed viscoelastic plate.The conditions for achieving a 3:1 frequency ratio between the first two orders of the system are analytically deduced.Notably,when the viscoelastic coefficient diminishes,the stability boundaries exhibit increased complexity,manifesting as the irregular W-shaped contours in the parameter space.Numerical examples comprehensively investigate the effects of viscoelasticity on both the stability region and the steady-state response under internal resonance conditions.Finally,the accuracy of the obtained results is validated through numerical computation.展开更多
This paper presents the extension of the global description approach of a discontinuous function, which is proposed in the previous paper, to a spectral domain decomposition method. This multi-domain spectral immersed...This paper presents the extension of the global description approach of a discontinuous function, which is proposed in the previous paper, to a spectral domain decomposition method. This multi-domain spectral immersed interlace method（IIM） divides the whole computation domain into the smooth and discontinuous parts. Fewer points on the smooth domains are used via taking advantage of the high accuracy property of the spectral method, but more points on the discontinuous domains are employed to enhance the resolution of the calculation. Two that the domain decomposition technique can placed around the discontinuity. The present reached, in spite of the enlarged computational discontinuous problems are tested to verify the present method. The results show reduce the error of the spectral IIM, especially when more collocation points are method is t：avorable for the reason that the same level of the accuracy can be domain.展开更多
A new simple fictitious domain method, the algebraic immersed interface and boundary (AIIB) method, is presented for elliptic equations with immersed interface conditions. This method allows jump conditions on immerse...A new simple fictitious domain method, the algebraic immersed interface and boundary (AIIB) method, is presented for elliptic equations with immersed interface conditions. This method allows jump conditions on immersed interfaces to be discretized with a good accuracy on a compact stencil. Auxiliary unknowns are created at existing grid locations to increase the degrees of freedom of the initial problem. These auxiliary unknowns allow imposing various constraints to the system on interfaces of complex shapes. For instance, the method is able to deal with immersed interfaces for elliptic equations with jump conditions on the solution or discontinuous coefficients with a second order of spatial accuracy. As the AIIB method acts on an algebraic level and only changes the problem matrix, no particular attention to the initial discretization is required. The method can be easily implemented in any structured grid code and can deal with immersed boundary problems too. Several validation problems are presented to demonstrate the interest and accuracy of the method.展开更多
A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equ...A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.展开更多
Bridges and tunnels are good solutions to transportation problems in large cities separated by large rivers. In bridge construction great success has been achieved in China, but large-sized immersed tube tunnel constr...Bridges and tunnels are good solutions to transportation problems in large cities separated by large rivers. In bridge construction great success has been achieved in China, but large-sized immersed tube tunnel construction is still new. Element immersing is an important process of immersed tube tunnel construction. The accuracy of tunnel element positioning directly determines the quality of tunnel construction. In order to study the behavior of elements during its lowering to the sea bed, the experiments carried out in the State Key Laboratory of Ocean Engineering of Shanghai Jiaotong University. In consideration of the construction experience abroad and by reference to published papers on the Oresund tunnel in Norway-Sweden and Tokyo Bay tunnel in Japan, an element model to an appropriate scale is developed. A concise description of the model experiment wave environments is carried out, and the feasibility of two immersing strategies is studied.展开更多
A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the ...A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the Reynolds number based on the hull length is 1.0x 105, An immersed boundary method based on the moving-least-squares reconstruction is used to handle the complex geometric boundaries. The adaptive mesh refinement is utilized to resolve the flows near the hull, The parallel scalabilities of the flow solver are tested on meshes with the number of cells varying from 50 million to 3.2 billion, The parallel solver reaches nearly linear scalability for the flows around the underwater vehicle model, The present simulation captures the essential features of the vortex structures near the hull and in the wake, Both of the time-averaged pressure coefficients and srreamwise velocity profiles obtained from the LES are consistent with the characteristics of the flows pass an appended axisymmetric body. The code efficiency and its correct predictions on flow features allow us to perform the full-scale simulations on tens of thousands of cores with billions of grid points for higher-Reynolds-number flows around the underwater vehicles.展开更多
Mooring,towing and installing tests on immersed tunnel caissons across the Yong River are carried out in basin to investigate the motions and line tensions.In the towing operation,the maximum tension exists in the cas...Mooring,towing and installing tests on immersed tunnel caissons across the Yong River are carried out in basin to investigate the motions and line tensions.In the towing operation,the maximum tension exists in the case that half caisson is exposed to waves and currents.The cooperation of operating manually on many lines at the same time is very difficult and important.The dynamic behavior of the caisson and the line tensions are affected more by the manual operations than by waves and currents.In the installing operation,the tensions of hoisting lines are much larger than the positioning lines.展开更多
The aim of this study is to develop a model of fluid and heat transfer in a biological tissue taking into account the exact structure of the related microvascular network,and to analyze the influence of structural cha...The aim of this study is to develop a model of fluid and heat transfer in a biological tissue taking into account the exact structure of the related microvascular network,and to analyze the influence of structural changes of such a network induced by diabetes.A cubic region representing local skin tissue is selected as the computational domain,which in turn includes two intravascular and extravascular sub-domains.To save computational resources,the capillary network is reduced to a 1D pipeline model and embedded into the extravascular region.On the basis of the immersed boundary method(IBM)strategy,fluid and heat fluxes across a capillary wall are distributed to the surrounding tissue nodes by a delta function.We consider both steady and periodic blood pressure conditions at the entrances of the capillary network.Under steady blood pressure conditions,both the interstitial fluid pressure and tissue temperature around the capillary network are larger than those in other places.When the periodic blood pressure condition is considered,tissue temperature tends to fluctuate with the same frequency of the forcing,but the related waveform displays a smaller amplitude and a certain time(phase)delay.When the connectivity of capillary network is diminished,the capacity of blood redistribution through the capillary network becomes weaker and a subset of the vessel branches lose blood flow,which further aggravates the amplitude attenuation and time delay of the skin temperature fluctuation.展开更多
Immersed tube tunnels are usually placed on soft soil layers in cross-sea tunnelling engineering.Owing to the influence of stratum conditions and slope design,the longitudinal distribution of substratum layers is gene...Immersed tube tunnels are usually placed on soft soil layers in cross-sea tunnelling engineering.Owing to the influence of stratum conditions and slope design,the longitudinal distribution of substratum layers is generally uneven.Thus,the inhomogeneous deformation of the element-joint becomes the key factor in the failure of the immersed tube tun-nel.Therefore,a corresponding calculation method for joint deformation is needed to explore the deformation law of immersed tube tunnels.By constructing a three-section immersed tube tunnel analysis model(TTM),the relationship between the two types of deformation of the immersed tube tunnel structure in a longitudinal nonuniform soft soil foundation is described,and the deformation characteristics of the immersed structure under different boundaries are discussed.Based on the mechanical behaviour of the joint and foundation,according to the Timoshenko beam on the Vlasov two-parameter foundation(VTM),considering the tidal cyclic load during the operation and maintenance period,an example analysis is given.Moreover,the deformation characteristics and development trend of the immersed tube tunnel under the influence of different soil layers are discussed.The obtained results have a certain guiding significance for the deformation calculation of immersed tube tunnels.展开更多
This paper introduces a new idea of controlling cavitation around a hydrofoil through a passive cavitation controller called artificial cavitation bubble generator （ACG）. Cyclic processes, namely, growth and implosi...This paper introduces a new idea of controlling cavitation around a hydrofoil through a passive cavitation controller called artificial cavitation bubble generator （ACG）. Cyclic processes, namely, growth and implosion of bubbles around an immersed body, are the main reasons for the destruction and erosion of the said body. This paper aims to create a condition in which the cavitation bubbles reach a steady-state situation and prevent the occurrence of the cyclic processes. For this purpose, the ACG is placed on the surface of an immersed body, in particular, the suction surface of a 2D hydrofoil. A simulation was performed with an implicit finite volume scheme based on a SIMPLE algorithm associated with the multiphase and cavitation model. The modified k-ε RNG turbulence model equipped with a modification of the turbulent viscosity was applied to overcome the turbulence closure problem. Numerical simulation of water flow over the hydrofoil equipped with the ACG shows that a low-pressure recirculation area is produced behind the ACG and artificially generates stationary cavitation bubbles. The location, shape, and size of this ACG are the crucial parameters in creating a proper control. Results show that the cavitation bubble is controlled well with a well-designed ACG.展开更多
We present a high performance modularly-built open-source software-OpenIFEM.OpenIFEM is a C++implementation of the modified immersed finite element method(mIFEM)to solve fluid-structure interaction(FSI)problems.This s...We present a high performance modularly-built open-source software-OpenIFEM.OpenIFEM is a C++implementation of the modified immersed finite element method(mIFEM)to solve fluid-structure interaction(FSI)problems.This software is modularly built to perform multiple tasks including fluid dynamics(incompressible and slightly compressible fluid models),linear and nonlinear solid mechanics,and fully coupled fluid-structure interactions.Most of open-source software packages are restricted to certain discretization methods;some are under-tested,under-documented,and lack modularity as well as extensibility.OpenIFEM is designed and built to include a set of generic classes for users to adapt so that any fluid and solid solvers can be coupled through the FSI algorithm.In addition,the package utilizes well-developed and tested libraries.It also comes with standard test cases that serve as software and algorithm validation.The software can be built on cross-platform,i.e.,Linux,Windows,and Mac OS,using CMake.Efficient parallelization is also implemented for high-performance computing for large-sized problems.OpenIFEM is documented using Doxygen and publicly available to download on GitHub.It is expected to benefit the future development of FSI algorithms and be applied to a variety of FSI applications.展开更多
The immersed tube tunnel section of the Shenzhen-Zhongshan Link exhibits complex geological conditions and high back sludge strength. The tunnel cushion adopts the gravel and flaky stone combined cushion. The major in...The immersed tube tunnel section of the Shenzhen-Zhongshan Link exhibits complex geological conditions and high back sludge strength. The tunnel cushion adopts the gravel and flaky stone combined cushion. The major influencing factors of the mechanical deformation characteristics of the gravel and flaky stone composite cushion are studied through a physical model experiment. The following results are reported.(1) The load–settlement curves of the flaky stone cushion become more compact with a dense increment under the design load. These curves can be regarded as nonlinear mechanical characteristics. The load–settlement curves of the gravel cushion and the gravel and flaky stone composite cushion exhibit the characteristics of a two-stage linear change.(2) The flatness of the top of flaky stone cushion considerably affects settlement and secant modulus. The flatness of the top of flaky stone should be ensured during construction.(3) Gradation and thickness exert no evident effect on the compressibility of a cushion. The preloading pressure caused by the construction height difference of the cushion materials plays an important role in improving the initial stiffness of a cushion and reducing initial settlement and overall settlement.(4) This study investigates the preloading under 30 kPa of the 0.7-m flaky stone and 1.0-m gravel combination cushion. It recommends the following secant modulus values: 48.89 MPa for the section of 0–30 kPa and 10.47 MPa for the section of 30–110 kPa.展开更多
This article introduces a numerical scheme on the basis of semi-implicit method for pressure-linked equations （SIMPLE） algorithm to simulate incompressible unsteady flows with fluid-structure interaction. The Navier...This article introduces a numerical scheme on the basis of semi-implicit method for pressure-linked equations （SIMPLE） algorithm to simulate incompressible unsteady flows with fluid-structure interaction. The Navier-Stokes equation is discretized spatially with collocated finite volume method and Eulerian implicit method in time domain. The hybrid method that combines immersed boundary method （IBM） and volume of fluid （VOF） method is used to deal with rigid body motion in fluid domain. The details of movement of immersed boundary （IB） and calculation of VOF are also described. This method can be easily applied to any existing finite-volume-based computational fluid dynamics （CFD） solver without complex operation, with which fluid flow interaction of arbitrarily complex geometry can be realized on a fixed mesh. The method is verified by low Reynolds number flows passing both stationary and oscillating cylinders. The drag and lift coefficients acquired by the study well accord with other published results, which indicate the reasonability of the proposed method.展开更多
A numerical study of vortex-induced rotations(VIRs) of an equivalent triangular cylinder, which is free to rotate in the azimuthal direction in a uniform flow, is presented. Based on an immersed boundary method, the n...A numerical study of vortex-induced rotations(VIRs) of an equivalent triangular cylinder, which is free to rotate in the azimuthal direction in a uniform flow, is presented. Based on an immersed boundary method, the numerical model is established, and is verified through the benchmark problem of flow past a freely rotating rectangular body.The computation is performed for a fixed reduced mass of m~*=2.0 and the structural stiffness and damping ratio are set to zero. The effects of Reynolds number(Re=25-180) on the characteristics of VIR are studied. It is found that the dynamic response of the triangular cylinder exhibits four distinct modes with increasing Re: a rest position,periodic rotational oscillation, random rotation and autorotation. For the rotational oscillation mode, the cylinder undergoes a periodic vibration around an equilibrium position with one side facing the incoming flow. Since the rotation effect, the outset of vortex shedding from cylinder shifts to a much lower Reynolds number. Further increase in Re leads to 2 P and P+S vortex shedding modes besides the typical 2 S pattern. Our simulation results also elucidate that the free rotation significantly changes the drag and lift forces. Inspired by these facts, the effect of free rotation on flow-induced vibration of a triangular cylinder in the in-line and transverse directions is investigated. The results show that when the translational vibration is coupled with rotation, the triangular cylinder presents a galloping response instead of vortex-induced vibration(VIV).展开更多
In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body...In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.展开更多
基金Supported by the National Natural Science Foundation of China(52201350,52201394,and 52271301)the Innovation Group Project of Southern Marine Science and Engineering Guangdong Laboratory(Zhuhai)(Grant No.SML2022008).
文摘To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced velocities and gap spacing to observe flow-induced vibrations(FIVs).The immersed boundary–lattice Boltzmann flux solver(IB–LBFS)was applied as a numerical solution method,allowing for straightforward application on a simple Cartesian mesh.The accuracy and rationality of this method have been verified through comparisons with previous numerical results,including studies on flow past three stationary circular cylinders arranged in a similar pattern and vortex-induced vibrations of a single cylinder across different reduced velocities.When examining the FIVs of three cylinders,numerical simulations were carried out across a range of reduced velocities(3.0≤Ur≤13.0)and gap spacing(L=3D,4D,and 5D).The observed vibration response included several regimes:the desynchronization regime,the initial branch,and the lower branch.Notably,the transverse amplitude peaked,and a double vortex street formed in the wake when the reduced velocity reached the lower branch.This arrangement of three cylinders proved advantageous for energy capture as the upstream cylinder’s vibration response mirrored that of an isolated cylinder,while the response of each downstream cylinder was significantly enhanced.Compared to a single cylinder,the vibration and flow characteristics of this system are markedly more complex.The maximum transverse amplitudes of the downstream cylinders are nearly identical and exceed those observed in a single-cylinder set-up.Depending on the gap spacing,the flow pattern varied:it was in-phase for L=3D,antiphase for L=4D,and exhibited vortex shedding for L=5D.The wake configuration mainly featured double vortex streets for L=3D and evolved into two pairs of double vortex streets for L=5D.Consequently,it well illustrates the coupling mechanism that dynamics characteristics and wake vortex change with gap spacing and reduced velocities.
基金partially supported by the National Natural Science Foundation of China(Nos.92271103,12202191)。
文摘The diffuse-interface immersed boundary method(IBM)possesses excellent capabilities for simulating flows around complex geometries and moving boundaries.In this method,the flow field is solved on a fixed Cartesian mesh,while the solid boundary is discretized into a series of Lagrangian points immersed in the flow field.The boundary condition is implemented by introducing a force term into the momentum equation,and the interaction between the immersed boundary and the fluid domain is achieved via an interpolation process.Over the past decades,the diffuse-interface IBM has gained popularity and spawned many variants,effectively handling a wide range of flow problems from isothermal to thermal flows,from laminar to turbulent flows,and from complex geometries to fluidstructure interaction scenarios.This paper first outlines the basic principles of the diffuse-interface IBM,then highlights recent advancements achieved by the authors’research group,and finally shows the method’s excellent numerical performance and wide applicability through several case studies involving complex moving boundary problems.
文摘Jet ventilation is widely used in the ventilation design of highway and railway tunnels as an important air supply method during tunnel operation and disaster periods.This ventilation method has also been applied for fire control in immersed tunnels.We conduct numerical simulations using computational fluid dynamics(CFD)to study positive ventilation in the upstream and reverse ventilation in the downstream(P-R)for an extra-wide immersed tunnel.The effects of fire source location and jet fan air velocity response strategy on the ceiling temperature decay,carbon monoxide(CO)distribution,and smoke exhaust efficiency were investigated for varying fire source locations.The results show that flames will be tilted to the side of the jet fan with a smaller air velocity.Additionally,the jet fan air velocity should be adjusted based on the relative distance between the fire source and the smoke vent.Among the studied scenarios,the most effective outcome was achieved when the air velocity was adjusted to 25 m/s on the side near the smoke vent.Also in this scenario,the phenomenon of smoke deposition was effectively mitigated and the average smoke exhaust efficiency reached 87%.Moreover,we found that the temperature decay of the tunnel follows an exponential decay law.The temperature decay rate is significantly higher on the side closest to the smoke vent compared to the farther side.This research provides a theoretical basis for smoke control strategies for fires that occur in immersed tunnels.
基金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.
基金supported by the National Natural Science Foundation of China(Grant Nos.52174084 and U23B20146)the Open Fund for State Key Laboratory of Water Resource Protection and Utilization in Coal Mining(Grant No.WPUKFJJ2022-07)。
文摘Safe operation of underground reservoirs in coal mines is crucial for the coordinated exploitation of coal and water resources in western China.Mine water infiltration significantly influences the stability of the coal pillar.Therefore,laboratory tests were systematically carried out on coal from the Daliuta Coal Mine in Northwest China.Samples were taken in the vertical and parallel bedding directions and soaked for 0 d,2 d,4 d,or 16 d.In this study,atomic absorption spectroscopy(AAS),X-ray diffraction(XRD),and scanning electron microscopy(SEM)were used to analyze the variations in the water absorption characteristics and corresponding internal structure of the coal.Uniaxial compression tests and synchronous acoustic emission(AE)monitoring revealed the sample failure process and mechanical properties of the samples.Finally,the time-and frequency-domain characteristics of the AE signal were comprehensively analyzed using fractal dimension,fast Fourier transform,and cluster analysis.The strength and elastic modulus demonstrate significant anisotropy with different bedding planes and reveal the existence of the optimum water content.Specifically,the sample strength increases after 2 d of immersion,with increments of 23.3% and 0.6% for the vertical and parallel bedding samples,and decreases after 16 d of immersion,with decreases of 29% and 45% for the vertical and parallel samples,respectively.Additionally,shear cracks account for over 60% during the damage development of the samples.The proportion of tensile cracks is higher for samples with longer immersion times and parallel bedding planes.This research provides a theoretical basis for stability evaluation and protection of coal pillars in underground reservoirs using the AE technique.
基金Project supported by the Science and Technology Support Plan for Youth Innovation of Colleges and Universities of Shandong Province of China(No.2023KJ215)the National Natural Science Foundation of China(Nos.12002142 and 52405274)the National Natural Science Foundation of Shanghai of China(No.ZR2023QE100)。
文摘In industrial applications,plate-like structures such as steel strips in continuous hot-dip galvanizing and papers under fan action are ubiquitous.The vibration issues that arise when these structures are in axial motion,and are influenced by fluids and thermal fields,have attracted significant attention from the academic community.This study focuses on the nonlinear dynamic behavior of axially transporting immersed viscoelastic plates with particular emphasis on internal resonance and speed-dependent tension.The governing equation and the related boundary conditions for the axially transporting viscoelastic immersed plate are derived with Hamilton's principle,prioritizing the impact of time-varying tension induced by speed perturbations.Based on the second-order Galerkin truncation,the governing equation is discretized into a system of second-order ordinary differential equations.The multi-scale method is used to analyze the stable steady-state response of the immersed viscoelastic plate.The conditions for achieving a 3:1 frequency ratio between the first two orders of the system are analytically deduced.Notably,when the viscoelastic coefficient diminishes,the stability boundaries exhibit increased complexity,manifesting as the irregular W-shaped contours in the parameter space.Numerical examples comprehensively investigate the effects of viscoelasticity on both the stability region and the steady-state response under internal resonance conditions.Finally,the accuracy of the obtained results is validated through numerical computation.
基金National Natural Science Foundation of China(51076006)
文摘This paper presents the extension of the global description approach of a discontinuous function, which is proposed in the previous paper, to a spectral domain decomposition method. This multi-domain spectral immersed interlace method（IIM） divides the whole computation domain into the smooth and discontinuous parts. Fewer points on the smooth domains are used via taking advantage of the high accuracy property of the spectral method, but more points on the discontinuous domains are employed to enhance the resolution of the calculation. Two that the domain decomposition technique can placed around the discontinuity. The present reached, in spite of the enlarged computational discontinuous problems are tested to verify the present method. The results show reduce the error of the spectral IIM, especially when more collocation points are method is t：avorable for the reason that the same level of the accuracy can be domain.
文摘A new simple fictitious domain method, the algebraic immersed interface and boundary (AIIB) method, is presented for elliptic equations with immersed interface conditions. This method allows jump conditions on immersed interfaces to be discretized with a good accuracy on a compact stencil. Auxiliary unknowns are created at existing grid locations to increase the degrees of freedom of the initial problem. These auxiliary unknowns allow imposing various constraints to the system on interfaces of complex shapes. For instance, the method is able to deal with immersed interfaces for elliptic equations with jump conditions on the solution or discontinuous coefficients with a second order of spatial accuracy. As the AIIB method acts on an algebraic level and only changes the problem matrix, no particular attention to the initial discretization is required. The method can be easily implemented in any structured grid code and can deal with immersed boundary problems too. Several validation problems are presented to demonstrate the interest and accuracy of the method.
基金Project supported by the National Natural Science Foundation of China(Nos.91752118,11672305,11232011,and 11572331)the Strategic Priority Research Program(No.XDB22040104)the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences(No.QYZDJ-SSWSYS002)
文摘A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.
文摘Bridges and tunnels are good solutions to transportation problems in large cities separated by large rivers. In bridge construction great success has been achieved in China, but large-sized immersed tube tunnel construction is still new. Element immersing is an important process of immersed tube tunnel construction. The accuracy of tunnel element positioning directly determines the quality of tunnel construction. In order to study the behavior of elements during its lowering to the sea bed, the experiments carried out in the State Key Laboratory of Ocean Engineering of Shanghai Jiaotong University. In consideration of the construction experience abroad and by reference to published papers on the Oresund tunnel in Norway-Sweden and Tokyo Bay tunnel in Japan, an element model to an appropriate scale is developed. A concise description of the model experiment wave environments is carried out, and the feasibility of two immersing strategies is studied.
基金supported by the National Natural Science Foundation of China (11302238, 11232011. and 11572331)support from the Strategic Priority Research Program (XDB22040104)+1 种基金the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences (QYZDJ-SSW-SYS002)the National Basic Research Program of China (973 Program 2013CB834100: Nonlinear science)
文摘A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the Reynolds number based on the hull length is 1.0x 105, An immersed boundary method based on the moving-least-squares reconstruction is used to handle the complex geometric boundaries. The adaptive mesh refinement is utilized to resolve the flows near the hull, The parallel scalabilities of the flow solver are tested on meshes with the number of cells varying from 50 million to 3.2 billion, The parallel solver reaches nearly linear scalability for the flows around the underwater vehicle model, The present simulation captures the essential features of the vortex structures near the hull and in the wake, Both of the time-averaged pressure coefficients and srreamwise velocity profiles obtained from the LES are consistent with the characteristics of the flows pass an appended axisymmetric body. The code efficiency and its correct predictions on flow features allow us to perform the full-scale simulations on tens of thousands of cores with billions of grid points for higher-Reynolds-number flows around the underwater vehicles.
基金the Science and Technology Committee of Shanghai (No. 07DZ22027)
文摘Mooring,towing and installing tests on immersed tunnel caissons across the Yong River are carried out in basin to investigate the motions and line tensions.In the towing operation,the maximum tension exists in the case that half caisson is exposed to waves and currents.The cooperation of operating manually on many lines at the same time is very difficult and important.The dynamic behavior of the caisson and the line tensions are affected more by the manual operations than by waves and currents.In the installing operation,the tensions of hoisting lines are much larger than the positioning lines.
基金This study was supported by National Natural Science Foundation of China(NSFC No.51576033)Dalian Innovative Funding of Science and Technology(2018J12SN076)NSFC No 11602053.
文摘The aim of this study is to develop a model of fluid and heat transfer in a biological tissue taking into account the exact structure of the related microvascular network,and to analyze the influence of structural changes of such a network induced by diabetes.A cubic region representing local skin tissue is selected as the computational domain,which in turn includes two intravascular and extravascular sub-domains.To save computational resources,the capillary network is reduced to a 1D pipeline model and embedded into the extravascular region.On the basis of the immersed boundary method(IBM)strategy,fluid and heat fluxes across a capillary wall are distributed to the surrounding tissue nodes by a delta function.We consider both steady and periodic blood pressure conditions at the entrances of the capillary network.Under steady blood pressure conditions,both the interstitial fluid pressure and tissue temperature around the capillary network are larger than those in other places.When the periodic blood pressure condition is considered,tissue temperature tends to fluctuate with the same frequency of the forcing,but the related waveform displays a smaller amplitude and a certain time(phase)delay.When the connectivity of capillary network is diminished,the capacity of blood redistribution through the capillary network becomes weaker and a subset of the vessel branches lose blood flow,which further aggravates the amplitude attenuation and time delay of the skin temperature fluctuation.
基金The work was financially supported by the Natural Science Foundation of Guangdong Province,China(Grant No.2022A1515011200)the State Key Laboratory for Geo-Mechanics and Deep Underground Engineering of China University of Mining&Technology(Grant No.SKLGDUEK2005)the Science and Technology Planning Project of Guangdong Province of China(Grant No.STKJ2021129).
文摘Immersed tube tunnels are usually placed on soft soil layers in cross-sea tunnelling engineering.Owing to the influence of stratum conditions and slope design,the longitudinal distribution of substratum layers is generally uneven.Thus,the inhomogeneous deformation of the element-joint becomes the key factor in the failure of the immersed tube tun-nel.Therefore,a corresponding calculation method for joint deformation is needed to explore the deformation law of immersed tube tunnels.By constructing a three-section immersed tube tunnel analysis model(TTM),the relationship between the two types of deformation of the immersed tube tunnel structure in a longitudinal nonuniform soft soil foundation is described,and the deformation characteristics of the immersed structure under different boundaries are discussed.Based on the mechanical behaviour of the joint and foundation,according to the Timoshenko beam on the Vlasov two-parameter foundation(VTM),considering the tidal cyclic load during the operation and maintenance period,an example analysis is given.Moreover,the deformation characteristics and development trend of the immersed tube tunnel under the influence of different soil layers are discussed.The obtained results have a certain guiding significance for the deformation calculation of immersed tube tunnels.
文摘This paper introduces a new idea of controlling cavitation around a hydrofoil through a passive cavitation controller called artificial cavitation bubble generator （ACG）. Cyclic processes, namely, growth and implosion of bubbles around an immersed body, are the main reasons for the destruction and erosion of the said body. This paper aims to create a condition in which the cavitation bubbles reach a steady-state situation and prevent the occurrence of the cyclic processes. For this purpose, the ACG is placed on the surface of an immersed body, in particular, the suction surface of a 2D hydrofoil. A simulation was performed with an implicit finite volume scheme based on a SIMPLE algorithm associated with the multiphase and cavitation model. The modified k-ε RNG turbulence model equipped with a modification of the turbulent viscosity was applied to overcome the turbulence closure problem. Numerical simulation of water flow over the hydrofoil equipped with the ACG shows that a low-pressure recirculation area is produced behind the ACG and artificially generates stationary cavitation bubbles. The location, shape, and size of this ACG are the crucial parameters in creating a proper control. Results show that the cavitation bubble is controlled well with a well-designed ACG.
文摘We present a high performance modularly-built open-source software-OpenIFEM.OpenIFEM is a C++implementation of the modified immersed finite element method(mIFEM)to solve fluid-structure interaction(FSI)problems.This software is modularly built to perform multiple tasks including fluid dynamics(incompressible and slightly compressible fluid models),linear and nonlinear solid mechanics,and fully coupled fluid-structure interactions.Most of open-source software packages are restricted to certain discretization methods;some are under-tested,under-documented,and lack modularity as well as extensibility.OpenIFEM is designed and built to include a set of generic classes for users to adapt so that any fluid and solid solvers can be coupled through the FSI algorithm.In addition,the package utilizes well-developed and tested libraries.It also comes with standard test cases that serve as software and algorithm validation.The software can be built on cross-platform,i.e.,Linux,Windows,and Mac OS,using CMake.Efficient parallelization is also implemented for high-performance computing for large-sized problems.OpenIFEM is documented using Doxygen and publicly available to download on GitHub.It is expected to benefit the future development of FSI algorithms and be applied to a variety of FSI applications.
基金supported by the National Key Research and Development Program of China(Nos.2018YFC0809600 and 2018YFC0809602)。
文摘The immersed tube tunnel section of the Shenzhen-Zhongshan Link exhibits complex geological conditions and high back sludge strength. The tunnel cushion adopts the gravel and flaky stone combined cushion. The major influencing factors of the mechanical deformation characteristics of the gravel and flaky stone composite cushion are studied through a physical model experiment. The following results are reported.(1) The load–settlement curves of the flaky stone cushion become more compact with a dense increment under the design load. These curves can be regarded as nonlinear mechanical characteristics. The load–settlement curves of the gravel cushion and the gravel and flaky stone composite cushion exhibit the characteristics of a two-stage linear change.(2) The flatness of the top of flaky stone cushion considerably affects settlement and secant modulus. The flatness of the top of flaky stone should be ensured during construction.(3) Gradation and thickness exert no evident effect on the compressibility of a cushion. The preloading pressure caused by the construction height difference of the cushion materials plays an important role in improving the initial stiffness of a cushion and reducing initial settlement and overall settlement.(4) This study investigates the preloading under 30 kPa of the 0.7-m flaky stone and 1.0-m gravel combination cushion. It recommends the following secant modulus values: 48.89 MPa for the section of 0–30 kPa and 10.47 MPa for the section of 30–110 kPa.
文摘This article introduces a numerical scheme on the basis of semi-implicit method for pressure-linked equations （SIMPLE） algorithm to simulate incompressible unsteady flows with fluid-structure interaction. The Navier-Stokes equation is discretized spatially with collocated finite volume method and Eulerian implicit method in time domain. The hybrid method that combines immersed boundary method （IBM） and volume of fluid （VOF） method is used to deal with rigid body motion in fluid domain. The details of movement of immersed boundary （IB） and calculation of VOF are also described. This method can be easily applied to any existing finite-volume-based computational fluid dynamics （CFD） solver without complex operation, with which fluid flow interaction of arbitrarily complex geometry can be realized on a fixed mesh. The method is verified by low Reynolds number flows passing both stationary and oscillating cylinders. The drag and lift coefficients acquired by the study well accord with other published results, which indicate the reasonability of the proposed method.
基金financially supported by the Fundamental Research Funds for the Central Universities(Grant Nos.2018B56414 and2019B12014)the National Natural Science Foundation of China(Grant No.51609077)
文摘A numerical study of vortex-induced rotations(VIRs) of an equivalent triangular cylinder, which is free to rotate in the azimuthal direction in a uniform flow, is presented. Based on an immersed boundary method, the numerical model is established, and is verified through the benchmark problem of flow past a freely rotating rectangular body.The computation is performed for a fixed reduced mass of m~*=2.0 and the structural stiffness and damping ratio are set to zero. The effects of Reynolds number(Re=25-180) on the characteristics of VIR are studied. It is found that the dynamic response of the triangular cylinder exhibits four distinct modes with increasing Re: a rest position,periodic rotational oscillation, random rotation and autorotation. For the rotational oscillation mode, the cylinder undergoes a periodic vibration around an equilibrium position with one side facing the incoming flow. Since the rotation effect, the outset of vortex shedding from cylinder shifts to a much lower Reynolds number. Further increase in Re leads to 2 P and P+S vortex shedding modes besides the typical 2 S pattern. Our simulation results also elucidate that the free rotation significantly changes the drag and lift forces. Inspired by these facts, the effect of free rotation on flow-induced vibration of a triangular cylinder in the in-line and transverse directions is investigated. The results show that when the translational vibration is coupled with rotation, the triangular cylinder presents a galloping response instead of vortex-induced vibration(VIV).
基金supported by the US ARO grants 49308-MA and 56349-MAthe US AFSOR grant FA9550-06-1-024+1 种基金he US NSF grant DMS-0911434the State Key Laboratory of Scientific and Engineering Computing of Chinese Academy of Sciences during a visit by Z.Li between July-August,2008.
文摘In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.
