A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C-unstructured non-orthogonal grid in the hor- izontal d...A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C-unstructured non-orthogonal grid in the hor- izontal direction and z-level grid in the vertical direction. The three-dimensional model is discretized by the semi-implicit finite-volume method, in that the free-surface and the vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly by an integral method. The partial cell method is used for resolving topography, which enables the model to better represent irregular topography. The model has been tested against analytical cases for wind and tidal oscillation circulation, and is applied to simulating the tidal flow in the Bohal Sea. The results are in good agreement both with the analytical solutions and measurement results.展开更多
This paper presents a parallel method for simulating real-time 3D deformable objects using the volume preservation mass-spring system method on tetrahedron meshes.In general,the conventional mass-spring system is mani...This paper presents a parallel method for simulating real-time 3D deformable objects using the volume preservation mass-spring system method on tetrahedron meshes.In general,the conventional mass-spring system is manipulated as a force-driven method because it is fast,simple to implement,and the parameters can be controlled.However,the springs in traditional mass-spring system can be excessively elongated which cause severe stability and robustness issues that lead to shape restoring,simulation blow-up,and huge volume loss of the deformable object.In addition,traditional method that uses a serial process of the central processing unit(CPU)to solve the system in every frame cannot handle the complex structure of deformable object in real-time.Therefore,the first order implicit constraint enforcement for a mass-spring model is utilized to achieve accurate visual realism of deformable objects with tough constraint error.In this paper,we applied the distance constraint and volume conservation constraints for each tetrahedron element to improve the stability of deformable object simulation using the mass-spring system and behave the same as its real-world counterparts.To reduce the computational complexity while ensuring stable simulation,we applied a method that utilizes OpenGL compute shader,a part of OpenGL Shading Language(GLSL)that executes on the graphic processing unit(GPU)to solve the numerical problems effectively.We applied the proposed methods to experimental volumetric models,and volume percentages of all objects are compared.The average volume percentages of all models during the simulation using the mass-spring system,distance constraint,and the volume constraint method were 68.21%,89.64%,and 98.70%,respectively.The proposed approaches are successfully applied to improve the stability of mass-spring system and the performance comparison from our experimental tests also shows that the GPU-based method is faster than CPU-based implementation for all cases.展开更多
This paper proposes a stable and efficient implicit block Lower-Upper Symmetric-Gauss-Seidel(LU-SGS)algorithm-based lattice Boltzmann flux solver(LBFS)for simulation of hypersonic flows.In this method,the finite volum...This paper proposes a stable and efficient implicit block Lower-Upper Symmetric-Gauss-Seidel(LU-SGS)algorithm-based lattice Boltzmann flux solver(LBFS)for simulation of hypersonic flows.In this method,the finite volume method(FVM)is applied to discretize the Navier-Stokes equations,and the LBFS is utilized to evaluate the numerical flux at the cell interface.In LBFS,the local solution of discrete velocity Boltzmann equation(DVBE)with the non-free parameter D1Q4 lattice Boltzmann model is adopted to reconstruct the inviscid flux across the cell interface,and the viscous flux is approximated by conventional smooth function approach.In order to improve the robustness and convergence rate of the simulation for hypersonic flows,especially for problems with complex geometry,the implicit block LU-SGS algorithm is introduced to solve resultant discrete governing equations.A double cone model at Mach number of Ma=9.86 is firstly simulated to validate the proposed scheme,and a hypersonic flight vehicle with wings and rudders at Mach number of Ma=5.56 is then calculated to extend the application in practical engineering problems.Numerical results show that the proposed scheme could offer a more accurate and effective prediction for hypersonic flows.展开更多
Based on rectangular partition and bilinear interpolation,we construct an alternating-direction implicit(ADI)finite volume element method,which combined the merits of finite volume element method and alternating direc...Based on rectangular partition and bilinear interpolation,we construct an alternating-direction implicit(ADI)finite volume element method,which combined the merits of finite volume element method and alternating direction implicit method to solve a viscous wave equation with variable coefficients.This paper presents a general procedure to construct the alternating-direction implicit finite volume element method and gives computational schemes.Optimal error estimate in L2 norm is obtained for the schemes.Compared with the finite volume element method of the same convergence order,our method is more effective in terms of running time with the increasing of the computing scale.Numerical experiments are presented to show the efficiency of our method and numerical results are provided to support our theoretical analysis.展开更多
Semi-implicit direct kinetics(SIDK)is an innovative method for the temporal discretization of neutronic equations proposed by J.Banfield.The key approximation of the SIDK method is to substitute a timeaveraged quantit...Semi-implicit direct kinetics(SIDK)is an innovative method for the temporal discretization of neutronic equations proposed by J.Banfield.The key approximation of the SIDK method is to substitute a timeaveraged quantity for the fission source term in the delayed neutron differential equations.Hence,these equations are decoupled from prompt neutron equations and an explicit analytical representation of precursor groups is obtained,which leads to a significant reduction in computational cost.As the fission source is not known in a time step,the original study suggested using a constant quantity pertaining to the previous time step for this purpose,and a reduction in the size of the time step was proposed to lessen the imposed errors.However,this remedy notably diminishes the main advantage of the SIDK method.We discerned that if the original method is properly introduced into the algorithm of the point-implicit solver along with some modifications,the mentioned drawbacks will be mitigated adequately.To test this idea,a novel multigroup,multi-dimensional diffusion code using the finitevolume method and a point-implicit solver is developed which works in both transient and steady states.In addition to the SIDK,two other kinetic methods,i.e.,direct kinetics and higher-order backward discretization,are programmed into the diffusion code for comparison with the proposed model.The final code is tested at different conditions of two well-known transient benchmark problems.Results indicate that while the accuracy of the improved SIDK is closely comparable with the best available kinetic methods,it reduces the total time required for computation by up to 24%.展开更多
In this paper, a two-dimensional nonlinear coupled Gray Scott system is simulated with a finite difference scheme and a finite volume technique. Pre and post-processing lead to a new solution called GSmFoam by underst...In this paper, a two-dimensional nonlinear coupled Gray Scott system is simulated with a finite difference scheme and a finite volume technique. Pre and post-processing lead to a new solution called GSmFoam by understandin<span>g geometry settings and mesh information. The concentration profile chan</span>ges over time, as does the intensity of the contour patterns. The OpenFoam solver gives you the confidence to compare the pattern result with efficient numerical algorithms on the Gray Scott model.展开更多
In this paper,a numerical simulation model of the flow field in a gearbox with an oil volume adjusting device is established for the first time to study its influence on the lubrication characteristics of a high-speed...In this paper,a numerical simulation model of the flow field in a gearbox with an oil volume adjusting device is established for the first time to study its influence on the lubrication characteristics of a high-speed electric multiple unit(EMU)gearbox.The moving particle semi-implicit(MPS)method is used to numerically simulate the internal flow field of the gearbox of the high-speed EMU under working conditions.The effects of the velocity of the high-speed EMU,the immersion depth,and the oil sump temperature on the power loss of the gears and the lubricant quantity of each bearing are studied and provide an effective tool for the quantitative evaluation of the lubrication characteristics of the gearbox.The lubrication characteristics of the gearbox under different working conditions are studied when the oil volume adjusting device is closed and opened.The results show that the oil volume adjusting device mainly changes the amount of lubricant stirred by the output gear by changing the flow rate of lubricant from the cavity pinion(Cavity P)to the cavity gear(Cavity G),and thus affects the power loss of gears and the lubricant quantity of each bearing.展开更多
We study numerical methods for level set like equations arising in image processing and curve evolution problems. Semi-implicit finite volume-element type schemes are constructed for the general level set like equati...We study numerical methods for level set like equations arising in image processing and curve evolution problems. Semi-implicit finite volume-element type schemes are constructed for the general level set like equation （image selective smoothing model） given by Alvarez et al. （Alvarez L, Lions P L, Morel J M. Image selective smoothing and edge detection by nonlinear diffusion II. SIAM J. Numer. Anal., 1992, 29： 845-866）. Through the reasonable semi-implicit discretization in time and co-volume method for space approximation, we give finite volume schemes, unconditionally stable in L∞ and W1＇2 （W1＇1） sense in isotropic （anisotropic） diffu- sion domain.展开更多
基金The Major State Basic Research Program of China under contract No. 2012CB417002the National Natural Science Foundation of China under contract Nos 50909065 and 51109039
文摘A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C-unstructured non-orthogonal grid in the hor- izontal direction and z-level grid in the vertical direction. The three-dimensional model is discretized by the semi-implicit finite-volume method, in that the free-surface and the vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly by an integral method. The partial cell method is used for resolving topography, which enables the model to better represent irregular topography. The model has been tested against analytical cases for wind and tidal oscillation circulation, and is applied to simulating the tidal flow in the Bohal Sea. The results are in good agreement both with the analytical solutions and measurement results.
基金This work was supported by the Basic Science Research Program through the National Research Foundation of Korea(NRF-2019R1F1A1062752)funded by the Ministry of Education+1 种基金was funded by BK21 FOUR(Fostering Outstanding Universities for Research)(No.:5199990914048)and was also supported by the Soonchunhyang University Research Fund.
文摘This paper presents a parallel method for simulating real-time 3D deformable objects using the volume preservation mass-spring system method on tetrahedron meshes.In general,the conventional mass-spring system is manipulated as a force-driven method because it is fast,simple to implement,and the parameters can be controlled.However,the springs in traditional mass-spring system can be excessively elongated which cause severe stability and robustness issues that lead to shape restoring,simulation blow-up,and huge volume loss of the deformable object.In addition,traditional method that uses a serial process of the central processing unit(CPU)to solve the system in every frame cannot handle the complex structure of deformable object in real-time.Therefore,the first order implicit constraint enforcement for a mass-spring model is utilized to achieve accurate visual realism of deformable objects with tough constraint error.In this paper,we applied the distance constraint and volume conservation constraints for each tetrahedron element to improve the stability of deformable object simulation using the mass-spring system and behave the same as its real-world counterparts.To reduce the computational complexity while ensuring stable simulation,we applied a method that utilizes OpenGL compute shader,a part of OpenGL Shading Language(GLSL)that executes on the graphic processing unit(GPU)to solve the numerical problems effectively.We applied the proposed methods to experimental volumetric models,and volume percentages of all objects are compared.The average volume percentages of all models during the simulation using the mass-spring system,distance constraint,and the volume constraint method were 68.21%,89.64%,and 98.70%,respectively.The proposed approaches are successfully applied to improve the stability of mass-spring system and the performance comparison from our experimental tests also shows that the GPU-based method is faster than CPU-based implementation for all cases.
文摘This paper proposes a stable and efficient implicit block Lower-Upper Symmetric-Gauss-Seidel(LU-SGS)algorithm-based lattice Boltzmann flux solver(LBFS)for simulation of hypersonic flows.In this method,the finite volume method(FVM)is applied to discretize the Navier-Stokes equations,and the LBFS is utilized to evaluate the numerical flux at the cell interface.In LBFS,the local solution of discrete velocity Boltzmann equation(DVBE)with the non-free parameter D1Q4 lattice Boltzmann model is adopted to reconstruct the inviscid flux across the cell interface,and the viscous flux is approximated by conventional smooth function approach.In order to improve the robustness and convergence rate of the simulation for hypersonic flows,especially for problems with complex geometry,the implicit block LU-SGS algorithm is introduced to solve resultant discrete governing equations.A double cone model at Mach number of Ma=9.86 is firstly simulated to validate the proposed scheme,and a hypersonic flight vehicle with wings and rudders at Mach number of Ma=5.56 is then calculated to extend the application in practical engineering problems.Numerical results show that the proposed scheme could offer a more accurate and effective prediction for hypersonic flows.
基金supported by the National Natural Science Foundation of China grants No.11971241.
文摘Based on rectangular partition and bilinear interpolation,we construct an alternating-direction implicit(ADI)finite volume element method,which combined the merits of finite volume element method and alternating direction implicit method to solve a viscous wave equation with variable coefficients.This paper presents a general procedure to construct the alternating-direction implicit finite volume element method and gives computational schemes.Optimal error estimate in L2 norm is obtained for the schemes.Compared with the finite volume element method of the same convergence order,our method is more effective in terms of running time with the increasing of the computing scale.Numerical experiments are presented to show the efficiency of our method and numerical results are provided to support our theoretical analysis.
文摘Semi-implicit direct kinetics(SIDK)is an innovative method for the temporal discretization of neutronic equations proposed by J.Banfield.The key approximation of the SIDK method is to substitute a timeaveraged quantity for the fission source term in the delayed neutron differential equations.Hence,these equations are decoupled from prompt neutron equations and an explicit analytical representation of precursor groups is obtained,which leads to a significant reduction in computational cost.As the fission source is not known in a time step,the original study suggested using a constant quantity pertaining to the previous time step for this purpose,and a reduction in the size of the time step was proposed to lessen the imposed errors.However,this remedy notably diminishes the main advantage of the SIDK method.We discerned that if the original method is properly introduced into the algorithm of the point-implicit solver along with some modifications,the mentioned drawbacks will be mitigated adequately.To test this idea,a novel multigroup,multi-dimensional diffusion code using the finitevolume method and a point-implicit solver is developed which works in both transient and steady states.In addition to the SIDK,two other kinetic methods,i.e.,direct kinetics and higher-order backward discretization,are programmed into the diffusion code for comparison with the proposed model.The final code is tested at different conditions of two well-known transient benchmark problems.Results indicate that while the accuracy of the improved SIDK is closely comparable with the best available kinetic methods,it reduces the total time required for computation by up to 24%.
文摘In this paper, a two-dimensional nonlinear coupled Gray Scott system is simulated with a finite difference scheme and a finite volume technique. Pre and post-processing lead to a new solution called GSmFoam by understandin<span>g geometry settings and mesh information. The concentration profile chan</span>ges over time, as does the intensity of the contour patterns. The OpenFoam solver gives you the confidence to compare the pattern result with efficient numerical algorithms on the Gray Scott model.
基金supported by the Natural Science Foundation of Sichuan Province,China(Nos.2022NSFSC0034 and 2022NSFSC1901)the National Railway Group Science and Technology Program(No.N2021J028)+1 种基金the Independent Research and Development Projects of State Key Laboratory of Heavy Duty AC Drive Electric Locomotive Systems Integration(No.R111720H01385)the Independent Research and Development Projects of State Key Laboratory of Traction Power(No.2022TPL-T02),China。
文摘In this paper,a numerical simulation model of the flow field in a gearbox with an oil volume adjusting device is established for the first time to study its influence on the lubrication characteristics of a high-speed electric multiple unit(EMU)gearbox.The moving particle semi-implicit(MPS)method is used to numerically simulate the internal flow field of the gearbox of the high-speed EMU under working conditions.The effects of the velocity of the high-speed EMU,the immersion depth,and the oil sump temperature on the power loss of the gears and the lubricant quantity of each bearing are studied and provide an effective tool for the quantitative evaluation of the lubrication characteristics of the gearbox.The lubrication characteristics of the gearbox under different working conditions are studied when the oil volume adjusting device is closed and opened.The results show that the oil volume adjusting device mainly changes the amount of lubricant stirred by the output gear by changing the flow rate of lubricant from the cavity pinion(Cavity P)to the cavity gear(Cavity G),and thus affects the power loss of gears and the lubricant quantity of each bearing.
文摘We study numerical methods for level set like equations arising in image processing and curve evolution problems. Semi-implicit finite volume-element type schemes are constructed for the general level set like equation （image selective smoothing model） given by Alvarez et al. （Alvarez L, Lions P L, Morel J M. Image selective smoothing and edge detection by nonlinear diffusion II. SIAM J. Numer. Anal., 1992, 29： 845-866）. Through the reasonable semi-implicit discretization in time and co-volume method for space approximation, we give finite volume schemes, unconditionally stable in L∞ and W1＇2 （W1＇1） sense in isotropic （anisotropic） diffu- sion domain.
