Search results

The main purpose of this paper is to find convenient methods to solve the differential‐algebraic equations which have great importance in various fields of science and engineering.

The paper applies a semi‐analytical approach, using both the homotopy analysis method (HAM) and the modified homotopy analysis method (MHAM) for finding the solution of linear and nonlinear DAEs.

The results show that the new modification can effectively reduce computational costs and accelerates the rapid convergence of the series solution.

Some high index DAEs are investigated to present a comparative study between the HAM and the MHAM.

The purpose of this paper is to present the methodology to the robust stability analysis of a vision‐based control loop in an uncalibrated environment. The type of uncertainties considered is the parametric uncertainties. The approach adopted in this paper utilizes quadratic Lyapunov function to determine the composite Jacobian matrix and ensures the robust stability using linear matrix inequality (LMI) optimization. The effectiveness of the proposed approach can be witnessed by applying it to two‐link robotic manipulator with the camera mounted on the end‐effector.

The objective of this research is the analysis of uncertain nonlinear system by representing it in differential‐algebraic form. By invoking the suitable system representation and Lyapunov analysis, the stability conditions are described in terms of linear matrix inequalities.

The proposed method is proved robust in the presence of parametric uncertainties.

Through a differential‐algebraic equation, LMI conditions are devised that ensure the stability of the uncertain system while providing an estimate of the domain of attraction based upon quadratic Lyapunov function.

The purpose of this paper is to propose the generalized integral transform technique (GITT) to the solution of convection-diffusion problems with nonlinear boundary conditions by employing the corresponding nonlinear eigenvalue problem in the construction of the expansion basis.

The original nonlinear boundary condition coefficients in the problem formulation are all incorporated into the adopted eigenvalue problem, which may be itself integral transformed through a representative linear auxiliary problem, yielding a nonlinear algebraic eigenvalue problem for the associated eigenvalues and eigenvectors, to be solved along with the transformed ordinary differential system. The nonlinear eigenvalues computation may also be accomplished by rewriting the corresponding transcendental equation as an ordinary differential system for the eigenvalues, which is then simultaneously solved with the transformed potentials.

An application on one-dimensional transient diffusion with nonlinear boundary condition coefficients is selected for illustrating some important computational aspects and the convergence behavior of the proposed eigenfunction expansions. For comparison purposes, an alternative solution with a linear eigenvalue problem basis is also presented and implemented.

This novel approach can be further extended to various classes of nonlinear convection-diffusion problems, either already solved by the GITT with a linear coefficients basis, or new challenging applications with more involved nonlinearities.

The purpose of this study is to further advance the multiple space/time subdomain framework with model reduction. Existing linear multistep (LMS) methods that are second-order time accurate, and useful for practical applications, have a significant limitation. They do not account for separable controllable numerical dissipation of the primary variables. Furthermore, they have little or no significant choices of altogether different algorithms that can be integrated in a single analysis to mitigate numerical oscillations that may occur. In lieu of such limitations, under the generalized single-step single-solve (GS4) umbrella, several of the deficiencies are circumvented.

The GS4 framework encompasses a wide variety of LMS schemes that are all second-order time accurate and offers controllable numerical dissipation. Unlike existing state-of-art, the present framework permits implicit–implicit and implicit–explicit coupling of algorithms via differential algebraic equations (DAE). As further advancement, this study embeds proper orthogonal decomposition (POD) to further reduce model sizes. This study also uses an iterative convergence check in acquiring sufficient snapshot data to adequately capture the physics to prescribed accuracy requirements. Simple linear/nonlinear transient numerical examples are presented to provide proof of concept.

The present DAE-GS4-POD framework has the flexibility of using different spatial methods and different time integration algorithms in altogether different subdomains in conjunction with the POD to advance and improve the computational efficiency.

The novelty of this paper is the addition of reduced order modeling features, how it applies to the previous DAE-GS4 framework and the improvement of the computational efficiency. The proposed framework/tool kit provides all the needed flexibility, robustness and adaptability for engineering computations.

The purpose of this paper is to describe a novel implementation of a multispatial method, multitime-scheme subdomain differential algebraic equation (DAE) framework allowing a mix of different space discretization methods and different time schemes by a robust generalized single step single solve (GS4) family of linear multistep (LMS) algorithms on a single body analysis for the first-order nonlinear transient systems.

This proposed method allows the coupling of different numerical methods, such as the finite element method and particle methods, and different implicit and/or explicit algorithms in each subdomain into a single analysis with the GS4 framework. The DAE, which constrains both space and time in multi-subdomain analysis, combined with the GS4 framework ensures the second-order time accuracy in all primary variables and Lagrange multiplier. With the appropriate GS4 parameters, the algorithmic temperature rate variable shift can be matched for all time steps using the DAE. The proposed method is used to solve various combinations of spatial methods and time schemes between subdomains in a single analysis of nonlinear first-order system problems.

The proposed method is capable of coupling different spatial methods for multiple subdomains and different implicit/explicit time integration schemes in the GS4 framework while sustaining second-order time accuracy.

Traditional approaches do not permit such robust and flexible coupling features. The proposed framework encompasses most of the LMS methods that are second-order time accurate and unconditionally stable.

Simple methods for the steady‐state analysis of a flow network are readily available, but the dynamic behavior of a large‐scale flow network is difficult to study due to the complex differential‐algebraic equation system resulting from its modeling. It is the aim of this paper to present two simple methods for the dynamic analysis of large‐scale flow networks and to demonstrate their use by examining the dynamics of a self‐similar branching tree network.

Two numerical projection methods are proposed for one‐dimensional dynamic analysis of large piping networks. Both are extensions of that suggested by Chorin for the nonlinear differential‐algebraic system resulting from the Navier‐Stokes equations. Each numerical algorithm is discussed and verified for turbulent flow in a nonlinear, self‐similar, branching tree network with constant friction factor for which an exact solution is available.

The dynamics of this network are calculated for more realistic friction factors and described as system parameters are varied. Self‐excited oscillations due to laminar‐turbulent transition are found for some parameter values and dynamic component behavior is observed in the network which is not observable in components apart from it.

It is shown that the dynamics of a flow network can exhibit unexpected behavior, reinforcing the need for simple methods to perform dynamic analysis.

This paper presents two numerical projection schemes for dynamic analysis of large‐scale flow networks to aid in their study and design.

This paper discusses algorithms for large displacement analysis of interconnected flexible and rigid multibody systems. Hydrostatic and hydrodynamic loads for systems being submerged in water are also considered. The systems may consist of cables and beams and may combine very flexible parts with rigid parts. Various ways of introducing structural joints are discussed. A special implementation of the Hilber‐Hughes‐Taylor time integration scheme for constrained non‐linear systems is outlined. The formulation is general and allows for displacements and rotational motion of unlimited size. Aspects concerning efficient solution of constrained dynamic problems are discussed. These capabilities have been implemented in a general purpose non‐linear finite element program. Applications involving static and dynamic analysis of a bi‐articulated tower and a floating tripod platform kept in place by three anchor lines are discussed.

The space discretization of eddy‐current problems in the magnetic vector potential formulation leads to a system of differential‐algebraic equations. They are typically time discretized by an implicit method. This requires the solution of large linear systems in the Newton iterations. The authors seek to speed up this procedure. In most relevant applications, several materials are non‐conducting and behave linearly, e.g. air and insulation materials. The corresponding matrix system parts remain constant but are repeatedly solved during Newton iterations and time‐stepping routines. The paper aims to exploit invariant matrix parts to accelerate the system solution.

Following the principle “reduce, reuse, recycle”, the paper proposes a Schur complement method to precompute a factorization of the linear parts. In 3D models this decomposition requires a regularization in non‐conductive regions. Therefore, the grad‐div regularization is revisited and tailored such that it takes anisotropies into account.

The reduced problem exhibits a decreased effective condition number. Thus, fewer preconditioned conjugate gradient iterations are necessary. Numerical examples show a decrease of the overall simulation time, if the step size is small enough. 3D simulations with large time step sizes might not benefit from this approach, because the better condition does not compensate for the computational costs of the direct solvers used for the Schur complement. The combination of the Schur approach with other more sophisticated preconditioners or multigrid solvers is subject to current research.

The Schur complement method is adapted for the eddy‐current problem. Therefore, a new partitioning approach into linear/non‐linear and static/dynamic domains is proposed. Furthermore, a new variant of the grad‐div gauging is introduced that allows for anisotropies and enables the Schur complement method in 3D.

The purpose of this paper is to show that the computation of time-periodic signals for coupled antenna-circuit problems can be substantially accelerated by means of the single shooting method. This allows an efficient analysis of nonlinearly loaded coupled loop antennas for near field communication (NFC) applications.

For the modelling of electrically small coupled field-circuit problems, the partial element equivalent circuit (PEEC) method shows to be very efficient. For analysing the circuit-like description of the coupled problem, this paper developed a generalised modified nodal analysis (MNA) and applied it to specific NFC problems.

It is shown that the periodic steady state (PSS) solution of the resulting differential-algebraic system can be computed very time efficiently by the single shooting method. A speedup of roughly 114 to conventional transient approaches can be achieved.

The proposed approach appears to be an efficient alternative for the computation of time PSS solutions for nonlinear circuit problems coupled with discretised conductive structures, where the homogeneous solution is not of interest.

The present paper explores the implementation and application of the shooting method for nonlinearly loaded coupled antenna-circuit problems based on the PEEC method and shows the efficiency of this approach.