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 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.

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.

A LEVEL SET METHOD FOR MICROSTRUCTURE DESIGN OF COMPOSITE MATERIALS

Based on a level set model and the homogenization theory, an optimization al- gorithm for ?nding the optimal con?guration of the microstructure with speci?ed properties is proposed, which extends current research on the level set method for structure topology opti- mization. The method proposed employs a level set model to implicitly describe the material interfaces of the microstructure and a Hamilton-Jacobi equation to continuously evolve the ma- terial interfaces until an optimal design is achieved. Meanwhile, the moving velocities of level set are obtained by conducting sensitivity analysis and gradient projection. Besides, how to handle the violated constraints is also discussed in the level set method for topological optimization, and a return-mapping algorithm is constructed. Numerical examples show that the method exhibits outstanding ?exibility of handling topological changes and ?delity of material interface represen- tation as compared with other conventional methods in literatures.

Topology and Shape Optimization of 2-D and 3-D Micro-ArchitecturedThermoelastic Metamaterials Using a Parametric Level Set Method

2-D and 3-D micro-architectured multiphase thermoelastic metamaterials are designed and analyzed using a parametric level set method for topology optimization and the finite element method.An asymptotic homogenization approach is employed to obtain the effective thermoelastic properties of the multiphase metamaterials.Theε-constraint multi-objective optimization method is adopted in the formulation.The coefficient of thermal expansion(CTE)and Poisson’s ratio(PR)are chosen as two objective functions,with the CTE optimized and the PR treated as a constraint.The optimization problems are solved by using the method of moving asymptotes.Effective isotropic and anisotropic CTEs and stiffness constants are obtained for the topologically optimized metamaterials with prescribed values of PR under the constraints of specified effective bulk modulus,volume fractions and material symmetry.Two solid materials along with one additional void phase are involved in each of the 2-D and 3-D optimal design examples.The numerical results reveal that the newly proposed approach can integrate shape and topology optimizations and lead to optimal microstructures with distinct topological boundaries.The current method can topologically optimize metamaterials with a positive,negative or zero CTE and a positive,negative or zero Poisson’s ratio.

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.

BSN:Photometric Light Curve Analysis of Two Contact Binary Systems LS Del and V997 Cyg

The light curve analyses and orbital period variations for two contact binary stars,LS Del and V997 Cyg,are presented in this work which was conducted in the frame of the Binary Systems of South and North project Ground-based photometric observations were performed at two observatories in France.We used the Transiting Exoplanet Survey Satellite(TESS)data for extracting times of minima and light curve analysis of the targe systems.The O-C diagram for both systems displays a parabolic trend.LS Del and V997 Cyg’s orbital periods are increasing at rates of dP/dt=7.20×10^(-08)days yr^(-1)and dP/dt=2.54×10^(-08)days yr^(-1),respectively Therefore,it can be concluded that mass is being transferred from the less massive star to the more massive component with a rate of dM/dt=-1.96×10^(-7)M_(⊙)yr^(-1)for the LS Del system,and dM/dt=-3.83×10^(-7)M_(⊙)yr-1for V997 Cyg.The parameters of a third possible object in the system are also considered.The PHysics Of Eclipsing BinariEs Python code was used to analyze the light curves.The light curve solutions needed a cold starspot due to the asymmetry in the LS Del system’s light curve maxima.The mass ratio fill-out factor,and star temperature all indicate that both systems are contact binary types in this investigation.Two methods were applied to estimate the absolute parameters of the systems:one method relied on the parallax of Gaia DR3,and the other used a P-M relationship.The positions of the systems are also depicted on the M-L,M-R q-L_(ratio),and logM_(tot)-logJ_(0)diagrams.We recommend that further observations and investigations be done on the existence of a fourth body in this system.

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法可捕捉非线性自由表面翻卷、射流等物理现象。基于以上方法模拟二维圆柱的常速入水砰击,并与试验和数值结果对比,验证本文数值方法的正确性。最后,数值模拟二维船型剖面和飞机剖面的入水砰击,分析了冲击压力、运动响应、压力分布和自由表面运动随入水角度的变化规律。同时,通过模拟可观察到复杂结构浸入水中后的流体分离、射流和小角度空气腔等典型砰击物理现象。

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.

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.

A two-dimensional numerical wave flume based on SA-MPLS method

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.

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.

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 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.

New normalized nonlocal hybrid level set method for image segmentation

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.

The shape optimization of the arterial graft design by level set methods

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.