The sealing performance of magnetic fluid is related to the magnetic fluid itself. Many factors can influence the magnetic field and the seal pressure differences of magnetic fluid seals, such as the sealing gap, the ...The sealing performance of magnetic fluid is related to the magnetic fluid itself. Many factors can influence the magnetic field and the seal pressure differences of magnetic fluid seals, such as the sealing gap, the shaft eccentricity, the shaft diameter, the volume of the magnetic fluid and the centrifugal force. These factors are analyzed by numerical computation. When the seal material and structure are the same, the seal pressure difference is directly proportional to the magnetic field intensity and the saturation magnetization of the magnetic fluid. The sealing performance of the magnetic fluid will reduce with the increase of the sealing gap and shaft eccentricity. The sealing performance will increase with the volume of the magnetic fluid and decrease with the increase of the shaft diameter taking gravity into account. The increase of the shaft diameter is the same as the reduction of the volume of the magnetic fluid. The magnetic fluid cross-section can change because of the centrifugal force. Some improvements can reducc the influence of the centrifugal force. The centrifugal force can be utilized to improve the sealing performance.展开更多
A finite different method is developed to predict the side force on a high speedtrain in a cross-wind at low yaw anglee. The k-εturbulence model with wallfunctions is employed. the solution algorithm is based on curv...A finite different method is developed to predict the side force on a high speedtrain in a cross-wind at low yaw anglee. The k-εturbulence model with wallfunctions is employed. the solution algorithm is based on curvilinearnonorthogonal coordinates,covariant velocity components ,and staggered gridarrangement. The convective fluxes are described by the Power tow Scheme.A highly deformed grid generated with an elliptic grid generator is used aroundthe comero of the cross-section of the train. The results obtained comparepositively with wind tunnel experinients.展开更多
In this paper a mixed finite element-characteristic mixed finite element method is discussed to simulate an incompressible miscible Darcy-Forchheimer problem.The flow equation is solved by a mixed finite element and t...In this paper a mixed finite element-characteristic mixed finite element method is discussed to simulate an incompressible miscible Darcy-Forchheimer problem.The flow equation is solved by a mixed finite element and the approximation accuracy of Darch-Forchheimer velocity is improved one order.The concentration equation is solved by the method of mixed finite element,where the convection is discretized along the characteristic direction and the diffusion is discretized by the zero-order mixed finite element method.The characteristics can confirm strong stability at sharp fronts and avoids numerical dispersion and nonphysical oscillation.In actual computations the characteristics adopts a large time step without any loss of accuracy.The scalar unknowns and its adjoint vector function are obtained simultaneously and the law of mass conservation holds in every element by the zero-order mixed finite element discretization of diffusion flux.In order to derive the optimal 3/2-order error estimate in L^(2) norm,a post-processing technique is included in the approximation to the scalar unknowns.Numerical experiments are illustrated finally to validate theoretical analysis and efficiency.This method can be used to solve such an important problem.展开更多
An energy method is proposed to investigate the critical transformation condition from a Taylor cone to a cone-jet. Based on the kinetic theorem, the system power allocation and the electrohydrodynamics stability are ...An energy method is proposed to investigate the critical transformation condition from a Taylor cone to a cone-jet. Based on the kinetic theorem, the system power allocation and the electrohydrodynamics stability are discussed. The numerical results indicate that the energy of the liquid cone tip experiences a maximum value during the transformation. With the proposed jetting energy, we give the critical transformation condition under which the derivative of jetting energy with respect to the surface area is greater than or equal to the energy required to form a unit of new liquid surface.展开更多
Using an improved computational fluid dynamics （CFD） method developed for highly unsteady three-dimensional flows, numerical simulations for oscillating flow cycles and detailed unsteady simulations of the flow and ...Using an improved computational fluid dynamics （CFD） method developed for highly unsteady three-dimensional flows, numerical simulations for oscillating flow cycles and detailed unsteady simulations of the flow and forces on the aortic vessels at the iliac bifurcation, for both healthy and diseased patients, are analyzed. Improvements in computational efficiency and acceleration in convergence are achieved by calculating both an unsteady pressure gradient which is due to fluid acceleration and a good global pressure field correction based on mass flow for the pressure Poisson equation. Applications of the enhanced method to oscillatory flow in curved pipes yield an order of magnitude increase in speed and efficiency, thus allowing the study of more complex flow problems such as flow through the mammalian abdominal aorta at the iliac arteries bifurcation. To analyze the large forces which can exist on stent graft of patients with abdominal aor- tic aneurysm （AAA） disease, a complete derivation of the force equations is presented. The accelerated numerical algorithm and the force equations derived are used to calculate flow and forces for two individuals whose geometry is obtained from CT data and whose respective blood pressure measurements are obtained experimentally. Although the use of endovascular stent grafts in diseased patients can alter vessel geometries, the physical characteristics of stents are still very different when compared to native blood vessels of healthy subjects. The geometry for the AAA stent graph patient studied in this investigation induced flows that resulted in large forces that are primarily caused by the blood pressure. These forces are also directly related to the flow cross-sectional area and the angle of the iliac arteries relative to the main descending aorta. Furthermore, the fluid flow is significantly disturbed in the diseased patient with large flow recirculation and stagnant regions which are not present for healthy subjects.展开更多
The numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science. The mathematical model is defined by an initialboundary nonlinear system of four partial differenti...The numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science. The mathematical model is defined by an initialboundary nonlinear system of four partial differential equations: an elliptic equation for electric potential, two convection-diffusion equations for electron concentration and hole concentration, and a heat conduction equation for temperature. The first equation is solved by the conservative block-centered method. The concentrations and temperature are computed by the block-centered upwind difference method on a changing mesh, where the block-centered method and upwind approximation are used to discretize the diffusion and convection, respectively. The computations on a changing mesh show very well the local special properties nearby the P-N junction. The upwind scheme is applied to approximate the convection, and numerical dispersion and nonphysical oscillation are avoided. The block-centered difference computes concentrations, temperature, and their adjoint vector functions simultaneously.The local conservation of mass, an important rule in the numerical simulation of a semiconductor device, is preserved during the computations. An optimal order convergence is obtained. Numerical examples are provided to show efficiency and application.展开更多
Based on the gradient-Hamiltonian decomposition （GHD） theory of vector fields, an algorithm （ called as GHD algorithm） is proposed in this paper. For the GHD algorithm, visual interpretations of the advantages in ...Based on the gradient-Hamiltonian decomposition （GHD） theory of vector fields, an algorithm （ called as GHD algorithm） is proposed in this paper. For the GHD algorithm, visual interpretations of the advantages in stability are given by using the eigenvalue curves. From the numerical results for linear decay systems, it reaches the conclusion that the GHD algorithm proposed in this paper has a better computational accuracy than other algorithms and presents a replication of long time qualitative properties of the underlying system.展开更多
A new type water-cooled heat dissipater for multiple high-power thyristors in explosion-proof shell used in coal mine was designed, and then, the numerical computation of the three-dimensional steady-state temperature...A new type water-cooled heat dissipater for multiple high-power thyristors in explosion-proof shell used in coal mine was designed, and then, the numerical computation of the three-dimensional steady-state temperature distributions under different working conditions for cooling core was conducted in order to understand in detail the heat transfer performance. Based on the computation results, the temperature differences and the maximum heat transfer rates were given. These results of the study on the heat dissipater lay a basis for optimising its structure design and guiding its operation.展开更多
基金Projects 10172031 supported by National Natural Science Foundation of China and HIT. MD2002.15 by the Multidiscipline Scientific Research Founda-tion of Harbin Institute of Technology
文摘The sealing performance of magnetic fluid is related to the magnetic fluid itself. Many factors can influence the magnetic field and the seal pressure differences of magnetic fluid seals, such as the sealing gap, the shaft eccentricity, the shaft diameter, the volume of the magnetic fluid and the centrifugal force. These factors are analyzed by numerical computation. When the seal material and structure are the same, the seal pressure difference is directly proportional to the magnetic field intensity and the saturation magnetization of the magnetic fluid. The sealing performance of the magnetic fluid will reduce with the increase of the sealing gap and shaft eccentricity. The sealing performance will increase with the volume of the magnetic fluid and decrease with the increase of the shaft diameter taking gravity into account. The increase of the shaft diameter is the same as the reduction of the volume of the magnetic fluid. The magnetic fluid cross-section can change because of the centrifugal force. Some improvements can reducc the influence of the centrifugal force. The centrifugal force can be utilized to improve the sealing performance.
文摘A finite different method is developed to predict the side force on a high speedtrain in a cross-wind at low yaw anglee. The k-εturbulence model with wallfunctions is employed. the solution algorithm is based on curvilinearnonorthogonal coordinates,covariant velocity components ,and staggered gridarrangement. The convective fluxes are described by the Power tow Scheme.A highly deformed grid generated with an elliptic grid generator is used aroundthe comero of the cross-section of the train. The results obtained comparepositively with wind tunnel experinients.
基金supported by the Natural ScienceFoundation of Shandong Province(ZR2021MA019)。
文摘In this paper a mixed finite element-characteristic mixed finite element method is discussed to simulate an incompressible miscible Darcy-Forchheimer problem.The flow equation is solved by a mixed finite element and the approximation accuracy of Darch-Forchheimer velocity is improved one order.The concentration equation is solved by the method of mixed finite element,where the convection is discretized along the characteristic direction and the diffusion is discretized by the zero-order mixed finite element method.The characteristics can confirm strong stability at sharp fronts and avoids numerical dispersion and nonphysical oscillation.In actual computations the characteristics adopts a large time step without any loss of accuracy.The scalar unknowns and its adjoint vector function are obtained simultaneously and the law of mass conservation holds in every element by the zero-order mixed finite element discretization of diffusion flux.In order to derive the optimal 3/2-order error estimate in L^(2) norm,a post-processing technique is included in the approximation to the scalar unknowns.Numerical experiments are illustrated finally to validate theoretical analysis and efficiency.This method can be used to solve such an important problem.
基金supported by the National Basic Research Program of China(Grant No.2013CB733004)
文摘An energy method is proposed to investigate the critical transformation condition from a Taylor cone to a cone-jet. Based on the kinetic theorem, the system power allocation and the electrohydrodynamics stability are discussed. The numerical results indicate that the energy of the liquid cone tip experiences a maximum value during the transformation. With the proposed jetting energy, we give the critical transformation condition under which the derivative of jetting energy with respect to the surface area is greater than or equal to the energy required to form a unit of new liquid surface.
文摘Using an improved computational fluid dynamics （CFD） method developed for highly unsteady three-dimensional flows, numerical simulations for oscillating flow cycles and detailed unsteady simulations of the flow and forces on the aortic vessels at the iliac bifurcation, for both healthy and diseased patients, are analyzed. Improvements in computational efficiency and acceleration in convergence are achieved by calculating both an unsteady pressure gradient which is due to fluid acceleration and a good global pressure field correction based on mass flow for the pressure Poisson equation. Applications of the enhanced method to oscillatory flow in curved pipes yield an order of magnitude increase in speed and efficiency, thus allowing the study of more complex flow problems such as flow through the mammalian abdominal aorta at the iliac arteries bifurcation. To analyze the large forces which can exist on stent graft of patients with abdominal aor- tic aneurysm （AAA） disease, a complete derivation of the force equations is presented. The accelerated numerical algorithm and the force equations derived are used to calculate flow and forces for two individuals whose geometry is obtained from CT data and whose respective blood pressure measurements are obtained experimentally. Although the use of endovascular stent grafts in diseased patients can alter vessel geometries, the physical characteristics of stents are still very different when compared to native blood vessels of healthy subjects. The geometry for the AAA stent graph patient studied in this investigation induced flows that resulted in large forces that are primarily caused by the blood pressure. These forces are also directly related to the flow cross-sectional area and the angle of the iliac arteries relative to the main descending aorta. Furthermore, the fluid flow is significantly disturbed in the diseased patient with large flow recirculation and stagnant regions which are not present for healthy subjects.
基金supported the Natural Science Foundation of Shandong Province(ZR2016AM08)Natural Science Foundation of Hunan Province(2018JJ2028)National Natural Science Foundation of China(11871312).
文摘The numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science. The mathematical model is defined by an initialboundary nonlinear system of four partial differential equations: an elliptic equation for electric potential, two convection-diffusion equations for electron concentration and hole concentration, and a heat conduction equation for temperature. The first equation is solved by the conservative block-centered method. The concentrations and temperature are computed by the block-centered upwind difference method on a changing mesh, where the block-centered method and upwind approximation are used to discretize the diffusion and convection, respectively. The computations on a changing mesh show very well the local special properties nearby the P-N junction. The upwind scheme is applied to approximate the convection, and numerical dispersion and nonphysical oscillation are avoided. The block-centered difference computes concentrations, temperature, and their adjoint vector functions simultaneously.The local conservation of mass, an important rule in the numerical simulation of a semiconductor device, is preserved during the computations. An optimal order convergence is obtained. Numerical examples are provided to show efficiency and application.
基金Sponsored by the National Natural Science Foundation of China (10572021)the Doctoral Programme Foundation of Institute of Higher Educationof China (20040007022)the Program for New Century Excellent Talents in University of Fujian Province
文摘Based on the gradient-Hamiltonian decomposition （GHD） theory of vector fields, an algorithm （ called as GHD algorithm） is proposed in this paper. For the GHD algorithm, visual interpretations of the advantages in stability are given by using the eigenvalue curves. From the numerical results for linear decay systems, it reaches the conclusion that the GHD algorithm proposed in this paper has a better computational accuracy than other algorithms and presents a replication of long time qualitative properties of the underlying system.
文摘A new type water-cooled heat dissipater for multiple high-power thyristors in explosion-proof shell used in coal mine was designed, and then, the numerical computation of the three-dimensional steady-state temperature distributions under different working conditions for cooling core was conducted in order to understand in detail the heat transfer performance. Based on the computation results, the temperature differences and the maximum heat transfer rates were given. These results of the study on the heat dissipater lay a basis for optimising its structure design and guiding its operation.
