From SOI to Abstract Electric-Thermal-1D Multiscale Modeling for First Order Thermal Effects

Self-heating occurs in integrated circuits, specially for SOI-based devices. Naturally, the heat distribution affects the circuit’s functionality. For reliable designs in SOI-chip technology, and other applications, the thermal aspects have to be addressed. Therefore we develop a model, which is based on distributed 1D and lumped 0D elements, and takes into account that heat is stored and slowly conducted between elements. The emerging coupled multiscale system of heat evolution and electric network consists of parabolic partial-differential (thermal part) and differential-algebraic equations (electric network part). For the thermal model, we verify properties as positivity and strict passivity. Since time scales differ largely, the coupled problem exhibits multirate potential.     

Structural analysis of electrical circuits including magnetoquasistatic devices

Modeling electric circuits that contain magnetoquasistatic (MQS) devices leads to a coupled system of differential-algebraic equations (DAEs). In our case, the MQS device is described by the eddy current problem being already discretized in space (via edge-elements). This yields a DAE with a properly stated leading term, which has to be solved in the time domain. We are interested in structural properties of this system, which are important for numerical integration. Applying a standard projection technique, we are able to deduce topological conditions such that the tractability index of the coupled problem does not exceed two. Although index-2, we can conclude that the numerical difficulties for this problem are not severe due to a linear dependency on index-2 variables.     

Convergence analysis of a partial differential algebraic system from coupling a semiconductor model to a circuit model

We are interested in circuit simulation including distributed semiconductor models. The circuit itself is modeled by the modified nodal analysis. The stationary drift diffusion equations are used to describe the semiconductors. The complete system is then a partial differential-algebraic system. We discretize it first in space with finite elements and the Scharfetter-Gummel discretization. The resulting semi-discrete system can be analyzed as a differential-algebraic equation with properly stated leading term. We present topological index one criteria. They coincide with previous results for the non-discretized partial differential-algebraic equation. For the time discretization we use standard BDF methods (implicit Gear formulas). Finally we derive a convergence estimate for the whole partial differential-algebraic system close to equilibrium.     

Perturbation index of linear partial differential-algebraic equations with a hyperbolic part

This paper deals with linear partial differential-algebraic equations (PDAEs) which have a hyperbolic part. If the spatial differential operator satisfies a Gårding-type inequality in a suitable function space setting, a perturbation index can be defined. Theoretical and practical examples are considered.     

On data-dependence of exponential stability and stability radii for linear time-varying differential-algebraic systems

This paper is addressed to some questions concerning the exponential stability and its robustness measure for linear time-varying differential-algebraic systems of index 1. First, the Bohl exponent theory that is well known for ordinary differential equations is extended to differential-algebraic equations. Then, it is investigated that how the Bohl exponent and the stability radii with respect to dynamic perturbations for a differential-algebraic system depend on the system data. The paper can be considered as a continued and complementary part to a recent paper on stability radii for time-varying differential-algebraic equations [N.H. Du, V.H. Linh, Stability radii for linear time-varying differential-algebraic equations with respect to dynamic perturbations, J. Differential Equations 230 (2006) 579-599].     

On coupled transversal and axial motions of two beams with a joint

In this paper we develop and analyze a mathematical model for combined axial and transverse motions of two Euler-Bernoulli beams coupled through a joint composed of two rigid bodies. The motivation for this problem comes from the need to accurately model damping and joints for the next generation of inflatable/rigidizable space structures. We assume Kelvin-Voigt damping in the two beams whose motions are coupled through a joint which includes an internal moment. The resulting equations of motion consist of four, second-order in time, partial differential equations, four second-order ordinary differential equations, and certain compatibility boundary conditions. The system is re-cast as an abstract second-order differential equation in an appropriate Hilbert space, consisting of function spaces describing the distributed beam deflections, and a finite-dimensional space that projects important features at the joint boundary. Semigroup theory is used to prove the system is well posed, and that with positive damping parameters the resulting semigroup is analytic and exponentially stable. The spectrum of the infinitesimal generator is characterized.     

A nonlinear truck model and its treatment as a multibody system

A planar vertical truck model with nonlinear suspension and its multibody system formulation are presented. The equations of motion of the model form a system of differential-algebraic equations (DAEs). All equations are given explicitly, including a complete set of parameter values, consistent initial values, and a sample road excitation. Thus the truck model allows various investigations of the specific DAE effects and represents a test problem for algorithms in control theory, mechanics of multibody systems, and numerical analysis. Several numerical tests show the properties of the model.     

Index-aware model order reduction for differential-algebraic equations

We introduce a model order reduction (MOR) procedure for differential-algebraic equations, which is based on the intrinsic differential equation contained in the starting system and on the remaining algebraic constraints. The decoupling procedure in differential and algebraic part is based on the projector and matrix chain which leads to the definition of tractability index. The differential part can be reduced by using any MOR method, we use Krylov-based projection methods to illustrate our approach. The reduction on the differential part induces a reduction on the algebraic part. In this paper, we present the method for index-1 differential-algebraic equations. We implement numerically this procedure and show numerical evidence of its validity.     

Wavelet-based adaptive grids for multirate partial differential-algebraic equations

The mathematical model of electric circuits yields systems of differential-algebraic equations (DAEs). In radio frequency applications, a multivariate model of oscillatory signals transforms the DAEs into a system of multirate partial differential-algebraic equations (MPDAEs). Considering quasiperiodic signals, an approach based on a method of characteristics yields efficient numerical schemes for the MPDAEs in time domain. If additionally digital signal structures occur, an adaptive grid is required to achieve the efficiency of the technique. We present a strategy applying a wavelet transformation to construct a mesh for resolving steep gradients in respective signals. Consequently, we employ finite difference methods to determine an approximative solution of characteristic systems in according grid points. Numerical simulations demonstrate the performance of the adaptive grid generation, where radio frequency signals with digital structures are resolved.     

微分代数方程去奇异化分析

利用去奇异化方法讨论了拟线性微分代数方程在奇点邻域内光滑解的性质.通过尺度参数的微分同胚变换,将拟线性微分代数方程转化为相应的常微分方程,从而构造出在孤立奇点邻域内的初始微分代数方程的光滑解,给出解存在的充分条件,并进一步讨论了解的性质.     

Simulating the motion of the leech: A biomechanical application of DAEs

In this paper a mathematical model is developed for the dynamical behaviour of a hydrostatic skeleton. The basic configuration is taken from the worm-like shape of the medicinal leech. It consists of a sequence of hexahedra with damped elastic springs as edges to model the various parts of the musculature. The system is stabilized by the constraint of constant volume either in the whole body or in prescribed compartments. We set up Lagrange's equations of motion with the Lagrange multipliers being the pressure values in the compartments. The equations of motion lead to a large differential-algebraic system which is solved by an application of semi-explicit numerical methods. Though the model has not yet been adapted to experimental data, first simulations with a simplified set of parameters show that it is capable of generating basic movements of the leech such as crawling and swimming. This revised version was published online in June 2006 with corrections to the Cover Date.     

A Joint DAE/PDE Model for Interconnected Electrical Networks

To model transmission lines effects in integrated circuits, we couple the network equations for the circuits with the telegrapher's equations for the transmission lines. This results in an initial/boundary value problem for a mixed system of Differential-Algebraic Equations (DAEs) and hyperbolic Partial Differential Equations (PDEs). By semidiscretization the system is transformed into differential-algebraic equations in time only. We apply this modeling approach to a CMOS ring oscillator, an oscillatory circuit with transmission lines as coupling units, and discuss the simulation results.     

Amplitude–shape approximation as an extension of separation of variables

Separation of variables is a well‐known technique for solving differential equations. However, it is seldom used in practical applications since it is impossible to carry out a separation of variables in most cases. In this paper, we propose the amplitude–shape approximation (ASA) which may be considered as an extension of the separation of variables method for ordinary differential equations. The main idea of the ASA is to write the solution as a product of an amplitude function and a shape function, both depending on time, and may be viewed as an incomplete separation of variables. In fact, it will be seen that such a separation exists naturally when the method of lines is used to solve certain classes of coupled partial differential equations. We derive new conditions which may be used to solve the shape equations directly and present a numerical algorithm for solving the resulting system of ordinary differential equations for the amplitude functions. Alternatively, we propose a numerical method, similar to the well‐established exponential time differencing method, for solving the shape equations. We consider stability conditions for the specific case corresponding to the explicit Euler method. We also consider a generalization of the method for solving systems of coupled partial differential equations. Finally, we consider the simple reaction diffusion equation and a numerical example from chemical kinetics to demonstrate the effectiveness of the method. The ASA results in far superior numerical results when the relative errors are compared to the separation of variables method. Furthermore, the method leads to a reduction in CPU time as compared to using the Rosenbrock semi‐implicit method for solving a stiff system of ordinary differential equations resulting from a method of lines solution of a coupled pair of partial differential equations. The present amplitude–shape method is a simplified version of previous ones due to the use of a linear approximation to the time dependence of the shape function. Copyright © 2007 John Wiley & Sons, Ltd.     

Mathematical modelling and analysis of temperature effects in MEMS

This article is concerned with the mathematical analysis of models for Micro-Electro-Mechanical Systems (MEMS). These models arise in the form of coupled partial differential equations with a moving boundary. Although MEMS devices are often operated in non-isothermal environments, temperature is often neglected in the mathematical investigations. In light of this finding the focus of our modelling is to incorporate temperature and the related material properties. We derive a full model for this coupling and discuss a simplified version as well. Lastly, we prove local well-posedness in time and also global well-posedness under additional assumptions on the model’s parameters.     

Global well-posedness of the stochastic 2D Boussinesq equations with partial viscosity

This paper deals with the stochastic 2D Boussinesq equations with partial viscosity. This is a coupled system of Navier-Stokes/Euler equations and the transport equation for temperature under additive noise. Global well-posedness result of this system under partial viscosity is proved by using classical energy estimates method.     

Global existence of weak solutions to a diffuse interface model for magnetic fluids

This article is devoted to the derivation and analysis of a system of partial differential equations modeling a diffuse interface flow of two Newtonian incompressible magnetic fluids. The system consists of the incompressible Navier–Stokes equations coupled with an evolutionary equation for the magnetization vector and the Cahn–Hilliard equations. We show global in time existence of weak solutions to the system using the time discretization method.     

Modeling of the Crosstalk Phenomenon by Bilateral Coupling of the Non-Stationary Maxwell Equations and the Circuit Equations

In this paper, we introduce a general modeling approach for the crosstalk phenomenon within electro-magnetic systems. We model the phenomenon by bilateral coupling of two sets of differential equations, namely non-stationary Maxwell equations and circuit equations. Considering these two sets of equations as input-output systems, the bilateral coupling approach connects the output of one system to the input of the other system and conversely, via coupling and re-coupling relations. The coupling of these two sets leads to a set of partial differential-algebraic equations modeling the crosstalk phenomenon. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamic optimization and event location in atmospheric chemistry

A model coupling differential equations and an optimization problem with equality and inequality constraints is described. The first order optimality conditions are coupled with the differential equations to form a differential-algebraic system. Inequality constraints for the optimization variables induce discontinuities in the first derivatives of the solution in time. Tracking techniques based on sensitivity analysis and second order multistep methods are proposed to locate the discontinuities. Application to the dynamics of organic atmospheric aerosol particles is given to illustrate the detection of the activation of constraints. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)