Output feedback pole placement problem is not solvable analytically for the plants having either number of inputs or number of outputs, more than two. There is even lack of a good numerical method which solves any gen...
详细信息
Output feedback pole placement problem is not solvable analytically for the plants having either number of inputs or number of outputs, more than two. There is even lack of a good numerical method which solves any general pole assignment problem. Due to the multi-linear nature of the pole placement problem there is a possibility to utilize multi-linear structure and arrive at a better numerical solution. This paper shows that it is possible to compute analytically the Jacobian and Hessian matrix in an easy manner and utilize them in an iterative numerical method to solve the pole-placement problem. Newton-Raphson method is used with analytical solution of Jacobian matrix to give an iterative solution of pole placement equations with better percentage success rate than the other methods quoted in literature.
A new algorithm to test percolation conditions for the solution of percolation problems on a lattice and continuum percolation for spaces of an arbitrary dimension has been proposed within the Newman-Ziff algorithm. T...
详细信息
A new algorithm to test percolation conditions for the solution of percolation problems on a lattice and continuum percolation for spaces of an arbitrary dimension has been proposed within the Newman-Ziff algorithm. The algorithm is based on the use of bitwise operators and does not reduce the efficiency of the operation of the Newman-Ziff algorithm as a whole. This algorithm makes it possible to verify the existence of both clusters touching boundaries at an arbitrary point and single-loop clusters continuously connecting the opposite boundaries in a percolating system with periodic boundary conditions. The existence of a cluster touching the boundaries of the system at an arbitrary point for each direction, the formation of a one-loop cluster, and the formation of a cluster with an arbitrary number of loops on a torus can be identified in one calculation by combining the proposed algorithm with the known approaches for the identification of the existence of a percolation cluster. The operation time of the proposed algorithm is linear in the number of objects in the system.
Fundamental properties of the angular-spatial symmetry of radiation fields in the uniform slab of a finite optical thickness are used for improvement of the numerical methods and algorithms of the classical radiative ...
详细信息
ISBN:
(纸本)9783642396403
Fundamental properties of the angular-spatial symmetry of radiation fields in the uniform slab of a finite optical thickness are used for improvement of the numerical methods and algorithms of the classical radiative transfer theory. A new notion of so called photometrical invariants is introduced. The basic boundary-value problem of the radiative transfer theory is reformulated in new terms for the subsequent simplification of algorithms of numerical modeling methods such as spherical harmonics, discrete ordinates, Gauss-Seidel and Case methods. This simplification leads to two-fold decrease of the ranks of linear algebraic equations with simultaneous reduction of numerical modeling intervals connected with angular and spatial variables.
Toughness in Ceramic Matrix Composites (CMCs) is achieved if crack deflection can occur at the fiber/matrix interface, preventing crack penetration into the fiber and enabling energy-dissipating fiber pullout. To inve...
详细信息
Toughness in Ceramic Matrix Composites (CMCs) is achieved if crack deflection can occur at the fiber/matrix interface, preventing crack penetration into the fiber and enabling energy-dissipating fiber pullout. To investigate toughening in nanoscale CMCs, direct atomistic models are used to study how matrix cracks behave as a function of the degree of interfacial bonding/sliding, as controlled by the density of C interstitial atoms, at the interface between carbon nanotubes (CNTs) and a diamond matrix. Under all interface conditions studied, incident matrix cracks do not penetrate into the nanotube. Under increased loading, weaker interfaces fail in shear while stronger interfaces do not fail and, instead, the CNT fails once the stress on the CNT reaches its tensile strength. An analytic shear lag model captures all of the micromechanical details as a function of loading and material parameters. Interface deflection versus fiber penetration is found to depend on the relative bond strengths of the interface and the CNT, with CNT failure occurring well below the prediction of the toughness-based continuum He-Hutchinson model. The shear lag model, in contrast, predicts the CNT failure point and shows that the nanoscale embrittlement transition occurs at an interface shear strength scaling as tau similar to epsilon(f),CNT sigma CNT, rather than tau s similar to sigma CNT typically prevailing for micron scale composites, where epsilon f,CNT and sigma CNT are the CNT failure strain and stress, respectively. Interface bonding also lowers the effective fracture strength in SWCNTs, due to formation of defects, but does not play a role in DWCNTs having interwall coupling, which are weaker than SWCNTs but less prone to damage in the outerwall. (C) 2013 Elsevier Ltd. All rights reserved.
Minimum distance to a solid wall is a commonly used parameter in turbulence closure formulations associated with the Reynolds Averaged form of the Navier Stokes Equations (RANS). This paper presents a new approach to ...
详细信息
Minimum distance to a solid wall is a commonly used parameter in turbulence closure formulations associated with the Reynolds Averaged form of the Navier Stokes Equations (RANS). This paper presents a new approach to efficiently compute the minimum distance between a set of points and a surface. The method is based on sphere voxelization, and uses fast integer arithmetic algorithms from the field of computer graphics. Using a simple test case where the number of points (N-p) and surface elements (N-b) can be independently specified, the present method is empirically estimated to be O((NpNb0.5)-N-0.8). An unstructured grid around an aircraft configuration (DLR-F6) is chosen as the test case for demonstration and validation. Multi-processor computations (up to 256 processors) are conducted to study efficiency and scalability. Encouraging results are obtained, with the sphere voxelization algorithm demonstrated to be more efficient than all of the alternate methods for computing minimum distances. However, a load imbalance does exist, which negatively impacts the scalability for large number of cores. A simple method for load re-balancing is formulated and tested, which results in significant improvements in both efficiency and scalability. (C) 2013 Elsevier Inc. All rights reserved.
A small-displacement, phenomenological 3D model for shape memory alloys is expanded for large deformations. The Eulerian rate type formulation from the large deformation theory is used, with an additive decomposition ...
详细信息
A small-displacement, phenomenological 3D model for shape memory alloys is expanded for large deformations. The Eulerian rate type formulation from the large deformation theory is used, with an additive decomposition of the stretching tensor. The expansion is performed by introducing, in the first stage, a thermodynamic frame for an isotropic material with logarithmic strains and scalar-valued internal state variables. The rate type characterisation of the elasticity is derived from the free energy. The relation is exactly integrable to define a free energy-based elastic relation. The formulation presented extends the previously presented ideas to possibly characterise changing elastic behaviours as a result of a phase transformation or other inelastic mechanisms. In the second part, a thermodynamic framework is used to derive the SMA constitutive equations by introducing large deformation counterparts in the small displacement theory for the Gibbs free energy, the transformation flow rule and the criteria for transformation. Then, a compatible integration algorithm is derived and implemented in Abaqus Standard FE software. The simulations presented demonstrate both a path-independent description of the elasticity and the model applicability to describe the unique properties of shape memory alloys. (C) 2012 Elsevier Ltd. All rights reserved.
Accelerating numerical algorithms for solving sparse linear systems on parallel architectures has attracted the attention of many researchers due to their applicability to many engineering and scientific problems. The...
详细信息
Accelerating numerical algorithms for solving sparse linear systems on parallel architectures has attracted the attention of many researchers due to their applicability to many engineering and scientific problems. The solution of sparse systems often dominates the overall execution time of such problems and is mainly solved by iterative methods. Preconditioners are used to accelerate the convergence rate of these solvers and reduce the total execution time. Sparse approximate inverse (SAI) preconditioners are a popular class of preconditioners designed to improve the condition number of large sparse matrices. We propose a GPU accelerated SAI preconditioning technique called GSAI, which parallelizes the computation of this preconditioner on NVIDIA graphic cards. The preconditioner is then used to enhance the convergence rate of the BiConjugate Gradient Stabilized (BiCGStab) iterative solver on the GPU. The SAI preconditioner is generated on average 28 and 23 times faster on the NVIDIA GTX480 and TESLA M2070 graphic cards, respectively, compared to ParaSails (a popular implementation of SAI preconditioners on CPU) single processor/core results. The proposed GSAI technique computes the SAI preconditioner in approximately the same time as ParaSails generates the same preconditioner on 16 AMD Opteron 252 processors.
This paper deals with the acoustic velocity field simulation generated by interaction of flow with moving bodies. Starting from the Ffowcs Williams and Hawkings equation, an analytical formulation of the acoustic velo...
详细信息
This paper deals with the acoustic velocity field simulation generated by interaction of flow with moving bodies. Starting from the Ffowcs Williams and Hawkings equation, an analytical formulation of the acoustic velocity is derived for sources in arbitrary motion. This makes the imposition of the boundary condition on a (rigid) scattering surface much more straightforward, as, if the traditional pressure formulation is used, then the pressure gradient must be calculated. Computational results for a pulsating sphere, dipole source, and a propeller case with subsonic tips verify this formulation.
A key component in enabling the application of model predictive control (MPC) in fields such as automotive, aerospace, and factory automation is the availability of low-complexity fast optimisation algorithms to solve...
详细信息
A key component in enabling the application of model predictive control (MPC) in fields such as automotive, aerospace, and factory automation is the availability of low-complexity fast optimisation algorithms to solve the MPC finite horizon optimal control problem in architectures with reduced computational capabilities. In this paper, we introduce a projection-free iterative optimisation algorithm and discuss its application to linear MPC. The algorithm, originally developed by Brand for non-negative quadratic programs, is based on a multiplicative update rule and it is shown to converge to a fixed point which is the optimum. An acceleration technique based on a projection-free line search is also introduced, to speed-up the convergence to the optimum. The algorithm is applied to MPC through the dual of the quadratic program (QP) formulated from the MPC finite time optimal control problem. We discuss how termination conditions with guaranteed degree of suboptimality can be enforced, and how the algorithm performance can be optimised by pre-computing the matrices in a parametric form. We show computational results of the algorithm in three common case studies and we compare such results with the results obtained by other available free and commercial QP solvers.
We present a new numerical method (based on the computation of deflating subspaces) for the gamma-iteration in H-infinity, control in the extended matrix pencil formulation. We introduce a permuted graph representatio...
详细信息
We present a new numerical method (based on the computation of deflating subspaces) for the gamma-iteration in H-infinity, control in the extended matrix pencil formulation. We introduce a permuted graph representation of these subspaces, which avoids the known difficulties that arise when the iteration is based on the solution of algebraic Riccati equations but at the same time makes use of the special symmetry structures that are present in the problems. We use this representation to perform both the deflation of spurious infinity eigenvalues of the even pencils and the implementation of the inverse-free sign iteration. We show that the new method returns accurate results and is applicable in many situations where conventional methods fail. (C) 2013 Elsevier Ltd. All rights reserved.
暂无评论