Institute
General information about the institute, such as our mission statement, organizational structure, staff directory, history, directions, etc.

See More
Research
Scientific profile with all research groups, topics, collaborations, as well as columns on research at the institute.

See More
News
News and press releases about the institute, as well as a media archive.

See More
- News Overview
- Press Releases
Events
Overview of all events around the institute, such as talks, seminars, lectures, workshops, conferences and public events.

See More
Publications
Overview of all scientific publications of the institute, as well as our preprint and software repositories.

See More
Career
Information on open positions at the institute, benefits of working with us, graduate school, and postdoctoral supervision.

See More

Events
Multiscale II: analysis and efficient simulation of multiscale problems.

MPI für Mathematik in den Naturwissenschaften Leipzig

G3 10 (Lecture hall)

conference

08.10.97 10.10.97

Multiscale II: analysis and efficient simulation of multiscale problems.

The goal of the workshop is to stimulate an exchange of ideas between analysts and computational experts, and to focus on new approaches towards a better understanding of the intrinsic difficulties in problems with multiple scales in space and time.
In particular the workshop will focus on the following topics:

Scientific problems with multiple scales in space and time: Phase transitions, molecular design/vibrations, microstructures.
New analytical tools: H-measures, semiclassical measures, two-scale Young measures, averaging techniques.
New computational methods: Essential molecular dynamics, quantum-classical coupling, symplectic/symmetric discretization schemes.
Interaction of analytical and computational approaches.

The workshop is organized within the DFG-Schwerpunktprogramm DANSE.

Program

09:10 - 10:00	Luc Tartar (CMU/Leipzig) H-measures and their variants When a sequence converges weakly to 0 in L2(Ω) but not strongly (Ω being a subset of RN), one says that it contains oscillations or concentration effects depending upon if the limit of its square has a N-dimensional density or not. One quantitative way to analyze these oscillations or concentration effects is to use H-measures (which I have introduced for a few applications and which Patrick GÉRARD has introduced independently for other ones); they are measures on Ω ⨯ SN-1 which have a quadratic microlocal character (and make more precise the quadratic theorem of Compensated Compactness Theory). For dealing with problems with one characteristic length, Patrick GÉRARD then introduced a variant, called semi-classical measures, living on Ω ⨯ RN (quite similar to the idea which I had proposed independently to use H-measures after adding one dimension). It was then shown by Pierre-Louis LIONS and Thierry PAUL that the semi-classical measures could be introduced using WIGNER transform. I then found with Patrick GÉRARD how to use correlations instead, and why these objects are not good enough when at least two scales are present.Some other variants of H-measures can be defined, adapted to different questions, and I will describe the advantages and defects of some of these variants.
10:10 - 11:00	Gilles Francfort (Paris) Energy density for the wave equation on thin domains In this joint work with P. Gérard, 3-d/2-d asymptotic analysis of the wave equation is undertaken for cylindrical domains with vanishing thickness. The limit equation is a 2-d wave equation but its energy density is not the weak-limit of its 3-d analogue. The computation of that limit is performed microlocally through the computation of the semi-classical measure associated to the space-time derivatives of the field.
11:40 - 12:30	Giovanni Alberti (Pisa) A new approach to variational problems with multiple scales Many problems in science involve the creation and interaction of microstructures on multiple spatial or temporal scales. Most mathematical approaches have addressed "infinitely fine" microstructures or the case of a single small scale. In this talk we propose a new idea to deal with multiple scales and illustrate it with some simple examples. This is joint work with S. Müller.
14:00 - 14:50	M. Dellnitz (Bayreuth) The approximation of complicated dynamical behavior Multilevel techniques are presented for the efficient numerical approximation of complicated dynamical behavior. Concretely we develop (adaptive) methods which allow to extract statistical information on the underlying dynamical system. This is done by an approximation of natural invariant measures as well as (almost) cyclic dynamical components. We discuss issues concerning the implementation (e.g. parallelization strategies) and indicate potential applications of these methods (e.g. to the computation of Lyapunov exponents). The results are illustrated by several numerical examples.
15:00 - 15:50	P. Deuflhard (ZIB) Essential molecular dynamics: Progress in a new algorithmic approach Recently, the author -- together with his co-authors Michael Dellnitz, Oliver Junge, and Christof Schütte -- developed a new algorithmic approach to molecular dynamics, which is based on the computation of almost invariant sets of Hamiltonian systems. In this approach, only well-conditioned short term subtrajectories in lieu of ill-posed long term trajectories (which are typically used in Monte Carlo simulations) are exploited. The aim is to directly compute averages of physical observables, comformations and conformational changes - informations that are actually desired by computational chemistry. Mathematically speaking, such informations come out of the computation of the invariant measures and sets (corresponding to the dominant eigenvalue 1 and infinite relaxation time) and almost invariant measures and sets (corresponding to eigenmodes for eigenvalues close to 1 and therefore finite, but "large" relaxation times). An adaptive multilevel box method or subdivision technique is presented, which helps to solve the arising stochastic eigenvalue problem "fast". The basic concepts of the new algorithm and recent speed-ups will be presented.

09:00 - 09:50	F. Bornemann (ZIB) Weak convergence methods and adiabatic results in classical and quantum mechanics We will show that weak convergence methods are conveniently suited to explicitly access singular limits of a certain family of mechanical systems with multiple time scales. This family turns out to be characterized by the existence of sufficiently many adiabatic invariants.The key step is the idenfication of the weak limits of all those quadratic quantities which carry important information for the limit system. This idenfication becomes possible by a weak convergence analogue of the Virial Theorem, resulting in certain matrix commutation relations.We will address natural mechanical systems with strong constraining potentials and the adiabatic theorem of quantum mechanics. These examples will show that our approach considerably extends the possibility of passing through resonances.
10:00 - 10:50	Christof Schütte (ZIB) Multiple scales and nonadiabatic effects in molecular quantum dynamics The full quantum dynamical (QD) description of molecules includes multiple scales in particular if particles of essentially different masses are included. This makes full QD simulations of larger molecules simply impossible. In real life molecular dynamics one deals with this problem by introducing simplified "mixed quantum-classical" models which describe the heavy particles of the molecule by the means of classical mechanics and only a small portion by the means of a wavefunction. This talk will present two of these simplified models. By using the techniques of the previous talk "Weak Convergence Methods and Adiabatic Results in Classical and Quantum Mechanics", it will be discussed under which circumstances these models gain a suitable approximation of the true QD solution. The discussion will be concentrated on the understanding of the nonadiabatic effects which may cause such simplified descriptions to fail and on the alternatives which may prevent such failure.
11:30 - 12:20	S. Reich (ZIB) Adiabatic invariants and their preservation under symplectic discretizations We will apply normal form theory to highly oscillatory Hamiltonian systems and discuss the preservation of adiabatic invariants over exponentially long periods of time. Next, the conservation of adiabatic invariants under a symplectic discretization will be considered. Backward error analysis will provide us the appropriate tool to show that, again, the invariants are preserved over exponentially long periods of time.
14:00 - 14:50	Carsten Carstensen (Kiel) Numerical analysis of the (relaxed) double well problem The mathematical modelling of phase transitions in alloys, or optimal design of composites, is concerned with non-convex variational problems, where a minimum may or may not be attained. Young's model example, the two-well or double-well problem (P), is considered in this talk. Convexification yields a relaxed problem (RP) related to (P). The loss of information we face in recasting (P) as (RP) is not substantial as long as we are interested in the global displacements, stress fields and Young measures (which describe oscillations creating microstructures); these variables defined by (P) can be computed from solutions of (RP).Typically, the problem (RP) is convex, but not strictly convex. Hence error estimates for finite element methods are more difficult to prove than for simpler uniformly convex problems. A priori and a posteriori error estimates are presented and illustrated in numerical examples which indicate that it is more reasonable to perform a numerical analysis of (RP) rather than of (P).

09:00 - 09:50	G. Schneider (Hannover) The amplitude equation for degenerate subcritical bifurcations in pattern formating systems We consider a model for a pattern formating system on the infinite line which leads to a subcritical bifurcation. For the degenerate case we use multiple scaling analysis to derive $A_T =a_0 A + a_1 A_{XX} + a_2 \|A\|^2 A+a_3 \|A\|^2 A_X + a_4 A^2 \overline{A}_X + a_5 \|A\|^4 A$ as an amplitude equation for the envelope A of the bifurcating pattern which is modulated slowly in time and in space. We show exact estimates between the approximations obtained via the amplitude equation and true solutions of the original system. Moreover, we show that every small solution of the original system develops in such a way that it can be described after a certain time by the solutions of the amplitude equation. The difficulty is to show the estimates on an $O(1/\delta^4)$-time scale in contrast to $O(1/\delta^4)$ for the classical Ginzburg-Landau equation, if $\delta$; is the order of the amplitude. This theory allows the description of modulated N-pulses in the original system.
10:00 - 10:50	Alexander Mielke (Hannover) Flow properties for Young-measure solutions of semilinear hyperbolic systems We consider a semilinear hyperbolic system of d equations of first order. We consider sequences of initial data which approximate a Young measure and study the Young-measure limit of the associated sequence of solutions. Thus we are able to define the notion of Cauchy problem and semigroups for Young measure solutions, at least for d≤2 by using compactness through compensation. For d≥3 the subclass of product measure solutions is considered, where each Cauchy problem has a unique solution. Continuous dependence of the associated semigroup on the initial data can be established in the Wasserstein topology and, under additional structural condition on the nonlinearity, also in the weak* topology. The reported research extends previous work of L. Tartar and is partly joint work with Florian Theil.
11:30 - 13:00	Round table discussion (volunteers welcome)

Multiscale II: analysis and efficient simulation of multiscale problems.

Program

Organizers

P. Deuflhard

P. Dolzmann

Stefan Müller