This paper concerns the sonic-supersonic structures of the transonic crossflow generated by the steady supersonic flow past an infinite cone of arbitrary cross section.Under the conical assumption,the three-dimensiona...This paper concerns the sonic-supersonic structures of the transonic crossflow generated by the steady supersonic flow past an infinite cone of arbitrary cross section.Under the conical assumption,the three-dimensional(3-D)steady Euler equations can be projected onto the unit sphere and the state of fluid can be characterized by the polar and azimuthal angles.Given a segment smooth curve as a conical-sonic line in the polar-azimuthal angle plane,we construct a classical conical-supersonic solution near the curve under some reasonable assumptions.To overcome the difficulty caused by the parabolic degeneracy,we apply the characteristic decomposition technique to transform the Euler equations into a new degenerate hyperbolic system in a partial hodograph plane.The singular terms are isolated from the highly nonlinear complicated system and then can be handled successfully.We establish a smooth local solution to the new system in a suitable weighted metric space and then express the solution in terms of the original variables.展开更多
Two high-order splitting schemes based on the idea of the operators splitting method are given. The three-dimensional advection-diffusion equation was split into several one-dimensional equations that were solved by t...Two high-order splitting schemes based on the idea of the operators splitting method are given. The three-dimensional advection-diffusion equation was split into several one-dimensional equations that were solved by these two schemes, only three computational grid points were needed in each direction but the accuracy reaches the spatial fourth-order. The third scheme proposed is based on the classical ADI scheme and the accuracy of the advection term of it can reach the spatial fourth-order. Finally, two typical numerical experiments show that the solutions of these three schemes compare well with that given by the analytical solution when the Peclet number is not bigger than 5.展开更多
In this work, we present a numerical model to solve the drift diffusion equations coupled with electromagnetic model, where all simulations codes are implemented using MATLAB code software. As first, we present a one-...In this work, we present a numerical model to solve the drift diffusion equations coupled with electromagnetic model, where all simulations codes are implemented using MATLAB code software. As first, we present a one-dimensional (1-D) PIN diode structure simulation achieved by solving the drift diffusion model (DDM). Backward Euler algorithm is used for the discretization of the proposed model. The aim is to accomplish time-domain integration. Also, finite different method (FDM) is considered to achieve space-Domain mesh. We introduced an iterative scheme to solve the obtained matrix systems, which combines the Gummel’s iteration with an efficient direct numerical UMFPACK method. The obtained solutions of the proposed algorithm provide the time and space distribution of the unknown functions like electrostatic potential and carrier’s concentration for the PIN diode. As second case, the finite-difference time-domain (FDTD) technique is adopted to analyze the entire 3-D structure of the stripline circuit including the lumped element PIN diode. The microwave circuit is located in an unbounded medium, requiring absorbing boundaries to avoid nonphysical reflections. Active device results were presented and show a good agreement with other reference. Electromagnetic results are qualitatively in agreement with other results obtained using SILVACO-TCAD.展开更多
In this work,di erent kinds of traveling wave solutions and uncategorized soliton wave solutions are obtained in a three dimensional(3-D)nonlinear evolution equations(NEEs)through the implementation of the modi ed ext...In this work,di erent kinds of traveling wave solutions and uncategorized soliton wave solutions are obtained in a three dimensional(3-D)nonlinear evolution equations(NEEs)through the implementation of the modi ed extended direct algebraic method.Bright-singular and dark-singular combo solitons,Jacobi's elliptic functions,Weierstrass elliptic functions,constant wave solutions and so on are attained beside their existing conditions.Physical interpretation of the solutions to the 3-D modi ed KdV-Zakharov-Kuznetsov equation are also given.展开更多
In this paper, we are first concerned with viscous approximations for the three-dimensional axisymmetric incompressible Euler equations. It is proved that the viscous approximations, which are the solutions of the cor...In this paper, we are first concerned with viscous approximations for the three-dimensional axisymmetric incompressible Euler equations. It is proved that the viscous approximations, which are the solutions of the corresponding Navier-Stokes equations, converge strongly in provided that they have strong convergence in the region away from the symmetry axis. This result has been proved by the authors for the approximate solutions generated by smoothing the initial data, with no restriction of the sign of the initial data. Then we discuss the decay rate for maximal vorticity function, which is established for both approximate solutions generated by smoothing the initial data and viscous approximations respectively. One sufficient condition to guarantee the strong convergence in the region away from the symmetry axis is given, and a decay rate for maximal vorticity function in the region away from the symmetry axis is obtained for non-negative initial vorticity.展开更多
In this paper,we establish the exponential convergence theory for the multipole and local expansions,shifting and translation operators for the Green's function of 3-dimensional Laplace equation in layered media.A...In this paper,we establish the exponential convergence theory for the multipole and local expansions,shifting and translation operators for the Green's function of 3-dimensional Laplace equation in layered media.An immediate application of the theory is to ensure the exponential convergence of the FMM which has been shown by the numerical results reported in[27].As the Green's function in layered media consists of free space and reaction field components and the theory for the free space components is well known,this paper will focus on the analysis for the reaction components.We first prove that the density functions in the integral representations of the reaction components are analytic and bounded in the right half complex wave number plane.Then,by using the Cagniard-de Hoop transform and contour deformations,estimates for the remainder terms of the truncated expansions are given,and,as a result,the exponential convergence for the expansions and translation operators is proven.展开更多
The generalized lattice Boltzmann equation(GLBE),with the addition of the standard Smagorinsky subgrid-stress(SGS) model,has been proved that it is more suitable for simulating high Reynolds number turbulent flows whe...The generalized lattice Boltzmann equation(GLBE),with the addition of the standard Smagorinsky subgrid-stress(SGS) model,has been proved that it is more suitable for simulating high Reynolds number turbulent flows when compared with the lattice BGK Boltzmann equation(LBGK).However,the computing efficiency of lattice Boltzmann method(LBM) is too low to make it for practical applications,unless using a massive parallel computing clusters facility.In this study,the massive parallel computing power from an inexpensive graphic processor unit(GPU) and a typical personal computer has been developed for improving the computing efficiency,more than 100 times.This developed three-dimensional(3-D) GLBE-SGS model,with the D3Q19 scheme for simplifying collision and streaming courses,has been successfully used to study 3-D rectangular cavity flows with Reynolds number up to 10000.展开更多
The principle and accuracy of 3-D coordinates acquisition using one single camera and the Aided Measuring Probe(AMP) are discussed in this paper. Using one single camera and one AMP which has several embedded targets ...The principle and accuracy of 3-D coordinates acquisition using one single camera and the Aided Measuring Probe(AMP) are discussed in this paper. Using one single camera and one AMP which has several embedded targets and one tip with known coordinates, the single camera′s orientation and location can be calculated. After orientation, the global coordinate system is obtained. During measurement, the camera is fixed firstly, then the AMP is held and the feature point is touched.The camera is triggered lastly. The position and orientation of the AMP are therefore calculated from the size and position of its image on the sensor. Since the tip point of AMP has known relation with the embedded targets, the feature point can be measured. Tests show that the accuracy of length measurement is 0.2 mm and accuracy for flatness measurement in XSY-plane is 0.1 mm.展开更多
Numerical analysis of three-dimensional(3-D)two-phase reacting flowfield in an annular combustor wity the dump diffuser is developed in arbitrary curvilinear coordi-nates.Combustor performances are estimated by the em...Numerical analysis of three-dimensional(3-D)two-phase reacting flowfield in an annular combustor wity the dump diffuser is developed in arbitrary curvilinear coordi-nates.Combustor performances are estimated by the em-pirical-analytical desing method.Ths influence of three inlet velocity profiles of the prediffuser and two operating conditions on combustor preformance and flow character-istic is predicted.展开更多
Thermal neutron albedo has been investigated for different thicknesses of mono-material and bi-material reflectors. An equation has been obtained for a bi-material reflector by considering the neutron diffusion equati...Thermal neutron albedo has been investigated for different thicknesses of mono-material and bi-material reflectors. An equation has been obtained for a bi-material reflector by considering the neutron diffusion equation. The bi-material reflector consists of binary combinations of water, graphite, lead, and polyethylene. An experimental measurement of thermal neutron albedo has also been conducted for mono-material and bi-material reflectors by using a^(241) Am–Be(5.2 Ci) neutron source and a BF3 detector. The maximum value of thermal neutron albedo was obtained for a polyethylene–water combination(0.95 ± 0.02).展开更多
In this article, we report the derivation of high accuracy finite difference method based on arithmetic average discretization for the solution of Un=F(x,u,u′)+∫K(x,s)ds , 0 x s < 1 subject to natural boundary co...In this article, we report the derivation of high accuracy finite difference method based on arithmetic average discretization for the solution of Un=F(x,u,u′)+∫K(x,s)ds , 0 x s < 1 subject to natural boundary conditions on a non-uniform mesh. The proposed variable mesh approximation is directly applicable to the integro-differential equation with singular coefficients. We need not require any special discretization to obtain the solution near the singular point. The convergence analysis of a difference scheme for the diffusion convection equation is briefly discussed. The presented variable mesh strategy is applicable when the internal grid points of the solution space are both even and odd in number as compared to the method discussed by authors in their previous work in which the internal grid points are strictly odd in number. The advantage of using this new variable mesh strategy is highlighted computationally.展开更多
A 3-D time-domain numerical coupled model for nonlinear waves acting on a ship in a harbor has been developed in the present study.The whole domain is divided into the inner domain and the outer domain.The inner domai...A 3-D time-domain numerical coupled model for nonlinear waves acting on a ship in a harbor has been developed in the present study.The whole domain is divided into the inner domain and the outer domain.The inner domain is the area around the ship,where the flow is expressed by the Laplace equation and numerically solved by the finite element method.The other area is the outer domain,where the flow is described by the higher-order Boussinesq equations and numerically solved by the finite difference method.The matching conditions on the interfaces between the inner domain and the outer domain,the procedure of coupled solution,the length of common domain and the mesh generation in the inner domain are discussed in detail.The other coupled model with the flow in the inner domain governed by the simplified linear Euler equations and relevant physical experiment are adopted to validate the present coupled model,and it is shown that the numerical results of the present model agree with the experimental data,so the present model can be used for the study on the effect of nonlinear waves acting on a fixed ship in a large area and provide a reference for the time-domain simulation of nonlinear wave forces on an arbitrary object in a large harbor and the 3-D district computation in the future.展开更多
Wave breaking plays an important role in wave-structure interaction. A novel control volume finite element method with adaptive unstructured meshes is employed here to study 3-D breaking waves. The numerical framework...Wave breaking plays an important role in wave-structure interaction. A novel control volume finite element method with adaptive unstructured meshes is employed here to study 3-D breaking waves. The numerical framework consists of a ＂volume of fluid＂ type method for the interface capturing and adaptive unstructured meshes to improve computational efficiency. The numerical model is validated against experimental measurements of breaking wave over a sloping beach and is then used to study the breaking wave impact on a vertical circular cylinder on a slope. Detailed complex interfacial structures during wave impact, such as plunging jet formation and splash-up are captured in the simulation, demonstrating the capability of the present method.展开更多
基金the two referees for very helpful comments and suggestions to improve the quality of the paper.This work was partially supported by the Natural Science Foundation of Zhejiang province of China(LY21A010017)the National Natural Science Foundation of China(12071106,12171130).
文摘This paper concerns the sonic-supersonic structures of the transonic crossflow generated by the steady supersonic flow past an infinite cone of arbitrary cross section.Under the conical assumption,the three-dimensional(3-D)steady Euler equations can be projected onto the unit sphere and the state of fluid can be characterized by the polar and azimuthal angles.Given a segment smooth curve as a conical-sonic line in the polar-azimuthal angle plane,we construct a classical conical-supersonic solution near the curve under some reasonable assumptions.To overcome the difficulty caused by the parabolic degeneracy,we apply the characteristic decomposition technique to transform the Euler equations into a new degenerate hyperbolic system in a partial hodograph plane.The singular terms are isolated from the highly nonlinear complicated system and then can be handled successfully.We establish a smooth local solution to the new system in a suitable weighted metric space and then express the solution in terms of the original variables.
文摘Two high-order splitting schemes based on the idea of the operators splitting method are given. The three-dimensional advection-diffusion equation was split into several one-dimensional equations that were solved by these two schemes, only three computational grid points were needed in each direction but the accuracy reaches the spatial fourth-order. The third scheme proposed is based on the classical ADI scheme and the accuracy of the advection term of it can reach the spatial fourth-order. Finally, two typical numerical experiments show that the solutions of these three schemes compare well with that given by the analytical solution when the Peclet number is not bigger than 5.
文摘In this work, we present a numerical model to solve the drift diffusion equations coupled with electromagnetic model, where all simulations codes are implemented using MATLAB code software. As first, we present a one-dimensional (1-D) PIN diode structure simulation achieved by solving the drift diffusion model (DDM). Backward Euler algorithm is used for the discretization of the proposed model. The aim is to accomplish time-domain integration. Also, finite different method (FDM) is considered to achieve space-Domain mesh. We introduced an iterative scheme to solve the obtained matrix systems, which combines the Gummel’s iteration with an efficient direct numerical UMFPACK method. The obtained solutions of the proposed algorithm provide the time and space distribution of the unknown functions like electrostatic potential and carrier’s concentration for the PIN diode. As second case, the finite-difference time-domain (FDTD) technique is adopted to analyze the entire 3-D structure of the stripline circuit including the lumped element PIN diode. The microwave circuit is located in an unbounded medium, requiring absorbing boundaries to avoid nonphysical reflections. Active device results were presented and show a good agreement with other reference. Electromagnetic results are qualitatively in agreement with other results obtained using SILVACO-TCAD.
文摘In this work,di erent kinds of traveling wave solutions and uncategorized soliton wave solutions are obtained in a three dimensional(3-D)nonlinear evolution equations(NEEs)through the implementation of the modi ed extended direct algebraic method.Bright-singular and dark-singular combo solitons,Jacobi's elliptic functions,Weierstrass elliptic functions,constant wave solutions and so on are attained beside their existing conditions.Physical interpretation of the solutions to the 3-D modi ed KdV-Zakharov-Kuznetsov equation are also given.
基金partially supported by National Natural Sciences Foundation of China (No.10101014)Beijing Natural Sciences Foundation+1 种基金the Key Project of NSFB-FBEC,by Grants from RGC of HKSAR CUHK4279/00P and CUHK4129/99Pthe generous hospitality and financial support of IMS of The Chinese University of Hong Kongpartially supported by Zheng Ge Ru Funds, Grants from RGC of HKSAR CUHK4279/00P and CUHK4129/99P
文摘In this paper, we are first concerned with viscous approximations for the three-dimensional axisymmetric incompressible Euler equations. It is proved that the viscous approximations, which are the solutions of the corresponding Navier-Stokes equations, converge strongly in provided that they have strong convergence in the region away from the symmetry axis. This result has been proved by the authors for the approximate solutions generated by smoothing the initial data, with no restriction of the sign of the initial data. Then we discuss the decay rate for maximal vorticity function, which is established for both approximate solutions generated by smoothing the initial data and viscous approximations respectively. One sufficient condition to guarantee the strong convergence in the region away from the symmetry axis is given, and a decay rate for maximal vorticity function in the region away from the symmetry axis is obtained for non-negative initial vorticity.
基金supported by the US National Science Foundation (Grant No.DMS-1950471)the US Army Research Office (Grant No.W911NF-17-1-0368)partially supported by NSFC (grant Nos.12201603 and 12022104)。
文摘In this paper,we establish the exponential convergence theory for the multipole and local expansions,shifting and translation operators for the Green's function of 3-dimensional Laplace equation in layered media.An immediate application of the theory is to ensure the exponential convergence of the FMM which has been shown by the numerical results reported in[27].As the Green's function in layered media consists of free space and reaction field components and the theory for the free space components is well known,this paper will focus on the analysis for the reaction components.We first prove that the density functions in the integral representations of the reaction components are analytic and bounded in the right half complex wave number plane.Then,by using the Cagniard-de Hoop transform and contour deformations,estimates for the remainder terms of the truncated expansions are given,and,as a result,the exponential convergence for the expansions and translation operators is proven.
基金supported by the Virginia Institute of Marine Science,College of William and Mary for the Study Environmentthe National Natural Science Foundation of China(Grant No.50679008)
文摘The generalized lattice Boltzmann equation(GLBE),with the addition of the standard Smagorinsky subgrid-stress(SGS) model,has been proved that it is more suitable for simulating high Reynolds number turbulent flows when compared with the lattice BGK Boltzmann equation(LBGK).However,the computing efficiency of lattice Boltzmann method(LBM) is too low to make it for practical applications,unless using a massive parallel computing clusters facility.In this study,the massive parallel computing power from an inexpensive graphic processor unit(GPU) and a typical personal computer has been developed for improving the computing efficiency,more than 100 times.This developed three-dimensional(3-D) GLBE-SGS model,with the D3Q19 scheme for simplifying collision and streaming courses,has been successfully used to study 3-D rectangular cavity flows with Reynolds number up to 10000.
文摘The principle and accuracy of 3-D coordinates acquisition using one single camera and the Aided Measuring Probe(AMP) are discussed in this paper. Using one single camera and one AMP which has several embedded targets and one tip with known coordinates, the single camera′s orientation and location can be calculated. After orientation, the global coordinate system is obtained. During measurement, the camera is fixed firstly, then the AMP is held and the feature point is touched.The camera is triggered lastly. The position and orientation of the AMP are therefore calculated from the size and position of its image on the sensor. Since the tip point of AMP has known relation with the embedded targets, the feature point can be measured. Tests show that the accuracy of length measurement is 0.2 mm and accuracy for flatness measurement in XSY-plane is 0.1 mm.
文摘Numerical analysis of three-dimensional(3-D)two-phase reacting flowfield in an annular combustor wity the dump diffuser is developed in arbitrary curvilinear coordi-nates.Combustor performances are estimated by the em-pirical-analytical desing method.Ths influence of three inlet velocity profiles of the prediffuser and two operating conditions on combustor preformance and flow character-istic is predicted.
文摘Thermal neutron albedo has been investigated for different thicknesses of mono-material and bi-material reflectors. An equation has been obtained for a bi-material reflector by considering the neutron diffusion equation. The bi-material reflector consists of binary combinations of water, graphite, lead, and polyethylene. An experimental measurement of thermal neutron albedo has also been conducted for mono-material and bi-material reflectors by using a^(241) Am–Be(5.2 Ci) neutron source and a BF3 detector. The maximum value of thermal neutron albedo was obtained for a polyethylene–water combination(0.95 ± 0.02).
文摘In this article, we report the derivation of high accuracy finite difference method based on arithmetic average discretization for the solution of Un=F(x,u,u′)+∫K(x,s)ds , 0 x s < 1 subject to natural boundary conditions on a non-uniform mesh. The proposed variable mesh approximation is directly applicable to the integro-differential equation with singular coefficients. We need not require any special discretization to obtain the solution near the singular point. The convergence analysis of a difference scheme for the diffusion convection equation is briefly discussed. The presented variable mesh strategy is applicable when the internal grid points of the solution space are both even and odd in number as compared to the method discussed by authors in their previous work in which the internal grid points are strictly odd in number. The advantage of using this new variable mesh strategy is highlighted computationally.
基金the National Natural Science Foundation of China(Grant Nos.59979002,50809008)the Hong Kong Research Grants Council(Grant No.HKU7171/06E)+2 种基金the China Postdoctoral Science Foundation(Grant No.20060400972)the Project of the Educational Department of Liaoning Province(Grant No.2005058)the Dalian Science and Technology Foundation(Grant No.2007J23JH027)
文摘A 3-D time-domain numerical coupled model for nonlinear waves acting on a ship in a harbor has been developed in the present study.The whole domain is divided into the inner domain and the outer domain.The inner domain is the area around the ship,where the flow is expressed by the Laplace equation and numerically solved by the finite element method.The other area is the outer domain,where the flow is described by the higher-order Boussinesq equations and numerically solved by the finite difference method.The matching conditions on the interfaces between the inner domain and the outer domain,the procedure of coupled solution,the length of common domain and the mesh generation in the inner domain are discussed in detail.The other coupled model with the flow in the inner domain governed by the simplified linear Euler equations and relevant physical experiment are adopted to validate the present coupled model,and it is shown that the numerical results of the present model agree with the experimental data,so the present model can be used for the study on the effect of nonlinear waves acting on a fixed ship in a large area and provide a reference for the time-domain simulation of nonlinear wave forces on an arbitrary object in a large harbor and the 3-D district computation in the future.
基金the financial support by the National Natural Science Foundation of China (Grant No. 51490673)the Open Awards of the State Key Laboratory of Coastal and Offshore Engineering+1 种基金funded by the EPSRC MEMPHIS multiphase Programme (Grant No. EP/K003976/1)funding from the European Union Seventh Framework Programme (FP7/20072013) under grant agreement No. 603663 for the research project PEARL (Preparing for Extreme and Rare events in coasta L regions)
文摘Wave breaking plays an important role in wave-structure interaction. A novel control volume finite element method with adaptive unstructured meshes is employed here to study 3-D breaking waves. The numerical framework consists of a ＂volume of fluid＂ type method for the interface capturing and adaptive unstructured meshes to improve computational efficiency. The numerical model is validated against experimental measurements of breaking wave over a sloping beach and is then used to study the breaking wave impact on a vertical circular cylinder on a slope. Detailed complex interfacial structures during wave impact, such as plunging jet formation and splash-up are captured in the simulation, demonstrating the capability of the present method.
