Two-level finite element approximation to stream function form of unsteady Navier-Stokes equations is studied.This algorithm involves solving one nonlinear system on a coarse grid and one linear problem on a fine grid...Two-level finite element approximation to stream function form of unsteady Navier-Stokes equations is studied.This algorithm involves solving one nonlinear system on a coarse grid and one linear problem on a fine grid.Moreover,the scaling between these two grid sizes is super-linear.Approximation,stability and convergence aspects of a fully discrete scheme are analyzed.At last a numrical example is given whose results show that the algorithm proposed in this paper is effcient.展开更多
The use of metal oxides has been extensively documented in the literature and applied in a variety of contexts,including but not limited to energy storage,chemical sensors,and biomedical applications.One of the most s...The use of metal oxides has been extensively documented in the literature and applied in a variety of contexts,including but not limited to energy storage,chemical sensors,and biomedical applications.One of the most significant applications of metal oxides is heterogeneous catalysis,which represents a pivotal technology in industrial production on a global scale.Catalysts serve as the primary enabling agents for chemical reactions,and among the plethora of catalysts,metal oxides including magnesium oxide(MgO),ceria(CeO_(2))and titania(TiO_(2)),have been identified to be particularly effective in catalyzing a variety of reactions[1].Theoretical calculations based on density functional theory(DFT)and a multitude of other quantum chemistry methods have proven invaluable in elucidating the mechanisms of metal-oxide-catalyzed reactions,thereby facilitating the design of high-performance catalysts[2].展开更多
With the continuous advancement in topology optimization and additive manufacturing(AM)technology,the capability to fabricate functionally graded materials and intricate cellular structures with spatially varying micr...With the continuous advancement in topology optimization and additive manufacturing(AM)technology,the capability to fabricate functionally graded materials and intricate cellular structures with spatially varying microstructures has grown significantly.However,a critical challenge is encountered in the design of these structures–the absence of robust interface connections between adjacent microstructures,potentially resulting in diminished efficiency or macroscopic failure.A Hybrid Level Set Method(HLSM)is proposed,specifically designed to enhance connectivity among non-uniform microstructures,contributing to the design of functionally graded cellular structures.The HLSM introduces a pioneering algorithm for effectively blending heterogeneous microstructure interfaces.Initially,an interpolation algorithm is presented to construct transition microstructures seamlessly connected on both sides.Subsequently,the algorithm enables the morphing of non-uniform unit cells to seamlessly adapt to interconnected adjacent microstructures.The method,seamlessly integrated into a multi-scale topology optimization framework using the level set method,exhibits its efficacy through numerical examples,showcasing its prowess in optimizing 2D and 3D functionally graded materials(FGM)and multi-scale topology optimization.In essence,the pressing issue of interface connections in complex structure design is not only addressed but also a robust methodology is introduced,substantiated by numerical evidence,advancing optimization capabilities in the realm of functionally graded materials and cellular structures.展开更多
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.展开更多
A two-level discretization method for eigenvalue problems is studied.Compared to the standard Galerkin finite element discretization technique performed on a fine gridthis method discretizes the eigenvalue problem on ...A two-level discretization method for eigenvalue problems is studied.Compared to the standard Galerkin finite element discretization technique performed on a fine gridthis method discretizes the eigenvalue problem on a coarse grid and obtains an improved eigenvector(eigenvalue) approximation by solving only a linear problem on the fine grid (or two linear problemsfor the case of eigenvalue approximation of nonsymmetric problems). The improved solution has theasymptotic accuracy of the Galerkin discretization solution. The link between the method and theiterated Galerkin method is established. Error estimates for the general nonsymmetric case arederived.展开更多
In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique ...In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.展开更多
A two-level stabilized finite element method for the Stokes eigenvalue problem based on the local Gauss integration is considered. This method involves solving a Stokes eigenvalue problem on a coarse mesh with mesh si...A two-level stabilized finite element method for the Stokes eigenvalue problem based on the local Gauss integration is considered. This method involves solving a Stokes eigenvalue problem on a coarse mesh with mesh size H and a Stokes problem on a fine mesh with mesh size h -- O（H2）, which can still maintain the asymptotically optimal accuracy. It provides an approximate solution with the convergence rate of the same order as the usual stabilized finite element solution, which involves solving a Stokes eigenvalue problem on a fine mesh with mesh size h. Hence, the two-level stabilized finite element method can save a large amount of computational time. Moreover, numerical tests confirm the theoretical results of the present method.展开更多
The potential energy curves of the ground state X2∑＋g of the fluorine molecule have been accurately reconstructed employing the Ryderg-Klein-Rees （RKR） method extrapolated by a Hulburt and Hirschfeler potential fu...The potential energy curves of the ground state X2∑＋g of the fluorine molecule have been accurately reconstructed employing the Ryderg-Klein-Rees （RKR） method extrapolated by a Hulburt and Hirschfeler potential function for longer internuclear distances. Solving the corresponding radial one-dimensional Schr？dinger equation of nuclear motion yields 22 bound vibrational levels above v=0. The comparison of these theoretical levels with the experimental data yields a mean absolute deviation of about 7.6 cm^-1 over the 23 levels. The highest vibrational level energy obtained using this method is 13308.16 cm？1 and the relative deviation compared with the experimental datum of 13408.49 cm^-1 is only 0.74%. The value from our method is much closer and more accurate than the value obtained by the quantum mechanical ab initio method by Bytautas. The reported agreement of the vibrational levels and dissociation energy with experiment is contingent upon the potential energy curve of the F2 ground state.展开更多
Based on domain decomposition, a parallel two-level finite element method for the stationary Navier-Stokes equations is proposed and analyzed. The basic idea of the method is first to solve the Navier-Stokes equations...Based on domain decomposition, a parallel two-level finite element method for the stationary Navier-Stokes equations is proposed and analyzed. The basic idea of the method is first to solve the Navier-Stokes equations on a coarse grid, then to solve the resulted residual equations in parallel on a fine grid. This method has low communication complexity. It can be implemented easily. By local a priori error estimate for finite element discretizations, error bounds of the approximate solution are derived. Numerical results are also given to illustrate the high efficiency of the method.展开更多
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 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.展开更多
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.展开更多
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.展开更多
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.展开更多
Segmentation of three-dimensional(3D) complicated structures is of great importance for many real applications.In this work we combine graph cut minimization method with a variant of the level set idea for 3D segmenta...Segmentation of three-dimensional(3D) complicated structures is of great importance for many real applications.In this work we combine graph cut minimization method with a variant of the level set idea for 3D segmentation based on the Mumford-Shah model.Compared with the traditional approach for solving the Euler-Lagrange equation we do not need to solve any partial differential equations.Instead,the minimum cut on a special designed graph need to be computed.The method is tested on data with complicated structures.It is rather stable with respect to initial value and the algorithm is nearly parameter free.Experiments show that it can solve large problems much faster than traditional approaches.展开更多
When describing the mechanical behavior of some engineering materials,such as composites,grains,biological materials and cellular solids,the Cosserat continuum theory has more powerful capabilities compared with the c...When describing the mechanical behavior of some engineering materials,such as composites,grains,biological materials and cellular solids,the Cosserat continuum theory has more powerful capabilities compared with the classical Cauchy elasticity since an additional local rotation of point and its counterpart(couple stress)are considered in the Cosserat elasticity to represent the material microscale effects.In this paper,a parameterized level set topology optimization method is developed based on the Cosserat elasticity for the minimum compliance problem of the Cosserat solids.The influence of material characteristic length and Cosserat shear modulus on the optimized structure is investigated in detail.It can be found that the microstructural constants in the Cosserat elasticity have a significant impact on the optimized topology configurations.In addition,the minimum feature size and the geometric complexity of the optimized structure can be controlled implicitly by adjusting the parameters of the characteristic length and Cosserat shear modulus easily.Furthermore,the optimized structure obtained by the developed Cosserat elasticity based parameterized level set method will degenerate to the result by using the classical Cauchy elasticity based parameterized level set method when the Cosserat shear modulus approaches zero.展开更多
A flexure hinge is a major component in designing compliant mechanisms that o ers unique possibilities in a wide range of application fields in which high positioning accuracy is required. Although various flexure hin...A flexure hinge is a major component in designing compliant mechanisms that o ers unique possibilities in a wide range of application fields in which high positioning accuracy is required. Although various flexure hinges with di erent configurations have been successively proposed, they are often designed based on designers' experiences and inspirations. This study presents a systematic method for topological optimization of flexure hinges by using the level set method. Optimization formulations are developed by considering the functional requirements and geometrical constraints of flexure hinges. The functional requirements are first constructed by maximizing the compliance in the desired direction while minimizing the compliances in the other directions. The weighting sum method is used to construct an objective function in which a self-adjust method is used to set the weighting factors. A constraint on the symmetry of the obtained configuration is developed. Several numerical examples are presented to demonstrate the validity of the proposed method. The obtained results reveal that the design of a flexure hinge starting from the topology level can yield more choices for compliant mechanism design and obtain better designs that achieve higher performance.展开更多
This article introduces a new normalized nonlocal hybrid level set method for image segmentation.Due to intensity overlapping,blurred edges with complex backgrounds,simple intensity and texture information,such kind o...This article introduces a new normalized nonlocal hybrid level set method for image segmentation.Due to intensity overlapping,blurred edges with complex backgrounds,simple intensity and texture information,such kind of image segmentation is still a challenging task.The proposed method uses both the region and boundary information to achieve accurate segmentation results.The region information can help to identify rough region of interest and prevent the boundary leakage problem.It makes use of normalized nonlocal comparisons between pairs of patches in each region,and a heuristic intensity model is proposed to suppress irrelevant strong edges and constrain the segmentation.The boundary information can help to detect the precise location of the target object,it makes use of the geodesic active contour model to obtain the target boundary.The corresponding variational segmentation problem is implemented by a level set formulation.We use an internal energy term for geometric active contours to penalize the deviation of the level set function from a signed distance function.At last,experimental results on synthetic images and real images are shown in the paper with promising results.展开更多
文摘Two-level finite element approximation to stream function form of unsteady Navier-Stokes equations is studied.This algorithm involves solving one nonlinear system on a coarse grid and one linear problem on a fine grid.Moreover,the scaling between these two grid sizes is super-linear.Approximation,stability and convergence aspects of a fully discrete scheme are analyzed.At last a numrical example is given whose results show that the algorithm proposed in this paper is effcient.
基金financial support from the National Key R&D Program of China(2021YFB3500700)the National Natural Science Foundation of China(22473042,22003016,and 92145302).
文摘The use of metal oxides has been extensively documented in the literature and applied in a variety of contexts,including but not limited to energy storage,chemical sensors,and biomedical applications.One of the most significant applications of metal oxides is heterogeneous catalysis,which represents a pivotal technology in industrial production on a global scale.Catalysts serve as the primary enabling agents for chemical reactions,and among the plethora of catalysts,metal oxides including magnesium oxide(MgO),ceria(CeO_(2))and titania(TiO_(2)),have been identified to be particularly effective in catalyzing a variety of reactions[1].Theoretical calculations based on density functional theory(DFT)and a multitude of other quantum chemistry methods have proven invaluable in elucidating the mechanisms of metal-oxide-catalyzed reactions,thereby facilitating the design of high-performance catalysts[2].
基金the National Key Research and Development Program of China(Grant Number 2021YFB1714600)the National Natural Science Foundation of China(Grant Number 52075195)the Fundamental Research Funds for the Central Universities,China through Program No.2172019kfyXJJS078.
文摘With the continuous advancement in topology optimization and additive manufacturing(AM)technology,the capability to fabricate functionally graded materials and intricate cellular structures with spatially varying microstructures has grown significantly.However,a critical challenge is encountered in the design of these structures–the absence of robust interface connections between adjacent microstructures,potentially resulting in diminished efficiency or macroscopic failure.A Hybrid Level Set Method(HLSM)is proposed,specifically designed to enhance connectivity among non-uniform microstructures,contributing to the design of functionally graded cellular structures.The HLSM introduces a pioneering algorithm for effectively blending heterogeneous microstructure interfaces.Initially,an interpolation algorithm is presented to construct transition microstructures seamlessly connected on both sides.Subsequently,the algorithm enables the morphing of non-uniform unit cells to seamlessly adapt to interconnected adjacent microstructures.The method,seamlessly integrated into a multi-scale topology optimization framework using the level set method,exhibits its efficacy through numerical examples,showcasing its prowess in optimizing 2D and 3D functionally graded materials(FGM)and multi-scale topology optimization.In essence,the pressing issue of interface connections in complex structure design is not only addressed but also a robust methodology is introduced,substantiated by numerical evidence,advancing optimization capabilities in the realm of functionally graded materials and cellular structures.
基金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.
文摘A two-level discretization method for eigenvalue problems is studied.Compared to the standard Galerkin finite element discretization technique performed on a fine gridthis method discretizes the eigenvalue problem on a coarse grid and obtains an improved eigenvector(eigenvalue) approximation by solving only a linear problem on the fine grid (or two linear problemsfor the case of eigenvalue approximation of nonsymmetric problems). The improved solution has theasymptotic accuracy of the Galerkin discretization solution. The link between the method and theiterated Galerkin method is established. Error estimates for the general nonsymmetric case arederived.
文摘In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.
基金Project supported by the National Natural Science Foundation of China(Nos.10901131,10971166, and 10961024)the National High Technology Research and Development Program of China (No.2009AA01A135)the Natural Science Foundation of Xinjiang Uygur Autonomous Region (No.2010211B04)
文摘A two-level stabilized finite element method for the Stokes eigenvalue problem based on the local Gauss integration is considered. This method involves solving a Stokes eigenvalue problem on a coarse mesh with mesh size H and a Stokes problem on a fine mesh with mesh size h -- O（H2）, which can still maintain the asymptotically optimal accuracy. It provides an approximate solution with the convergence rate of the same order as the usual stabilized finite element solution, which involves solving a Stokes eigenvalue problem on a fine mesh with mesh size h. Hence, the two-level stabilized finite element method can save a large amount of computational time. Moreover, numerical tests confirm the theoretical results of the present method.
基金This work was supported by the National Natural Science Foundation of China （No.20273066）.
文摘The potential energy curves of the ground state X2∑＋g of the fluorine molecule have been accurately reconstructed employing the Ryderg-Klein-Rees （RKR） method extrapolated by a Hulburt and Hirschfeler potential function for longer internuclear distances. Solving the corresponding radial one-dimensional Schr？dinger equation of nuclear motion yields 22 bound vibrational levels above v=0. The comparison of these theoretical levels with the experimental data yields a mean absolute deviation of about 7.6 cm^-1 over the 23 levels. The highest vibrational level energy obtained using this method is 13308.16 cm？1 and the relative deviation compared with the experimental datum of 13408.49 cm^-1 is only 0.74%. The value from our method is much closer and more accurate than the value obtained by the quantum mechanical ab initio method by Bytautas. The reported agreement of the vibrational levels and dissociation energy with experiment is contingent upon the potential energy curve of the F2 ground state.
基金Project supported by the National Natural Science Foundation of China(No.11001061)the Science and Technology Foundation of Guizhou Province of China(No.[2008]2123)
文摘Based on domain decomposition, a parallel two-level finite element method for the stationary Navier-Stokes equations is proposed and analyzed. The basic idea of the method is first to solve the Navier-Stokes equations on a coarse grid, then to solve the resulted residual equations in parallel on a fine grid. This method has low communication complexity. It can be implemented easily. By local a priori error estimate for finite element discretizations, error bounds of the approximate solution are derived. Numerical results are also given to illustrate the high efficiency of the method.
基金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.
基金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 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.
基金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.
基金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.
基金support from the Centre for Integrated Petroleum Research(CIPR),University of Bergen, Norway,and Singapore MOE Grant T207B2202NRF2007IDMIDM002-010
文摘Segmentation of three-dimensional(3D) complicated structures is of great importance for many real applications.In this work we combine graph cut minimization method with a variant of the level set idea for 3D segmentation based on the Mumford-Shah model.Compared with the traditional approach for solving the Euler-Lagrange equation we do not need to solve any partial differential equations.Instead,the minimum cut on a special designed graph need to be computed.The method is tested on data with complicated structures.It is rather stable with respect to initial value and the algorithm is nearly parameter free.Experiments show that it can solve large problems much faster than traditional approaches.
基金This work was supported by the National Natural Science Foundation of China(Grants 12072242,11772237,and 11472196)the Hubei Provincial Natural Science Foundation(Grant 2020CFB816)the Fundamental Research Funds for the Central Universities(Grant 2042018kf0016).
文摘When describing the mechanical behavior of some engineering materials,such as composites,grains,biological materials and cellular solids,the Cosserat continuum theory has more powerful capabilities compared with the classical Cauchy elasticity since an additional local rotation of point and its counterpart(couple stress)are considered in the Cosserat elasticity to represent the material microscale effects.In this paper,a parameterized level set topology optimization method is developed based on the Cosserat elasticity for the minimum compliance problem of the Cosserat solids.The influence of material characteristic length and Cosserat shear modulus on the optimized structure is investigated in detail.It can be found that the microstructural constants in the Cosserat elasticity have a significant impact on the optimized topology configurations.In addition,the minimum feature size and the geometric complexity of the optimized structure can be controlled implicitly by adjusting the parameters of the characteristic length and Cosserat shear modulus easily.Furthermore,the optimized structure obtained by the developed Cosserat elasticity based parameterized level set method will degenerate to the result by using the classical Cauchy elasticity based parameterized level set method when the Cosserat shear modulus approaches zero.
基金Supported by National Natural Science Foundation of China(Grant Nos.51605166,51820105007)Fundamental Research Funds for the Central Universities of China
文摘A flexure hinge is a major component in designing compliant mechanisms that o ers unique possibilities in a wide range of application fields in which high positioning accuracy is required. Although various flexure hinges with di erent configurations have been successively proposed, they are often designed based on designers' experiences and inspirations. This study presents a systematic method for topological optimization of flexure hinges by using the level set method. Optimization formulations are developed by considering the functional requirements and geometrical constraints of flexure hinges. The functional requirements are first constructed by maximizing the compliance in the desired direction while minimizing the compliances in the other directions. The weighting sum method is used to construct an objective function in which a self-adjust method is used to set the weighting factors. A constraint on the symmetry of the obtained configuration is developed. Several numerical examples are presented to demonstrate the validity of the proposed method. The obtained results reveal that the design of a flexure hinge starting from the topology level can yield more choices for compliant mechanism design and obtain better designs that achieve higher performance.
基金supported in part by the National Natural Science Foundation of China(11626214,11571309)the General Research Project of Zhejiang Provincial Department of Education(Y201635378)+3 种基金the Zhejiang Provincial Natural Science Foundation of China(LY17F020011)J.Peng is supported by the National Natural Science Foundation of China(11771160)the Research Promotion Program of Huaqiao University(ZQN-PY411)Natural Science Foundation of Fujian Province(2015J01254)
文摘This article introduces a new normalized nonlocal hybrid level set method for image segmentation.Due to intensity overlapping,blurred edges with complex backgrounds,simple intensity and texture information,such kind of image segmentation is still a challenging task.The proposed method uses both the region and boundary information to achieve accurate segmentation results.The region information can help to identify rough region of interest and prevent the boundary leakage problem.It makes use of normalized nonlocal comparisons between pairs of patches in each region,and a heuristic intensity model is proposed to suppress irrelevant strong edges and constrain the segmentation.The boundary information can help to detect the precise location of the target object,it makes use of the geodesic active contour model to obtain the target boundary.The corresponding variational segmentation problem is implemented by a level set formulation.We use an internal energy term for geometric active contours to penalize the deviation of the level set function from a signed distance function.At last,experimental results on synthetic images and real images are shown in the paper with promising results.
