24,132 research outputs found
A Numerical Slow Manifold Approach to Model Reduction for Optimal Control of Multiple Time Scale ODE
Time scale separation is a natural property of many control systems that can
be ex- ploited, theoretically and numerically. We present a numerical scheme to
solve optimal control problems with considerable time scale separation that is
based on a model reduction approach that does not need the system to be
explicitly stated in singularly perturbed form. We present examples that
highlight the advantages and disadvantages of the method
Asymptotology of Chemical Reaction Networks
The concept of the limiting step is extended to the asymptotology of
multiscale reaction networks. Complete theory for linear networks with well
separated reaction rate constants is developed. We present algorithms for
explicit approximations of eigenvalues and eigenvectors of kinetic matrix.
Accuracy of estimates is proven. Performance of the algorithms is demonstrated
on simple examples. Application of algorithms to nonlinear systems is
discussed.Comment: 23 pages, 8 figures, 84 refs, Corrected Journal Versio
Empirical balanced truncation of nonlinear systems
Novel constructions of empirical controllability and observability gramians
for nonlinear systems for subsequent use in a balanced truncation style of
model reduction are proposed. The new gramians are based on a generalisation of
the fundamental solution for a Linear Time-Varying system. Relationships
between the given gramians for nonlinear systems and the standard gramians for
both Linear Time-Invariant and Linear Time-Varying systems are established as
well as relationships to prior constructions proposed for empirical gramians.
Application of the new gramians is illustrated through a sample test-system.Comment: LaTeX, 11 pages, 2 figure
Balanced truncation of perturbative representations of nonlinear systems
The paper presents a novel approach for a balanced truncation style of model reduction of a perturbative representation of a nonlinear system. Empirical controllability and observability gramians for nonlinear systems are employed to define a projection matrix. However, the projection matrix is applied to the perturbative representation of the system rather than directly to the exact nonlinear system. This is to achieve the required increase in efficiency desired of a reduced-order model. Application of the new method is illustrated through a sample test-system. The technique will be compared to the standard approach for reducing a perturbative representation of a nonlinear system
WavePacket: A Matlab package for numerical quantum dynamics. II: Open quantum systems, optimal control, and model reduction
WavePacket is an open-source program package for numeric simulations in
quantum dynamics. It can solve time-independent or time-dependent linear
Schr\"odinger and Liouville-von Neumann-equations in one or more dimensions.
Also coupled equations can be treated, which allows, e.g., to simulate
molecular quantum dynamics beyond the Born-Oppenheimer approximation.
Optionally accounting for the interaction with external electric fields within
the semi-classical dipole approximation, WavePacket can be used to simulate
experiments involving tailored light pulses in photo-induced physics or
chemistry. Being highly versatile and offering visualization of quantum
dynamics 'on the fly', WavePacket is well suited for teaching or research
projects in atomic, molecular and optical physics as well as in physical or
theoretical chemistry. Building on the previous Part I which dealt with closed
quantum systems and discrete variable representations, the present Part II
focuses on the dynamics of open quantum systems, with Lindblad operators
modeling dissipation and dephasing. This part also describes the WavePacket
function for optimal control of quantum dynamics, building on rapid
monotonically convergent iteration methods. Furthermore, two different
approaches to dimension reduction implemented in WavePacket are documented
here. In the first one, a balancing transformation based on the concepts of
controllability and observability Gramians is used to identify states that are
neither well controllable nor well observable. Those states are either
truncated or averaged out. In the other approach, the H2-error for a given
reduced dimensionality is minimized by H2 optimal model reduction techniques,
utilizing a bilinear iterative rational Krylov algorithm
- âŠ