Structural Health Monitoring(SHM)systems have become a crucial tool for the operational management of long tunnels.For immersed tunnels exposed to both traffic loads and the effects of the marine environment,efficient...Structural Health Monitoring(SHM)systems have become a crucial tool for the operational management of long tunnels.For immersed tunnels exposed to both traffic loads and the effects of the marine environment,efficiently identifying abnormal conditions from the extensive unannotated SHM data presents a significant challenge.This study proposed amodel-based approach for anomaly detection and conducted validation and comparative analysis of two distinct temporal predictive models using SHM data from a real immersed tunnel.Firstly,a dynamic predictive model-based anomaly detectionmethod is proposed,which utilizes a rolling time window for modeling to achieve dynamic prediction.Leveraging the assumption of temporal data similarity,an interval prediction value deviation was employed to determine the abnormality of the data.Subsequently,dynamic predictive models were constructed based on the Autoregressive Integrated Moving Average(ARIMA)and Long Short-Term Memory(LSTM)models.The hyperparameters of these models were optimized and selected using monitoring data from the immersed tunnel,yielding viable static and dynamic predictive models.Finally,the models were applied within the same segment of SHM data,to validate the effectiveness of the anomaly detection approach based on dynamic predictive modeling.A detailed comparative analysis discusses the discrepancies in temporal anomaly detection between the ARIMA-and LSTM-based models.The results demonstrated that the dynamic predictive modelbased anomaly detection approach was effective for dealing with unannotated SHM data.In a comparison between ARIMA and LSTM,it was found that ARIMA demonstrated higher modeling efficiency,rendering it suitable for short-term predictions.In contrast,the LSTM model exhibited greater capacity to capture long-term performance trends and enhanced early warning capabilities,thereby resulting in superior overall performance.展开更多
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.展开更多
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 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.展开更多
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,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of fini...In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of finite element (FE) formulation and the simple structure of Cartesian grids,the IFE discretization is used in this paper.Two-grid schemes are formulated to linearize the FE equations.It is theoretically and numerically illustrated that the coarse space can be selected as coarse as H =O(h^1/4)(or H =O(h^1/8)),and the asymptotically optimal approximation can be achieved as the nonlinear schemes.As a result,we can settle a great majority of nonlinear equations as easy as linearized problems.In order to estimate the present two-grid algorithms,we derive the optimal error estimates of the IFE solution in the L^p norm.Numerical experiments are given to verify the theorems and indicate that the present two-grid algorithms can greatly improve the computing efficiency.展开更多
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.展开更多
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.展开更多
A micro-layer model is proposed to account for the lubrication effect of liquid layer near collisions of immersed particles at moderate particle Reynolds number.This new model is to allow determination of the pressure...A micro-layer model is proposed to account for the lubrication effect of liquid layer near collisions of immersed particles at moderate particle Reynolds number.This new model is to allow determination of the pressure profile within the micro-layer including the fluid inertia and viscosity.Then a correction based on the micro-layer model is applied to unsteady 3-D direct simulation of a particle approaching another one.The simulation is based on a modified immersed boundary method with direct force scheme.The quantitative agreement between numerical and experimental results validates the model presented in the study.The simulation results show that the fluid is squeezed prior to contact.When a particle approaches a flat wall or another particle,the lubrication force,obtained by integrating the pressure profile over the particle surface,is increased and prevents the particle from approaching.The model predicts that the velocity of approaching particle starts to decrease when separation distance of particles is less than 0.1dp,where dp is the particle diameter.展开更多
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.展开更多
Immersed tunnel is an important part of the Hong Kong–Zhuhai–Macao Bridge(HZMB) project. In immersed tunnel floating, translation which includes straight and transverse movements is the main working mode. To decide ...Immersed tunnel is an important part of the Hong Kong–Zhuhai–Macao Bridge(HZMB) project. In immersed tunnel floating, translation which includes straight and transverse movements is the main working mode. To decide the magnitude and direction of the towing force for each tug, a particle swarm-based translation control method is presented for non-power immersed tunnel element. A sort of linear weighted logarithmic function is exploited to avoid weak subgoals. In simulation, the particle swarm-based control method is evaluated and compared with traditional empirical method in the case of the HZMB project. Simulation results show that the presented method delivers performance improvement in terms of the enhanced surplus towing force.展开更多
Currently, many studies on the local discontinuous Galerkin method focus on the Cartesian grid with low computational e ciency and poor adaptability to complex shapes. A new immersed boundary method is presented, and ...Currently, many studies on the local discontinuous Galerkin method focus on the Cartesian grid with low computational e ciency and poor adaptability to complex shapes. A new immersed boundary method is presented, and this method employs the adaptive Cartesian grid to improve the adaptability to complex shapes and the immersed boundary to increase computational e ciency. The new immersed boundary method employs different boundary cells(the physical cell and ghost cell) to impose the boundary condition and the reconstruction algorithm of the ghost cell is the key for this method. The classical model elliptic equation is used to test the method. This method is tested and analyzed from the viewpoints of boundary cell type, error distribution and accuracy. The numerical result shows that the presented method has low error and a good rate of the convergence and works well in complex geometries. The method has good prospect for practical application research of the numerical calculation research.展开更多
A method combining the immersed boundary technique and a multirelaxation-time(MRT) lattice Boltzmann flux solver(LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic cylinders...A method combining the immersed boundary technique and a multirelaxation-time(MRT) lattice Boltzmann flux solver(LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic cylinders and NACA 0012 Airfoil. The method uses a simple Cartesian mesh to simulate flows past immersed complicated bodies. With the Chapman-Enskog expansion analysis, a transform is performed between the Navier-Stokes and lattice Boltzmann equations(LBEs). The LBFS is used to discretize the macroscopic differential equations with a finite volume method and evaluate the interface fluxes through local reconstruction of the lattice Boltzmann solution.The immersed boundary technique is used to correct the intermediate velocity around the solid boundary to satisfy the no-slip boundary condition. Agreement of simulation results with the data found in the literature shows reliability of the proposed method in simulating laminar flows on a Cartesian mesh.展开更多
The Immersed Interface Method (IIM) is derived to solve the corresponding Fokker-Planck equation of Brownian motion with pure dry friction, which is one of the simplest models of piecewise-smooth stochastic systems. T...The Immersed Interface Method (IIM) is derived to solve the corresponding Fokker-Planck equation of Brownian motion with pure dry friction, which is one of the simplest models of piecewise-smooth stochastic systems. The IIM is capable of treating a discontinuity in the drift of Fokker-Planck equation and it is readily extended to the dry and viscous friction model. Analytic results of the considered model are used to confirm the effectiveness and design accuracy of the scheme.展开更多
We present a cut-cell method for the simulation of 2D incompressible flows past obstacles.It consists in using the MAC scheme on cartesian grids and imposing Dirchlet boundary conditions for the velocity field on the ...We present a cut-cell method for the simulation of 2D incompressible flows past obstacles.It consists in using the MAC scheme on cartesian grids and imposing Dirchlet boundary conditions for the velocity field on the boundary of solid structures following the Shortley-Weller formulation.In order to ensure local conservation properties,viscous and convecting terms are discretized in a finite volume way.The scheme is second order implicit in time for the linear part,the linear systems are solved by the use of the capacitance matrix method for non-moving obstacles.Numerical results of flows around an impulsively started circular cylinder are presented which confirm the efficiency of the method,for Reynolds numbers 1000 and 3000.An example of flows around a moving rigid body at Reynolds number 800 is also shown,a solver using the PETSc-Library has been prefered in this context to solve the linear systems.展开更多
The electrochemical quartz crystal impedance system (EQCIS) has been used for the study of a partially immersed An electrode in 0.2 mol/L NaCIO4. aqueous solution. The influences of the immersed area and height of t...The electrochemical quartz crystal impedance system (EQCIS) has been used for the study of a partially immersed An electrode in 0.2 mol/L NaCIO4. aqueous solution. The influences of the immersed area and height of the electrode on the EQCIS responses were evaluated, showing the highest response sensitivity to liquid loading at the center of the piezoelectric quartz crystal electrode. The increase in the immersed height of the An electrode at oxygen reduction potentials during potential cycling was measured by this technique.展开更多
In this paper, the transportation mode of tube during the immersed tunnel engineering in Yangtze River incorporated, mathematical models of the combination of the tube and four barges are established when it is towed ...In this paper, the transportation mode of tube during the immersed tunnel engineering in Yangtze River incorporated, mathematical models of the combination of the tube and four barges are established when it is towed by a tugboat through an elastic rope and pushed by two pushers. Based on evaluating the hydrodynamic derivatives by regressive formula, the simulations and analysis of the whole maneuvering motion course are accomplished. They include sway and yaw motion simulation of the combination and prediction of maneuvering motion such as turning circles of different rudder angle under the condition of being pushed by two pushers. The coincidence of the prediction with the results of model test is good.展开更多
A new and computationally efficient version of the immersed boundary method,which is combined with the coarse-graining method,is introduced for modeling inextensible filaments immersed in low-Reynolds number flows.Thi...A new and computationally efficient version of the immersed boundary method,which is combined with the coarse-graining method,is introduced for modeling inextensible filaments immersed in low-Reynolds number flows.This is used to represent actin biopolymers,which are constituent elements of the cytoskeleton,a complex network-like structure that plays a fundamental role in shape morphology.An extension of the traditional immersed boundary method to include a stochastic stress tensor is also proposed in order to model the thermal fluctuations in the fluid at smaller scales.By way of validation,the response of a single,massless,inextensible semiflexible filament immersed in a thermally fluctuating fluid is obtained using the suggested numerical scheme and the resulting time-averaged contraction of the filament is compared to the theoretical value obtained from the worm-like chain model.展开更多
Most algorithms of the immersed boundary method originated by Peskin are explicit when it comes to the computation of the elastic forces exerted by the immersed boundary to the fluid. A drawback of such an explicit ap...Most algorithms of the immersed boundary method originated by Peskin are explicit when it comes to the computation of the elastic forces exerted by the immersed boundary to the fluid. A drawback of such an explicit approach is a severe restriction on the time step size for maintaining numerical stability. An implicit immersed boundary method in two dimensions using the lattice Boltzmann approach has been proposed. This paper reports an extension of the method to three dimensions and its application to simulation of a massive flexible sheet interacting with an incompressible viscous flow.展开更多
基金supported by the Research and Development Center of Transport Industry of New Generation of Artificial Intelligence Technology(Grant No.202202H)the National Key R&D Program of China(Grant No.2019YFB1600702)the National Natural Science Foundation of China(Grant Nos.51978600&51808336).
文摘Structural Health Monitoring(SHM)systems have become a crucial tool for the operational management of long tunnels.For immersed tunnels exposed to both traffic loads and the effects of the marine environment,efficiently identifying abnormal conditions from the extensive unannotated SHM data presents a significant challenge.This study proposed amodel-based approach for anomaly detection and conducted validation and comparative analysis of two distinct temporal predictive models using SHM data from a real immersed tunnel.Firstly,a dynamic predictive model-based anomaly detectionmethod is proposed,which utilizes a rolling time window for modeling to achieve dynamic prediction.Leveraging the assumption of temporal data similarity,an interval prediction value deviation was employed to determine the abnormality of the data.Subsequently,dynamic predictive models were constructed based on the Autoregressive Integrated Moving Average(ARIMA)and Long Short-Term Memory(LSTM)models.The hyperparameters of these models were optimized and selected using monitoring data from the immersed tunnel,yielding viable static and dynamic predictive models.Finally,the models were applied within the same segment of SHM data,to validate the effectiveness of the anomaly detection approach based on dynamic predictive modeling.A detailed comparative analysis discusses the discrepancies in temporal anomaly detection between the ARIMA-and LSTM-based models.The results demonstrated that the dynamic predictive modelbased anomaly detection approach was effective for dealing with unannotated SHM data.In a comparison between ARIMA and LSTM,it was found that ARIMA demonstrated higher modeling efficiency,rendering it suitable for short-term predictions.In contrast,the LSTM model exhibited greater capacity to capture long-term performance trends and enhanced early warning capabilities,thereby resulting in superior overall performance.
文摘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.
文摘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.
基金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.
基金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).
基金Project supported by the National Natural Science Foundation of China(Nos.11671157 and11826212)
文摘In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of finite element (FE) formulation and the simple structure of Cartesian grids,the IFE discretization is used in this paper.Two-grid schemes are formulated to linearize the FE equations.It is theoretically and numerically illustrated that the coarse space can be selected as coarse as H =O(h^1/4)(or H =O(h^1/8)),and the asymptotically optimal approximation can be achieved as the nonlinear schemes.As a result,we can settle a great majority of nonlinear equations as easy as linearized problems.In order to estimate the present two-grid algorithms,we derive the optimal error estimates of the IFE solution in the L^p norm.Numerical experiments are given to verify the theorems and indicate that the present two-grid algorithms can greatly improve the computing efficiency.
基金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.
文摘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.
文摘A micro-layer model is proposed to account for the lubrication effect of liquid layer near collisions of immersed particles at moderate particle Reynolds number.This new model is to allow determination of the pressure profile within the micro-layer including the fluid inertia and viscosity.Then a correction based on the micro-layer model is applied to unsteady 3-D direct simulation of a particle approaching another one.The simulation is based on a modified immersed boundary method with direct force scheme.The quantitative agreement between numerical and experimental results validates the model presented in the study.The simulation results show that the fluid is squeezed prior to contact.When a particle approaches a flat wall or another particle,the lubrication force,obtained by integrating the pressure profile over the particle surface,is increased and prevents the particle from approaching.The model predicts that the velocity of approaching particle starts to decrease when separation distance of particles is less than 0.1dp,where dp is the particle diameter.
基金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.
基金financially supported by the Ministry of Education of Humanities and Social Science Project(Grant Nos.15YJC630145 and 15YJC630059)the Natural Science Foundation of Shanghai Science and Technology Committee(Grant No.15ZR1420200)
文摘Immersed tunnel is an important part of the Hong Kong–Zhuhai–Macao Bridge(HZMB) project. In immersed tunnel floating, translation which includes straight and transverse movements is the main working mode. To decide the magnitude and direction of the towing force for each tug, a particle swarm-based translation control method is presented for non-power immersed tunnel element. A sort of linear weighted logarithmic function is exploited to avoid weak subgoals. In simulation, the particle swarm-based control method is evaluated and compared with traditional empirical method in the case of the HZMB project. Simulation results show that the presented method delivers performance improvement in terms of the enhanced surplus towing force.
基金Supported by National Natural Science Foundation of China(Grant No.51405375)National Key Basic Research and Development Program of China(973 Program,Grant No.2011CB706606)
文摘Currently, many studies on the local discontinuous Galerkin method focus on the Cartesian grid with low computational e ciency and poor adaptability to complex shapes. A new immersed boundary method is presented, and this method employs the adaptive Cartesian grid to improve the adaptability to complex shapes and the immersed boundary to increase computational e ciency. The new immersed boundary method employs different boundary cells(the physical cell and ghost cell) to impose the boundary condition and the reconstruction algorithm of the ghost cell is the key for this method. The classical model elliptic equation is used to test the method. This method is tested and analyzed from the viewpoints of boundary cell type, error distribution and accuracy. The numerical result shows that the presented method has low error and a good rate of the convergence and works well in complex geometries. The method has good prospect for practical application research of the numerical calculation research.
文摘A method combining the immersed boundary technique and a multirelaxation-time(MRT) lattice Boltzmann flux solver(LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic cylinders and NACA 0012 Airfoil. The method uses a simple Cartesian mesh to simulate flows past immersed complicated bodies. With the Chapman-Enskog expansion analysis, a transform is performed between the Navier-Stokes and lattice Boltzmann equations(LBEs). The LBFS is used to discretize the macroscopic differential equations with a finite volume method and evaluate the interface fluxes through local reconstruction of the lattice Boltzmann solution.The immersed boundary technique is used to correct the intermediate velocity around the solid boundary to satisfy the no-slip boundary condition. Agreement of simulation results with the data found in the literature shows reliability of the proposed method in simulating laminar flows on a Cartesian mesh.
文摘The Immersed Interface Method (IIM) is derived to solve the corresponding Fokker-Planck equation of Brownian motion with pure dry friction, which is one of the simplest models of piecewise-smooth stochastic systems. The IIM is capable of treating a discontinuity in the drift of Fokker-Planck equation and it is readily extended to the dry and viscous friction model. Analytic results of the considered model are used to confirm the effectiveness and design accuracy of the scheme.
文摘We present a cut-cell method for the simulation of 2D incompressible flows past obstacles.It consists in using the MAC scheme on cartesian grids and imposing Dirchlet boundary conditions for the velocity field on the boundary of solid structures following the Shortley-Weller formulation.In order to ensure local conservation properties,viscous and convecting terms are discretized in a finite volume way.The scheme is second order implicit in time for the linear part,the linear systems are solved by the use of the capacitance matrix method for non-moving obstacles.Numerical results of flows around an impulsively started circular cylinder are presented which confirm the efficiency of the method,for Reynolds numbers 1000 and 3000.An example of flows around a moving rigid body at Reynolds number 800 is also shown,a solver using the PETSc-Library has been prefered in this context to solve the linear systems.
基金This work was supported by the National Natural Science Foundation of China the Science and Technology Foundation of Hunan P
文摘The electrochemical quartz crystal impedance system (EQCIS) has been used for the study of a partially immersed An electrode in 0.2 mol/L NaCIO4. aqueous solution. The influences of the immersed area and height of the electrode on the EQCIS responses were evaluated, showing the highest response sensitivity to liquid loading at the center of the piezoelectric quartz crystal electrode. The increase in the immersed height of the An electrode at oxygen reduction potentials during potential cycling was measured by this technique.
文摘In this paper, the transportation mode of tube during the immersed tunnel engineering in Yangtze River incorporated, mathematical models of the combination of the tube and four barges are established when it is towed by a tugboat through an elastic rope and pushed by two pushers. Based on evaluating the hydrodynamic derivatives by regressive formula, the simulations and analysis of the whole maneuvering motion course are accomplished. They include sway and yaw motion simulation of the combination and prediction of maneuvering motion such as turning circles of different rudder angle under the condition of being pushed by two pushers. The coincidence of the prediction with the results of model test is good.
文摘A new and computationally efficient version of the immersed boundary method,which is combined with the coarse-graining method,is introduced for modeling inextensible filaments immersed in low-Reynolds number flows.This is used to represent actin biopolymers,which are constituent elements of the cytoskeleton,a complex network-like structure that plays a fundamental role in shape morphology.An extension of the traditional immersed boundary method to include a stochastic stress tensor is also proposed in order to model the thermal fluctuations in the fluid at smaller scales.By way of validation,the response of a single,massless,inextensible semiflexible filament immersed in a thermally fluctuating fluid is obtained using the suggested numerical scheme and the resulting time-averaged contraction of the filament is compared to the theoretical value obtained from the worm-like chain model.
基金supported by the US National Science Foundation (DMS-0713718)
文摘Most algorithms of the immersed boundary method originated by Peskin are explicit when it comes to the computation of the elastic forces exerted by the immersed boundary to the fluid. A drawback of such an explicit approach is a severe restriction on the time step size for maintaining numerical stability. An implicit immersed boundary method in two dimensions using the lattice Boltzmann approach has been proposed. This paper reports an extension of the method to three dimensions and its application to simulation of a massive flexible sheet interacting with an incompressible viscous flow.