The basic theory and effect of the new farming method of ＂Fenlong＂ cultivation which has been included in the main extension technology of Ministry of Agriculture of the People＇s Republic of China is fully illustra...The basic theory and effect of the new farming method of ＂Fenlong＂ cultivation which has been included in the main extension technology of Ministry of Agriculture of the People＇s Republic of China is fully illustrated for the first time, and it is the fourth set （generation） of farming modes and methods following manpower, animal and mechanical （tractor） farming. It follows the natural law to achieve soil activation, water saving, oxygen increase, warming and desalination through the active use of natural resources like soil, rainfall and solar energy, thereby promoting a new round of natural agricultural production and quality improvement and water con- servation, which has crop yield increase by 10%-30%, quality improvement of 5%, natural precipitation retaining increase by100%. The characteristics and mechanism are the use of spiral drill for one-time completion of the land preparation by drilling vertically to 30-50 cm of soil layer through high speed peeling. After instant high temperature and many fierce impacts, mechanical frictions, it could achieve the multiplication of the number of loose soil, soil physical modification and expansion of the soil nutrients, reservoirs, oxygen, microorganisms （＂Four pools＂）. Fenlong cultivation can give birth to new farming culture and civilization, and it can achieve the physical ＂desalinized＂ transformation and utilization of saline soil. The formation of Fenlong green farming technology system makes it possible to invent the farming tools of ＂serf-propelled Fenlong machinery＂ that has got the patent, and it is the method for farmland （dry land, paddy field） Fenlong cultivation, saline-alkali soil smash-ridging cultivation and for the abundance of grass ecology on degraded grassland. The application of Fenlong ＂4＋1＂ （arable, saline-alkali soil, grasslands, Sponge City＋rivers） green development in China can achieve the ＂double safety＂ of food and living space.展开更多
The goal of the arterial graft design problem is to find an optimal graft built on an occluded artery, which can be mathematically modeled by a fluid based shape optimization problem. The smoothness of the graft is on...The goal of the arterial graft design problem is to find an optimal graft built on an occluded artery, which can be mathematically modeled by a fluid based shape optimization problem. The smoothness of the graft is one of the important aspects in the arterial graft design problem since it affects the flow of the blood significantly. As an attractive design tool for this problem, level set methods are quite efficient for obtaining better shape of the graft. In this paper, a cubic spline level set method and a radial basis function level set method are designed to solve the arterial graft design problem. In both approaches, the shape of the arterial graft is implicitly tracked by the zero-level contour of a level set function and a high level of smoothness of the graft is achieved. Numerical results show the efficiency of the algorithms in the arterial graft design.展开更多
Some basic problems on the level set methods were discussed, such as the method used to preserve the distance junction , the existence and uniqueness of solution for the level set equations. The main contribution is t...Some basic problems on the level set methods were discussed, such as the method used to preserve the distance junction , the existence and uniqueness of solution for the level set equations. The main contribution is to prove that in a neighborhood of the initial zero level set, the level set equations with the restriction of the distance function have a unique solution, which must be the signed distance function with respect to the evolving surface. Some skillful approaches were used: Noticing that any solution for the original equation was a distance function, the original level set equations were transformed into a simpler alternative form. Moreover, since the new system was not a classical one, the system was transformed into an ordinary one, for which the implicit function method was adopted.展开更多
This paper presents an improved level set method for topology optimization of geometrically nonlinear structures accounting for the effect of thermo-mechanical couplings.It derives a new expression for element couplin...This paper presents an improved level set method for topology optimization of geometrically nonlinear structures accounting for the effect of thermo-mechanical couplings.It derives a new expression for element coupling stress resulting from the combination of mechanical and thermal loading,using geometric nonlinear finite element analysis.A topological model is then developed to minimize compliance while meeting displacement and frequency constraints to fulfill design requirements of structural members.Since the conventional Lagrange multiplier search method is unable to handle convergence instability arising from large deformation,a novel Lagrange multiplier search method is proposed.Additionally,the proposed method can be extended to multi-constrained geometrically nonlinear topology optimization,accommodating multiple physical field couplings.展开更多
This article studies numerically a familiar important phenomenon in spray combustion which is deformation and breakup of liquid drops in gas flow. The SIMPLER method is used to solve the two-dimensional （2D） unstead...This article studies numerically a familiar important phenomenon in spray combustion which is deformation and breakup of liquid drops in gas flow. The SIMPLER method is used to solve the two-dimensional （2D） unsteady axisymmetric Navier-Stokes equations for both the drop and the ambient gas flow. The level set method is applied to capturing the liquid/gas interface. Through calculation are obtained four typical breakup modes--oscillation, bag breakup, sheet stripping breakup and shear breakup governed by four non-dimensional numbers which are gas Weber number （Weg）, liquid Reynolds number （Rel）, gas Reynolds number （Reg） and density ratio （γ）. Their effects upon each mode are analyzed. The results indicate that among the four numbers, Weg is of the highest importance with Rel, Reg and γfollowing up. By widening the range of the density ratio up to 1 000, the breakup mode is discovered to be so complicated that a new one called multimode breakup mode turns up. This mode contains the shearing breakup and piercing breakup, which successively happen. The calculation results agree well with what is observed from the experiments.展开更多
A new program is developed for gas-liquid two-phase mold filling simulation in casting. The gas fluid, the superheated liquid metal and the liquid metal containing solid grains are assumed to be governed by Navier-Sto...A new program is developed for gas-liquid two-phase mold filling simulation in casting. The gas fluid, the superheated liquid metal and the liquid metal containing solid grains are assumed to be governed by Navier-Stokes equations and solved through Projection method. The Level set method is used to track the gas-liquid interface boundary. In order to demonstrate the correctness of this new program for simulation of gas-liquid two-phase mold filling in casting, a benchmark filling experiment is simulated （this benchmark test is designed by XU and the filling process is recorded by a 16-mm film camera）. The simulated results agree very well with the experimental results, showing that this new program can be used to properly predicate the gas-liquid two-phase mold filling simulation in casting.展开更多
Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-materia...Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.展开更多
Combining the vector level set model,the shape sensitivity analysis theory with the gradient projection technique,a level set method for topology optimization with multi-constraints and multi-materials is presented in...Combining the vector level set model,the shape sensitivity analysis theory with the gradient projection technique,a level set method for topology optimization with multi-constraints and multi-materials is presented in this paper.The method implicitly describes structural material in- terfaces by the vector level set and achieves the optimal shape and topology through the continuous evolution of the material interfaces in the structure.In order to increase computational efficiency for a fast convergence,an appropriate nonlinear speed mapping is established in the tangential space of the active constraints.Meanwhile,in order to overcome the numerical instability of general topology opti- mization problems,the regularization with the mean curvature flow is utilized to maintain the interface smoothness during the optimization process.The numerical examples demonstrate that the approach possesses a good flexibility in handling topological changes and gives an interface representation in a high fidelity,compared with other methods based on explicit boundary variations in the literature.展开更多
A non-isothermal injection molding process for a non-Newtonian viscous pseudoplastic fluid is simulated.A conservative interface capturing technique and the flow field solving method are coupled to perform a dynamic s...A non-isothermal injection molding process for a non-Newtonian viscous pseudoplastic fluid is simulated.A conservative interface capturing technique and the flow field solving method are coupled to perform a dynamic simulation.The validity of the numerical method is verified by a benchmark problem.The melt interface evolution versus time is captured and the physical quantities such as temperature,velocity and pressure at each time step are obtained with corresponding analysis.A"frozen skin"layer with the thickness increasing versus time during the injection process is found.The fact that the"frozen skin"layer can be reduced by increasing the injection velocity is numerically verified.The fountain flow phenomenon near the melt interface is also captured.Moreover,comparisons with the non-isothermal Newtonian case show that the curvatures of the interface arcs and the pressure contours near the horizontal mid-line of the cavity for the non-Newtonian pseudoplastic case is larger than that for the Newtonian case.The velocity profiles are different at different positions for the non-Newtonian pseudoplastic case,while in the case of Newtonian flow the velocity profiles are parabolic and almost the same at different positions.展开更多
Pore structure of porous media, including pore size and topology, is rather complex. In immiscible twophase displacement process, the capillary force affected by pore size dominates the two-phase flow in the porous me...Pore structure of porous media, including pore size and topology, is rather complex. In immiscible twophase displacement process, the capillary force affected by pore size dominates the two-phase flow in the porous media, affecting displacement results. Direct observation of the flow patterns in the porous media is difficult, and therefore knowledge about the two-phase displacement flow is insufficient. In this paper, a two-dimensional(2D) pore structure was extracted from a sandstone sample, and the flow process that CO_2 displaces resident brine in the extracted pore structure was simulated using the Navier eStokes equation combined with the conservative level set method. The simulation results reveal that the pore throat is a crucial factor for determining CO_2 displacement process in the porous media. The two-phase meniscuses in each pore throat were in a self-adjusting process. In the displacement process,CO_2 preferentially broke through the maximum pore throat. Before breaking through the maximum pore throat, the pressure of CO_2 continually increased, and the curvature and position of two-phase interfaces in the other pore throats adjusted accordingly. Once the maximum pore throat was broken through by the CO_2, the capillary force in the other pore throats released accordingly; subsequently, the interfaces withdrew under the effect of capillary fore, preparing for breaking through the next pore throat.Therefore, the two-phase displacement in CO_2 injection is accompanied by the breaking through and adjusting of the two-phase interfaces.展开更多
A level set method of non-uniform grids is used to simulate the whole evolution of a cavitation bubble, including its growth, collapse and rebound near a rigid wall. Single-phase Navier-Stokes equation in the liquid r...A level set method of non-uniform grids is used to simulate the whole evolution of a cavitation bubble, including its growth, collapse and rebound near a rigid wall. Single-phase Navier-Stokes equation in the liquid region is solved by MAC projection algorithm combined with second-order ENO scheme for the advection terms. The moving inter-face is captured by the level set function, and the interface velocity is resolved by ＂one-side＂ velocity extension from the liquid region to the bubble region, complementing the second-order weighted least squares method across the interface and projection inside bubble. The use of non-uniform grid overcomes the difficulty caused by the large computational domain and very small bubble size. The computation is very stable without suffering from large flow-field gradients, and the results are in good agreements with other studies. The bubble interface kinematics, dynamics and its effect on the wall are highlighted, which shows that the code can effectively capture the ＂shock wave＂-like pressure and velocity at jet impact, toroidal bubble, and complicated pressure structure with peak, plateau and valley in the later stage of bubble oscillating.展开更多
Target tracking is one of the most important issues in computer vision and has been applied in many fields of science, engineering and industry. Because of the occlusion during tracking, typical approaches with single...Target tracking is one of the most important issues in computer vision and has been applied in many fields of science, engineering and industry. Because of the occlusion during tracking, typical approaches with single classifier learn much of occluding background information which results in the decrease of tracking performance, and eventually lead to the failure of the tracking algorithm. This paper presents a new correlative classifiers approach to address the above problem. Our idea is to derive a group of correlative classifiers based on sample set method. Then we propose strategy to establish the classifiers and to query the suitable classifiers for the next frame tracking. In order to deal with nonlinear problem, particle filter is adopted and integrated with sample set method. For choosing the target from candidate particles, we define a similarity measurement between particles and sample set. The proposed sample set method includes the following steps. First, we cropped positive samples set around the target and negative samples set far away from the target. Second, we extracted average Haar-like feature from these samples and calculate their statistical characteristic which represents the target model. Third, we define the similarity measurement based on the statistical characteristic of these two sets to judge the similarity between candidate particles and target model. Finally, we choose the largest similarity score particle as the target in the new frame. A number of experiments show the robustness and efficiency of the proposed approach when compared with other state-of-the-art trackers.展开更多
A spatially adaptive (SA) two-dimensional (2-D) numerical wave flume is presented based on the quadtree mesh system,in which a new multiple particle level set (MPLS) method is proposed to solve the problem of interfac...A spatially adaptive (SA) two-dimensional (2-D) numerical wave flume is presented based on the quadtree mesh system,in which a new multiple particle level set (MPLS) method is proposed to solve the problem of interface tracking,in which common intersection may be traversed by multiple interfaces.By using the adaptive mesh technique and the MPLS method,mesh resolution is updated automatically with time according to flow characteristics in the modeling process with higher resolution around the free surface and the solid boundary and lower resolution in less important area.The model has good performance in saving computer memory and CPU time and is validated by computational examples of small amplitude wave,second-order Stokes wave and cnoidal wave.Computational results also indicate that standing wave and wave overtopping are also reasonably simulated by the model.展开更多
This paper proposes a multiscale isogeometric topology optimization(ITO)method where the configuration and layout of microstructures are optimized simultaneously.At micro scale,a shape deformation method is presented ...This paper proposes a multiscale isogeometric topology optimization(ITO)method where the configuration and layout of microstructures are optimized simultaneously.At micro scale,a shape deformation method is presented to transform a prototype microstructure(PM)for obtaining a series of graded microstructures(GMs),where microstructural skeleton based on the level set framework is applied to retain more topology features and improve the connectability.For the macro scale calculation,the effective mechanical properties can be estimated by means of the numerical homogenization method.By adopting identical non-uniform rational basis splines(NURBS)as basis functions for both parameterized level set model and isogeometric calculation model,the isogeometric analysis(IGA)is integrated into the level set method,which contributes to improving the accuracy and efficiency.Numerical examples demonstrate that,the proposed method is effective in improving the performance and manufacturability.展开更多
Segmenting a complex 3D surface model into some visually meaningful sub-parts is one of the fundamental problems in digital geometry processing. In this paper, a novel segmentation approach of point-sampled surfaces i...Segmenting a complex 3D surface model into some visually meaningful sub-parts is one of the fundamental problems in digital geometry processing. In this paper, a novel segmentation approach of point-sampled surfaces is proposed, which is based on the level set evolution scheme. To segment the model so as to align the patch boundaries with high curvature zones, the driven speed function for the zero level set inside narrow band is defined by the extended curvature field, which approaches zero speed as the propagating front approaches high curvature zone. The effectiveness of the proposed approach is demonstrated by our ex- perimental results. Furthermore, two applications of model segmentation are illustrated, such as piecewise parameterization and local editing for point-sampled geometry.展开更多
The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated includ...The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated including spherical, oblate and spherical. The bubble shape and drag coefficient were compared with experimental results. It is observed that the simulated results show good conformity to experimental results over a wide range of Reynolds number. In addition, the detailed flow field based on the reference coordinate system moving with the bubble is obtained, and the relationship among flow field, bubble shape and velocity is discussed.展开更多
Knowledge-Based Engineering （KBE） is introduced into the ship structural design in this paper. From the implementation of KBE, the design solutions for both Rules Design Method （RDM） and Interpolation Design Meth...Knowledge-Based Engineering （KBE） is introduced into the ship structural design in this paper. From the implementation of KBE, the design solutions for both Rules Design Method （RDM） and Interpolation Design Method （IDM） are generated. The corresponding Finite Element （FE） models are generated. Topological design of the longitudinal structures is studied where the Gaussian Process （GP） is employed to build the surrogate model for FE analysis. Multi-objective optimization methods inspired by Pareto Front are used to reduce the design tank weight and outer surface area simultaneously. Additionally, an enhanced Level Set Method （LSM） which employs implicit algorithm is applied to the topological design of typical bracket plate which is used extensively in ship structures. Two different sets of boundary conditions are considered. The proposed methods show satisfactory efficiency and accuracy.展开更多
Based on a level set model, a topology optimization method has been suggestedrecently. It uses a level set to express the moving structural boundary, which can flexibly handlecomplex topological changes. By combining ...Based on a level set model, a topology optimization method has been suggestedrecently. It uses a level set to express the moving structural boundary, which can flexibly handlecomplex topological changes. By combining vector level set models with gradient projectiontechnology, the level set method for topological optimization is extended to a topologicaloptimization problem with multi-constraints, multi-materials and multi-load cases. Meanwhile, anappropriate nonlinear speed, mapping is established in the tangential space of the activeconstraints for a fast convergence. Then the method is applied to structure designs, mechanism andmaterial designs by a number of benchmark examples. Finally, in order to further improvecomputational efficiency and overcome the difficulty that the level set method cannot generate newmaterial interfaces during the optimization process, the topological derivative analysis isincorporated into the level set method for topological optimization, and a topological derivativeand level set algorithm for topological optimization is proposed.展开更多
基金Supported by the National Key Technology R&D Program of China(2014BAD06B05)the Major Project of Science and Technology of Guangxi(2017AA22015)~~
文摘The basic theory and effect of the new farming method of ＂Fenlong＂ cultivation which has been included in the main extension technology of Ministry of Agriculture of the People＇s Republic of China is fully illustrated for the first time, and it is the fourth set （generation） of farming modes and methods following manpower, animal and mechanical （tractor） farming. It follows the natural law to achieve soil activation, water saving, oxygen increase, warming and desalination through the active use of natural resources like soil, rainfall and solar energy, thereby promoting a new round of natural agricultural production and quality improvement and water con- servation, which has crop yield increase by 10%-30%, quality improvement of 5%, natural precipitation retaining increase by100%. The characteristics and mechanism are the use of spiral drill for one-time completion of the land preparation by drilling vertically to 30-50 cm of soil layer through high speed peeling. After instant high temperature and many fierce impacts, mechanical frictions, it could achieve the multiplication of the number of loose soil, soil physical modification and expansion of the soil nutrients, reservoirs, oxygen, microorganisms （＂Four pools＂）. Fenlong cultivation can give birth to new farming culture and civilization, and it can achieve the physical ＂desalinized＂ transformation and utilization of saline soil. The formation of Fenlong green farming technology system makes it possible to invent the farming tools of ＂serf-propelled Fenlong machinery＂ that has got the patent, and it is the method for farmland （dry land, paddy field） Fenlong cultivation, saline-alkali soil smash-ridging cultivation and for the abundance of grass ecology on degraded grassland. The application of Fenlong ＂4＋1＂ （arable, saline-alkali soil, grasslands, Sponge City＋rivers） green development in China can achieve the ＂double safety＂ of food and living space.
基金Supported by National Foundation of Natural Science(11471092)Natural Science Foundation of Zhejiang Province(LZ13A010003)Foundation of Zhejiang Educational Committee(Y201121891)
文摘The goal of the arterial graft design problem is to find an optimal graft built on an occluded artery, which can be mathematically modeled by a fluid based shape optimization problem. The smoothness of the graft is one of the important aspects in the arterial graft design problem since it affects the flow of the blood significantly. As an attractive design tool for this problem, level set methods are quite efficient for obtaining better shape of the graft. In this paper, a cubic spline level set method and a radial basis function level set method are designed to solve the arterial graft design problem. In both approaches, the shape of the arterial graft is implicitly tracked by the zero-level contour of a level set function and a high level of smoothness of the graft is achieved. Numerical results show the efficiency of the algorithms in the arterial graft design.
基金the National Natural Science Foundation of China (6001161942, 60203003)
文摘Some basic problems on the level set methods were discussed, such as the method used to preserve the distance junction , the existence and uniqueness of solution for the level set equations. The main contribution is to prove that in a neighborhood of the initial zero level set, the level set equations with the restriction of the distance function have a unique solution, which must be the signed distance function with respect to the evolving surface. Some skillful approaches were used: Noticing that any solution for the original equation was a distance function, the original level set equations were transformed into a simpler alternative form. Moreover, since the new system was not a classical one, the system was transformed into an ordinary one, for which the implicit function method was adopted.
基金supported by grants from the National Natural Science Foundation of China (51478130)the Guangzhou Municipal Education Bureau’s Scientific Research Project, China (2024312217)+1 种基金the China Scholarship Council (201808440070)the 111 Project of China (D21021).
文摘This paper presents an improved level set method for topology optimization of geometrically nonlinear structures accounting for the effect of thermo-mechanical couplings.It derives a new expression for element coupling stress resulting from the combination of mechanical and thermal loading,using geometric nonlinear finite element analysis.A topological model is then developed to minimize compliance while meeting displacement and frequency constraints to fulfill design requirements of structural members.Since the conventional Lagrange multiplier search method is unable to handle convergence instability arising from large deformation,a novel Lagrange multiplier search method is proposed.Additionally,the proposed method can be extended to multi-constrained geometrically nonlinear topology optimization,accommodating multiple physical field couplings.
文摘This article studies numerically a familiar important phenomenon in spray combustion which is deformation and breakup of liquid drops in gas flow. The SIMPLER method is used to solve the two-dimensional （2D） unsteady axisymmetric Navier-Stokes equations for both the drop and the ambient gas flow. The level set method is applied to capturing the liquid/gas interface. Through calculation are obtained four typical breakup modes--oscillation, bag breakup, sheet stripping breakup and shear breakup governed by four non-dimensional numbers which are gas Weber number （Weg）, liquid Reynolds number （Rel）, gas Reynolds number （Reg） and density ratio （γ）. Their effects upon each mode are analyzed. The results indicate that among the four numbers, Weg is of the highest importance with Rel, Reg and γfollowing up. By widening the range of the density ratio up to 1 000, the breakup mode is discovered to be so complicated that a new one called multimode breakup mode turns up. This mode contains the shearing breakup and piercing breakup, which successively happen. The calculation results agree well with what is observed from the experiments.
基金Projects(51304145,51301118,51304152)supported by the National Natural Science Foundation of ChinaProject(2013JQ7016)supported by the Natural Science Foundation of Shannxi Province,China+1 种基金Project(2013T002)supported by the Science Foundation of Taiyuan University of Technology,ChinaProject(2013JK0904)supported by Shannxi Provincial Education Department,China
文摘A new program is developed for gas-liquid two-phase mold filling simulation in casting. The gas fluid, the superheated liquid metal and the liquid metal containing solid grains are assumed to be governed by Navier-Stokes equations and solved through Projection method. The Level set method is used to track the gas-liquid interface boundary. In order to demonstrate the correctness of this new program for simulation of gas-liquid two-phase mold filling in casting, a benchmark filling experiment is simulated （this benchmark test is designed by XU and the filling process is recorded by a 16-mm film camera）. The simulated results agree very well with the experimental results, showing that this new program can be used to properly predicate the gas-liquid two-phase mold filling simulation in casting.
基金Supported by the National Natural Science Foundation of China(10476011)~~
文摘Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.
基金The project supported by the National Natural Science Foundation of China (59805001,10332010) and Key Science and Technology Research Project of Ministry of Education of China (No.104060)
文摘Combining the vector level set model,the shape sensitivity analysis theory with the gradient projection technique,a level set method for topology optimization with multi-constraints and multi-materials is presented in this paper.The method implicitly describes structural material in- terfaces by the vector level set and achieves the optimal shape and topology through the continuous evolution of the material interfaces in the structure.In order to increase computational efficiency for a fast convergence,an appropriate nonlinear speed mapping is established in the tangential space of the active constraints.Meanwhile,in order to overcome the numerical instability of general topology opti- mization problems,the regularization with the mean curvature flow is utilized to maintain the interface smoothness during the optimization process.The numerical examples demonstrate that the approach possesses a good flexibility in handling topological changes and gives an interface representation in a high fidelity,compared with other methods based on explicit boundary variations in the literature.
基金Supported by the National Natural Science Foundation of China（10871159） the National Basic Research Program of China（2005CB321704）
文摘A non-isothermal injection molding process for a non-Newtonian viscous pseudoplastic fluid is simulated.A conservative interface capturing technique and the flow field solving method are coupled to perform a dynamic simulation.The validity of the numerical method is verified by a benchmark problem.The melt interface evolution versus time is captured and the physical quantities such as temperature,velocity and pressure at each time step are obtained with corresponding analysis.A"frozen skin"layer with the thickness increasing versus time during the injection process is found.The fact that the"frozen skin"layer can be reduced by increasing the injection velocity is numerically verified.The fountain flow phenomenon near the melt interface is also captured.Moreover,comparisons with the non-isothermal Newtonian case show that the curvatures of the interface arcs and the pressure contours near the horizontal mid-line of the cavity for the non-Newtonian pseudoplastic case is larger than that for the Newtonian case.The velocity profiles are different at different positions for the non-Newtonian pseudoplastic case,while in the case of Newtonian flow the velocity profiles are parabolic and almost the same at different positions.
基金funded by Key Laboratory of Coal-based CO_2 Capture and Geological Storage,Jiangsu Province,ChinaUS Advanced Coal Technology Consortium(No.2013 DFB60140-08)
文摘Pore structure of porous media, including pore size and topology, is rather complex. In immiscible twophase displacement process, the capillary force affected by pore size dominates the two-phase flow in the porous media, affecting displacement results. Direct observation of the flow patterns in the porous media is difficult, and therefore knowledge about the two-phase displacement flow is insufficient. In this paper, a two-dimensional(2D) pore structure was extracted from a sandstone sample, and the flow process that CO_2 displaces resident brine in the extracted pore structure was simulated using the Navier eStokes equation combined with the conservative level set method. The simulation results reveal that the pore throat is a crucial factor for determining CO_2 displacement process in the porous media. The two-phase meniscuses in each pore throat were in a self-adjusting process. In the displacement process,CO_2 preferentially broke through the maximum pore throat. Before breaking through the maximum pore throat, the pressure of CO_2 continually increased, and the curvature and position of two-phase interfaces in the other pore throats adjusted accordingly. Once the maximum pore throat was broken through by the CO_2, the capillary force in the other pore throats released accordingly; subsequently, the interfaces withdrew under the effect of capillary fore, preparing for breaking through the next pore throat.Therefore, the two-phase displacement in CO_2 injection is accompanied by the breaking through and adjusting of the two-phase interfaces.
基金the National Natural Science Foundation of China(10272032 and 10672043).
文摘A level set method of non-uniform grids is used to simulate the whole evolution of a cavitation bubble, including its growth, collapse and rebound near a rigid wall. Single-phase Navier-Stokes equation in the liquid region is solved by MAC projection algorithm combined with second-order ENO scheme for the advection terms. The moving inter-face is captured by the level set function, and the interface velocity is resolved by ＂one-side＂ velocity extension from the liquid region to the bubble region, complementing the second-order weighted least squares method across the interface and projection inside bubble. The use of non-uniform grid overcomes the difficulty caused by the large computational domain and very small bubble size. The computation is very stable without suffering from large flow-field gradients, and the results are in good agreements with other studies. The bubble interface kinematics, dynamics and its effect on the wall are highlighted, which shows that the code can effectively capture the ＂shock wave＂-like pressure and velocity at jet impact, toroidal bubble, and complicated pressure structure with peak, plateau and valley in the later stage of bubble oscillating.
基金supported by the National Science Foundation of China(61472289)National Key Research and Development Project(2016YFC0106305)The Key Technology R&D Program of Hubei Provence(2014BAA153)
文摘Target tracking is one of the most important issues in computer vision and has been applied in many fields of science, engineering and industry. Because of the occlusion during tracking, typical approaches with single classifier learn much of occluding background information which results in the decrease of tracking performance, and eventually lead to the failure of the tracking algorithm. This paper presents a new correlative classifiers approach to address the above problem. Our idea is to derive a group of correlative classifiers based on sample set method. Then we propose strategy to establish the classifiers and to query the suitable classifiers for the next frame tracking. In order to deal with nonlinear problem, particle filter is adopted and integrated with sample set method. For choosing the target from candidate particles, we define a similarity measurement between particles and sample set. The proposed sample set method includes the following steps. First, we cropped positive samples set around the target and negative samples set far away from the target. Second, we extracted average Haar-like feature from these samples and calculate their statistical characteristic which represents the target model. Third, we define the similarity measurement based on the statistical characteristic of these two sets to judge the similarity between candidate particles and target model. Finally, we choose the largest similarity score particle as the target in the new frame. A number of experiments show the robustness and efficiency of the proposed approach when compared with other state-of-the-art trackers.
基金The Innovative Research Groups of the National Natural Science Foundation of China under contract No.51021004the National Natural Science Foundation for Youth of China under contract No. 51109018+2 种基金the Open Foundation of Water & Sediment Science and Water Hazard Prevention Hunan Provincial Key Laboratory under contract No. 2011SS05the Open Foundation of Port,Coastal and offshore Engineering Hunan Provincial Key Discipline under contract No. 20110815001the Open Foundation of State Key Laboratory of Hydraulic Engineering Simulation and Safety under contract No.HSSKLTJU-201208.
文摘A spatially adaptive (SA) two-dimensional (2-D) numerical wave flume is presented based on the quadtree mesh system,in which a new multiple particle level set (MPLS) method is proposed to solve the problem of interface tracking,in which common intersection may be traversed by multiple interfaces.By using the adaptive mesh technique and the MPLS method,mesh resolution is updated automatically with time according to flow characteristics in the modeling process with higher resolution around the free surface and the solid boundary and lower resolution in less important area.The model has good performance in saving computer memory and CPU time and is validated by computational examples of small amplitude wave,second-order Stokes wave and cnoidal wave.Computational results also indicate that standing wave and wave overtopping are also reasonably simulated by the model.
基金National Key R&D Program of China(2018YFB1700803,2018YFB1700804).
文摘This paper proposes a multiscale isogeometric topology optimization(ITO)method where the configuration and layout of microstructures are optimized simultaneously.At micro scale,a shape deformation method is presented to transform a prototype microstructure(PM)for obtaining a series of graded microstructures(GMs),where microstructural skeleton based on the level set framework is applied to retain more topology features and improve the connectability.For the macro scale calculation,the effective mechanical properties can be estimated by means of the numerical homogenization method.By adopting identical non-uniform rational basis splines(NURBS)as basis functions for both parameterized level set model and isogeometric calculation model,the isogeometric analysis(IGA)is integrated into the level set method,which contributes to improving the accuracy and efficiency.Numerical examples demonstrate that,the proposed method is effective in improving the performance and manufacturability.
基金Project supported by the National Basic Research Program (973) of China (No. 2002CB312101)the National Natural Science Foundation of China (Nos. 60503056, 60373036, 60333010)the Education Department of Zhejiang Province, China (No. 20060797)
文摘Segmenting a complex 3D surface model into some visually meaningful sub-parts is one of the fundamental problems in digital geometry processing. In this paper, a novel segmentation approach of point-sampled surfaces is proposed, which is based on the level set evolution scheme. To segment the model so as to align the patch boundaries with high curvature zones, the driven speed function for the zero level set inside narrow band is defined by the extended curvature field, which approaches zero speed as the propagating front approaches high curvature zone. The effectiveness of the proposed approach is demonstrated by our ex- perimental results. Furthermore, two applications of model segmentation are illustrated, such as piecewise parameterization and local editing for point-sampled geometry.
基金Project(21406141)supported by the National Natural Science Foundation of ChinaProject(20141078)supported by the Scientific Research Starting Foundation for Doctors of Liaoning Province,China+1 种基金Project(L2014060)supported by the Foundation of Department of Education of Liaoning Province,ChinaProject(157B21)supported by the Scientific Research Starting Foundation for Doctors of Shenyang Aerospace University,China
文摘The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated including spherical, oblate and spherical. The bubble shape and drag coefficient were compared with experimental results. It is observed that the simulated results show good conformity to experimental results over a wide range of Reynolds number. In addition, the detailed flow field based on the reference coordinate system moving with the bubble is obtained, and the relationship among flow field, bubble shape and velocity is discussed.
基金financially supported by the Project of Ministry of Education and Finance of China(Grant Nos.200512 and 201335)the Project of the State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University(Grant No.GKZD010053-10)
文摘Knowledge-Based Engineering （KBE） is introduced into the ship structural design in this paper. From the implementation of KBE, the design solutions for both Rules Design Method （RDM） and Interpolation Design Method （IDM） are generated. The corresponding Finite Element （FE） models are generated. Topological design of the longitudinal structures is studied where the Gaussian Process （GP） is employed to build the surrogate model for FE analysis. Multi-objective optimization methods inspired by Pareto Front are used to reduce the design tank weight and outer surface area simultaneously. Additionally, an enhanced Level Set Method （LSM） which employs implicit algorithm is applied to the topological design of typical bracket plate which is used extensively in ship structures. Two different sets of boundary conditions are considered. The proposed methods show satisfactory efficiency and accuracy.
基金This project is supported by National Natural Science Foundation of China(No.598005001, No.10332010) and Key Science and Technology Research Project of Ministry of Education (No.104060).
文摘Based on a level set model, a topology optimization method has been suggestedrecently. It uses a level set to express the moving structural boundary, which can flexibly handlecomplex topological changes. By combining vector level set models with gradient projectiontechnology, the level set method for topological optimization is extended to a topologicaloptimization problem with multi-constraints, multi-materials and multi-load cases. Meanwhile, anappropriate nonlinear speed, mapping is established in the tangential space of the activeconstraints for a fast convergence. Then the method is applied to structure designs, mechanism andmaterial designs by a number of benchmark examples. Finally, in order to further improvecomputational efficiency and overcome the difficulty that the level set method cannot generate newmaterial interfaces during the optimization process, the topological derivative analysis isincorporated into the level set method for topological optimization, and a topological derivativeand level set algorithm for topological optimization is proposed.
