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.展开更多
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.展开更多
We propose a novel type of nonlinear solver acceleration for systems of nonlinear partial differential equations(PDEs)that is based on online/adaptive learning.It is applied in the context of multiphase flow in porous...We propose a novel type of nonlinear solver acceleration for systems of nonlinear partial differential equations(PDEs)that is based on online/adaptive learning.It is applied in the context of multiphase flow in porous media.The proposed method rely on four pillars:(i)dimensionless numbers as input parameters for the machine learning model,(ii)simplified numerical model(two-dimensional)for the offline training,(iii)dynamic control of a nonlinear solver tuning parameter(numerical relaxation),(iv)and online learning for time real-improvement of the machine learning model.This strategy decreases the number of nonlinear iterations by dynamically modifying a single global parameter,the relaxation factor,and by adaptively learning the attributes of each numerical model on-the-run.Furthermore,this work performs a sensitivity study in the dimensionless parameters(machine learning features),assess the efficacy of various machine learning models,demonstrate a decrease in nonlinear iterations using our method in more intricate,realistic three-dimensional models,and fully couple a machine learning model into an open-source multiphase flow simulator achieving up to 85%reduction in computational time.展开更多
以Aspen Open Solver接口集中的非线性代数方程组(NLA)部分作为研究对象,在对接口集进行系统地分析之后,利用AspenTech提供的接口代码将分别基于梯度和非基于梯度的四种求解算法嵌入生成solver组件,并实现用Aspen Plus调用该solver组件...以Aspen Open Solver接口集中的非线性代数方程组(NLA)部分作为研究对象,在对接口集进行系统地分析之后,利用AspenTech提供的接口代码将分别基于梯度和非基于梯度的四种求解算法嵌入生成solver组件,并实现用Aspen Plus调用该solver组件观察各种算法嵌入的结果。展开更多
In this paper,a wave generating approach for long-crest irregular waves in a numerical tank by our in-house solver naoe-FOAM-SJTU is presented.The naoe-FOAM-SJTU solver is developed using an open source tool kit,Open ...In this paper,a wave generating approach for long-crest irregular waves in a numerical tank by our in-house solver naoe-FOAM-SJTU is presented.The naoe-FOAM-SJTU solver is developed using an open source tool kit,Open FOAM.Reynolds-averaged Navier？Stokes（RANS） equations are chosen as governing equations and the volume of fluid（VOF） is employed to capture the two phases interface.Incoming wave group is generated by imposing the boundary conditions of the tank inlet.A spectrum based correction procedure is developed to make the measured spectrum approaching to the target spectrum.This procedure can automatically adjust the wave generation signal based on the measured wave elevation by wave height probe in numerical wave tank.After 3 to 4 iterations,the measured spectrum agrees well with the target one.In order to validate this method,several wave spectra are chosen and validated in the numerical wave tank,with comparison between the final measured and target spectra.In order to investigate a practical situation,a modified Wigley hull is placed in the wave tank with incoming irregular waves.The wave-induced heave and pitch motions are treated by Fourier analysis to obtain motion responses,showing good agreements with the measurements.展开更多
A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the gov...A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.展开更多
The elliptic mild slope equation is used to simulate linear wave propagation over variable sea bed topography with mild slopes. The governing equation is discretized by the finite difference method. Based on the BI-CG...The elliptic mild slope equation is used to simulate linear wave propagation over variable sea bed topography with mild slopes. The governing equation is discretized by the finite difference method. Based on the BI-CGSTAB technique, an attractive variant bf BI-Conjugate Gradients (BI-CG) method, the obtained linear algebraic system of equations is solved. Numerical experiments show that the BI-CGSTAB method is efficient for solving the elliptic mild slope equation. The results obtained by the BI-CGSTAB-Based method are much the same as those obtained by other authors with different solution methods, but the convergence rate is much faster than that of other methods.展开更多
Simulation of solitary wave run-up on a vertical circular cylinder is carried out in a viscous numerical wave tank developed based on the open source codes Open FOAM. An incompressible two-phase flow solver naoe-FOAM-...Simulation of solitary wave run-up on a vertical circular cylinder is carried out in a viscous numerical wave tank developed based on the open source codes Open FOAM. An incompressible two-phase flow solver naoe-FOAM-SJTU is used to solve the Reynolds-Averaged Navier–Stokes（RANS） equations with the SST k ？？ turbulence model. The PISO algorithm is utilized for the pressure-velocity coupling. The air-water interface is captured via Volume of Fluid（VOF） technique. The present numerical model is validated by simulating the solitary wave run-up and reflected against a vertical wall, and solitary wave run-up on a vertical circular cylinder. Comparisons between numerical results and available experimental data show satisfactory agreement. Furthermore, simulations are carried out to study the solitary wave run-up on the cylinder with different incident wave height H and different cylinder radius a. The relationships of the wave run-up height with the incident wave height H, cylinder radius a are analyzed. The evolutions of the scattering free surface and vortex shedding are also presented to give a better understanding of the process of nonlinear wave-cylinder interaction.展开更多
A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver is built with elastic waves （HLLCE） for one-dimensional elastic-plastic flows with a hypo- elastic constitutive model and the von Mises＇ yielding cr...A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver is built with elastic waves （HLLCE） for one-dimensional elastic-plastic flows with a hypo- elastic constitutive model and the von Mises＇ yielding criterion. Based on the HLLCE, a third-order cell-centered Lagrangian scheme is built for one-dimensional elastic-plastic problems. A number of numerical experiments are carried out. The numerical results show that the proposed third-order scheme achieves the desired order of accuracy. The third-order scheme is used to the numerical solution of the problems with elastic shock waves and elastic rarefaction waves. The numerical results are compared with a reference solution and the results obtained by other authors. The comparison shows that the pre- sented high-order scheme is convergent, stable, and essentially non-oscillatory. Moreover, the HLLCE is more efficient than the two-rarefaction Riemann solver with elastic waves （TRRSE）展开更多
In this paper,the finite difference weighted essentially non-oscillatory (WENO) scheme is incorporated into the recently developed four kinds of lattice Boltzmann flux solver (LBFS) to simulate compressible flows,incl...In this paper,the finite difference weighted essentially non-oscillatory (WENO) scheme is incorporated into the recently developed four kinds of lattice Boltzmann flux solver (LBFS) to simulate compressible flows,including inviscid LBFS Ⅰ,viscous LBFS Ⅱ,hybrid LBFS Ⅲ and hybrid LBFS Ⅳ.Hybrid LBFS can automatically realize the switch between inviscid LBFS Ⅰ and viscous LBFS Ⅱ through introducing a switch function.The resultant hybrid WENO-LBFS scheme absorbs the advantages of WENO scheme and hybrid LBFS.We investigate the performance of WENO scheme based on four kinds of LBFS systematically.Numerical results indicate that the devopled hybrid WENO-LBFS scheme has high accuracy,high resolution and no oscillations.It can not only accurately calculate smooth solutions,but also can effectively capture contact discontinuities and strong shock waves.展开更多
In this paper,a hybrid Lattice Boltzmann Flux Solver(LBFS)with an improved switch function is proposed for simulation of integrated hypersonic fluid-thermal-structural problems.In the solver,the macroscopic Navier–St...In this paper,a hybrid Lattice Boltzmann Flux Solver(LBFS)with an improved switch function is proposed for simulation of integrated hypersonic fluid-thermal-structural problems.In the solver,the macroscopic Navier–Stokes equations and structural heat transfer equation are discretized by the finite volume method,and the numerical fluxes at the cell interface are reconstructed by the local solution of the Boltzmann equation.To compute the numerical fluxes,two equilibrium distribution functions are introduced.One is the D1Q4 discrete velocity model for calculating the inviscid flux across the cell interface of Navier–Stokes equations,and the other is the D2Q4 model for evaluating the flux of structural energy equation.In this work,a new dual thermal resistance model is proposed to calculate the thermal properties at the fluid–solid interface.The accuracy and stability of the present hybrid solver are validated by simulating several numerical examples,including the fluid-thermal-structural problem of cylindrical leading edge.Numerical results show that the present solver can accurately predict the thermal properties of hypersonic fluid-thermalstructural problems and has the great potential for solving fluid-thermal-structural problems of long-endurance high-speed vehicles.展开更多
Inspired by the recently proposed Legendre orthogonal polynomial representation for imaginary-time Green s functions G（τ）,we develop an alternate and superior representation for G（τ） and implement it in the hybr...Inspired by the recently proposed Legendre orthogonal polynomial representation for imaginary-time Green s functions G（τ）,we develop an alternate and superior representation for G（τ） and implement it in the hybridization expansion continuous-time quantum Monte Carlo impurity solver.This representation is based on the kernel polynomial method,which introduces some integral kernel functions to filter the numerical fluctuations caused by the explicit truncations of polynomial expansion series and can improve the computational precision significantly.As an illustration of the new representation,we re-examine the imaginary-time Green＇s functions of the single-band Hubbard model in the framework of dynamical mean-field theory.The calculated results suggest that with carefully chosen integral kernel functions,whether the system is metallic or insulating,the Gibbs oscillations found in the previous Legendre orthogonal polynomial representation have been vastly suppressed and remarkable corrections to the measured Green＇s functions have been obtained.展开更多
基金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.
基金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.
基金MUFFINS,MUltiphase Flow-induced Fluid-flexible structure InteractioN in Subsea applications(EP/P033180/1)the PREMIERE programme grant(EP/T000414/1)SMARTRES,Smart assessment,management and optimization of urban geothermal resources(NE/X005607/1).
文摘We propose a novel type of nonlinear solver acceleration for systems of nonlinear partial differential equations(PDEs)that is based on online/adaptive learning.It is applied in the context of multiphase flow in porous media.The proposed method rely on four pillars:(i)dimensionless numbers as input parameters for the machine learning model,(ii)simplified numerical model(two-dimensional)for the offline training,(iii)dynamic control of a nonlinear solver tuning parameter(numerical relaxation),(iv)and online learning for time real-improvement of the machine learning model.This strategy decreases the number of nonlinear iterations by dynamically modifying a single global parameter,the relaxation factor,and by adaptively learning the attributes of each numerical model on-the-run.Furthermore,this work performs a sensitivity study in the dimensionless parameters(machine learning features),assess the efficacy of various machine learning models,demonstrate a decrease in nonlinear iterations using our method in more intricate,realistic three-dimensional models,and fully couple a machine learning model into an open-source multiphase flow simulator achieving up to 85%reduction in computational time.
基金financially supported by National Natural Science Foundation of China(Grant Nos.51379125,51411130131,11432009,and 51490675)the Chang Jiang Scholars Program(Grant No.T2014099)+3 种基金the Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China(Grant No.2016-23)the Foundation of State key Laboratory of Ocean Engineering(Grant No.GKZD010065)the Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(Grant No.2013022)center for HPC at Shanghai Jiao Tong University,and Lloyd’s Register Foundation(LRF)
文摘In this paper,a wave generating approach for long-crest irregular waves in a numerical tank by our in-house solver naoe-FOAM-SJTU is presented.The naoe-FOAM-SJTU solver is developed using an open source tool kit,Open FOAM.Reynolds-averaged Navier？Stokes（RANS） equations are chosen as governing equations and the volume of fluid（VOF） is employed to capture the two phases interface.Incoming wave group is generated by imposing the boundary conditions of the tank inlet.A spectrum based correction procedure is developed to make the measured spectrum approaching to the target spectrum.This procedure can automatically adjust the wave generation signal based on the measured wave elevation by wave height probe in numerical wave tank.After 3 to 4 iterations,the measured spectrum agrees well with the target one.In order to validate this method,several wave spectra are chosen and validated in the numerical wave tank,with comparison between the final measured and target spectra.In order to investigate a practical situation,a modified Wigley hull is placed in the wave tank with incoming irregular waves.The wave-induced heave and pitch motions are treated by Fourier analysis to obtain motion responses,showing good agreements with the measurements.
基金Supported by the National Natural Science Foundation of China(11272153)
文摘A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.
基金National Natural Science Foundation of China under grant No.59839330China Postdoctoral Science Foundation
文摘The elliptic mild slope equation is used to simulate linear wave propagation over variable sea bed topography with mild slopes. The governing equation is discretized by the finite difference method. Based on the BI-CGSTAB technique, an attractive variant bf BI-Conjugate Gradients (BI-CG) method, the obtained linear algebraic system of equations is solved. Numerical experiments show that the BI-CGSTAB method is efficient for solving the elliptic mild slope equation. The results obtained by the BI-CGSTAB-Based method are much the same as those obtained by other authors with different solution methods, but the convergence rate is much faster than that of other methods.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51379125,51411130131,and 11432009)Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(Grant No.2013022)the National Basic Research Program of China(973 Program,Grant No.2013CB036103)
文摘Simulation of solitary wave run-up on a vertical circular cylinder is carried out in a viscous numerical wave tank developed based on the open source codes Open FOAM. An incompressible two-phase flow solver naoe-FOAM-SJTU is used to solve the Reynolds-Averaged Navier–Stokes（RANS） equations with the SST k ？？ turbulence model. The PISO algorithm is utilized for the pressure-velocity coupling. The air-water interface is captured via Volume of Fluid（VOF） technique. The present numerical model is validated by simulating the solitary wave run-up and reflected against a vertical wall, and solitary wave run-up on a vertical circular cylinder. Comparisons between numerical results and available experimental data show satisfactory agreement. Furthermore, simulations are carried out to study the solitary wave run-up on the cylinder with different incident wave height H and different cylinder radius a. The relationships of the wave run-up height with the incident wave height H, cylinder radius a are analyzed. The evolutions of the scattering free surface and vortex shedding are also presented to give a better understanding of the process of nonlinear wave-cylinder interaction.
基金Project supported by the National Natural Science Foundation of China(Nos.11172050 and11672047)the Science and Technology Foundation of China Academy of Engineering Physics(No.2013A0202011)
文摘A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver is built with elastic waves （HLLCE） for one-dimensional elastic-plastic flows with a hypo- elastic constitutive model and the von Mises＇ yielding criterion. Based on the HLLCE, a third-order cell-centered Lagrangian scheme is built for one-dimensional elastic-plastic problems. A number of numerical experiments are carried out. The numerical results show that the proposed third-order scheme achieves the desired order of accuracy. The third-order scheme is used to the numerical solution of the problems with elastic shock waves and elastic rarefaction waves. The numerical results are compared with a reference solution and the results obtained by other authors. The comparison shows that the pre- sented high-order scheme is convergent, stable, and essentially non-oscillatory. Moreover, the HLLCE is more efficient than the two-rarefaction Riemann solver with elastic waves （TRRSE）
基金This study was supported by the National Natural Science Foundation of China(Grants 11372168,11772179).
文摘In this paper,the finite difference weighted essentially non-oscillatory (WENO) scheme is incorporated into the recently developed four kinds of lattice Boltzmann flux solver (LBFS) to simulate compressible flows,including inviscid LBFS Ⅰ,viscous LBFS Ⅱ,hybrid LBFS Ⅲ and hybrid LBFS Ⅳ.Hybrid LBFS can automatically realize the switch between inviscid LBFS Ⅰ and viscous LBFS Ⅱ through introducing a switch function.The resultant hybrid WENO-LBFS scheme absorbs the advantages of WENO scheme and hybrid LBFS.We investigate the performance of WENO scheme based on four kinds of LBFS systematically.Numerical results indicate that the devopled hybrid WENO-LBFS scheme has high accuracy,high resolution and no oscillations.It can not only accurately calculate smooth solutions,but also can effectively capture contact discontinuities and strong shock waves.
基金co-supported by the Postgraduate Research&Practice Innovation Program of Jiangsu Province of China(No.KYCX17_0235)Nanjing University of Aeronautics and Astronautics Ph D Short-term Visiting Scholar Project(No.180602DF01)National Numerical Wind Tunnel Project(Nos.NNW2018-ZT3B08 and NNW2019-ZT7B30)。
文摘In this paper,a hybrid Lattice Boltzmann Flux Solver(LBFS)with an improved switch function is proposed for simulation of integrated hypersonic fluid-thermal-structural problems.In the solver,the macroscopic Navier–Stokes equations and structural heat transfer equation are discretized by the finite volume method,and the numerical fluxes at the cell interface are reconstructed by the local solution of the Boltzmann equation.To compute the numerical fluxes,two equilibrium distribution functions are introduced.One is the D1Q4 discrete velocity model for calculating the inviscid flux across the cell interface of Navier–Stokes equations,and the other is the D2Q4 model for evaluating the flux of structural energy equation.In this work,a new dual thermal resistance model is proposed to calculate the thermal properties at the fluid–solid interface.The accuracy and stability of the present hybrid solver are validated by simulating several numerical examples,including the fluid-thermal-structural problem of cylindrical leading edge.Numerical results show that the present solver can accurately predict the thermal properties of hypersonic fluid-thermalstructural problems and has the great potential for solving fluid-thermal-structural problems of long-endurance high-speed vehicles.
基金Project supported by the National Natural Science Foundation of China(Grant No.11504340)
文摘Inspired by the recently proposed Legendre orthogonal polynomial representation for imaginary-time Green s functions G（τ）,we develop an alternate and superior representation for G（τ） and implement it in the hybridization expansion continuous-time quantum Monte Carlo impurity solver.This representation is based on the kernel polynomial method,which introduces some integral kernel functions to filter the numerical fluctuations caused by the explicit truncations of polynomial expansion series and can improve the computational precision significantly.As an illustration of the new representation,we re-examine the imaginary-time Green＇s functions of the single-band Hubbard model in the framework of dynamical mean-field theory.The calculated results suggest that with carefully chosen integral kernel functions,whether the system is metallic or insulating,the Gibbs oscillations found in the previous Legendre orthogonal polynomial representation have been vastly suppressed and remarkable corrections to the measured Green＇s functions have been obtained.
