32,506 research outputs found
Nonlinear analysis of spacecraft thermal models
We study the differential equations of lumped-parameter models of spacecraft
thermal control. Firstly, we consider a satellite model consisting of two
isothermal parts (nodes): an outer part that absorbs heat from the environment
as radiation of various types and radiates heat as a black-body, and an inner
part that just dissipates heat at a constant rate. The resulting system of two
nonlinear ordinary differential equations for the satellite's temperatures is
analyzed with various methods, which prove that the temperatures approach a
steady state if the heat input is constant, whereas they approach a limit cycle
if it varies periodically. Secondly, we generalize those methods to study a
many-node thermal model of a spacecraft: this model also has a stable steady
state under constant heat inputs that becomes a limit cycle if the inputs vary
periodically. Finally, we propose new numerical analyses of spacecraft thermal
models based on our results, to complement the analyses normally carried out
with commercial software packages.Comment: 29 pages, 4 figure
Index Reduction for Differential-Algebraic Equations with Mixed Matrices
Differential-algebraic equations (DAEs) are widely used for modeling of
dynamical systems. The difficulty in solving numerically a DAE is measured by
its differentiation index. For highly accurate simulation of dynamical systems,
it is important to convert high-index DAEs into low-index DAEs. Most of
existing simulation software packages for dynamical systems are equipped with
an index-reduction algorithm given by Mattsson and S\"{o}derlind.
Unfortunately, this algorithm fails if there are numerical cancellations.
These numerical cancellations are often caused by accurate constants in
structural equations. Distinguishing those accurate constants from generic
parameters that represent physical quantities, Murota and Iri introduced the
notion of a mixed matrix as a mathematical tool for faithful model description
in structural approach to systems analysis. For DAEs described with the use of
mixed matrices, efficient algorithms to compute the index have been developed
by exploiting matroid theory.
This paper presents an index-reduction algorithm for linear DAEs whose
coefficient matrices are mixed matrices, i.e., linear DAEs containing physical
quantities as parameters. Our algorithm detects numerical cancellations between
accurate constants, and transforms a DAE into an equivalent DAE to which
Mattsson--S\"{o}derlind's index-reduction algorithm is applicable. Our
algorithm is based on the combinatorial relaxation approach, which is a
framework to solve a linear algebraic problem by iteratively relaxing it into
an efficiently solvable combinatorial optimization problem. The algorithm does
not rely on symbolic manipulations but on fast combinatorial algorithms on
graphs and matroids. Furthermore, we provide an improved algorithm under an
assumption based on dimensional analysis of dynamical systems.Comment: A preliminary version of this paper is to appear in Proceedings of
the Eighth SIAM Workshop on Combinatorial Scientific Computing, Bergen,
Norway, June 201
On dual Schur domain decomposition method for linear first-order transient problems
This paper addresses some numerical and theoretical aspects of dual Schur
domain decomposition methods for linear first-order transient partial
differential equations. In this work, we consider the trapezoidal family of
schemes for integrating the ordinary differential equations (ODEs) for each
subdomain and present four different coupling methods, corresponding to
different algebraic constraints, for enforcing kinematic continuity on the
interface between the subdomains.
Method 1 (d-continuity) is based on the conventional approach using
continuity of the primary variable and we show that this method is unstable for
a lot of commonly used time integrators including the mid-point rule. To
alleviate this difficulty, we propose a new Method 2 (Modified d-continuity)
and prove its stability for coupling all time integrators in the trapezoidal
family (except the forward Euler). Method 3 (v-continuity) is based on
enforcing the continuity of the time derivative of the primary variable.
However, this constraint introduces a drift in the primary variable on the
interface. We present Method 4 (Baumgarte stabilized) which uses Baumgarte
stabilization to limit this drift and we derive bounds for the stabilization
parameter to ensure stability.
Our stability analysis is based on the ``energy'' method, and one of the main
contributions of this paper is the extension of the energy method (which was
previously introduced in the context of numerical methods for ODEs) to assess
the stability of numerical formulations for index-2 differential-algebraic
equations (DAEs).Comment: 22 Figures, 49 pages (double spacing using amsart
Recommended from our members
Acoustic scattering in waveguides that are discontinuous in geometry and material property
The scattering of acoustic waves at the discontinuity between two ducts of different heights is considered. At least one of the ducts is bounded by a membrane and, thus, the underlying eigenproblem is non-Sturm–Liouville. A mode-matching procedure, based on an appropriate orthogonality relation, reduces the problem to that of truncating and solving an infinite system
of linear equations. The distribution of power between the fluid regions and the membrane(s) is analysed. Further, it is shown that a fundamental
property of the truncated system is that the expression for power balance is always satisfied
- …