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.展开更多
For eliminating the zero-order image in digital holography, a new method using the differential of the hologram intensity instead of the hologram itself for numerical reconstruction is proposed. This method is based o...For eliminating the zero-order image in digital holography, a new method using the differential of the hologram intensity instead of the hologram itself for numerical reconstruction is proposed. This method is based on digital image processing. By analyzing the spatial spectrum of the off-axis digital hologram, it theoretically proves that the zero-order image can be effectively eliminated by differential before reconstruction. Then, the detected hologram is processed in the program with differential and reconstruction. Both the theoretical analysis and digital reconstruction results show that it can effectively eliminate the large bright spot in the center of the reconstructed image caused by the zero-order image, improve the image quality significantly, and render a better contrast of the reconstructed image. This method is very simple and convenient due to no superfluous optical elements and requiring only one time record.展开更多
In this research, we have improved a relaxation method for triangular meshes intended for finite element fluid simulations which contain discrete element particles. The triangle edges are treated as springs which rela...In this research, we have improved a relaxation method for triangular meshes intended for finite element fluid simulations which contain discrete element particles. The triangle edges are treated as springs which relax their lengths towards a “better” force equilibrium where the triangles are closer to equilateral shape. The actual kernel is an improved zero order integrator which is able to follow reconfigurations of the particles faster than earlier methods. The improved relaxation allows larger timesteps in the flow simulation and leads to more stable, faster mesh reconfigurations for fast moving particles in the flow. Additionally, this demonstrates how integrators of the same order zero can nevertheless have different convergence speeds towards展开更多
Shell structures have increasingly widespread applications in biomedical ultrasound fields such as contrast agents and drug delivery,which requires the precise prediction of the acoustic radiation force under various ...Shell structures have increasingly widespread applications in biomedical ultrasound fields such as contrast agents and drug delivery,which requires the precise prediction of the acoustic radiation force under various circumstances to improve the system efficiency.The acoustic radiation force exerted by a zero-order quasi-Bessel-Gauss beam on an elastic spherical shell near an impedance boundary is theoretically and numerically studied in this study.By means of the finite series method and the image theory,a zero-order quasi-Bessel-Gauss beam is expanded in terms of spherical harmonic functions,and the exact solution of the acoustic radiation force is derived based on the acoustic scattering theory.The acoustic radiation force function,which represents the radiation force per unit energy density and per unit cross-sectional surface,is especially investigated.Some simulated results for a polymethyl methacrylate shell and an aluminum shell are provided to illustrate the behavior of acoustic radiation force in this case.The simulated results show the oscillatory property and the negative radiation force caused by the impedance boundary.An appropriate relative thickness of the shell can generate sharp peaks for a polymethyl methacrylate shell.Strong radiation force can be obtained at small half-cone angles and the beam waist only affects the results at high frequencies.Considering that the quasi-Bessel-Gauss beam possesses both the energy focusing property and the non-diffracting advantage,this study is expected to be useful in the development of acoustic tweezers,contrast agent micro-shells,and drug delivery applications.展开更多
With derivative spectrophotometric method,the complexes of lanthanides with 2-(5-bromo-2-pyridylazo)-5-diethylaminophen- ol(5-Br-PADAP)in the presence of octylphenol poly(ethyleneg- lycol)ether(TX-100)were studied.It ...With derivative spectrophotometric method,the complexes of lanthanides with 2-(5-bromo-2-pyridylazo)-5-diethylaminophen- ol(5-Br-PADAP)in the presence of octylphenol poly(ethyleneg- lycol)ether(TX-100)were studied.It is found that the maximum absorption of fourth-order derivative spectra of the neodymium complex by 4f electron transitions is at 579(+)nm and 582(-)nm with molar derivative absorptivities of 5.2×10~3 1.mol^(-1).cm^(-1). The maximum absorption of the zero-order derivative spectra of the complexes for neodymium and praseodymium is at 580hm,the molar absorptivities are 1.47×10~5 l.mol^(-1).cm^(-1)for Nd and 1.26 ×10~5 l.mol^(-1).cm^(-1)for Pr.The component ratio of the complex is Nd:5-Br-PADAP=1:4.Beer's law is obeyed in the range of 0-3.75μg of Nd and Pr in 25ml of solution.The method has been used for the simultaneous determination of Nd and Pr in four synthetic samples and the results obtained are satisfactory.展开更多
The aim of this study is to create a fast and stable iterative technique for numerical solution of a quasi-linear elliptic pressure equation. We developed a modified version of the Anderson acceleration(AA)algorithm t...The aim of this study is to create a fast and stable iterative technique for numerical solution of a quasi-linear elliptic pressure equation. We developed a modified version of the Anderson acceleration(AA)algorithm to fixed-point(FP) iteration method. It computes the approximation to the solutions at each iteration based on the history of vectors in extended space, which includes the vector of unknowns, the discrete form of the operator, and the equation's right-hand side. Several constraints are applied to AA algorithm, including a limitation of the time step variation during the iteration process, which allows switching to the base FP iterations to maintain convergence. Compared to the base FP algorithm, the improved version of the AA algorithm enables a reliable and rapid convergence of the iterative solution for the quasi-linear elliptic pressure equation describing the flow of particle-laden yield-stress fluids in a narrow channel during hydraulic fracturing, a key technology for stimulating hydrocarbon-bearing reservoirs. In particular, the proposed AA algorithm allows for faster computations and resolution of unyielding zones in hydraulic fractures that cannot be calculated using the FP algorithm. The quasi-linear elliptic pressure equation under consideration describes various physical processes, such as the displacement of fluids with viscoplastic rheology in a narrow cylindrical annulus during well cementing,the displacement of cross-linked gel in a proppant pack filling hydraulic fractures during the early stage of well production(fracture flowback), and multiphase filtration in a rock formation. We estimate computational complexity of the developed algorithm as compared to Jacobian-based algorithms and show that the performance of the former one is higher in modelling of flows of viscoplastic fluids. We believe that the developed algorithm is a useful numerical tool that can be implemented in commercial simulators to obtain fast and converged solutions to the non-linear problems described above.展开更多
The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially i...The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.展开更多
以Aspen Open Solver接口集中的非线性代数方程组(NLA)部分作为研究对象,在对接口集进行系统地分析之后,利用AspenTech提供的接口代码将分别基于梯度和非基于梯度的四种求解算法嵌入生成solver组件,并实现用Aspen Plus调用该solver组件...以Aspen Open Solver接口集中的非线性代数方程组(NLA)部分作为研究对象,在对接口集进行系统地分析之后,利用AspenTech提供的接口代码将分别基于梯度和非基于梯度的四种求解算法嵌入生成solver组件,并实现用Aspen Plus调用该solver组件观察各种算法嵌入的结果。展开更多
基金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.
基金The Natural Science Foundation of Jiangsu Province (No.BK2006102)the National Natural Science Foundation of China(No.10772086)
文摘For eliminating the zero-order image in digital holography, a new method using the differential of the hologram intensity instead of the hologram itself for numerical reconstruction is proposed. This method is based on digital image processing. By analyzing the spatial spectrum of the off-axis digital hologram, it theoretically proves that the zero-order image can be effectively eliminated by differential before reconstruction. Then, the detected hologram is processed in the program with differential and reconstruction. Both the theoretical analysis and digital reconstruction results show that it can effectively eliminate the large bright spot in the center of the reconstructed image caused by the zero-order image, improve the image quality significantly, and render a better contrast of the reconstructed image. This method is very simple and convenient due to no superfluous optical elements and requiring only one time record.
文摘In this research, we have improved a relaxation method for triangular meshes intended for finite element fluid simulations which contain discrete element particles. The triangle edges are treated as springs which relax their lengths towards a “better” force equilibrium where the triangles are closer to equilateral shape. The actual kernel is an improved zero order integrator which is able to follow reconfigurations of the particles faster than earlier methods. The improved relaxation allows larger timesteps in the flow simulation and leads to more stable, faster mesh reconfigurations for fast moving particles in the flow. Additionally, this demonstrates how integrators of the same order zero can nevertheless have different convergence speeds towards
基金Project supported by the National Natural Science Foundation of China(Grant Nos.81527901,11604361,and 91630309)。
文摘Shell structures have increasingly widespread applications in biomedical ultrasound fields such as contrast agents and drug delivery,which requires the precise prediction of the acoustic radiation force under various circumstances to improve the system efficiency.The acoustic radiation force exerted by a zero-order quasi-Bessel-Gauss beam on an elastic spherical shell near an impedance boundary is theoretically and numerically studied in this study.By means of the finite series method and the image theory,a zero-order quasi-Bessel-Gauss beam is expanded in terms of spherical harmonic functions,and the exact solution of the acoustic radiation force is derived based on the acoustic scattering theory.The acoustic radiation force function,which represents the radiation force per unit energy density and per unit cross-sectional surface,is especially investigated.Some simulated results for a polymethyl methacrylate shell and an aluminum shell are provided to illustrate the behavior of acoustic radiation force in this case.The simulated results show the oscillatory property and the negative radiation force caused by the impedance boundary.An appropriate relative thickness of the shell can generate sharp peaks for a polymethyl methacrylate shell.Strong radiation force can be obtained at small half-cone angles and the beam waist only affects the results at high frequencies.Considering that the quasi-Bessel-Gauss beam possesses both the energy focusing property and the non-diffracting advantage,this study is expected to be useful in the development of acoustic tweezers,contrast agent micro-shells,and drug delivery applications.
文摘With derivative spectrophotometric method,the complexes of lanthanides with 2-(5-bromo-2-pyridylazo)-5-diethylaminophen- ol(5-Br-PADAP)in the presence of octylphenol poly(ethyleneg- lycol)ether(TX-100)were studied.It is found that the maximum absorption of fourth-order derivative spectra of the neodymium complex by 4f electron transitions is at 579(+)nm and 582(-)nm with molar derivative absorptivities of 5.2×10~3 1.mol^(-1).cm^(-1). The maximum absorption of the zero-order derivative spectra of the complexes for neodymium and praseodymium is at 580hm,the molar absorptivities are 1.47×10~5 l.mol^(-1).cm^(-1)for Nd and 1.26 ×10~5 l.mol^(-1).cm^(-1)for Pr.The component ratio of the complex is Nd:5-Br-PADAP=1:4.Beer's law is obeyed in the range of 0-3.75μg of Nd and Pr in 25ml of solution.The method has been used for the simultaneous determination of Nd and Pr in four synthetic samples and the results obtained are satisfactory.
基金partial financial support from Gazpromneft Science and Technology Center。
文摘The aim of this study is to create a fast and stable iterative technique for numerical solution of a quasi-linear elliptic pressure equation. We developed a modified version of the Anderson acceleration(AA)algorithm to fixed-point(FP) iteration method. It computes the approximation to the solutions at each iteration based on the history of vectors in extended space, which includes the vector of unknowns, the discrete form of the operator, and the equation's right-hand side. Several constraints are applied to AA algorithm, including a limitation of the time step variation during the iteration process, which allows switching to the base FP iterations to maintain convergence. Compared to the base FP algorithm, the improved version of the AA algorithm enables a reliable and rapid convergence of the iterative solution for the quasi-linear elliptic pressure equation describing the flow of particle-laden yield-stress fluids in a narrow channel during hydraulic fracturing, a key technology for stimulating hydrocarbon-bearing reservoirs. In particular, the proposed AA algorithm allows for faster computations and resolution of unyielding zones in hydraulic fractures that cannot be calculated using the FP algorithm. The quasi-linear elliptic pressure equation under consideration describes various physical processes, such as the displacement of fluids with viscoplastic rheology in a narrow cylindrical annulus during well cementing,the displacement of cross-linked gel in a proppant pack filling hydraulic fractures during the early stage of well production(fracture flowback), and multiphase filtration in a rock formation. We estimate computational complexity of the developed algorithm as compared to Jacobian-based algorithms and show that the performance of the former one is higher in modelling of flows of viscoplastic fluids. We believe that the developed algorithm is a useful numerical tool that can be implemented in commercial simulators to obtain fast and converged solutions to the non-linear problems described above.
基金the National Supercomputer Center in Tianjin for their patient assistance in providing the compilation environment.We thank the editor,Huajian Yao,for handling the manuscript and Mingming Li and another anonymous reviewer for their constructive comments.The research leading to these results has received funding from National Natural Science Foundation of China projects(Grant Nos.92355302 and 42121005)Taishan Scholar projects(Grant No.tspd20210305)others(Grant Nos.XDB0710000,L2324203,XK2023DXC001,LSKJ202204400,and ZR2021ZD09).
文摘The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.
