Optimal combination of spatial basis functions for the model reduction of nonlinear distributed parameter systems

Spectral methods are among the most extensively used techniques for model reduction of distributed parameter systems in various fields, including fluid dynamics, quantum mechanics, heat conduction, and weather prediction. However, the model dimension is not minimized for a given desired accuracy because of general spatial basis functions. New spatial basis functions are obtained by linear combination of general spatial basis functions in spectral method, whereas the basis function transformation matrix is derived from straightforward optimization techniques. After the expansion and truncation of spatial basis functions, the present spatial basis functions can provide a lower dimensional and more precise ordinary differential equation system to approximate the dynamics of the systems. The numerical example shows the feasibility and effectiveness of the optimal combination of spectral basis functions for model reduction of nonlinear distributed parameter systems.     

二阶耗散动力系统的降维对解长期行为的误差估计

基于非线性动力学理论,对一类高维二阶耗散自治动力系统的降维及其对解的长期行为的影响进行了理论分析．该分析将方程的解投影到控制方程的线性算子的特征向量所张成的完备空间中,并在相空间中引入一距离的概念,方便地解决了缩减后系统与原始系统解之间的误差或距离的描述．基于此距离定义,首先,分析了由于高阶模态的截取对解的长期行为的影响,并推导出了相应的误差估计,该估计表明由于降维对系统长期行为的影响不仅与系统的高阶子空间中的固有频率和阻尼比乘积的最小值有关,并且与高阶子空间中的某一最大固有频率有关．然后,将一般的模态截取视为对原系统的解的一个扰动,对一些文献中由于降维程度的不同而造成解的拓扑性质发生变化的现象进行了定性的解释．     

Locomotion of Articulated Bodies in a Perfect Fluid

This paper is concerned with modeling the dynamics of N articulated solid bodies submerged in an ideal fluid. The model is used to analyze the locomotion of aquatic animals due to the coupling between their shape changes and the fluid dynamics in their environment. The equations of motion are obtained by making use of a two-stage reduction process which leads to significant mathematical and computational simplifications. The first reduction exploits particle relabeling symmetry: that is, the symmetry associated with the conservation of circulation for ideal, incompressible fluids. As a result, the equations of motion for the submerged solid bodies can be formulated without explicitly incorporating the fluid variables. This reduction by the fluid variables is a key difference with earlier methods, and it is appropriate since one is mainly interested in the location of the bodies, not the fluid particles. The second reduction is associated with the invariance of the dynamics under superimposed rigid motions. This invariance corresponds to the conservation of total momentum of the solid-fluid system. Due to this symmetry, the net locomotion of the solid system is realized as the sum of geometric and dynamic phases over the shape space consisting of allowable relative motions, or deformations, of the solids. In particular, reconstruction equations that govern the net locomotion at zero momentum, that is, the geometric phases, are obtained. As an illustrative example, a planar three-link mechanism is shown to propel and steer itself at zero momentum by periodically changing its shape. Two solutions are presented: one corresponds to a hydrodynamically decoupled mechanism and one is based on accurately computing the added inertias using a boundary element method. The hydrodynamically decoupled model produces smaller net motion than the more accurate model, indicating that it is important to consider the hydrodynamic interaction of the links.     

The effects of body composition on setpoint based weight loss

The setpoint model of human metabolic adaption to diet-induced weight loss has been introduced and examined considering a constant energy density in the body weight lost. Body composition analysis indicates that the ratio of fat and nonfat lost is constant for an individual but is dependent on the initial percent body fat, producing an energy density varying from individual to individual. In this study, body composition and the setpoint model are used to examine weight loss dynamics. Comparison is made to a Harris-Benedict (HB) based model using a nominal energy density (3500 calories/lb). Conclusions: the weight dynamics for a given reduction in calorie intake are completely determined by the initial percentage of body fat. The HB model is predicted to accurately track the weight loss for the very obese who experience little metabolic adaption and whose energy density is close to the nominal value. However, for subjects with decreased initial body fat percentage, the HB model first underestimates then overestimates the weight loss. The crossover time and the maximum overestimate is dependent on the initial body fat percent and the percent dietary calorie reduction. A characteristic time and an energy density ratio are defined and used to calculate a maximum rate of weight loss for a given calorie reduction.     

Dynamics of a ratio‐dependent stage‐structured predator‐prey model with delay

In this paper, we investigate the dynamics of a time‐delay ratio‐dependent predator‐prey model with stage structure for the predator. This predator‐prey system conforms to the realistically biological environment. The existence and stability of the positive equilibrium are thoroughly analyzed, and the sufficient and necessary conditions for the stability and instability of the positive equilibrium are obtained for the case without delay. Then, the influence of delay on the dynamics of the system is investigated using the geometric criterion developed by Beretta and Kuang. 26 We show that the positive steady state can be destabilized through a Hopf bifurcation and there exist stability switches under some conditions. The formulas determining the direction and the stability of Hopf bifurcations are explicitly derived by using the center manifold reduction and normal form theory. Finally, some numerical simulations are performed to illustrate and expand our theoretical results.     

Stabilization of nonlinear systems in compound critical cases

In this paper, we study the stabilization of nonlinear systems in critical cases by using the center manifold reduction technique. Three degenerate cases are considered, wherein the linearized model of the system has two zero eigenvalues, one zero eigenvalue and a pair of nonzero pure imaginary eigenvalues, or two distinct pairs of nonzero pure imaginary eigenvalues; while the remaining eigenvalues are stable. Using a local nonlinear mapping (normal form reduction) and Liapunov stability criteria, one can obtain the stability conditions for the degenerate reduced models in terms of the original system dynamics. The stabilizing control laws, in linear and/or nonlinear feedback forms, are then designed for both linearly controllable and linearly uncontrollable cases. The normal form transformations obtained in this paper have been verified by using code MACSYMA.     

Reduced-order model for laminar vortex-induced vibration of a rigid circular cylinder with an internal nonlinear absorber

The nonlinear interaction of a laminar flow and a sprung rigid circular cylinder results in vortex-induced vibration (VIV) of the cylinder. Passive suppression of the VIV by attaching an internal nonlinear vibration absorber that acts, in essence, as a nonlinear energy sink (NES) to the cylinder has been observed in finite-element computations involving thousands of degrees of freedom (DOF). A single-DOF self-excited oscillator is developed to approximate the limit-cycle oscillation (LCO) of the cylinder undergoing VIV. This self-excited oscillator models the interaction of the flow and the cylinder. Then, a two-DOF reduced-order model for the system with the internal NES is constructed by coupling the single-DOF NES to the single-DOF self-excited oscillator. Hence, the complicated high-dimensional system of flow-cylinder-NES involving thousands of DOF is reduced to a two-DOF model. The two targeted energy transfer mechanisms responsible for passive VIV suppression that are observed in the finite-element computations are fully reproduced using the two-DOF reduced-order model. This reduction of the dynamics to an easily tractable low-dimensional reduced-order model facilitates the approximate analysis of the underlying dynamics. Moreover, the underlying assumptions of the order reduction, and the parameter ranges of validity of the reduced-order model are formulated and systematically studied.     

Modelling of Heat Transfer and Vehicle Dynamics for Thermal Load Reduction by Hypersonic Flight Optimization

Heat input reduction by appropriate, optimal trajectory control is considered for the range cruise and the return-to-base cruise of a hypersonic vehicle propelled by a turbo/ram jet engines combination. A mathematical model is developed for describing the unsteady heat transfer through the thermal protection system. This model is coupled to the model of the dynamics of the vehicle. An efficient optimization technique is applied for constructing a solution for the two cruise problems. The results show that significant heat input reductions can be achieved with only a small penalty in fuel consumption.     

A MODEL OF THE ENZOOTIOLOGY OF LYME DISEASE IN THE ATLANTIC NORTHEAST OF THE UNITED STATES

A mathematical model is presented for the dynamics of the rate of infection of the Lyme disease vector tick Ixodes dammini (Acari: Ixodidae) by the spirochete Borrelia burgdorferi, in the Atlantic Northeast of the United States. According to this model, moderate reductions in the abundance of white-tailed deer Odocoileus virginianus may either decrease or increase the spirochete infection rate in ticks, provided the deer are not reservoir hosts for Lyme disease. Expressions for the basic reproductive rate of the disease are computed analytically for special cases, and it is shown that as the basic reproductive rate increases, a proportional reduction in the tick population produces a smaller proportional reduction in the infection rate, so that vector control is less effective far above the threshold. The model also shows that control of the mouse reservoir hosts Peromyscus leucopus could reduce the infection rate if the survivorship of juvenile stages of ticks were reduced as a consequence. If the survivorship of juvenile stages does not decline as the rodent population is reduced, then rodent reduction can increase the spirochete infection rate in the ticks.     

Varying Workpiece Dynamics in Milling Stability Analysis

Dynamic stability of a milling process with varying workpiece dynamics is investigated. The milling tool moves along the workpiece with a prescribed feed rate, whereby the contact point shifts. Furthermore, the workpiece dynamics is affected by material removal. The resultant varying workpiece dynamics is taken into account by parametric model order reduction including modal truncation. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analog to digital conversion in physical measurements

There exist measuring devices where an analog input is converted into a digital output. Such converters can have a nonlinear internal dynamics. We show how measurements with such converting devices can be understood using concepts from symbolic dynamics. Our approach is based on a nonlinear one-to-one mapping between the analog input and the digital output of the device. We analyze the Bernoulli shift and the tent map which are realized in specific analog/digital (A/D) converters. Furthermore, we discuss the sources of errors that are inevitable in physical realizations of such systems and suggest methods for error reduction.     

Model order reduction of finite element models: improved component mode synthesis

The finite element (FE) approach constitutes an essential methodology when modelling the elastic properties of structures in various research disciplines such as structural mechanics, engine dynamics and so on. Because of increased accuracy requirements, the FE method results in discretized models, which are described by higher order ordinary differential equations, or, in FE terms, by a large number of degrees of freedom (DoF). In this regard, the application of an additional methodology, referred to as the model order reduction (MOR) or DoF condensation, is rather compulsory. Herein, a reduced dimension set of ordinary differential equations is generated, i.e. the initially large number of DoF is condensed, while aiming to keep the dynamics of the original model as intact as possible. In the commercially available FE software tools, the static and the component mode syntheses (CMS) are the only available integrated condensation methods. The latter represents the state of the art generating well-correlated reduced order models (ROMs), which can be further utilized for FE or multi-body systems simulations. Taking into consideration the information loss of the CMS, which is introduced by its part-static nature, the improved CMS (ICMS) method is proposed. Here the algorithmic scheme of the standard CMS is adopted, which is qualitatively improved by adequately considering the advantageous characteristics of another MOR approach, the so-called improved reduction system method. The ICMS results in better correlated reduced order models in comparison to all the aforementioned methods, while preserving the required structural properties of the original FE model.     

Approximate reduction of non-linear discrete models with two time scales

The aim of this work is to present a general class of nonlinear discrete time models with two time scales whose dynamics is susceptible of being approached by means of a reduced system. The reduction process is included in the so-called approximate aggregation of variables methods which consist of describing the dynamics of a complex system involving many coupled variables through the dynamics of a reduced system formulated in terms of a few global variables. For the time unit of the discrete system we use that of the slow dynamics and assume that fast dynamics acts a large number of times during it. After introducing a general two-time scales nonlinear discrete model we present its reduced accompanying model and the relationships between them. The main result proves that certain asymptotic behaviours, hyperbolic asymptotically stable (A.S.) periodic solutions, to the aggregated system entail that to the original system.     

Analysis of a Class of Symmetric Equilibrium Configurations for a Territorial Model

<正>Motivated by an animal territoriality model,we consider a centroidal Voronoi tessellation algorithm from a dynamical systems perspective.In doing so,we discuss the stability of an aligned equilibrium configuration for a rectangular domain that exhibits interesting symmetry properties.We also demonstrate the procedure for performing a center manifold reduction on the system to extract a set of coordinates which capture the long term dynamics when the system is close to a bifurcation.Bifurcations of the system restricted to the center manifold are then classified and compared to numerical results.Although we analyze a specific set-up,these methods can in principle be applied to any bifurcation point of any equilibrium for any domain.     

On a special type of stability of differential equations for induction machines with double squirrel-cage rotor

This paper puts forward a new mathematical model describing the dynamics of an induction machine with double squirrel cage rotor. The limit load problem for induction machines is considered. The nonlocal reduction method is applied to estimate the limit load. Estimates for the stability domain of the system of differential equations for inductions machines in question are obtained.     

Aggregation and emergence in systems of ordinary differential equations

The aim of this article is to present aggregation methods for a system of ordinary differential equations (ODE's) involving two time scales. The system of ODE's is composed of the sum of fast parts and a perturbation. The fast dynamics are assumed to be conservative. The corresponding first integrals define a few global variables. Aggregation corresponds to the reduction of the dimension of the dynamical system which is replaced by an aggregated system governing the global variables at the slow time scale. The centre manifold theorem is used in order to get the slow reduced model as a Taylor expansion of a small parameter. We particularly look for the conditions necessary to get emerging properties in the aggregated model with respect to the nonaggregated one. We define two different types of emergences, functional and dynamical. Functional emergence corresponds to different functions for the two dynamics, aggregated and nonaggregated. Dynamical emergence means that both dynamics are qualitatively different. We also present averaging methods for aggregation when the fast system converges towards a stable limit cycle.     

Model reduction on inertial manifolds for N–S equations approached by multilevel finite element method

Approximate Inertial Manifolds (AIMs) is approached by multilevel finite element method, which can be referred to as a Post-processed nonlinear Galerkin finite element method, and is applied to the model reduction for fluid dynamics, a typical kind of nonlinear continuous dynamic system from viewpoint of nonlinear dynamics. By this method, each unknown variable, namely, velocity and pressure, is divided into two components, that is the large eddy and small eddy components. The interaction between large eddy and small eddy components, which is negligible if standard Galerkin algorithm is used to approach the original governing equations, is considered essentially by AIMs, and consequently a coarse grid finite element space and a fine grid incremental finite element space are introduced to approach the two components. As an example, the flow field of incompressible flows around airfoil is simulated numerically and discussed, and velocity and pressure distributions of the flow field are obtained accurately. The results show that there exists less essential degrees-of-freedom which can dominate the dynamic behaviors of the discretized system in comparison with the traditional methods, and large computing time can be saved by this efficient method. In a sense, the small eddy component can be captured by AIMs with fewer grids, and an accurate result can also be obtained.     

Low-dimensional I/O modeling of Distributed Parameter Systems

Many model reduction techniques take a semi-discretization of the original PDE model as starting point and aim then at an accurate approximation of its input/output map. In this contribution, we discuss the direct discretization of the i/o map of the original infinite-dimensional system. First, the input and output signals are discretized in space and time, second, the system dynamics are approximated in form of the underlying evolution operator, leading to an approximated i/o map with matrix representation. The discretization framework, corresponding error estimations, a SVD-based system reduction method and a numerical application in an optimization problem are presented for the example of a linear time-invariant heat equation. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

基于卷积神经网络气动力降阶模型的翼型优化方法

针对非线性大扰动翼型气动力优化问题,提出了基于卷积神经网络气动力降阶模型的优化方法.该方法用不同形状参数下翼型的气动力数据作为训练信号,训练卷积神经网络翼型气动力降阶模型.采用该气动力降阶模型,以最大升阻比为目标,对翼型进行优化,结果表明该方法可用于大扰动下翼型气动力的预测和优化.该文同时还讨论了池化法和径向基法的训练...