Search results

The purpose of this paper is to present novel search formulations in gradient‐type methods for prediction of boundary heat flux distribution in two‐dimensional nonlinear heat conduction problems.

The performance of gradient‐type methods is strongly contingent upon the effective determination of the search direction. Based on the definition of this parameter, gradient‐based methods such as steepest descent, various versions of both conjugate gradient and quasi‐Newton can be distinguished. By introducing new search techniques, several examples in the presence of noise in data are studied and discussed to verify the accuracy and efficiency of the present strategies.

The verification of the proposed methods for recovering time and space varying heat flux. The performance of the proposed methods via comparisons with the classical methods involved in its derivation.

The innovation of the present method is to use a hybridization of a conjugate gradient and a quasi‐Newton method to determine the search directions in gradient‐based approaches.

The present paper is directed towards elasto‐plastic large deformation analysis of thin shells based on the concept of degenerated solids. The main aspect of the paper is the derivation of an efficient computational strategy placing emphasis on consistent elasto‐plastic tangent moduli and stress integration with the radial return method under the restriction of ‘zero normal stress condition’ in thickness direction. The advantageous performance of the standard Newton iteration using a consistent tangent stiffness matrix is compared to the classical scheme with an iteration matrix based on the infinitesimal elasto‐plastic constitutive tensor. Several numerical examples also demonstrate the effectiveness of the standard Newton iteration with respect to modified and quasi‐Newton methods like BFGS and others.

The present study aims to perform inverse analysis of a conjugate heat transfer problem including conduction and forced convection via the quasi-Newton method. The inverse analysis is defined for a heat source that is surrounded by a solid medium which is exposed to a free stream in external flow.

The objective of the inverse design problem is finding temperature distribution of the heat source as thermal boundary condition to establish a prescribed temperature along the interface of solid body and fluid. This problem is a simplified version of thermal-based ice protection systems in which the formation of ice is avoided by maintaining the interface of fluid and solid at a specified temperature.

The effects of the different pertinent parameters such as Reynolds number, interface temperature and thermal conductivity ratio of fluid and solid mediums are analyzed.

This paper fulfils the analysis to study how thermal based anti-icing system can be used with different heat source shapes.

In a model resulting from Maxwell's equations with a constitutive law using Preisach operators for incorporating magnetization hysteresis, this paper aims at identifying the hysteresis operator, i.e. the Preisach weight function, from indirect measurements.

Dealing with a nonlinear inverse problem, one has to apply iterative methods for its numerical solution. For this purpose several approaches are proposed based on fixed point or Newton type ideas. In the latter case, one has to take into account nondifferentiability of the hysteresis operator. This is done by using differentiable substitutes or quasi‐Newton methods.

Numerical tests with synthetic data show that fixed point methods based on fitting after a full forward sweep (alternating iteration) and Newton type iterations using the hysteresis centerline or commutation curve exhibit a satisfactory convergence behavior, while fixed point iterations based on subdividing the time interval (Kaczmarz) suffer from instability problems and quasi Newton iterations (Broyden) are too slow in some cases.

Application of the proposed methods to measured data will be the subject of future research work.

The proposed methodologies allow to determine material parameters in hysteresis models from indirect measurements.

Taking into account the full PDE model, one can expect to get accurate and reliable results in this model identification problem. Especially the use of Newton type methods – taking into account nondifferentiability – is new in this context.

In this paper, with the goal of reducing the fuel consumption of UAV, the engine performance optimization is studied and on the basis of aircraft/engine integrated control, the minimum fuel consumption optimization method of engine given thrust is proposed. In the case of keeping the given thrust of the engine unchanged, the main fuel flow of the engine without being connected to the afterburner is optimally controlled so as to minimize the fuel consumption.

In this study, the reference model real-time optimization control method is adopted. The engine reference model uses a nonlinear real-time mathematical model of a certain engine component method. The quasi-Newton method is adopted in the optimization algorithm. According to the optimization variable nozzle area, the turbine drop-pressure ratio corresponding to the optimized nozzle area is calculated, which is superimposed with the difference of the drop-pressure ratio of the conventional control plan and output to the conventional nozzle controller of the engine. The nozzle area is controlled by the conventional nozzle controller.

The engine real-time minimum fuel consumption optimization control method studied in this study can significantly reduce the engine fuel consumption rate under a given thrust. At the work point, this is a low-altitude large Mach work point, which is relatively close to the edge of the flight envelope. Before turning on the optimization controller, the fuel consumption is 0.8124 kg/s. After turning on the optimization controller, you can see that the fuel supply has decreased by about 4%. At this time, the speed of the high-pressure rotor is about 94% and the temperature after the turbine can remain stable all the time.

The optimal control method of minimum fuel consumption for the given thrust of UAV is proposed in this paper and the optimal control is carried out for the nozzle area of the engine. At the same time, a method is proposed to indirectly control the nozzle area by changing the turbine pressure ratio. The relevant UAV and its power plant designers and developers may consider the results of this study to reach a feasible solution to reduce the fuel consumption of UAV.

Fuel consumption optimization can save fuel consumption during aircraft cruising, increase the economy of commercial aircraft and improve the combat radius of military aircraft. With the increasingly wide application of UAVs in military and civilian fields, the demand for energy-saving and emission reduction will promote the UAV industry to improve the awareness of environmental protection and reduce the cost of UAV use and operation.

This paper aims to improve the alignment accuracy of large components in aircraft assembly and an evaluation algorithm, which is based on manufacture accuracy and coordination accuracy, is proposed.

With relative deviations of manufacturing feature points and coordinate feature points, an evaluation function of assembly error is constructed. Then the optimization model of large aircraft digital alignment is established to minimize the synthesis assembly error with tolerance requirements, which consist of three-dimensional (3D) tolerance of manufacturing feature points and relative tolerance between coordination feature points. The non-linear constrained optimization problem is solved by Lagrange multiplier method and quasi-Newton method with its initial value provided by the singular value decomposition method.

The optimized postures of large components are obtained, which makes the tolerance of both manufacturing and coordination requirements be met. Concurrently, the synthesis assembly error is minimized. Compared to the result of the singular value decomposition method, the algorithm is validated in three typical cases with practical data.

The proposed method has been used in several aircraft assembly projects and gained a good effect.

This paper proposes a method to optimize the manufacturing and coordination accuracy with tolerance constraints when the postures of several components are adjusted at the same time. The results of this paper will help to improve the quality of component assemblies.

This paper presents a numerical study of inverse parameter identification problems in fracture mechanics. Inverse methodology is applied to the detection of subsurface cracks and to the study of propagating cracks. The procedure for detecting subsurface cracks combines the finite element method with a sequential quadratic programming algorithm to solve for the unknown geometric parameters associated with the internal flaw. The procedure utilizes finite element substructuring capabilities in order to minimize the processing and solution time for practical problems. The finite element method and non‐linear optimization are also used in determining the direction a crack will propagate in a heterogeneous planar domain. This procedure involves determining the direction that produces the maximum strain energy release for a given increment of crack growth. The procedure is applied to several numerical examples. The results of these numerical studies coincide with theoretical predictions and experimentally observed crack behavior.

Bayesian cubature (BC) has emerged to be one of most competitive approach for estimating the multi-dimensional integral especially when the integrand is expensive to evaluate, and alternative acquisition functions, such as the Posterior Variance Contribution (PVC) function, have been developed for adaptive experiment design of the integration points. However, those sequential design strategies also prevent BC from being implemented in a parallel scheme. Therefore, this paper aims at developing a parallelized adaptive BC method to further improve the computational efficiency.

By theoretically examining the multimodal behavior of the PVC function, it is concluded that the multiple local maxima all have important contribution to the integration accuracy as can be selected as design points, providing a practical way for parallelization of the adaptive BC. Inspired by the above finding, four multimodal optimization algorithms, including one newly developed in this work, are then introduced for finding multiple local maxima of the PVC function in one run, and further for parallel implementation of the adaptive BC.

The superiority of the parallel schemes and the performance of the four multimodal optimization algorithms are then demonstrated and compared with the k-means clustering method by using two numerical benchmarks and two engineering examples.

Multimodal behavior of acquisition function for BC is comprehensively investigated. All the local maxima of the acquisition function contribute to adaptive BC accuracy. Parallelization of adaptive BC is realized with four multimodal optimization methods.

A simple numerical algorithm is developed for the implementation of the harmonic balance method to analyse periodic responses of a general dynamic system having geometrical non‐linearities of the quadratic and cubic types. The resulting non‐linear algebraic equations which are not explicitly determined are solved by non‐linear equation routines available in most mathematical libraries. Various non‐linear responses, such as the combinational resonances of a hinged‐clamped beam, the non‐linear effect on degenerate vibration modes of a square plate and the non‐linear oscillation of thin rings, are presented to demonstrate the versatility of the algorithm.

This paper examines several modifications to Newton's method for the numerical solution of the nonlinear Poisson equation which describes the electrostatic potential distribution in a semiconductor device. Two methods for a more efficient solution of the equation when the device is a Metal‐Oxide‐Semiconductor Field Effect Transistor are proposed. Their extension to the solution of the fully coupled system of equations is also discussed. The modifications to Newton's method are also compared numerically.