Refine
Year of publication
Document Type
- Preprint (52)
- Article (45)
- Monograph/Edited Volume (21)
- Doctoral Thesis (1)
- Postprint (1)
Keywords
- aerosol size distribution (2)
- inversion (2)
- Ill-posed problem (1)
- MHD-equations (1)
- Magnetfeld-Satellit (1)
- Magnetic field measurements (1)
- Magnetische Feldmessungen (1)
- Magnetometer-Kalibrierung (1)
- Multiwavelength LIDAR (1)
- Planetary Rings (1)
Institute
- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (120) (remove)
The ill-posed inversion of multiwavelength lidar data by a hybrid method of variable projection
(1999)
Based on the data of the Magion2 subsatellite of the Intercosmos24 satellite, an example of small-scale irregularities of the electron concentration with linear dimensions l ~ 100-300 m in the polar ion- osphere of the morning sector under field-aligned currents at altitudes of 1800-2030 km during the main phase of the magnetic storm of June 13, 1990 is presented. The dependence of the spectral index of the above small-scale irregularities on latitude is determined for the first time. Certain mechanisms of the generation of these small-scale irregularities are also qualitatively discussed.
We analyze the variability in the x-ray lightcurves of the black hole candidate Cygnus X-1 by linear and nonlinear time series analysis methods. While a linear model describes the overall second order properties of the observed data well, surrogate data analysis reveals a significant deviation from linearity. We discuss the relation between shot noise models usually applied to analyze these data and linear stochastic autoregressive models. We debate statistical and interpretational issues of surrogate data testing for the present context. Finally, we suggest a combination of tools from linear and nonlinear time series analysis methods as a procedure to test the predictions of astrophysical models on observed data.
Laser beam melt ablation - a contact-free machining process - offers several advantages compared to conventional processing mechanisms: there exists no tool wear and even extremely hard or brittle materials can be processed. During ablation the workpiece is molten by a CO2-laser beam, this melt is then driven out by the impulse of a process gas. The idea behind laser ablation is rather simple, but it has a major limitation in practical applications: with increasing ablation rates surface quality of the workpiece processed declines rapidly. At high ablation rates, depending on the process parameters different periodic-like structures can be observed on the ablated surface. These structures show a dependence on the line energy, which has been identified as a fundamental control parameter. In dependence on this parameter several regimes with different behaviours of the process have been separated. These regimes are distinguishable as well in the surfaces obtained as in the signals gained by the measurement of the process emissions. Further aim is to identify the different modes of the system and reach a deeper understanding of the dynamics of the molten material in order to understand the formation of these surface structures. With this it should be possible to influence the system in the direction of avoiding structure formation even at high ablation rates. Relying on the results on-line monitoring and control of the process should be studied.
We analyse the X-ray light curves of compact objects using linear and nonlinear time series analysis methods. A Power Density Spectrum (PDS) describes the overall second order properties of the observed data well. To look beyond we propose the nonlinear Q-statistic to detect an asymmetry of the time series. This allows us to find relevant time scales. This method even grants a subclassification of the known states of X-ray sources.
Acoustic emission signals generated during high speed cutting of steel are investigated. The data are represen ted in time-folded form. Several methods from linear and nonlinear data analysis based on time- and frequency- domain are applied to the data and reveal signatures of the observed acoustic emission signal. These investiga tions are necessary for modeling the cutting process by means of differential equations.
We investigate the relationship between precipitation and runoff data from a small forested catchment in the Harz mountains (Germany). For this purpose, we develop a conceptual model including memory effects to predict the runoff signal using the precipitation data as input. An enhanced variant of the model also includes air temperature as input variable. We show in terms of correlation functions that this model describes main dynamical properties of the runoff, especially the delay between rain event and runoff response as the annual persistence in the runoff data.
In den letzten 2 Jahrzehnten des 20. Jahrhunderts hat sich mit der rasanten Entwicklung der Nichtlinearen Wissenschaften ein weiterer Umbruch vollzogen, der eine ausgepraegte Nachhaltigkeit in Wissenschaft und Technik ebenso wie in der Gesellschaft erwarten laesst. Die Nichtlinearen Wissenschaften werden auch als Nichtlineare Dynamik, Wissenschaft Komplexer Systeme oder etwas eingegrenzt Chaostheorie bezeichnet.
The 300 year record of the yearly sunspot numbers and numerically generated trajectory of the solar inertial motion (SIM) were subjects of a synchronization analysis. Phase synchronization of the sunspot cycle and a fast component of the SIM have been found and confirmed with statistical significance in three epochs (1727-1757, 1802-1832 and 1863-1922) of the entire 1700-1997 record. This result can be considered as a quantitative support for the hypothesis that there is a weak interaction of gravity and solar activity.
The structure of time series and letter sequences is investigated using the concepts of entropy and complexity. First conditional entropy and transinformation are introduced and several generalizations are discussed. Further several measures of complexity are introduced and discussed. The capability of these concepts to describe the structure of time series and letter sequences generated by nonlinear maps, data series from meteorology, astrophysics, cardiology, cognitive psychology and finance is investigated. The relation between the complexity and the predictability of informational strings is discussed. The relation between local order and the predictability of time series is investigated.
Observational data of natural systems, as measured in medical measurements are typically quite different from those obtained in laboratories. Due to the peculiarities of these data, wellknown characteristics, such as power spectra or fractal dimension, often do not provide a suitable description. To study such data, we present here some measures of complexity, which are basing on symbolic dynamics. Firstly, a motivation for using symbolic dynamics and measures of complexity in data analysis based on the logistic map is given and next, two applications to medical data are shown. We demonstrate that symbolic dynamics is a useful tool for the risk assessment of patients after myocardial infarction as well as for the evaluation of th e architecture of human cancellous bone.
Charged dust grains in circumplanetary environments experience, beyond various deterministic forces, also stochastic perturbations caused, by fluctuations of the magnetic field, the charge of the grains, by chaotic rotation of aspherical grains, etc. Here we investigate the dynamics of a dust population in a circular orbit around a planet which is perturbed by a stochastic planetary magnetic field B', modeled by an isotropically Gaussian white noise. The resulting perturbation equations give rise to a modified diffusion of the inclinations i and eccentricities e. The diffusion coefficient is found to be D proportional to w^2 O /n^2 , where the gyrofrequency, the Kepler frequency, and the synodic frequency are denoted by w , O, and n, respectively. This behavior has been checked against numerical simulations. We have chosen dust grains (1 m in radius) ejected from Jupiter's satellite Europa in circular equatorial orbits around Jupiter and integrated numerically their trajectories over their typical lifetimes (100 years). The particles were exposed to a Gaussian fluctuating magnetic field B' with the same statistical properties as in the analytical treatment. These simulations have confirmed the analytical results. The theoretical studies showed the statistical properties of B' to be of decisive importance. To estimate them, we analyzed the magnetic field data obtained by the Galileo spacecraft magnetometer at Jupiter and found almost Gaussian fluctuations of about 5% of the mean field and exponentially decaying correlations. This results in a diffusion of orbital inclinations and eccentricities of the dust grains of about ten percent over the lifetime of the particles. For smaller dusty motes or for close-in particles (e.g., in Jovian gossamer rings) stochastics might well dominate the dynamics.
In the last decade, there has been an increasing interest in compensating thermally induced errors to improve the manufacturing accuracy of modular tool systems. These modular tool systems are interfaces between spindle and workpiece and consist of several complicatedly formed parts. Their thermal behavior is dominated by nonlinearities, delay and hysteresis effects even in tools with simpler geometry and it is difficult to describe it theoretically. Due to the dominant nonlinear nature of this behavior the so far used linear regression between the temperatures and the displacements is insufficient. Therefore, in this study we test the hypothesis whether we can reliably predict such thermal displacements via nonlinear temperature-displacement regression functions. These functions are estimated firstly from learning measurements using the alternating conditional expectation (ACE) algorithm and then tested on independent data sets. First, we analyze data that were generated by a finite element spindle model. We find that our approach is a powerful tool to describe the relation between temperatures and displacements for simulated data. Next, we analyze the temperature-displacement relationship in a silent real experimental setup, where the tool system is thermally forced. Again, the ACE-algorithm is powerful to estimate the deformation with high precision. The corresponding errors obtained by using the nonlinear regression approach are 10-fold lower in comparison to multiple linear regression analysis. Finally, we investigate the thermal behavior of a modular tool system in a working milling machine and get again promising results. The thermally induced errors can be estimated with 1-2${mu m}$ accuracy using this nonlinear regression analysis. Therefore, this approach seems to be very useful for the development of new modular tool systems.
In the last decade, there has been an increasing interest in compensating thermally induced errors to improve the manufacturing accuracy of modular tool systems. These modular tool systems are interfaces between spindle and workpiece and consist of several complicatedly formed parts. Their thermal behavior is dominated by nonlinearities, delay and hysteresis effects even in tools with simpler geometry and it is difficult to describe it theoretically. Due to the dominant nonlinear nature of this behavior the so far used linear regression between the temperatures and the displacements is insufficient. Therefore, in this study we test the hypothesis whether we can reliably predict such thermal displacements via nonlinear temperature-displacement regression functions. These functions are estimated firstly from learning measurements using the alternating conditional expectation (ACE) algorithm and then tested on independent data sets. First, we analyze data that were generated by a finite element spindle model. We find that our approach is a powerful tool to describe the relation between temperatures and displacements for simulated data. Next, we analyze the temperature-displacement relationship in a silent real experimental setup, where the tool system is thermally forced. Again, the ACE-algorithm is powerful to estimate the deformation with high precision. The corresponding errors obtained by using the nonlinear regression approach are 10-fold lower in comparison to multiple linear regression analysis. Finally, we investigate the thermal behavior of a modular tool system in a working milling machine and get again promising results. The thermally inducedaccuracy using this nonlinear regression analysis. Therefore, this approach seems to be very useful for the development of new modular tool systems. errors can be estimated with 1-2 micrometer
Tätigkeitsbericht 1994-2000
(2004)
Das Interdisziplinäre Zentrum für Nichtlineare Dynamik an der Universität Potsdam verbindet theoretisch-methodische Untersuchungen in Mathematik und theoretischer Physik mit einer Vielzahl anderer Wissenschaften und zielt auf eine fruchtbare Wechselwirkung zwischen Theorie und Experiment. Unter Einbezug von Instituten und Großforschungseinrichtungen, die insbesondere im Potsdamer Raum angesiedelt sind, soll sich ein überregional bedeutender Schwerpunkt entwickeln, wie er an keiner anderen deutschen Universität in gleicher Weise interdisziplinär angelegt ist.
CHAMP (CHAllenging Minisatellite Payload) is a German small satellite mission to study the earth's gravity field, magnetic field and upper atmosphere. Thanks to the good condition of the satellite so far, the planned 5 years mission is extended to year 2009. The satellite provides continuously a large quantity of measurement data for the purpose of Earth study. The measurements of the magnetic field are undertaken by two Fluxgate Magnetometers (vector magnetometer) and one Overhauser Magnetometer (scalar magnetometer) flown on CHAMP. In order to ensure the quality of the data during the whole mission, the calibration of the magnetometers has to be performed routinely in orbit. The scalar magnetometer serves as the magnetic reference and its readings are compared with the readings of the vector magnetometer. The readings of the vector magnetometer are corrected by the parameters that are derived from this comparison, which is called the scalar calibration. In the routine processing, these calibration parameters are updated every 15 days by means of scalar calibration. There are also magnetic effects coming from the satellite which disturb the measurements. Most of them have been characterized during tests before launch. Among them are the remanent magnetization of the spacecraft and fields generated by currents. They are all considered to be constant over the mission life. The 8 years of operation experience allow us to investigate the long-term behaviors of the magnetometers and the satellite systems. According to the investigation, it was found that for example the scale factors of the FGM show obvious long-term changes which can be described by logarithmic functions. The other parameters (offsets and angles between the three components) can be considered constant. If these continuous parameters are applied for the FGM data processing, the disagreement between the OVM and the FGM readings is limited to \pm1nT over the whole mission. This demonstrates, the magnetometers on CHAMP exhibit a very good stability. However, the daily correction of the parameter Z component offset of the FGM improves the agreement between the magnetometers markedly. The Z component offset plays a very important role for the data quality. It exhibits a linear relationship with the standard deviation of the disagreement between the OVM and the FGM readings. After Z offset correction, the errors are limited to \pm0.5nT (equivalent to a standard deviation of 0.2nT). We improved the corrections of the spacecraft field which are not taken into account in the routine processing. Such disturbance field, e.g. from the power supply system of the satellite, show some systematic errors in the FGM data and are misinterpreted in 9-parameter calibration, which brings false local time related variation of the calibration parameters. These corrections are made by applying a mathematical model to the measured currents. This non-linear model is derived from an inversion technique. If the disturbance field of the satellite body are fully corrected, the standard deviation of scalar error \triangle B remains about 0.1nT. Additionally, in order to keep the OVM readings a reliable standard, the imperfect coefficients of the torquer current correction for the OVM are redetermined by solving a minimization problem. The temporal variation of the spacecraft remanent field is investigated. It was found that the average magnetic moment of the magneto-torquers reflects well the moment of the satellite. This allows for a continuous correction of the spacecraft field. The reasons for the possible unknown systemic error are discussed in this thesis. Particularly, both temperature uncertainties and time errors have influence on the FGM data. Based on the results of this thesis the data processing of future magnetic missions can be designed in an improved way. In particular, the upcoming ESA mission Swarm can take advantage of our findings and provide all the auxiliary measurements needed for a proper recovery of the ambient magnetic field.
Aus dem Inhalt: 1. Einführung 2. Motivation für die nichtlineare Dynamik 3. Logistische Abbildung (Parabel-Abbildung) 4. Lorenz-Gleichungen 5. Fraktale Selbstähnlichkeit 6. Die Brownsche Bewegung 7. Stöße & Billards 8. Körper mit gravitativer Wechselwirkung 9. Glossar 10. Turbo-Pascal-Texte 11. IDL-Texte 12. Reduce-Texte
The statistical analysis of the variations of the dayly-mean frequency of the maximum ionospheric electron density foF2 is performed in connection with the occurrence of (more than 60) earthquakes with magnitudes M > 6.0, depths h < 80 km and distances from the vertical sounding station R < 1000 km. For the study, data of the Tokyo sounding station are used, which were registered every hour in the years 1957-1990. It is shown that, on the average, foF2 decreases before the earthquakes. One day before the shock the decrease amounts to about 5 %. The statistical reliability of this phenomenon is obtained to be better than 0.95. Further, the variations of the occurrence probability of the turbulization of the F-layer (F spread) are investigated for (more than 260) earthquakes with M > 5.5, h < 80 km, R < 1000 km. For the analysis, data of the Japanese station Akita from 1969-1990 are used, which were obtained every hour. It is found that before the earthquakes the occurrence probability of F spread decreases. In the week before the event, the decrease has values of more than 10 %. The statistical reliability of this phenomenon is also larger than 0.95. Examining the seismo-ionospheric effects, here periods of time with weak heliogeomagnetic disturbances are considered, the Wolf number is less than 100 and the index ∑ Kp is smaller than 30.
In the present work, phenomena in the ionosphere are studied, which are connected with earthquakes (16 events) having a depth of less than 50 km and a magnitude M larger than 4. Analysed are night-time Es-spread effects using data of the vertical sounding station Petropavlovsk- Kanchatsky (φ=53.0°, λ=158.7°) from May 2004 until August 2004 registered every 15 minutes. It is found that the maximum distance of the earthquake from the sounding station, where pre-seismic phenomena are yet observable, depends on the magnitude of the earthquake. Further it is shown that 1-2 days before the earthquakes, in the premidnight hours, the appearance of Es-spread increases. The reliability of this increase amounts to 0.95.
A model of the generation of pulses of local electric fields with characteristic time scales of 1–10 minutes is considered for atmospheric conditions above fracture regions of earthquakes. In the model, it is proposed that aerosols, increased ionization velocity and upstreaming air flows occur at night-time conditions. The pulses of local electric fields cause respective pulses of infrared emissions. But infrared emissions with time scales of 1–10 minutes were not observed up to now experimentally. The authors think, that the considered non-stationary field and radiation effects might be a new-type of applicable earthquake indicators and ask to perform special earth-based and satellite observations of the night-time atmosphere in seismoactive fracture regions.
We numerically investigate nonlinear asymmetric square patterns in a horizontal convection layer with up-down reflection symmetry. As a novel feature we find the patterns to appear via the skewed varicose instability of rolls. The time-independent nonlinear state is generated by two unstable checkerboard (symmetric square) patterns and their nonlinear interaction. As the bouyancy forces increase, the interacting modes give rise to bifurcations leading to a periodic alternation between a nonequilateral hexagonal pattern and the square pattern or to different kinds of standing oscillations.
In this paper an analysis of the excitation conditions of mirror waves is done, which propagate parallel to an external magnetic field. There are found analytical expressions for the dispersion relations of the waves in case of different plasma conditions. These relations may be used in future to develop the nonlinear theory of mirror waves. In comparison with former analytical works, in the study the inuence of the magnetic field and nite temperatures of the ions parallel to the magnetic field are taken into account. Application is done for the earth's magnetosheath.
Basing on recent solar models, the excitation of ion-acoustic turbulence in the weaklycollisional, fully and partially-ionized regions of the solar atmosphere is investigated. Within the frame of hydrodynamics, conditions are found under which the heating of the plasma by ion-acoustic type waves is more effective than the Joule heating. Taking into account wave and Joule heating effects, a nonlinear differential equation is derived, which describes the evolution of nonlinear ion-acoustic waves in the collisional plasma.
A numerical MHD model is developed to investigate acceleration and heating of both thermal and auroral plasma. This is done for magnetospheric flux tubes in which intensive field aligned currents flow. To give each of these tubes, the empirical Tsyganenko model of the magnetospheric field is used. The parameters of the background plasma outside the flux tube as well as the strength of the electric field of magnetospheric convection are given. Performing the numerical calculations, the distributions of the plasma densities, velocities, temperatures, parallel electric field and current, and of the coefficients of thermal conductivity are obtained in a self-consistent way. It is found that EIC turbulence develops effectively in the thermal plasma. The parallel electric field develops under the action of the anomalous resistivity. This electric field accelerates both the thermal and the auroral plasma. The thermal turbulent plasma is also subjected to an intensive heating. The increase of the plasma of the Earth's ionosphere. Besides, studying the growth and dispersion properties of oblique ion cyclotron waves excited in a drifting magnetized plasma, it is shown that under non-stationary conditions such waves may reveal the properties of bursts of polarized transverse electromagnetic waves at frequencies near the patron gyrofrequency.
This paper deals with the Mie scattering kernels for multi-spectral data. The kernels may be represented in form of power series. Furthermore, the singular-value spectrum and the degree of ill-posedness in dependence on the refractive index of the particles are numerically approximated. A special hybrid regularization technique allows us to determine via inversion the particle distributions of different types.
Contents: 1 Introduction 2 Experiment 3 Data 4 Symbolic dynamics 4.1 Symbolic dynamics as a tool for data analysis 4.2 2-symbols coding 4.3 3-symbols coding 5 Measures of complexity 5.1 Word statistics 5.2 Shannon entropy 6 Testing for stationarity 6.1 Stationarity 6.2 Time series of cycle durations 6.3 Chi-square test 7 Control parameters in the production of rhythms 8 Analysis of relative phases 9 Discussion 10 Outlook
A numerical bifurcation analysis of the electrically driven plane sheet pinch is presented. The electrical conductivity varies across the sheet such as to allow instability of the quiescent basic state at some critical Hartmann number. The most unstable perturbation is the two-dimensional tearing mode. Restricting the whole problem to two spatial dimensions, this mode is followed up to a time-asymptotic steady state, which proves to be sensitive to three-dimensional perturbations even close to the point where the primary instability sets in. A comprehensive three-dimensional stability analysis of the two-dimensional steady tearing-mode state is performed by varying parameters of the sheet pinch. The instability with respect to three-dimensional perturbations is suppressed by a sufficiently strong magnetic field in the invariant direction of the equilibrium. For a special choice of the system parameters, the unstably perturbed state is followed up in its nonlinear evolution and is found to approach a three-dimensional steady state.
We investigate numerically the appearance of heteroclinic behavior in a three-dimensional, buoyancy-driven fluid layer with stress-free top and bottom boundaries, a square horizontal periodicity with a small aspect ratio, and rotation at low to moderate rates about a vertical axis. The Prandtl number is 6.8. If the rotation is not too slow, the skewed-varicose instability leads from stationary rolls to a stationary mixed-mode solution, which in turn loses stability to a heteroclinic cycle formed by unstable roll states and connections between them. The unstable eigenvectors of these roll states are also of the skewed-varicose or mixed-mode type and in some parameter regions skewed-varicose like shearing oscillations as well as square patterns are involved in the cycle. Always present weak noise leads to irregular horizontal translations of the convection pattern and makes the dynamics chaotic, which is verified by calculating Lyapunov exponents. In the nonrotating case, the primary rolls lose, depending on the aspect ratio, stability to traveling waves or a stationary square pattern. We also study the symmetries of the solutions at the intermittent fixed points in the heteroclinic cycle.
The dynamics of tail-like current sheets under the influence of small-scale plasma turbulence
(1999)
A 2D-magnetohydrodynamic model of current-sheet dynamics caused by anomalous electrical resistivity as result of small-scale plasma turbulence is proposed. The anomalous resistivity is assumed to be proportional to the square of the gradient of the magnetic pressure as may be valid for instance in the case of lower-hybrid-drift turbulence. The initial resistivity pulse is given. Then the temporal and spatial evolution of the magnetic and electric fields, plasma density, pressure, convection and resistivity are considered. The motion of the induced electric field is discussed as indicator of the plasma disturbances. The obtained results found using much improved numerical methods show a magnetic field evolution with x-line formation and plasma acceleration. Besides, in the current sheet, three types of magnetohydrodynamic waves occur, fast magnetoacoustic waves of compression and rarefaction as well as slow magnetoacoustic waves.
The nonlinear interaction of waves excited by the modified two-stream instability (Farley-Buneman instability) is considered. It is found that, during the linear stage of wave growth, the enhanced pressure of the high-frequency part of the waves locally generates a ponderomotive force. This force acts on the plasma particles and redistributes them. Thus an additional electrostatic polarization field occurs, which influences the low-frequency part of the waves. Then, the low-frequency waves also cause a redistribution of the high-frequency waves. In the paper, a self-consistent system of equations is obtained, which describes the nonlinear interaction of the waves. It is shown that the considered mechanism of wave interaction causes a nonlinear stabilization of the high-frequency waves’ growth and a formation of local density structures of the charged particles. The density modifications of the charged particles during the non-linear stage of wave growth and the possible interval of aspect angles of the high-frequency waves are estimated.
The ill-posed inversion of multiwavelength lidar data by a hybrid method of variable projection
(1999)
The ill-posed problem of aerosol distribution determination from a small number of backscatter and extinction lidar measurements was solved successfully via a hybrid method by a variable dimension of projection with B-Splines. Numerical simulation results with noisy data at different measurement situations show that it is possible to derive a reconstruction of the aerosol distribution only with 4 measurements.
Contents: 1 Introduction 2 Formation and destruction of sporadic E-layers 3 Temporal variations of parameters of sporadic E-layers during earthquake preparation 3.1 Temporal variations of fbEs with time-scales of a few hours 3.2 Study of fbEs variations with characteristic time-scales of 0.5-3 hours 3.3 Variations of the parameters of sporadic E-layers with characteristic time-scales of 15-45 minutes 3.4 Sporadic E-layer variations with characteristic time-scales of 2-15 minutes 4 On the spatial scales of sporadic E-layer disturbances related to seismic activity 5 Complex experimental researches of the ionosphere, electromagnetic noise and the geomagnetic field 5.1 Ionospheric and electromagnetic phenomena of the Kayraccum earthquake in 1985 5.2 Comparison of anomalies with characteristic time-scales of 2-3 hours for ionospheric E- and F-layers, and temporal behaviour of electromagnetic noise emission intensity 5.3 Night airglow emissions in the E-region before earthquakes and sporadic E-layer variations 6 Physical models of lithosphere-ionosphere links 6.1 Lithosphere-ionosphere links due to AGW 6.2 Electromagnetic models for the lithosphere-ionosphere coupling 6.3 Sporadic E-layers as current generators 7 Discussion and conclusion
In single photon emission computed tomography (SPECT) one is interested in reconstructing the activity distribution f of some radiopharmaceutical. The data gathered suffer from attenuation due to the tissue density µ. Each imaged slice incorporates noisy sample values of the nonlinear attenuated Radon transform (formular at this place in the original abstract) Traditional theory for SPECT reconstruction treats µ as a known parameter. In practical applications, however, µ is not known, but either crudely estimated, determined in costly additional measurements or plainly neglected. We demonstrate that an approximation of both f and µ from SPECT data alone is feasible, leading to quantitatively more accurate SPECT images. The result is based on nonlinear Tikhonov regularization techniques for parameter estimation problems in differential equations combined with Gauss-Newton-CG minimization.
The aim of this paper is to describe an efficient strategy for descritizing ill-posed linear operator equations of the first kind: we consider Tikhonov-Phillips-regularization χ^δ α = (a * a + α I)^-1 A * y ^δ with a finite dimensional approximation A n instead of A. We propose a sparse matrix structure which still leads to optimal convergences rates but requires substantially less scalar products for computing A n compared with standard methods.
The determination of the atmospheric aerosol size distribution is an inverse illposed problem. The shape and the material composition of the air-carried particles are two substantial model parameters. Present evaluation algorithms only used an approximation with spherical homogeneous particles. In this paper we propose a new numerically efficient recursive algorithm for inhomogeneous multilayered coated and absorbing particles. Numerical results of real existing particles show that the influence of the two parameters on the model is very important and therefore cannot be ignored.
We have numerically studied the bifurcations and transition to chaos in a two-dimensional fluid for varying values of the Reynolds number. These investigations have been motivated by experiments in fluids, where an array of vortices was driven by an electromotive force. In these experiments, successive changes leading to a complex motion of the vortices, due to increased forcing, have been explored [Tabeling, Perrin, and Fauve, J. Fluid Mech. 213, 511 (1990)]. We model this experiment by means of two-dimensional Navier-Stokes equations with a special external forcing, driving a linear chain of eight counter-rotating vortices, imposing stress-free boundary conditions in the vertical direction and periodic boundary conditions in the horizontal direction. As the strength of the forcing or the Reynolds number is raised, the original stationary vortex array becomes unstable and a complex sequence of bifurcations is observed. Several steady states and periodic branches and a period doubling cascade appear on the route to chaos. For increasing values of the Reynolds number, shear flow develops, for which the spatial scale is large compared to the scale of the forcing. Furthermore, we have investigated the influence of the aspect ratio of the container as well as the effect of no-slip boundary conditions at the top and bottom, on the bifurcation scenario.
This paper deals with the electrical conductivity problem in geophysics. It is formulated as an elliptic boundary value problem of second order for a large class of bounded and unbounded domains. A special boundary condition, the so called "Complete Electrode Model", is used. Poincaré inequalities are formulated and proved in the context of weighted Sobolev spaces, leading to existence and uniqueness statements for the boundary value problem. In addition, a parameter-to-solution operator arising from the inverse conductivity problem in medicine (EIT) and geophysics is investigated mathematically and is shown to be smooth and analytic.
The stability of the quiescent ground state of an incompressible, viscous and electrically conducting fluid sheet, bounded by stress-free parallel planes and driven by an external electric field tangential to the boundaries, is studied numerically. The electrical conductivity varies as cosh–2(x1/a), where x1 is the cross-sheet coordinate and a is the half width of a current layer centered about the midplane of the sheet. For a <~ 0.4L, where L is the distance between the boundary planes, the ground state is unstable to disturbances whose wavelengths parallel to the sheet lie between lower and upper bounds depending on the value of a and on the Hartmann number. Asymmetry of the configuration with respect to the midplane of the sheet, modelled by the addition of an externally imposed constant magnetic field to a symmetric equilibrium field, acts as a stabilizing factor.
The bifurcations in a three-dimensional incompressible, electrically conducting fluid with an external forcing of the Roberts type have been studied numerically. The corresponding flow can serve as a model for the convection in the outer core of the Earth and is realized in an ongoing laboratory experiment aimed at demonstrating a dynamo effect. The symmetry group of the problem has been determined and special attention has been paid to symmetry breaking by the bifurcations. The nonmagnetic, steady Roberts flow loses stability to a steady magnetic state, which in turn is subject to secondary bifurcations. The secondary solution branches have been traced until they end up in chaotic states.
Two-dimensional bouyancy-driven convection in a horizontal fluid layer with stress-free boundary conditions at top and bottom and periodic boundary conditions in the horizontal direction is investigated by means of numerical simulation and bifurcation-analysis techniques. As the bouyancy forces increase, the primary stationary and symmetric convection rolls undergo successive Hopf bifurcations, bifurcations to traveling waves, and phase lockings. We pay attention to symmetry breaking and its connection with the generation of large-scale horizontal flows. Calculations of Lyapunov exponents indicate that at a Rayleigh number of 2.3×105 no temporal chaos is reached yet, but the system moves nonchaotically on a 4-torus in phase space.
We investigate the cognitive control in polyrhythmic hand movements as a model paradigm for bimanual coordination. Using a symbolic coding of the recorded time series, we demonstrate the existence of qualitative transitions induced by experimental manipulation of the tempo. A nonlinear model with delayed feedback control is proposed, which accounts for these dynamical transitions in terms of bifurcations resulting from variation of the external control parameter. Furthermore, it is shown that transitions can also be observed due to fluctuations in the timing control level. We conclude that the complexity of coordinated bimanual movements results from interactions between nonlinear control mechanisms with delayed feedback and stochastic timing components.
The stability of the quiescent ground state of an incompressible viscous fluid sheet bounded by two parallel planes, with an electrical conductivity varying across the sheet, and driven by an external electric field tangential to the boundaries is considered. It is demonstrated that irrespective of the conductivity profile, as magnetic and kinetic Reynolds numbers (based on the Alfvén velocity) are raised from small values, two-dimensional perturbations become unstable first.
Three-dimensional bouyancy-driven convection in a horizontal fluid layer with stress-free boundary conditions at the top and bottom and periodic boundary conditions in the horizontal directions is investigated by means of numerical simulation and bifurcation-analysis techniques. The aspect ratio is fixed to a value of 2√2 and the Prandtl number to a value of 6.8. Two-dimensional convection rolls are found to be stable up to a Rayleigh number of 17 950, where a Hopf bifurcation leads to traveling waves. These are stable up to a Rayleigh number of 30 000, where a secondary Hopf bifurcation generates modulated traveling waves. We pay particular attention to the symmetries of the solutions and symmetry breaking by the bifurcations.
The ill-posed problem of aerosol size distribution determination from a small number of backscatter and extinction measurements was solved successfully with a mollifier method which is advantageous since the ill-posed part is performed on exactly given quantities, the points r where n(r) is evaluated may be freely selected. A new twodimensional model for the troposphere is proposed.
We have studied the bifurcation structure of the incompressible two-dimensional Navier-Stokes equations with a special external forcing driving an array of 8×8 counterrotating vortices. The study has been motivated by recent experiments with thin layers of electrolytes showing, among other things, the formation of large-scale spatial patterns. As the strength of the forcing or the Reynolds number is raised the original stationary vortex array becomes unstable and a complex sequence of bifurcations is observed. The bifurcations lead to several periodic branches, torus and chaotic solutions, and other stationary solutions. Most remarkable is the appearance of solutions characterized by structures on spatial scales large compared to the scale of the forcing. We also characterize the different dynamic regimes by means of tracers injected into the fluid. Stretching rates and Hausdorff dimensions of convected line elements are calculated to quantify the mixing process. It turns out that for time-periodic velocity fields the mixing can be very effective.
The dynamics of noisy bistable systems is analyzed by means of Lyapunov exponents and measures of complexity. We consider both the classical Kramers problem with additive white noise and the case when the barrier fluctuates due to additional external colored noise. In case of additive noise we calculate the Lyapunov exponents and all measures of complexity analytically as functions of the noise intensity resp. the mean escape time. For the problem of fluctuating barrier the usual description of the dynamics with the mean escape time is not sufficient. The application of the concept of measures of complexity allows to describe the structures of motion in more detail. Most complexity measures sign the value of correlation time at which the phenomenon of resonant activation occurs with an extremum.
The usage of nonlinear Galerkin methods for the numerical solution of partial differential equations is demonstrated by treating an example. We desribe the implementation of a nonlinear Galerkin method based on an approximate inertial manifold for the 3D magnetohydrodynamic equations and compare its efficiency with the linear Galerkin approximation. Special bifurcation points, time-averaged values of energy and enstrophy as well as Kaplan-Yorke dimensions are calculated for both schemes in order to estimate the number of modes necessary to correctly describe the behavior of the exact solutions.
Anhand eines paradigmatischen Modellbeispiels werden die Konsequenzen der Koexistenz vieler Attraktoren auf die globale Dynamik schwach dissipativer Systeme studiert. Es wird gezeigt, dass diese Systeme eine sehr reichhaltige Dynamik besitzen und extrem sensitiv gegenüber Störungen in den Anfangsbedingungen sind. Diese Systeme zeichnen sich durch eine extrem hohe Flexibilität ihres Verhaltens aus.