A Hybrid Level Set Optimization Design Method of Functionally Graded Cellular Structures Considering Connectivity

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.

Topology Optimization for Steady-State Navier-Stokes Flow Based on Parameterized Level Set Based Method

In this paper,we consider solving the topology optimization for steady-state incompressibleNavier-Stokes problems via a new topology optimization method called parameterized level set method,which can maintain a relatively smooth level set function with a local optimality condition.The objective of topology optimization is to􀀀nd an optimal con􀀀guration of theuid and solid materials that minimizes power dissipation under a prescribeduid volume fraction constraint.An arti􀀀cial friction force is added to the Navier-Stokes equations to apply the no-slip boundary condition.Although a great deal of work has been carried out for topology optimization ofuidow in recent years,there are few researches on the topology optimization ofuidow with physical body forces.To simulate theuidow in reality,the constant body force(e.g.,gravity)is considered in this paper.Several 2D numerical examples are presented to discuss the relationships between the proposed method with Reynolds number and initial design,and demonstrate the feasibility and superiority of the proposed method in dealing with unstructuredmesh problems.Three 3D numerical examples demonstrate the proposedmethod is feasible in three-dimensional.

A Low Mach Number IMEX Flux Splitting for the Level Set Ghost Fluid Method

Considering droplet phenomena at low Mach numbers,large differences in the magnitude of the occurring characteristic waves are presented.As acoustic phenomena often play a minor role in such applications,classical explicit schemes which resolve these waves suffer from a very restrictive timestep restriction.In this work,a novel scheme based on a specific level set ghost fluid method and an implicit-explicit(IMEX)flux splitting is proposed to overcome this timestep restriction.A fully implicit narrow band around the sharp phase interface is combined with a splitting of the convective and acoustic phenomena away from the interface.In this part of the domain,the IMEX Runge-Kutta time discretization and the high order discontinuous Galerkin spectral element method are applied to achieve high accuracies in the bulk phases.It is shown that for low Mach numbers a significant gain in computational time can be achieved compared to a fully explicit method.Applica-tions to typical droplet dynamic phenomena validate the proposed method and illustrate its capabilities.

Exponential Time Differencing Method for a Reaction-Diffusion System with Free Boundary

For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.

An Evidence-Based CoCoSo Framework with Double Hierarchy Linguistic Data for Viable Selection of Hydrogen Storage Methods

Hydrogen is the new age alternative energy source to combat energy demand and climate change.Storage of hydrogen is vital for a nation’s growth.Works of literature provide different methods for storing the produced hydrogen,and the rational selection of a viable method is crucial for promoting sustainability and green practices.Typically,hydrogen storage is associated with diverse sustainable and circular economy(SCE)criteria.As a result,the authors consider the situation a multi-criteria decision-making(MCDM)problem.Studies infer that previous models for hydrogen storage method(HSM)selection(i)do not consider preferences in the natural language form;(ii)weights of experts are not methodically determined;(iii)hesitation of experts during criteria weight assessment is not effectively explored;and(iv)three-stage solution of a suitable selection of HSM is unexplored.Driven by these gaps,in this paper,authors put forward a new integrated framework,which considers double hierarchy linguistic information for rating,criteria importance through inter-criteria correlation(CRITIC)for expert weight calculation,evidence-based Bayesian method for criteria weight estimation,and combined compromise solution(CoCoSo)for ranking HSMs.The applicability of the developed framework is testified by using a case example of HSM selection in India.Sensitivity and comparative analysis reveal the merits and limitations of the developed framework.

A Simple Embedding Method for the Laplace-Beltrami Eigenvalue Problem onImplicit Surfaces

We propose a simple embedding method for computing the eigenvalues and eigenfunctions of the Laplace-Beltrami operator on implicit surfaces.The approach follows an embedding approach for solving the surface eikonal equation.We replace the differential operator on the interface with a typical Cartesian differential operator in the surface neighborhood.Our proposed algorithm is easy to implement and efficient.We will give some two-and three-dimensional numerical examples to demonstrate the effectiveness of our proposed approach.

All-electron basis sets for H to Xe specific for ZORA calculations:Applications in atoms and molecules

A segmented basis set of quadruple zeta valence quality plus polarization functions(QZP)for H through Xe was developed to be used in conjunction with the ZORA Hamiltonian.This set was augmented with diffuse functions to describe electrons farther away from the nuclei adequately.Using the ZORA-CCSD(T)/QZP-ZORA theoretical model,atomic ionization energies and bond lengths,harmonic vibrational frequencies,and atomization energies of some molecules were calculated.The addition of core-valence corrections has been shown to improve the agreement between theoretical and experimental results for molecular properties.For atomization energies,a similar observation emerges when considering spin-orbit couplings.With the augmented QZP-ZORA set,static mean dipole polarizabilities of a set of atoms were calculated and compared with previously published recommended and experimental values.Performance evaluations of the ZORA and Douglas–Kroll–Hess Hamiltonians were made for each property studied.

Approximation Maintenance When Adding a Conditional Value in Set-Valued Ordered Decision Systems

The integration of set-valued ordered rough set models and incremental learning signify a progressive advancement of conventional rough set theory, with the objective of tackling the heterogeneity and ongoing transformations in information systems. In set-valued ordered decision systems, when changes occur in the attribute value domain, such as adding conditional values, it may result in changes in the preference relation between objects, indirectly leading to changes in approximations. In this paper, we effectively addressed the issue of updating approximations that arose from adding conditional values in set-valued ordered decision systems. Firstly, we classified the research objects into two categories: objects with changes in conditional values and objects without changes, and then conducted theoretical studies on updating approximations for these two categories, presenting approximation update theories for adding conditional values. Subsequently, we presented incremental algorithms corresponding to approximation update theories. We demonstrated the feasibility of the proposed incremental update method with numerical examples and showed that our incremental algorithm outperformed the static algorithm. Ultimately, by comparing experimental results on different datasets, it is evident that the incremental algorithm efficiently reduced processing time. In conclusion, this study offered a promising strategy to address the challenges of set-valued ordered decision systems in dynamic environments.

基于梯度增量level-set方法的复杂物体入水砰击数值模拟

复杂结构物自由表面砰击问题,如空投鱼雷、水空无人航行器潜水和高速船出入水,一直是海洋工程领域的研究热点。本文采用梯度增量level-set和虚拟网格法数值模拟复杂结构物的入水砰击,采用时间半隐式有限差分法求解不可压缩Navier-Stokes方程,虚拟网格法可重构复杂结构物界面插值模块并施加无滑移边界条件,梯度增量level-set法可捕捉非线性自由表面翻卷、射流等物理现象。基于以上方法模拟二维圆柱的常速入水砰击,并与试验和数值结果对比,验证本文数值方法的正确性。最后,数值模拟二维船型剖面和飞机剖面的入水砰击,分析了冲击压力、运动响应、压力分布和自由表面运动随入水角度的变化规律。同时,通过模拟可观察到复杂结构浸入水中后的流体分离、射流和小角度空气腔等典型砰击物理现象。

Fenlong Cultivation-the Fourth Set of Farming Methods Invented in China

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.

Topological and Shape Optimization of Flexure Hinges for Designing Compliant Mechanisms Using the Level Set Method

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.

A LEVEL SET METHOD FOR STRUCTURAL TOPOLOGY OPTIMIZATION WITH MULTI-CONSTRAINTS AND MULTI-MATERIALS

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.

Level set method for numerical simulation of a cavitation bubble,its growth, collapse and rebound near a rigid wall

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.

Structural Topology Design of Container Ship Based on Knowledge-Based Engineering and Level Set Method

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.

Simulation of Non-isothermal Injection Molding for a Non-Newtonian Fluid by Level Set Method

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.

Fully implicational methods for approximate reasoning based on interval-valued fuzzy sets

The aim of this paper is to discuss the approximate rea- soning problems with interval-valued fuzzy environments based on the fully implicational idea. First, this paper constructs a class of interval-valued fuzzy implications by means of a type of impli- cations and a parameter on the unit interval, then uses them to establish fully implicational reasoning methods for interval-valued fuzzy modus ponens （IFMP） and interval-valued fuzzy modus tel- lens （IFMT） problems. At the same time the reversibility properties of these methods are analyzed and the reversible conditions are given. It is shown that the existing unified forms of α-triple I （the abbreviation of triple implications） methods for FMP and FMT can be seen as the particular cases of our methods for IFMP and IFMT.

LEVEL SET METHOD FOR TOPOLOGICAL OPTIMIZATION APPLYING TO STRUCTURE, MECHANISM AND MATERIAL DESIGNS

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.

A Simple Weighted Integration Method for Calculating Surface Tension Force to Suppress Parasitic Flow in the Level Set Approach

Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded interface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.

Numerical simulation of single bubble rising in shear-thinning fluids by level set method

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.

A Simple Weighted Integration Method for Calculating Surface Tension Force to Suppress Parasitic Flow in the Level Set Approach

Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded in- terface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.