An adaptive method for high-resolution topology design

For the purpose of achieving high-resolution optimal solutions this paper proposes a nodal design variablebased adaptive method for topology optimization of continuum structures. The analysis mesh-independent density field, interpolated by the nodal design variables at a given set of density points, is adaptively refined/coarsened accord- ing to a criterion regarding the gray-scale measure of local regions. New density points are added into the gray regions and redundant ones are removed from the regions occupied by purely solid/void phases for decreasing the number of de- sign variables. A penalization factor adaptivity technique is employed-to prevent premature convergence of the optimiza- tion iterations. Such an adaptive scheme not only improves the structural boundary description quality, but also allows for sufficient further topological evolution of the structural layout in higher adaptivity levels and thus essentially enables high-resolution solutions. Moreover, compared with the case with uniformly and finely distributed density points, the proposed adaptive method can achieve a higher numerical efficiency of the optimization process.

Expediting carbon dots synthesis by the active adaptive method with machine learning and applications in dental diagnosis and treatment

Synthesis of functional nanostructures with the least number of tests is paramount towards the propelling materials development. However, the synthesis method containing multivariable leads to high uncertainty, exhaustive attempts, and exorbitant manpower costs. Machine learning (ML) burgeons and provokes an interest in rationally designing and synthesizing materials. Here, we collect the dataset of nano-functional materials carbon dots (CDs) on synthetic parameters and optical properties. ML is applied to assist the synthesis process to enhance photoluminescence quantum yield (QY) by building the methodology named active adaptive method (AAM), including the model selection, max points screen, and experimental verification. An interactive iteration strategy is the first time considered in AAM with the constant acquisition of the furnished data by itself to perfect the model. CDs exhibit a strong red emission with QY up to 23.3% and enhancement of around 200% compared with the pristine value obtained through the AAM guidance. Furthermore, the guided CDs are applied as metal ions probes for Co^(2+) and Fe^(3+), with a concentration range of 0–120 and 0–150 µM, and their detection limits are 1.17 and 0.06 µM. Moreover, we also apply CDs for dental diagnosis and treatment using excellent optical ability. It can effectively detect early caries and treat mineralization combined with gel. The study shows that the error of experiment verification gradually decreases and QY improves double with the effective feedback loops by AAM, suggesting the great potential of utilizing ML to guide the synthesis of novel materials. Finally, the code is open-source and provided to be referenced for further investigation on the novel inorganic material prediction.

An adaptive finite-difference method for seismic traveltime modeling based on 3D eikonal equation

3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.

Coupled-generalized nonlinear Schr¨odinger equations solved by adaptive step-size methods in interaction picture

We extend two adaptive step-size methods for solving two-dimensional or multi-dimensional generalized nonlinear Schr ¨odinger equation(GNLSE): one is the conservation quantity error adaptive step-control method(RK4IP-CQE), and the other is the local error adaptive step-control method(RK4IP-LEM). The methods are developed in the vector form of fourthorder Runge–Kutta iterative scheme in the interaction picture by converting a vector equation in frequency domain. By simulating the supercontinuum generated from the high birefringence photonic crystal fiber, the calculation accuracies and the efficiencies of the two adaptive step-size methods are discussed. The simulation results show that the two methods have the same global average error, while RK4IP-LEM spends more time than RK4IP-CQE. The decrease of huge calculation time is due to the differences in the convergences of the relative photon number error and the approximated local error between these two adaptive step-size algorithms.

A novel adaptive harmonic balance method with an asymptotic harmonic selection

The harmonic balance method(HBM)is one of the most widely used methods in solving nonlinear vibration problems,and its accuracy and computational efficiency largely depend on the number of the harmonics selected.The adaptive harmonic balance(AHB)method is an improved HBM method.This paper presents a modified AHB method with the asymptotic harmonic selection(AHS)procedure.This new harmonic selection procedure selects harmonics from the frequency spectra of nonlinear terms instead of estimating the contribution of each harmonic to the whole nonlinear response,by which the additional calculation is avoided.A modified continuation method is proposed to deal with the variable size of nonlinear algebraic equations at different values of path parameters,and then all solution branches of the amplitude-frequency response are obtained.Numerical experiments are carried out to verify the performance of the AHB-AHS method.Five typical nonlinear dynamic equations with different types of nonlinearities and excitations are chosen as the illustrative examples.Compared with the classical HBM and Runge-Kutta methods,the proposed AHB-AHS method is of higher accuracy and better convergence.The AHB-AHS method proposed in this paper has the potential to investigate the nonlinear vibrations of complex high-dimensional nonlinear systems.

Sensitivity Analysis of Electromagnetic Scattering from Dielectric Targets with Polynomial Chaos Expansion and Method of Moments

In this paper,an adaptive polynomial chaos expansion method(PCE)based on the method of moments(MoM)is proposed to construct surrogate models for electromagnetic scattering and further sensitivity analysis.The MoM is applied to accurately solve the electric field integral equation(EFIE)of electromagnetic scattering from homogeneous dielectric targets.Within the bistatic radar cross section(RCS)as the research object,the adaptive PCE algorithm is devoted to selecting the appropriate order to construct the multivariate surrogate model.The corresponding sensitivity results are given by the further derivative operation,which is compared with those of the finite difference method(FDM).Several examples are provided to demonstrate the effectiveness of the proposed algorithm for sensitivity analysis of electromagnetic scattering from homogeneous dielectric targets.

Phase-field modeling of dendritic growth under forced flow based on adaptive finite element method

A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.

Global Optimization Method Using SLE and Adaptive RBF Based on Fuzzy Clustering

High fidelity analysis models,which are beneficial to improving the design quality,have been more and more widely utilized in the modern engineering design optimization problems.However,the high fidelity analysis models are so computationally expensive that the time required in design optimization is usually unacceptable.In order to improve the efficiency of optimization involving high fidelity analysis models,the optimization efficiency can be upgraded through applying surrogates to approximate the computationally expensive models,which can greately reduce the computation time.An efficient heuristic global optimization method using adaptive radial basis function（RBF） based on fuzzy clustering（ARFC） is proposed.In this method,a novel algorithm of maximin Latin hypercube design using successive local enumeration（SLE） is employed to obtain sample points with good performance in both space-filling and projective uniformity properties,which does a great deal of good to metamodels accuracy.RBF method is adopted for constructing the metamodels,and with the increasing the number of sample points the approximation accuracy of RBF is gradually enhanced.The fuzzy c-means clustering method is applied to identify the reduced attractive regions in the original design space.The numerical benchmark examples are used for validating the performance of ARFC.The results demonstrates that for most application examples the global optima are effectively obtained and comparison with adaptive response surface method（ARSM） proves that the proposed method can intuitively capture promising design regions and can efficiently identify the global or near-global design optimum.This method improves the efficiency and global convergence of the optimization problems,and gives a new optimization strategy for engineering design optimization problems involving computationally expensive models.

Adaptive element free Galerkin method applied to analysis of earthquake induced liquefaction

An automatically adaptive element free method is presented to analyze the seismic response of liquefiable soils. The method is based on the element free Galerkin method （EFGM） and the fission procedure that is part of h-refinement, indicated by error estimation. In the proposed method, a posteriori error estimate procedure that depends on the energy norm of stress and the T-Belytschko （TB） stress recovery scheme is incorporated. The effective cyclic elasto-plastic constitutive model is used to describe the nonlinear behavior of the saturated soil. The governing equations are established by u-p formulation. The proposed method can effectively avoid the volumetric locking due to large deformation that usually occurs in numerical computations using the finite element method （FEM）. The efficiency of the proposed method is demonstrated by evaluating the seismic response of an embankment and comparing it to results obtained through FEM. It is shown that the proposed method provides an accurate seismic analysis of saturated soil that includes the effects of liquefaction .

Key Techniques and Applications of Adaptive Growth Method for Stiffener Layout Design of Plates and Shells

The application of the adaptive growth method is limited because several key techniques during the design process need manual intervention of designers. Key techniques of the method including the ground structure construction and seed selection are studied, so as to make it possible to improve the effectiveness and applicability of the adaptive growth method in stiffener layout design optimization of plates and shells. Three schemes of ground structures, which are comprised by different shell elements and beam elements, are proposed. It is found that the main stiffener layouts resulted from different ground structures are almost the same, but the ground structure comprised by 8-nodes shell elements and both 3-nodes and 2-nodes beam elements can result in clearest stiffener layout, and has good adaptability and low computational cost. An automatic seed selection approach is proposed, which is based on such selection rules that the seeds should be positioned on where the structural strain energy is great for the minimum compliance problem, and satisfy the dispersancy requirement. The adaptive growth method with the suggested key techniques is integrated into an ANSYS-based program, which provides a design tool for the stiffener layout design optimization of plates and shells. Typical design examples, including plate and shell structures to achieve minimum compliance and maximum bulking stability are illustrated. In addition, as a practical mechanical structural design example, the stiffener layout of an inlet structure for a large-scale electrostatic precipitator is also demonstrated. The design results show that the adaptive growth method integrated with the suggested key techniques can effectively and flexibly deal with stiffener layout design problem for plates and shells with complex geometrical shape and loading conditions to achieve various design objectives, thus it provides a new solution method for engineering structural topology design optimization.

Adaptive hp finite element method for fluorescence molecular tomography with simplified spherical harmonics approximation

Recently,the simplified spherical harmonics equations(SP)model has at tracted much att entionin modeling the light propagation in small tissue ggeometriesat visible and near-infrared wave-leng ths.In this paper,we report an eficient numerical method for fluorescence moleeular tom-ography(FMT)that combines the advantage of SP model and adaptive hp finite elementmethod(hp-FEM).For purposes of comparison,hp-FEM and h-FEM are,respectively applied tothe reconstruction pro cess with diffusion approximation and SPs model.Simulation experiments on a 3D digital mouse atlas and physical experiments on a phantom are designed to evaluate thereconstruction methods in terms of the location and the reconstructed fluorescent yield.Theexperimental results demonstrate that hp-FEM with SPy model,yield more accurate results thanh-FEM with difusion approximation model does.The phantom experiments show the potentialand feasibility of the proposed approach in FMT applications.

Adaptive Finite Element Method Based on Optimal Error Estimates for Linear Elliptic Problems on Nonconvex Polygonal Domains

The subject of this work is to propose adaptive finite element methods based on an optimal maximum norm error control estimate.Using estimators of the local regularity of the unknown exact solution derived from computed approximate solutions,the proposed procedures are analyzed in detail for a non-trivial class of corner problems and shown to be efficient in the sense that they generate the correct type of refinement and lead to the desired control under consideration.

Adaptive split-step Fourier method for simulating ultrashort laser pulse propagation in photonic crystal fibres

In this paper, the generalized nonlinear Schrodinger equation （GNLSE） is solved by an adaptive split-step Fourier method （ASSFM）. It is found that ASSFM must be used to solve GNLSE to ensure precision when the soliton selffrequency shift is remarkable and the photonic crystal fibre （PCF） parameters vary with the frequency considerably. The precision of numerical simulation by using ASSFM is higher than that by using split-step Fourier method in the process of laser pulse propagation in PCFs due to the fact that the variation of fibre parameters with the peak frequency in the pulse spectrum can be taken into account fully.

Prediction of Rural Residents’ Consumption Expenditure Based on Lasso and Adaptive Lasso Methods

When the variable of model is large, the Lasso method and the Adaptive Lasso method can effectively select variables. This paper prediction the rural residents’ consumption expenditure in China, based on respectively using the Lasso method and the Adaptive Lasso method. The results showed that both can effectively and accurately choose the appropriate variable, but the Adaptive Lasso method is better than the Lasso method in prediction accuracy and prediction error. It shows that in variable selection and parameter estimation, Adaptive Lasso method is better than the Lasso method.

METHOD FOR ADAPTIVE MESH GENERATION BASED ON GEOMETRICAL FEATURES OF 3D SOLID

In order to provide a guidance to specify the element size dynamically during adaptive finite element mesh generation, adaptive criteria are firstly defined according to the relationships between the geometrical features and the elements of 3D solid. Various modes based on different datum geometrical elements, such as vertex, curve, surface, and so on, are then designed for generating local refined mesh. With the guidance of the defmed criteria, different modes are automatically selected to apply on the appropriate datum objects to program the element size in the local special areas. As a result, the control information of element size is successfully programmed covering the entire domain based on the geometrical features of 3D solid. A new algorithm based on Delatmay triangulation is then developed for generating 3D adaptive finite element mesh, in which the element size is dynamically specified to catch the geometrical features and suitable tetrahedron facets are selected to locate interior nodes continuously. As a result, adaptive mesh with good-quality elements is generated. Examples show that the proposed method can be successfully applied to adaptive finite element mesh automatic generation based on the geometrical features of 3D solid.

The Adaptive Wavelet Collocation Method and Its Application in Front Simulation

The adaptive wavelet collocation method （AWCM） is a variable grid technology for solving partial differential equations （PDEs） with high singularities. Based on interpolating wavelets, the AWCM adapts the grid so that a higher resolution is automatically attributed to domain regions with high singularities. Accuracy problems with the AWCM have been reported in the literature, and in this paper problems of efficiency with the AWCM are discussed in detail through a simple one-dimensional （1D） nonlinear advection equation whose analytic solution is easily obtained. A simple and efficient implementation of the AWCM is investigated. Through studying the maximum errors at the moment of frontogenesis of the 1D nonlinear advection equation with different initial values and a comparison with the finite difference method （FDM） on a uniform grid, the AWCM shows good potential for modeling the front efficiently. The AWCM is also applied to a two-dimensional （2D） unbalanced frontogenesis model in its first attempt at numerical simulation of a meteorological front. Some important characteristics about the model are revealed by the new scheme.

ADAPTIVE FINITE ELEMENT METHOD BASED ON OPTIMAL ERROR ESTIMATES FOR LINEAR ELLIPTIC PROBLEMS ON CONCAVE CORNER DOMAINS

An adaptive finite element procedure designed for specific computational goals is presented,using mesh refinement strategies based on optimal or nearly optimal a priori error estimates for the finite element method and using estimators of the local regularity of the unknown exact solution derived from computed approximate solutions.The proposed procedure is analyzed in detail for a non-trivial class of corner problems and shown to be efficient in the sense that the method can generate the correct type of refinements and lead to the desired control under consideration.

Quasi-optimal complexity of adaptive finite element method for linear elasticity problems in two dimensions

This paper introduces an adaptive finite element method （AFEM） using the newest vertex bisection and marking exclusively according to the error estimator without special treatment of oscillation. By the combination of the global lower bound and the localized upper bound of the posteriori error estimator, perturbation of oscillation, and cardinality of the marked element set, it is proved that the AFEM is quasi-optimal for linear elasticity problems in two dimensions, and this conclusion is verified by the numerical examples.

An Adaptive Finite Element Method Based on Optimal Error Estimates for Linear Elliptic Problems

The subject of the work is to propose a series of papers about adaptive finite element methods based on optimal error control estimate. This paper is the third part in a series of papers on adaptive finite element methods based on optimal error estimates for linear elliptic problems on the concave corner domains. In the preceding two papers (part 1:Adaptive finite element method based on optimal error estimate for linear elliptic problems on concave corner domain; part 2:Adaptive finite element method based on optimal error estimate for linear elliptic problems on nonconvex polygonal domains), we presented adaptive finite element methods based on the energy norm and the maximum norm. In this paper, an important result is presented and analyzed. The algorithm for error control in the energy norm and maximum norm in part 1 and part 2 in this series of papers is based on this result.

PARALLEL ADAPTIVELY MODIFIED CHARACTERISTIC BASIS FUNCTION METHOD BASED ON STATIC LOAD BALANCE

Characteristic Basis Function Method (CBFM) is a novel approach for analyzing the ElectroMagnetic (EM) scattering from electrically large objects. Based on dividing the studied object into small blocks, the CBFM is suitable for parallel computing. In this paper, a static load balance parallel method is presented by combining Message Passing Interface (MPI) with Adaptively Modified CBFM (AMCBFM). In this method, the object geometry is partitioned into distinct blocks, and the serial number of blocks is sent to related nodes according to a certain rule. Every node only needs to calculate the information on local blocks. The obtained results confirm the accuracy and efficiency of the proposed method in speeding up solving large electrical scale problems.