A problem very often arising in applications is presented: finding the minimal root of an equation with the objective function being multiextremal and nondifferentiable. Applications from the field of electronic measu...
详细信息
A problem very often arising in applications is presented: finding the minimal root of an equation with the objective function being multiextremal and nondifferentiable. Applications from the field of electronic measurements are given. Three methods based on global optimization ideas are introduced for solving this problem. The first one uses an a priori estimate of the global Lipschitz constant. The second method adaptively estimates the global Lipschitz constant. The third algorithm adaptively estimates local Lipschitz constants during the search. All the methods either find the minimal root or determine the global minimizers (in the case when the equation under consideration has no roots). Sufficient convergence conditions of the new methods to the desired solution are established. numerical results including wide experiments with test functions, stability study, and a real-life applied problem are also presented.
A detailed survey is presented of the literature on attitude representation dating from the early work of Euler and Hamilton to recent publications in fields such as navigation and control. The scope is limited to the...
详细信息
A detailed survey is presented of the literature on attitude representation dating from the early work of Euler and Hamilton to recent publications in fields such as navigation and control. The scope is limited to the development of the aircraft kinematic transformation equations in terms of four different attitude representations, including the well-known Enter angles, the Euler-axis rotation parameters, the direction cosines, and the Euler-Rodrigues quaternion. The emphasis is directed at the application of the quaternion formulation to aircraft kinematics. Results are presented that reinforce observations that the quaternion formulation, typically implemented to eliminate singularities associated with the Enter angle formulation, is far superior to the other commonly used formulations based on computational efficiency alone. A development of quaternion constraints necessary to independently constrain roll, pitch, yaw, bank angle, elevation angle, and/or azimuth angle is presented. For verification of simulation codes, a general closed-form solution to the quaternion formulation, for the case of constant rotation, is also presented. Additionally, a discussion is provided of numerical integration methods and numerical errors for the quaternion formulation. This discussion is especially important for simulations that may still utilize a common error reduction scheme originally developed for analog computers.
The discrete-time optimal projection equations for systems with white parameters are strengthened. For the class of minimal ms (mean square) stabilizing compensators the strengthened discrete-time optimal projection e...
详细信息
ISBN:
(纸本)085296708X
The discrete-time optimal projection equations for systems with white parameters are strengthened. For the class of minimal ms (mean square) stabilizing compensators the strengthened discrete-time optimal projection equations are proved to be equivalent to first-order necessary optimality conditions for optimal reduced-order dynamic compensation of systems with white parameters. The conventional discrete-time optimal projection equations are proved to be weaker. As a result solutions of the conventional discrete-time optimal projection equations may not correspond to optimal reduced-order compensators. To compute optimal reduced-order compensators two numerical algorithms are proposed. One is a homotopy algorithm and one is based on iteration of the strengthened discrete-time optimal projection equations. The latter algorithm is a generalization of the algorithm that solves the full-order problem, which in turn is a generalization of the algorithm that solves the two Riccati equations of full-order LQG control through iteration. Therefore the efficiency of these three types of algorithms is comparable. It is demonstrated that, despite the strengthening of the optimal projection equations, the optimal reduced-order compensation problem, in general, may posses multiple extrema.
Diffraction problems in analytical treatments and numeric have been constant topics for a very long time and there are a lot of different analytical expressions as well as numerical algorithms. The modern approaches a...
详细信息
ISBN:
(纸本)0819430382
Diffraction problems in analytical treatments and numeric have been constant topics for a very long time and there are a lot of different analytical expressions as well as numerical algorithms. The modern approaches are concentrated in the area of particle sizing, fibber optical sizing and loss evaluation in laser physics and techniques. Some problems between the fibber optics dimensioning and particles of cylindrical shape dimensioning are chosen and based on literature;some ameliorations in numerical approaches are made. The evolved expressions can be used for different laser wavelengths, different dimensions of disturbance objects tin scattering and diffraction phenomena) as well as different indices of refractions where the particles tin a very general way) are stratified. Some existing programs for scattering are discussed and the new computer architecture is presented here.
The design of an algorithm for the numerical solution of a generalized quadratic algebraic matrix Riccati equation is presented. The approach is based on probability-1 homotopy methods. The algorithm is illustrated wi...
详细信息
The design of an algorithm for the numerical solution of a generalized quadratic algebraic matrix Riccati equation is presented. The approach is based on probability-1 homotopy methods. The algorithm is illustrated with numerical examples.
Some numerical solutions of acoustic propagation problems using linearized Euler equations are studied. The two-dimensional Euler equations are linearized around a known stationary mean flow. The computed solution is ...
详细信息
Some numerical solutions of acoustic propagation problems using linearized Euler equations are studied. The two-dimensional Euler equations are linearized around a known stationary mean flow. The computed solution is obtained by using is dispersion-relation-preserving scheme in space, combined with a fourth-order Runge-Kutta algorithm in time. This numerical integration leads to very good results in terms of accuracy, stability, and low storage. The implementation of source terms in these equations is studied very carefully in various configurations, inasmuch as the final goal is to improve and to validate the stochastic noise generation and radiation model. In this approach, the turbulent velocity field is modeled by a sum of random Fourier modes through a source term in the linearized Euler equations to predict the noise from subsonic flows. The radiation of a point source in a subsonic and a supersonic uniform mean flow is investigated. The numerical estimates are shown to be in excellent agreement with the analytical solutions. Then, the emphasis is on the ability of the method to describe correctly the multipolar structure of aeroacoustic sources. The radiation of dipolar and guadrupolar extended sources is, thus, studied. Next, a typical problem in jet noise is considered with the propagation of acoustic waves in a sheared mean flow. The numerical solution compares favorably with ray tracing. Finally, a nonlinear formulation of Euler's equations is solved to limit the growth of instability waves excited by the acoustic source terms.
We consider the changing face of computing and the evolving interface of numerical algorithms from subroutines to: spreadsheets, high-level scripting languages, visual programming, Web browsers. Through examples we co...
详细信息
We consider the changing face of computing and the evolving interface of numerical algorithms from subroutines to: spreadsheets, high-level scripting languages, visual programming, Web browsers. Through examples we consider the movement from the subroutine language interface to sophisticated computational environments and trace how NAG libraries and products have evolved along this path. (C) 2000 IMACS. Published by Elsevier Science B.V. All rights reserved.
作者:
Yang, KSInha Univ
Dept Engn Mech Nam Ku Inchon 402751 South Korea
Instability and transition in an obstructed channel how are investigated using direct numerical simulation. The flow geometry under consideration is a plane channel with two-dimensional thin obstacles (baffles) mounte...
详细信息
Instability and transition in an obstructed channel how are investigated using direct numerical simulation. The flow geometry under consideration is a plane channel with two-dimensional thin obstacles (baffles) mounted symmetrically in the vertical direction and periodically in the streamwise direction. The flow is steady and symmetric at low Reynolds numbers. Above a certain Reynolds number, the flow undergoes a Hopf bifurcation, leading to a solution periodic in time. At high Reynolds numbers the unsteady flow exhibits a space-time symmetry-breaking bifurcation. A secondary instability is also observed at high Reynolds numbers, which is believed to be responsible for a subsequent chaotic breakdown of the flow. To study the secondary instability, we take the periodic solution, which results from the Hopf bifurcation, as the basic flow to be disturbed with small-amplitude random noise. Depending on the Reynolds number the basic Row becomes unstable to three-dimensional disturbances, which results in a chaotic how The numerical results obtained are consistent with experimental findings currently available.
This paper. deals with the influence of martensitic transformation latent heat on the superelastic behaviour modelling. Exothermic and endothermic effects are responsible for a strong evolution on the temperature fiel...
详细信息
This paper. deals with the influence of martensitic transformation latent heat on the superelastic behaviour modelling. Exothermic and endothermic effects are responsible for a strong evolution on the temperature field inside the material that modifies its mechanical response. This is responsible for a strain late effect that is taken into account by introducing a coupling equation between the production rate of martensite and the temperature change, into a micro macro modelling of the superthermoelastic behaviour with the assumption that the temperature field remains uniform but different to the test temperature imposed. numerical simulations so obtained show a good agreement with experimental results performed on Cu-based superelastic alloys. (C) 2000 Published by Elsevier Science Ltd.
This paper deals with the simulation of the mechanical response and texture evolution of cubic crystals and polycrystals for a rate-independent elastic-plastic constitutive law. No viscous effects are considered. An a...
详细信息
This paper deals with the simulation of the mechanical response and texture evolution of cubic crystals and polycrystals for a rate-independent elastic-plastic constitutive law. No viscous effects are considered. An algorithm is introduced to treat the difficult case of multi-surface plasticity. This algorithm allows the computation of the mechanical response of a single crystal. The corresponding yield surface is made of the intersection of several hyper-planes in the stress space. The problem of the multiplicity of the slip systems is solved thanks to a pseudo-inversion method. Self and latent hardening are taken into account. In order to compute the response of a polycrystal, a Taylor homogenization scheme is used. The stress-strain response of single crystals and polycrystals is computed for various loading cases. The texture evolution predicted for compression, plane strain compression and simple shear are compared with the results given by a visco-plastic polycrystalline model. (C) 2000 Elsevier Science Ltd. All rights reserved.
暂无评论