2D multi-scale hybrid optimization method for geophysical inversion and its application

Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm （SAGA） and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A＇nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.

THREE-DIMENSIONAL TRAJECTORY OPTIMIZATION WITH DIRECT METHOD

The principle of direct method used in optimal control problem is introduced. Details of applying this method to flight trajectory generation are presented including calculation of velocity and controls histories. And capabilities of flight and propulsion systems are considered also. Combined with digital terrain map technique, the direct method is applied to the three dimensional trajectory optimization for low altitude penetration, and simplex algorithm is used to solve the parameters in optimization. For the small number of parameters, the trajectory can be optimized in real time on board.

GLOBAL CONVERGENCE OF A CLASS OF OPTIMALLY CONDITIONED SSVM METHODS

This paper explores the convergence of a class of optimally conditioned self scaling variable metric (OCSSVM) methods for unconstrained optimization. We show that this class of methods with Wolfe line search are globally convergent for general convex functions.

PARETO FRONT CAPTURE USING DETERMINISTIC OPTIMIZATION METHODS IN MULTI-CRITERION AERODYNAMIC DESIGN

Deterministic optimization methods are combined with the Pareto front concept to solve multi-criterion design problems. The algorithm and the numerical implementation are applied to aerodynamic designs. Evolutionary algorithms （EAs） and the Pareto front concept are used to solve practical design problems in industry for its robustness in capturing convex, concave, discrete or discontinuous Pareto fronts of multi-objective optimization problems. However, the process is time-consuming. Therefore, deterministic optimization methods are introduced to capture the Pareto front, and the types of the captured Pareto front are explained. Numerical experiments show that the deterministic optimization method is a good alternative to EAs for capturing any convex and some concave Pareto fronts in multi-criterion aerodynamic optimization problems due to its efficiency.

Methods for extraction of microorganism DNA from glacier surface snow

In order to thoroughly investigate the diversity of glacier microorganisms, four DNA extraction methods with differem lysis pat- terns were tested and two screened methods （the Bosshard-Bano method and the Zhou method） were optimized for the most ef- fective form of the filter membrane （cut vs. uncut）, the DNA extraction method, and the precipitation method. The two optimized methods were then compared with the commercial Mo-Bio DNA extraction kit, and the results showed that the kit was generally suitable for extraction of microorganism DNA fi＇om glacier surface snow. Procedurally, it was found that a modified Boss- hard-Bano method （i.e., cutting the filter membrane into pieces, using a specific lysis pattern [lysozyme （5 mg/mL）-protease K （ 1 mg/mL）-CTAB （ 1%）-SDS （ 1%）], performing the extraction only once by chloroform-isoamyl alcohol （24： 1）, and conducting DNA precipitation by pure ethanol） was also an effective and less expensive method for extraction of microorganism DNA from glacier surface snow.

A NEW DERIVATIVE FREE OPTIMIZATION METHOD BASED ON CONIC INTERPOLATION MODEL

In this paper, a new derivative free trust region method is developed based on the conic interpolation model for the unconstrained optimization. The conic interpolation model is built by means of the quadratic model function, the collinear scaling formula, quadratic approximation and interpolation. All the parameters in this model are determined by objective function interpolation condition. A new derivative free method is developed based upon this model and the global convergence of this new method is proved without any information on gradient.

An improved optimal elemental method for updating finite element models

The optimal matrix method and optimal elemental method used to update finite element models may not provide accurate results.This situation occurs when the test modal model is incomplete,as is often the case in practice.An improved optimal elemental method is presented that defines a new objective function,and as a byproduct,circumvents the need for mass normalized modal shapes,which are also not readily available in practice.To solve the group of nonlinear equations created by the improved optimal method,the Lagrange multiplier method and Matlab function fmincon are employed.To deal with actual complex structures, the float-encoding genetic algorithm(FGA)is introduced to enhance the capability of the improved method.Two examples,a 7- degree of freedom(DOF)mass-spring system and a 53-DOF planar frame,respectively,are updated using the improved method. The example results demonstrate the advantages of the improved method over existing optimal methods,and show that the genetic algorithm is an effective way to update the models used for actual complex structures.

PARAMETER IDENTIFICATION BY OPTIMIZATION METHOD FOR A POLLUTION PROBLEM IN POROUS MEDIA

In the present work, we investigate the inverse problem of reconstructing the parameter of an integro-differential parabolic equation, which comes from pollution problems in porous media, when the final observation is given. We use the optimal control framework to establish both the existence and necessary condition of the minimizer for the cost func- tional. Furthermore, we prove the stability and the local uniqueness of the minimizer. Some numerical results will be presented and discussed.

Fuzzy pattern recognition method for assessing groundwater vulnerability to pollution in the Zhangji area

Based on the widely used DRASTIC method, a fuzzy pattern recognition and optimization method was proposed and applied to the fissured-karstic aquifer of Zhangji area for assessing groundwater vulnerability to pollution. The result is compared with DRASTIC method. It is shown that by taking the fuzziness into consideration, the fuzzy pattern recognition and optimization method reflects more efficiently the fuzzy nature of the groundwater vulnerability to pollution and is more applicable in reality.

The ICM method with objective function transformed by variable discrete condition for continuum structure

ICM （Independent Continuous Mapping） method can solve topological optimization problems with the minimized weight as the objective and subjected to displacement constraints. To get a clearer topological configuration, by introducing the discrete condition of topological variables and integrating with the original objective, an optimal model with multi-objectives is formulated to make the topological variables approach 0 or 1 as near as possible, and the model reduces the effect of deleting rate on the result. The image-filtering method is employed to eliminate the checkerboard patterns and mesh dependence that occurred in the topology optimization of a continuum structure. The computational efficiency is enhanced through selecting quasi-active displacement constraints and a design region. Numerical examples indicate that this algorithm is robust and practicable, though the number of iterations is slightly increased with respect to the original algorithm.

A New Rational-based Optimal Design Strategy of Ship Structure Based on Multi-level Analysis and Super-element Modeling Method

A new multi-level analysis method of introducing the super-element modeling method, derived from the multi-level analysis method first proposed by O. F. Hughes, has been proposed in this paper to solve the problem of high time cost in adopting a rational-based optimal design method for ship structural design. Furthermore,the method was verified by its effective application in optimization of the mid-ship section of a container ship. A full 3-D FEM model of a ship,suffering static and quasi-static loads, was used as the analyzing object for evaluating the structural performance of the mid-ship module, including static strength and buckling performance. Research results reveal that this new method could substantially reduce the computational cost of the rational-based optimization problem without decreasing its accuracy, which increases the feasibility and economic efficiency of using a rational-based optimal design method in ship structural design.

Development of hybrid optimization algorithm for structures furnished with seismic damper devices using the particle swarm optimization method and gravitational search algorithm

Previous studies about optimizing earthquake structural energy dissipation systems indicated that most existing techniques employ merely one or a few parameters as design variables in the optimization process,and thereby are only applicable only to simple,single,or multiple degree-of-freedom structures.The current approaches to optimization procedures take a specific damper with its properties and observe the effect of applying time history data to the building;however,there are many different dampers and isolators that can be used.Furthermore,there is a lack of studies regarding the optimum location for various viscous and wall dampers.The main aim of this study is hybridization of the particle swarm optimization(PSO) and gravitational search algorithm(GSA) to optimize the performance of earthquake energy dissipation systems(i.e.,damper devices) simultaneously with optimizing the characteristics of the structure.Four types of structural dampers device are considered in this study:(ⅰ) variable stiffness bracing(VSB) system,(ⅱ) rubber wall damper(RWD),(ⅲ) nonlinear conical spring bracing(NCSB) device,(iv) and multi-action stiffener(MAS) device.Since many parameters may affect the design of seismic resistant structures,this study proposes a hybrid of PSO and GSA to develop a hybrid,multi-objective optimization method to resolve the aforementioned problems.The characteristics of the above-mentioned damper devices as well as the section size for structural beams and columns are considered as variables for development of the PSO-GSA optimization algorithm to minimize structural seismic response in terms of nodal displacement(in three directions) as well as plastic hinge formation in structural members simultaneously with the weight of the structure.After that,the optimization algorithm is implemented to identify the best position of the damper device in the structural frame to have the maximum effect and minimize the seismic structure response.To examine the performance of the proposed PSO-GSA optimization method,it has been applied to a three-story reinforced structure equipped with a seismic damper device.The results revealed that the method successfully optimized the earthquake energy dissipation systems and reduced the effects of earthquakes on structures,which significantly increase the building’s stability and safety during seismic excitation.The analysis results showed a reduction in the seismic response of the structure regarding the formation of plastic hinges in structural members as well as the displacement of each story to approximately 99.63%,60.5%,79.13% and 57.42% for the VSB device,RWD,NCSB device,and MAS device,respectively.This shows that using the PSO-GSA optimization algorithm and optimized damper devices in the structure resulted in no structural damage due to earthquake vibration.

Inverse displacement analysis of the general six degree-of-freedom serial robot based on optimization method

The paper presents a new solution of inverse displacement analysis of the general six degree-of-freedom serial robot.The inverse displacement analysis of the general serial robot is transformed into a minimization problem and then the optimization method is adopted to solve the nonlinear least squares problem with the analytic form of new Jacobian matrix.In this way,joint variables of the general serial robot can be searched out quickly under the desired precision when positions of the three non-collinear end effector points are given.Compared with the general Newton iterative method,the proposed algorithm can search out the solution when the robot is at the singular configuration and the initial configuration used in the optimization method may also be the singular configuration.So the convergence domain is bigger than that of the general Newton iterative method.Another advantage of the proposed algorithm is that positions of the three non-collinear end effector points are usually much easier to be measured than the orientation of the end effector.The inverse displacement analysis of the general 6R(six-revolute-joint) serial robot is illustrated as an example and the simulation results verify the efficiency of the proposed algorithm.Because the three non-collinear points can be selected at random,the method can be applied to any other types of serial robots.

Preconcentration of rare earth elements from Iranian monazite ore by spiral separator using multi-response optimization method

The present work dealt with the preconcentration of rare earth elements in Saghand ore(Yazd province,Iran)which was achieved by Humphrey spiral using orthogonal optimization method after scrubbing the sample at 45%solid pulp density for 30 min.The pulp was diluted and was fed to a Humphrey spiral for upgrading.The process parameters considered were feed size,feed solids and feed rate,and Taguchi’s L9(34)orthogonal array(OA)was selected for optimization of the process.The results show that the feed rate and feed size were more significant than the other operation parameters of the process.It was also found that under optimal conditions,the concentrate grade of rare earth elements increased from2860 10 6to 6050 10 6and recovery reached to 58%.

Six-DOF trajectory optimization for reusable launch vehicles via Gauss pseudospectral method

To be close to the practical flight process and increase the precision of optimal trajectory, a six-degree-offreedom（6-DOF） trajectory is optimized for the reusable launch vehicle（RLV） using the Gauss pseudospectral method（GPM）. Different from the traditional trajectory optimization problem which generally considers the RLV as a point mass, the coupling between translational dynamics and rotational dynamics is taken into account. An optimization problem is formulated to minimize a performance index subject to 6-DOF equations of motion, including translational and rotational dynamics. A two-step optimal strategy is then introduced to reduce the large calculations caused by multiple variables and convergence confinement in 6-DOF trajectory optimization. The simulation results demonstrate that the 6-DOF trajectory optimal strategy for RLV is feasible.

Shape optimization of axisymmetric solids with the finite cell method using a fixed grid

In this work, a design procedure extending the B-spline based finite cell method into shape optimization is developed for axisymmetric solids involving the centrifugal force effect. We first replace the traditional conforming mesh in the finite element method with structured cells that are fixed during the whole design process with a view to avoid the sophisticated re-meshing and eventual mesh distortion.Then, B-spline shape functions are further implemented to yield a high-order continuity field along the cell boundary in stress analysis. By means of the implicit description of the shape boundary, stress sensitivity is analytically derived with respect to shape design variables. Finally, we illustrate the efficiency and accuracy of the proposed protocol by several numerical test cases as well as a whole design procedure carried out on an aeronautic turbine disk.

Research on Design Method of the Full Form Ship with Minimum Thrust Deduction Factor

In the preliminary design stage of the full form ships, in order to obtain a hull form with low resistance and maximum propulsion efficiency, an optimization design program for a full form ship with the minimum thrust deduction factor has been developed, which combined the potential flow theory and boundary layer theory with the optimization technique. In the optimization process, the Sequential Unconstrained Minimization Technique（SUMT） interior point method of Nonlinear Programming（NLP） was proposed with the minimum thrust deduction factor as the objective function. An appropriate displacement is a basic constraint condition, and the boundary layer separation is an additional one. The parameters of the hull form modification function are used as design variables. At last, the numerical optimization example for lines of after-body of 50000 DWT product oil tanker was provided, which indicated that the propulsion efficiency was improved distinctly by this optimal design method.

An optimization method for passive muzzle arc control devices in augmented railguns

In this paper,a passive muzzle arc control device(PMACD)of the augmented railguns is studied.By discussing its performance at different numbers of extra rails,a parameter optimization model is proposed.Through the calculation model,it is found that the PMACD works well in the simple railgun,which refers to the gun that there is only one pair of rails in the inner bore.The PMACD may decrease the simple railgun’s armature peak current and muzzle arc,but affect its muzzle velocity not much.However,in the augmented railguns it has different characteristics.If the parameters of the PMACD are not selected suitable.It may increase the armature peak current and muzzle arc,but greatly decrease the velocity.The reason for this problem is that the extra rails generate a strong magnetic field in front of the armature,which induces a large current to change the armature current.It is also found that when the resistance and inductance parameters of the PMACD satisfy with the optimization formula,the PMACD can also play a good role in arc suppression in the augmented railguns.Experiments of an augmented railgun with a stainless steel PMACD are carried out to verify this optimization method.Results show that the muzzle arc is obviously controlled.This work may provide a reference for the design of the muzzle arc control device.

Design and optimization of a SiC thermal emitter/absorber composed of periodic microstructures based on a non-linear method

Spectral and directional control of thermal emission based on excitation of confined electromagnetic resonant modes paves a viable way for the design and construction of microscale thermal emitters/absorbers. In this paper, we present numerical simulation results of the thermal radiative properties of a silicon carbide（Si C） thermal emitter/absorber composed of periodic microstructures. We illustrate different electromagnetic resonant modes which can be excited with the structure,such as surface phonon polaritons, magnetic polaritons and photonic crystal modes, and the process of radiation spectrum optimization based on a non-linear optimization algorithm. We show that the spectral and directional control of thermal emission/absorption can be efficiently achieved by adjusting the geometrical parameters of the structure. Moreover, the optimized spectrum is insensitive to 3% dimension modification.

A SUPERLINEARLY CONVERGENT SPLITTING FEASIBLE SEQUENTIAL QUADRATIC OPTIMIZATION METHOD FOR TWO-BLOCK LARGE-SCALE SMOOTH OPTIMIZATION

This paper discusses the two-block large-scale nonconvex optimization problem with general linear constraints.Based on the ideas of splitting and sequential quadratic optimization(SQO),a new feasible descent method for the discussed problem is proposed.First,we consider the problem of quadratic optimal(QO)approximation associated with the current feasible iteration point,and we split the QO into two small-scale QOs which can be solved in parallel.Second,a feasible descent direction for the problem is obtained and a new SQO-type method is proposed,namely,splitting feasible SQO(SF-SQO)method.Moreover,under suitable conditions,we analyse the global convergence,strong convergence and rate of superlinear convergence of the SF-SQO method.Finally,preliminary numerical experiments regarding the economic dispatch of a power system are carried out,and these show that the SF-SQO method is promising.